Skip to main content

EXTL3 could serve as a potential biomarker of prognosis and immunotherapy for prostate cancer and its potential mechanisms

Abstract

Background

Exostosin like glycosyltransferase 3 (EXTL3) had been reported to be associated with immune deficiency and play prognostic roles in various cancers. However, little is known about the associations between EXTL3 and prostate cancer (PCa). Hence, this article was designed to clarify their associations.

Methods

All original data were downloaded from The Cancer Genome Atlas (TCGA) database. Gene set enrichment analysis (GSEA) and CellMiner database was utilized, respectively, to identify EXTL3-related signaling pathways and drugs. We explored the relationships between EXTL3 expression and immunity to further evaluate the involvement of EXTL3 in response to immunotherapies. LncRNA/RBP/EXTL3 mRNA networks were also identified for its potential mechanism.

Results

Compared with normal prostate samples, EXTL3 was poorly expressed in PCa samples not only in mRNA expression levels, but also in protein expression levels, with worse overall survival (P < 0.05) and this gene could be an independent prognostic biomarker for PCa (both P < 0.05). EXTL3 was revealed to be markedly linked with seven signaling pathways in PCa by GSEA, including calcium, chemokine, ERBB, JAK STAT, MAPK, WNT, oxidative phosphorylation pathways. EXTL3 expression was also revealed to be significantly associated with MSI, immune cells, immune checkpoint molecules, tumor microenvironment and immune cells infiltration. We further predicted immune responses of EXTL3 gene to immunotherapies by TIDE database and the IMvigor210 cohort. A total of six LncRNA/RBP/EXTL3 mRNA networks were eventually identified for its potential mechanisms.

Conclusions

EXTL3 could serve as a potential biomarker of prognosis and immunotherapy for PCa and six LncRNA/RBP/EXTL3 mRNA networks were also identified for its potential mechanisms.

Introduction

Prostate cancer (PCa), as the most common male malignancy, ranks the first leading cause of newly estimated cases with 248,530 patients and the second leading cause of death with 34,130 cases in men, in the United States, 2021 [1]. PCa is asymptomatic in the early stage and patients usually seek for doctors when symptoms of hematuria or dysuria occur [2]. There are several risk factors commonly known for PCa including family history, age, ethnicity, lifestyle, exposure (occupational or environmental) and so on [3]. Since 1990s, prostate-specific antigen (PSA) was applied into clinical use and early detection rates of PCa had improved significantly [4]. However, several factors were also reported to be able to affect serum PSA levels, containing age, body mass index (BMI), obesity, hypertriglyceridemia, smoking and so on [5]. Currently, various therapeutic methods are available to PCa, including androgen deprivation therapy, radiotherapy, radical prostatectomy, active surveillance, immunotherapy and so on [6]. Therein, androgen deprivation therapy remained the major therapeutic option for PCa [7]. Although these advances in early diagnosis and further treatment of PCa, the 5-year survival rates in distant stages were merely 30.2% and patients were prone to relapse [8, 9]. Hence, there were broad prospects for exploring novel biomarkers to detect the occurrence, progress and relapse of PCa.

Exostosin like glycosyltransferase 3 (EXTL3), also known as RPR; BOTV; REGR; EXTR1; ISDNA; EXTL1L, was located at chromosome 8p21.1 and belonged to the EXT family, encoding glycosyltransferases to biosynthesize heparan sulfate and taking part in the signal transduction of islet-derived proteins, which played roles in stimulating the growth of islet β-cell, gut immune defenses, tissue regeneration and so on [10]. As reported, Volpi et al. shed light on that mutations in EXTL3 could lead to skeletal dysplasia, developmental delay, and severe immune deficiency [11]. Mizuno et al. found that EXTL3 was also able to enhance NF-κB activity mediated by TNF-α [12]. As for the relationships between EXTL3 and tumors, Karibe et al. suggested that EXTL3 promoter methylation could reduce EXTL3 expression and then result in the loss of heparan sulfate in colorectal cancer [13]. Hara et al. shed light on the ubiquitous expression of regenerating gene (REG) Iα receptor EXTL3 in gastric cancer vessel cells and tumor cells [14]. Zhang et al. revealed the potential receptors of EXTL3 for Reg3A in gastrointestinal cancer [15]. Taken together, above-mentioned researches indicated the important roles of EXTL3 in human diseases.

However, little was known about the associations between EXTL3 and PCa. Therefore, we comprehensively analyzed the prognostic roles of EXTL3 in PCa, based on The Cancer Genome Atlas (TCGA), The Human Protein Atlas (HPA) databases. Sangerbox website was utilized to assess the relationships between EXTL3 expression and immunity and Tumor Immune Dysfunction and Exclusion (TIDE) dataset was applied to predict immune responses. We also used starBase v2.0 to identify LncRNA/RNA binding proteins (RBP)/EXTL3 axes for its potential mechanisms. This study offered insights to provide a novel therapeutic target and evidences of EXTL3 in anticancer immunotherapy for PCa.

Materials and methods

Data download and EXTL3 gene expression

All original data in this article were downloaded from TCGA (TCGA-PRAD, http://cancergenome.nih.gov/) database, containing 499 PCa tumor tissues and 52 adjacent normal prostate tissues. When TCGA-PRAD samples were short of key clinical information, they shall be excluded. Then, we extracted the EXTL3 single gene matrix and corresponding clinical data for further analyses. TIMER database (https://cistrome.shinyapps.io/timer/) was utilized to reveal the EXTL3 mRNA expression in pan-cancer [16] and HPA database (https://www.proteinatlas.org/) was applied to show the EXTL3 protein expression in PCa tumor tissues and adjacent normal prostate tissues [17].

R software data processing

All analyses in this article were performed by R software version 4.1.1 (https://www.r-project.org/). The “limma” R package was utilized to display the EXTL3 mRNA distribution in PCa tumor tissues and adjacent normal prostate tissues, with the threshold of adjusted P-value (FDR) < 0.05 and |log2-fold change|≥ 1. K–M survival analysis was performed to see the survival differences in high-EXTL3 and low-EXTL3 subgroups, classified by its median expression. A total of 7 clinical factors (Gleason’s score, age, lymph nodes, staged T, cancer status, PSA value, staged N) and EXTL3 expression were enrolled in this article to conduct univariate/multivariate Cox regression analyses to find independent prognostic factors for PCa, with the cutoff criteria of P < 0.05 [18]. The “rms” R package was utilized to establish a nomogram based on Gleason’s score, age, lymph nodes, staged T, cancer status, PSA value, staged N, recurrence and EXTL3. C-index, 5-year ROC and calibration curves were carried out to assess its performance [19].

To predict the immune responses of EXTL3 to immunotherapies, TIDE database (http://tide.dfci.harvard.edu/) was utilized to calculate TIDE score, T cell dysfunction score, immune exclusion score in high-EXTL3 and low-EXTL3 subgroups, classified by its median expression. When patients had a higher TIDE score, they shall have a higher chance of immune exclusion, suggesting a less possibility of benefiting from immunotherapy [20]. We also selected the IMvigor210 cohort to predict the immune responses of EXTL3 to anti-PD-L1 treatment atezolizumab in different EXTL3 subgroups [21]. P < 0.05 were regarded as statistical significance. CellMiner database (https://discover.nci.nih.gov/cellminer/) was applied to reveal the associations between EXTL3 expression and drug sensitivity IC50 values, with the threshold of correlation coefficient  ≥ 0.3 and P < 0.01 [22].

Gene set enrichment analysis (GSEA) software and Sangerbox tools

GSEA software version 4.0.0 was downloaded from GSEA website (http://www.gsea-msigdb.org/gsea/index.jsp) and gene set of “c2.cp.kegg.v7.1.symbols.gmt” was obtained from Molecular Signatures Database (MSigDB). We utilized GSEA software to identify EXTL3-relevant signaling pathways, according to EXTL3 expression in TCGA-PRAD, with the cutoff criteria of |normalized enrichment score (NES)|> 1.5 as well as nominal P < 0.05 [23].

As previously published articles described [24, 25], we utilized the single gene pan-cancer tools in Sangerbox version 2.0 website (http://www.sangerbox.com/tool) to evaluate the relationships between EXTL3 expression based on TCGA-PRAD and immune cells, microsatellite instability (MSI), immune checkpoint molecules, tumor neoantigen burden (TNB), tumor microenvironment, tumor mutational burden (TMB), immune cells infiltration, with the cutoff criteria of P < 0.001. STRING database (https://www.string-db.org/) was also searched by us to display EXTL3-involved protein–protein interaction (PPI) networks [26].

StarBase v2.0 identified EXTL3-related mechanisms

As previously published articles described [27,28,29], we utilized the StarBase v2.0 website (https://starbase.sysu.edu.cn/) to identify EXTL3-related mechanisms of LncRNA/RBP/mRNA networks. We first searched EXTL3-targeted RBPs with the cutoff criteria of pan-cancer ≥ 10 as well as medium stringency (≥ 2). Then, the selected RBP-targeted LncRNAs were sought once again with the cutoff criteria of pan-cancer ≥ 10, hub LncRNAs in TCGA PRAD (LncRNAs were differently expressed in TCGA-PRAD with threshold of |log2-fold change|≥ 1, FDR < 0.05, and were significantly associated with survival with threshold of P < 0.05), LncRNAs positively correlated with EXTL3 in TCGA-PRAD (correlation coefficient ≥ 0.3 and P < 0.001), medium stringency (≥ 2). Cytoscape 3.6.1 software was finally utilized by us to display the LncRNA/RBP/EXTL3 mRNA networks.

Results

EXTL3 gene expression and survival prognosis in PCa

Figure 1A details EXTL3 mRNA expressions in pan-caner. It was poorly expressed in cancers including BLCA, BRCA, KIRC, PRAD and highly expressed in cancers containing CHOL, COAD, ESCA, LIHC, LUSC, STAD (all P < 0.05). Figure 1B exclusively shows the EXTL3 mRNA expression was poorly expressed in TCGA PCa samples, compared with normal prostate samples (P = 1.696e− 05). Survival analysis indicated that patients in high-EXTL3 group shall have a longer overall survival (OS) than patients in low-EXTL3 group (P = 0.004; Fig. 1C). Similar to EXTL3 mRNA expression, EXTL3 protein expression was poorly expressed in PCa tissues, compared with normal prostate tissues by CAB025387 antibody in HPA database (Fig. 1D, E). All of these indicated the anticarcinogenic role of EXTL3 in PCa.

Fig. 1
figure 1

EXTL3 gene expression and survival prognosis in PCa: A EXTL3 mRNA expressions in pan-cancer from TCGA; B EXTL3 mRNA expressions in PCa from TCGA; C survival analysis of EXTL3 in PCa from TCGA; immunohistochemistry of EXTL3 protein expression in D normal prostate tissues and E PCa tissues from HPA database; *P < 0.05; **P < 0.01;***P < 0.001

Prognostic values of EXTL3 in PCa from TCGA and nomogram establishment

Univariate Cox regression analysis showed that EXTL3 (HR = 0.201, 95%CI = 0.051–0.789, P = 0.021) and PSA value (HR = 1.214, 95%CI = 1.031–1.429, P = 0.020) were both prognostic factors in PCa from TCGA. Multivariate Cox regression analysis indicated that EXTL3 (HR = 0.119, 95%CI = 0.022–0.657, P = 0.015) was an independent prognostic biomarker for PCa from TCGA (Table 1). Univariate and multivariate Cox regression analyses revealed the independent prognostic ability of EXTL3 in PCa. To further predict PCa patients’ OS possibilities, we established a nomogram based on Gleason’s score, age, lymph nodes, staged T, cancer status, PSA value, staged N, recurrence and EXTL3, with the help of the “rms” R package (Fig. 2A). C-index as well as the area under the curve (AUC) value of 5-year ROC curve of this nomogram were 0.921, 0.896, respectively (Fig. 2B). 5-year calibration curves confirmed the good performance of this nomogram (Fig. 2C).

Table 1 Univariate and multivariate Cox analysis based on EXTL3 and clinicopathologic characteristics of OS for PCa
Fig. 2
figure 2

Nomogram establishment for PCa from TCGA: A nomogram; B 5-year ROC curve; C 5-year calibration curve

EXTL3-related signaling pathways in PCa from TCGA

Based on the gene set of “c2.cp.kegg.v7.1.symbols.gmt” obtained from MSigDB, we utilized the GSEA version 4.0.0 software to identify EXTL3-relevant signaling pathways in PCa. As shown in Fig. 3 and Table 2, EXTL3 was significantly associated with calcium pathway (NES = 2.511; nominal P < 0.001), chemokine pathway (NES = 2.284; nominal P < 0.001), ERBB pathway (NES = 2.303; nominal P < 0.001), JAK STAT pathway (NES = 2.281; nominal P < 0.001), MAPK pathway (NES = 2.335; nominal P < 0.001), WNT pathway (NES = 2.416; nominal P < 0.001), oxidative phosphorylation pathway (NES = − 2.251; nominal P < 0.001) in PCa. All of these indicated the potential pathways related to EXTL3 in PCa from TCGA.

Fig. 3
figure 3

EXTL3-related signaling pathways: A calcium pathway; B chemokine pathway; C ERBB pathway; D JAK STAT pathway; E MAPK pathway; F WNT pathway; G oxidative phosphorylation pathway; H all of these seven pathways

Table 2 Gene set enrichment analysis results

Correlations between EXTL3 expression and MSI, TNB, TMB, immunity

We utilized the Sangerbox version 2.0 website single gene pan-cancer tools to evaluate the relationships between EXTL3 expression and immune cells, microsatellite instability (MSI), immune checkpoint molecules, tumor neoantigen burden (TNB), tumor microenvironment, tumor mutational burden (TMB), immune cells infiltration. EXTL3-involved PPI networks are displayed in Fig. 4A. MSI was significantly related to EXTL3 expression (P = 0.00055) in TCGA-PRAD (Fig. 4B). TNB (P = 0.66) and TMB (P = 0.44) did not have significant associations with EXTL3 expression in TCGA-PRAD (Fig. 4C, D).

Fig. 4
figure 4

Correlations between EXTL3 expression and A PPI networks; B MSI; C TNB; D TMB

Four aspects of immunity were analyzed, including immune cells, immune checkpoint molecules, tumor microenvironment, immune cells infiltration. Figure 5A presents the correlations between EXTL3 expression and immune cells infiltration levels. We could find that EXTL3 expression was significantly linked with infiltration levels of dendritic cells, B cells, neutrophil cells, CD8+T cells, macrophage cells, CD4+T cells (all P < 0.001). As for tumor microenvironment, EXTL3 expression was also remarkably correlated with ESTIMATEScore, ImmuneScore, StromalScore (all P < 0.001; Fig. 5B). In terms of the correlations between EXTL3 expression and immune cells, immune checkpoint molecules, EXTL3 expression was markedly involved with CD274, BTLA, CTLA4, CD27, CD244, CD276, activated B cell, Activated CD4 T cell, activated dendritic cell, central memory CD4 T cell, MDSC, type 2 T helper cell and so on (all P < 0.05; Fig. 5C, D). All of these indicated significant correlations between EXTL3 expression and immunity.

Fig. 5
figure 5

Correlations between EXTL3 expression and A immune cells infiltration; B tumor microenvironment; C immune checkpoint molecules; D immune cells; *P < 0.05; **P < 0.01;***P < 0.001

The involvement of EXTL3 in response to immunotherapies

Based on previously published articles and above-mentioned results, EXTL3 had significant associations with immunity. Moreover, the distribution of EXTL3 in pan-cancer immune subtypes is detailed in Fig. 6A (P = 0.001). We utilized TIDE database to predict the involvement of EXTL3 in response to immunotherapies by calculating TIDE score, T cell dysfunction score, immune exclusion score in different EXTL3 subgroups. As presented in Fig. 6B–D, low-EXTL3 subgroup was associated with higher TIDE, T cell dysfunction, immune exclusion scores, suggesting a less possibility of benefiting from immunotherapy of low-EXTL3 subgroup patients. We also selected the IMvigor210 cohort to predict the immune responses of EXTL3 to anti-PD-L1 treatment atezolizumab in different EXTL3 subgroups. Figure 6E indicates that patients in high-EXTL3 subgroup were more sensitive to anti-PD-L1 treatment atezolizumab.

Fig. 6
figure 6

Immune responses of EXTL3 gene to immunotherapies: A distribution of EXTL3 in pan-cancer immune subtypes; B TIDE score in different EXTL3 subgroups; C T cell dysfunction score in different EXTL3 subgroups; D immune exclusion score in different EXTL3 subgroups; E immune responses of EXTL3 to anti-PD-L1 treatment atezolizumab in the IMvigor210 cohort. **P < 0.01;***P < 0.001

EXTL3-related drugs from CellMiner database

We applied CellMiner database to reveal the associations between EXTL3 expression and drug sensitivity IC50 values, with the threshold of correlation coefficient ≥ 0.3 and P < 0.01. As detailed in Fig. 7, EXTL3 expression was dramatically correlated with AFP464 (correlation coefficient = − 0.355; P = 0.005), AP-26113 (correlation coefficient = − 0.361; P = 0.005), dexrazoxane (correlation coefficient = − 0.341; P = 0.008), LDK-378 (correlation coefficient = − 0.355; P = 0.005), palbociclib (correlation coefficient = − 0.356; P = 0.005), rapamycin (correlation coefficient = 0.336; P = 0.009). All of these indicated that five negatively correlated drugs (AFP464, AP-26113, dexrazoxane, LDK-378, palbociclib) had the potential to prevent the progress of PCa, while rapamycin might promote this disease.

Fig. 7
figure 7

Associations between EXTL3 expression and A AFP464; B AP-26113; C Dexrazoxane; D LDK-378; E palbociclib; F rapamycin; in the CellMiner database

EXTL3-related mechanism of LncRNA/RBP/EXTL3 mRNA networks

We utilized the StarBase v2.0 website to identify EXTL3-related mechanisms of LncRNA/RBP/EXTL3 mRNA networks. The detailed flowchart is summarized in Fig. 8A. A total of 18 potential RBPs were identified after first searching RBP-targeted LncRNAs. However, only six LncRNA/RBP/EXTL3 mRNA networks (FENDRR/DGCR8/EXTL3 axis, FENDRR/HNRNPA1/EXTL3 axis, FENDRR/HNRNPC/EXTL3 axis, FENDRR/RBFOX2/EXTL3 axis, FENDRR/SRSF1/EXTL3 axis, FENDRR/U2AF2/EXTL3 axis) were remained after second searching a selected RBP-targeted LncRNAs, based on the results of Venn diagrams (Fig. 8B–G). Cytoscape 3.6.1 software was finally utilized by us to display these six LncRNA/RBP/EXTL3 mRNA networks (Fig. 8H). A total of six RBP-involved mechanisms of LncRNA/RBP/EXTL3 mRNA networks were eventually identified for EXTL3-related mechanisms of PCa in this article.

Fig. 8
figure 8

EXTL3-related mechanism of LncRNA/RBP/EXTL3 mRNA networks: A flowchart; B Venn diagrams of FENDRR/DGCR8/EXTL3 axis; C Venn diagrams of FENDRR/HNRNPA1/EXTL3 axis; D Venn diagrams of FENDRR/HNRNPC/EXTL3 axis; E Venn diagrams of FENDRR/RBFOX2/EXTL3 axis; F Venn diagrams of FENDRR/SRSF1/EXTL3 axis; G Venn diagrams of FENDRR/U2AF2/EXTL3 axis; H LncRNA/RBP/EXTL3 mRNA networks by Cytoscape 3.6.1

Discussion

PCa is an indolent cancer, characterized by its insidious onset, low drug sensitivity, high possibilities of metastasis or relapse, associated with a 5-year survival of less than 30% [30]. As we know, PCa patients’ morbidity and mortality in the United States and Europe were ranked first and second in the male population [31]. In developing countries, its incidence is also increasing year by year [32]. Moreover, as we had described previously, EXTL3 belonged to the EXT family, playing vital roles in predicting various cancers’ prognosis and immune deficiency [13, 14, 33]. However, associations between EXTL3 and PCa had seldom been explored. Hence, in this article, we comprehensively analyzed the prognostic roles of EXTL3 in PCa and assessed the relationships between EXTL3 expression and immunity to predict immune responses. We also identified LncRNA/RBP/EXTL3 axes for its potential mechanisms. Our findings were expected to provide a novel treatment target and evidence of EXTL3 in anticancer immunotherapy for PCa.

Our results shed light on that EXTL3 was poorly expressed in PCa samples not only in mRNA expression levels, but also in protein expression levels, compared with normal prostate samples. Moreover, PCa patients in high-EXTL3 groups shall have a longer OS than patients in poorly expressed groups and this gene could be an independent prognostic biomarker for PCa. Consistent with previously published articles, Zhang et al. revealed the potential receptors of EXTL3 for Reg3A in gastrointestinal cancer [15]. EXTL3 promoter methylation could reduce EXTL3 expression, resulting in the loss of heparan sulfate in colorectal cancer [13] and the ubiquitous expression of regenerating gene (REG) Iα receptor EXTL3 was found in gastric cancer vessel cells and tumor cells [14]. All of these studies indicated the prognostic value of EXTL3 in cancer.

As for TCGA data mining, nomogram was often utilized to intuitively display the relationships between clinical factors and patients’ survival [34, 35]. To further predict PCa patients’ OS possibilities, we established a nomogram in this article based on Gleason’s score, age, lymph nodes, staged T, cancer status, PSA value, staged N, recurrence and EXTL3, with good performance. For single gene, GSEA was commonly applied by us to underline the potential signaling pathways related to this gene [36, 37]. According to our research, GSEA also identified several EXTL3-relevant signaling pathways, including the calcium pathway, chemokine pathway, ERBB pathway, JAK STAT pathway, MAPK pathway, WNT pathway, oxidative phosphorylation pathway. As reported by Tam et al., the IL-6R/JAK/STAT3 signaling pathway activation could lead to the progress of hormone-refractory PCa [38]. Wang et al. reported that GADD45B could promote PCa’s chemosensitivity and cell apoptosis by MAPK signaling pathway [39]. Zhao et al. concluded that AGAP2-AS1/miR-628-5p/FOXP2 axis was able to facilitate the growth of PCa cells by means of WNT signaling pathway [40]. As for oxidative phosphorylation pathways, Xiao et al. revealed that phenethyl isothiocyanate might trigger reactive oxygen species (ROS)-mediated of PCa cells’ death via inhibiting oxidative phosphorylation [41].

MSI, TNB and TMB had been regarded as vital factors influencing cancer patients’ prognosis and they could also be as characterized as biomarkers of immunotherapies [42,43,44]. Based on our results, MSI was significantly related to EXTL3 expression in TCGA-PRAD, while TMB and TNB did not have significant associations with EXTL3 expression. To further seek EXTL3-related drugs in PCa, CellMiner database was utilized and this database was designed for facilitating the study and selection of anticancer drugs [45, 46]. Our results showed that EXTL3 expression was negatively correlated with AFP464, AP-26113, dexrazoxane, LDK-378, palbociclib, having the potential to prevent the progress of PCa. Moreover, EXTL3 expression was positively correlated with rapamycin, having the potential to promote this disease.

Currently, there had been various types of immunotherapies applied in the treatment of cancer, including cancer vaccines, oncolytic virus therapies, adoptive cell transfer, immune checkpoint inhibitors, cytokine therapies, and so on [47]. Despite the successful application of these immunotherapies, merely a part of patients with cancer could benefit from them [48]. Hence, there was also an urgent need to better understand tumor microenvironment. As reported, the interactions among tumor and immune cells in tumor microenvironment shall result in tumorigenesis, metastasis and relapse [49]. Moreover, immune cell infiltration, as an important part of tumor microenvironment, was significantly linked with tumor progression and could predict immune responses [50]. In our article, four aspects of immunity were analyzed, including immune cells, immune checkpoint molecules, tumor microenvironment, immune cells infiltration levels. As for immune cells infiltration levels, we could find that EXTL3 expression was significantly linked with infiltration levels of dendritic cells, B cells, neutrophil cells, CD8 + T cells, macrophage cells, CD4 + T cells. In terms of tumor microenvironment, EXTL3 expression was also remarkably correlated with ESTIMATEScore, ImmuneScore, StromalScore. For the correlations between EXTL3 expression and immune cells, immune checkpoint molecules, EXTL3 expression was markedly involved with CD274, BTLA, CTLA4, CD27, CD244,CD276, activated B cell, activated CD4 T cell, activated dendritic cell, central memory CD4 T cell, MDSC, Type 2 T helper cell and so on. All of these indicated significant associations between EXTL3 with immunity in PCa.

Although significant responses to immunotherapies had been showed in other metastatic cancers including melanoma, renal cancer or lung cancer, only a small number of PCa sufferers showed responses to immunotherapies, due to complex tumor environment interactions among immune and malignant PCa cells. Currently, immunotherapies for PCa included immune checkpoint inhibitors and vaccine-based therapies [51]. As reported, the only approved PCa immunotherapy was sipuleucel-T vaccine therapy. Moreover, Ipilimumab or olaparib also showed survival benefits in PCa [52]. All of these indicated the urgent need to identify biomarkers for PCa immunotherapy. In this article, we further predicted immune responses of EXTL3 gene to immunotherapies, by means of TIDE database and the IMvigor210 cohort [53, 54]. Based on our results, the low-EXTL3 subgroup was associated with higher TIDE, T cell dysfunction, immune exclusion scores, suggesting a lesser possibility of these patients benefiting from immunotherapy. We also selected the IMvigor210 cohort to predict the immune responses of EXTL3 to anti-PD-L1 treatment atezolizumab in different EXTL3 subgroups, indicating that patients in high-EXTL3 subgroup were more sensitive to anti-PD-L1 treatment atezolizumab. All in all, patients with high EXTL3 expression were prone to benefit from immunotherapy.

Based on above-mentioned, we could easily find that EXTL3 served as an anti-oncogene in PCa. In order to further reveal EXTL3-related mechanisms, we referred to previously published articles as well as the functions of LncRNAs or RBPs [55,56,57]. Finally, the RBP-involved mechanisms of LncRNA/RBP/mRNA networks were revealed. As previously reported, Zou et al. found that the LINC00324/HuR (RBP)/FAM83B axis could promote the gastric cancer cells’ proliferation [28]. Yu et al. revealed that circ_0003258 could drive prostate cancer metastasis via both a competing endogenous RNA (ceRNA) mechanism of circ_0003258/miR-653-5p/ARHGAP5 and a RBP-involved mechanism of circ_0003258/IGF2BP3 (RBP)/HDAC4 [58]. In this article, the RBP-involved mechanisms of LncRNA/RBP/EXTL3 mRNA networks were identified by us with the help of StarBase v2.0 website. Therein, a total of six LncRNA/RBP/EXTL3 mRNA networks containing the FENDRR/DGCR8/EXTL3 axis, the FENDRR/HNRNPA1/EXTL3 axis, the FENDRR/HNRNPC/EXTL3 axis, the FENDRR/RBFOX2/EXTL3 axis, the FENDRR/SRSF1/EXTL3 axis, the FENDRR/U2AF2/EXTL3 axis, were eventually identified in this article for EXTL3-related mechanisms for PCa.

Several limitations should not be ignored too. Firstly, all analyses were based on bioinformatics analysis, without experimental verification. Further in vivo or in vitro experiments were required to validate the expression of EXTL3, EXTL3-related signaling pathways and its related mechanisms of LncRNA/RBP/EXTL3 mRNA networks. Secondly, most of these analyses were based on EXTL3 mRNA expressions. This might not be consistent with its protein expressions. Thirdly, other information (such as prostate cancer family history, immunotherapy, and so on) was limited in TCGA database. Hence, we currently had difficulties in analyzing these factors. Finally, no targeted EXTL3 therapy was currently available. Hence, we could not directly obtain clinical information and immune responses of this targeted therapy, affecting our results to some certain extent.

Conclusions

In summary, EXTL3 could serve as an anti-oncogene in PCa and it was found to be markedly linked with seven signaling pathways in PCa by GSEA, including calcium, chemokine, ERBB, JAK STAT, MAPK, WNT, oxidative phosphorylation pathways. EXTL3 expression was also revealed to be significantly associated with MSI, immune cells, immune checkpoint molecules, tumor microenvironment and immune cells infiltration. We further predicted the involvement of EXTL3 in response to immunotherapies by TIDE database and the IMvigor210 cohort. A total of six LncRNA/RBP/EXTL3 mRNA networks including the FENDRR/DGCR8/EXTL3 axis, the FENDRR/HNRNPA1/EXTL3 axis, the FENDRR/HNRNPC/EXTL3 axis, the FENDRR/RBFOX2/EXTL3 axis, the FENDRR/SRSF1/EXTL3 axis, the FENDRR/U2AF2/EXTL3 axis, were eventually identified in this article for its potential mechanisms.

Availability of data and materials

The RNA-sequencing data and corresponding clinical information were downloaded from the Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/).

Abbreviations

EXTL3:

Exostosin like glycosyltransferase 3

PCa:

Prostate cancer

TCGA:

The Cancer Genome Atlas

GSEA:

Gene set enrichment analysis

OS:

Overall survival

BMI:

Body mass index

REG:

Regenerating gene

HPA:

The Human Protein Atlas

TIDE:

Tumor immune dysfunction and exclusion

RBP:

RNA binding proteins

FDR:

Adjusted P-value

MSigDB:

Molecular signatures database

NES:

Normalized enrichment score

TNB:

Tumor neoantigen burden

TMB:

Tumor mutational burden

MSI:

Microsatellite instability

PPI:

Protein–protein interaction

BLCA:

Bladder urothelial carcinoma

BRCA:

Breast invasive carcinoma

CHOL:

Cholangiocarcinoma

COAD:

Colon adenocarcinoma

ESCA:

Esophageal carcinoma

HNSC:

Head and neck squamous cell carcinoma

KIRC:

Kidney renal clear cell carcinoma

LIHC:

Liver hepatocellular carcinoma

LUSC:

Lung squamous cell carcinoma

PRAD:

Prostate adenocarcinoma

SKCM:

Skin cutaneous melanoma

STAD:

Stomach adenocarcinoma

References

  1. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics. CA. 2021;71(1):7–33.

    PubMed  Google Scholar 

  2. Miller KD, Siegel RL, Lin CC, Mariotto AB, Kramer JL, Rowland JH, Stein KD, Alteri R, Jemal A. Cancer treatment and survivorship statistics. CA. 2016;66(4):271–89.

    PubMed  Google Scholar 

  3. Crocetto F, di Zazzo E, Buonerba C, Aveta A, Pandolfo SD, Barone B, Trama F, Caputo VF, Scafuri L, Ferro M, et al. Kaempferol, myricetin and fisetin in prostate and bladder cancer: a systematic review of the literature. Nutrients. 2021;13(11):3750.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Center MM, Jemal A, Lortet-Tieulent J, Ward E, Ferlay J, Brawley O, Bray F. International variation in prostate cancer incidence and mortality rates. Eur Urol. 2012;61(6):1079–92.

    Article  PubMed  Google Scholar 

  5. Tarantino G, Crocetto F, Vito CD, Martino R, Pandolfo SD, Creta M, Aveta A, Buonerba C, Imbimbo C. Clinical factors affecting prostate-specific antigen levels in prostate cancer patients undergoing radical prostatectomy: a retrospective study. Future Sci OA. 2021;7(3):643.

    Article  CAS  Google Scholar 

  6. Malouff TD, Stross WC, Seneviratne DS, Waddle MR, May BC, Buskirk SJ, Tzou KS. Current use of stereotactic body radiation therapy for low and intermediate risk prostate cancer: a national cancer database analysis. Prostate Cancer Prostatic Dis. 2020;23(2):349–55.

    Article  PubMed  Google Scholar 

  7. Di Zazzo E, Galasso G, Giovannelli P, Di Donato M, Di Santi A, Cernera G, Rossi V, Abbondanza C, Moncharmont B, Sinisi AA, et al. Prostate cancer stem cells: the role of androgen and estrogen receptors. Oncotarget. 2016;7(1):193–208.

    Article  PubMed  Google Scholar 

  8. Reynolds T. “User-friendly” cancer statistics review features new format. J Natl Cancer Inst. 1992;84(19):1470–1.

    Article  CAS  PubMed  Google Scholar 

  9. Zelefsky MJ, Eastham JA, Cronin AM, Fuks Z, Zhang Z, Yamada Y, Vickers A, Scardino PT. Metastasis after radical prostatectomy or external beam radiotherapy for patients with clinically localized prostate cancer: a comparison of clinical cohorts adjusted for case mix. J Clin Oncol. 2010;28(9):1508–13.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Yamada S. Specific functions of Exostosin-like 3 (EXTL3) gene products. Cell Mol Biol Lett. 2020;25:39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Volpi S, Yamazaki Y, Brauer PM, van Rooijen E, Hayashida A, Slavotinek A, Sun Kuehn H, Di Rocco M, Rivolta C, Bortolomai I, et al. EXTL3 mutations cause skeletal dysplasia, immune deficiency, and developmental delay. J Exp Med. 2017;214(3):623–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Mizuno K, Irie S, Sato TA. Overexpression of EXTL3/EXTR1 enhances NF-kappaB activity induced by TNF-alpha. Cell Signal. 2001;13(2):125–30.

    Article  CAS  PubMed  Google Scholar 

  13. Karibe T, Fukui H, Sekikawa A, Shiratori K, Fujimori T. EXTL3 promoter methylation down-regulates EXTL3 and heparan sulphate expression in mucinous colorectal cancers. J Pathol. 2008;216(1):32–42.

    Article  CAS  PubMed  Google Scholar 

  14. Hara K, Fukui H, Sun C, Kitayama Y, Eda H, Yamasaki T, Kondo T, Tomita T, Oshima T, Watari J, et al. Effect of REG Iα protein on angiogenesis in gastric cancer tissues. Oncol Rep. 2015;33(5):2183–9.

    Article  CAS  PubMed  Google Scholar 

  15. Zhang MY, Wang J, Guo J. Role of regenerating islet-derived protein 3A in gastrointestinal cancer. Front Oncol. 2019;9:1449.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Can Res. 2017;77(21):e108–10.

    Article  CAS  Google Scholar 

  17. Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson Å, Kampf C, Sjöstedt E, Asplund A, et al. Proteomics tissue-based map of the human proteome. Science. 2015;347(6220):1260419.

    Article  PubMed  CAS  Google Scholar 

  18. Wu H, Zhang J. Decreased expression of TFAP2B in endometrial cancer predicts poor prognosis: a study based on TCGA data. Gynecol Oncol. 2018;149(3):592–7.

    Article  CAS  PubMed  Google Scholar 

  19. Lin WP, Xing KL, Fu JC, Ling YH, Li SH, Yu WS, Zhang YF, Zhong C, Wang JH, Chen ZY, et al. Development and validation of a model including distinct vascular patterns to estimate survival in hepatocellular carcinoma. JAMA Netw Open. 2021;4(9):e2125055.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Cheng Q, Duan W, He S, Li C, Cao H, Liu K, Ye W, Yuan B, Xia Z. Multi-omics data integration analysis of an immune-related gene signature in lgg patients with epilepsy. Front Cell Dev Biol. 2021;9:686909.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, Kadel EE III, Koeppen H, Astarita JL, Cubas R, et al. TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature. 2018;554(7693):544–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Reinhold WC, Sunshine M, Liu H, Varma S, Kohn KW, Morris J, Doroshow J, Pommier Y. Cell Miner: a web-based suite of genomic and pharmacologic tools to explore transcript and drug patterns in the NCI-60 cell line set. Can Res. 2012;72(14):3499–511.

    Article  CAS  Google Scholar 

  23. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61.

    Article  CAS  PubMed  Google Scholar 

  24. Liu S, Wang Y, Miao C, Xing Q, Wang Z. High expression of CDCA7 predicts poor prognosis for clear cell renal cell carcinoma and explores its associations with immunity. Cancer Cell Int. 2021;21(1):140.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wang G, Qu F, Liu S, Zhou J, Wang Y. Nucleolar protein NOP2 could serve as a potential prognostic predictor for clear cell renal cell carcinoma. Bioengineered. 2021;12(1):4841–55.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607-d613.

    Article  CAS  PubMed  Google Scholar 

  27. Zhang Y, Huang Z, Sheng F, Yin Z. MYC upregulated LINC00319 promotes human acute myeloid leukemia (AML) cells growth through stabilizing SIRT6. Biochem Biophys Res Commun. 2019;509(1):314–21.

    Article  CAS  PubMed  Google Scholar 

  28. Zou Z, Ma T, He X, Zhou J, Ma H, Xie M, Liu Y, Lu D, Di S, Zhang Z. Long intergenic non-coding RNA 00324 promotes gastric cancer cell proliferation via binding with HuR and stabilizing FAM83B expression. Cell Death Dis. 2018;9(7):717.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Wang A, Bao Y, Wu Z, Zhao T, Wang D, Shi J, Liu B, Sun S, Yang F, Wang L, et al. Long noncoding RNA EGFR-AS1 promotes cell growth and metastasis via affecting HuR mediated mRNA stability of EGFR in renal cancer. Cell Death Dis. 2019;10(3):154.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Barry MJ, Simmons LH. Prevention of prostate cancer morbidity and mortality: primary prevention and early detection. Med Clin N Am. 2017;101(4):787–806.

    Article  PubMed  Google Scholar 

  31. Pezaro C, Woo HH, Davis ID. Prostate cancer: measuring PSA. Intern Med J. 2014;44(5):433–40.

    Article  CAS  PubMed  Google Scholar 

  32. Zhang J, Chen Z, Wu J, Liu K. An intelligent decision-making support system for the detection and staging of prostate cancer in developing countries. Comput Math Methods Med. 2020;2020:5363549.

    PubMed  PubMed Central  Google Scholar 

  33. Van Hul W, Wuyts W, Hendrickx J, Speleman F, Wauters J, De Boulle K, Van Roy N, Bossuyt P, Willems PJ. Identification of a third EXT-like gene (EXTL3) belonging to the EXT gene family. Genomics. 1998;47(2):230–7.

    Article  PubMed  Google Scholar 

  34. He L, Wang X, Jin Y, Xu W, Lyu J, Guan Y, Wu J, Han S, Liu G. A prognostic nomogram for predicting overall survival in pediatric Wilms tumor based on an autophagy-related gene signature. Comb Chem High Throughput Screen. 2021. https://doi.org/10.2174/1386207324666210826143727.

    Article  PubMed  Google Scholar 

  35. Ye D, Qu J, Wang J, Li G, Sun B, Xu Q. A new clinical nomogram from the TCGA database to predict the prognosis of hepatocellular carcinoma. Front Oncol. 2021;11:698980.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Chen R, Cao J, Jiang W, Wang S, Cheng J. Upregulated expression of CYBRD1 predicts poor prognosis of patients with ovarian cancer. J Oncol. 2021;2021:7548406.

    PubMed  PubMed Central  Google Scholar 

  37. Zhu J, Luo C, Zhao J, Zhu X, Lin K, Bu F, Yu Z, Zou F, Zhu Z. Expression of LOX suggests poor prognosis in gastric cancer. Front Med. 2021;8:718986.

    Article  Google Scholar 

  38. Tam L, McGlynn LM, Traynor P, Mukherjee R, Bartlett JM, Edwards J. Expression levels of the JAK/STAT pathway in the transition from hormone-sensitive to hormone-refractory prostate cancer. Br J Cancer. 2007;97(3):378–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Wang Q, Wu W, Gao Z, Li K, Peng S, Fan H, Xie Z, Guo Z, Huang H. GADD45B is a potential diagnostic and therapeutic target gene in chemotherapy-resistant prostate cancer. Front Cell Dev Biol. 2021;9:716501.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Zhao X, Liu Y, Luo C, Zuo Y. AGAP2-AS1/miR-628-5p/FOXP2 feedback loop facilitates the growth of prostate cancer via activating WNT pathway. Carcinogenesis. 2021. https://doi.org/10.1093/carcin/bgab062.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Xiao D, Powolny AA, Moura MB, Kelley EE, Bommareddy A, Kim SH, Hahm ER, Normolle D, Van Houten B, Singh SV. Phenethyl isothiocyanate inhibits oxidative phosphorylation to trigger reactive oxygen species-mediated death of human prostate cancer cells. J Biol Chem. 2010;285(34):26558–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Miao Y, Wang J, Li Q, Quan W, Wang Y, Li C, Wu J, Mi D. Prognostic value and immunological role of PDCD1 gene in pan-cancer. Int Immunopharmacol. 2020;89(Pt B): 107080.

    Article  CAS  PubMed  Google Scholar 

  43. Cui G, Wang C, Lin Z, Feng X, Wei M, Miao Z, Sun Z, Wei F. Prognostic and immunological role of Ras-related protein Rap1b in pan-cancer. Bioengineered. 2021;12(1):4828–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Hu J, Qiu D, Yu A, Hu J, Deng H, Li H, Yi Z, Chen J, Zu X. YTHDF1 is a potential pan-cancer biomarker for prognosis and immunotherapy. Front Oncol. 2021;11:607224.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Shankavaram UT, Varma S, Kane D, Sunshine M, Chary KK, Reinhold WC, Pommier Y, Weinstein JN. Cell miner: a relational database and query tool for the NCI-60 cancer cell lines. BMC Genomics. 2009;10:277.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Zhang X, Huang T, Li Y, Qiu H. Upregulation of THBS1 is related to immunity and chemotherapy resistance in gastric cancer. Int J Gen Med. 2021;14:4945–57.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Zhang Y, Zhang Z. The history and advances in cancer immunotherapy: understanding the characteristics of tumor-infiltrating immune cells and their therapeutic implications. Cell Mol Immunol. 2020;17(8):807–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Darvin P, Toor SM, Sasidharan Nair V, Elkord E. Immune checkpoint inhibitors: recent progress and potential biomarkers. Exp Mol Med. 2018;50(12):1–11.

    Article  PubMed  CAS  Google Scholar 

  49. Wang Z, Song K, Zhao W, Zhao Z. Dendritic cells in tumor microenvironment promoted the neuropathic pain via paracrine inflammatory and growth factors. Bioengineered. 2020;11(1):661–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Topalian SL, Drake CG, Pardoll DM. Immune checkpoint blockade: a common denominator approach to cancer therapy. Cancer Cell. 2015;27(4):450–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Fay EK, Graff JN. Immunotherapy in prostate cancer. Cancers. 2020. https://doi.org/10.3390/cancers12071752.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Mitsogiannis I, Tzelves L, Dellis A, Issa H, Papatsoris A, Moussa M. Prostate cancer immunotherapy. Expert Opin Biol Ther. 2022;22(5):577–90.

    Article  CAS  PubMed  Google Scholar 

  53. Zhang B, Wu Q, Li B, Wang D, Wang L, Zhou YL. m(6)A regulator-mediated methylation modification patterns and tumor microenvironment infiltration characterization in gastric cancer. Mol Cancer. 2020;19(1):53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Xu F, Huang X, Li Y, Chen Y, Lin L. m(6)A-related lncRNAs are potential biomarkers for predicting prognoses and immune responses in patients with LUAD. Mol Ther Nucleic Acids. 2021;24:780–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Pereira B, Billaud M, Almeida R. RNA-binding proteins in cancer: old players and new actors. Trends Cancer. 2017;3(7):506–28.

    Article  CAS  PubMed  Google Scholar 

  56. Jiang MC, Ni JJ, Cui WY, Wang BY, Zhuo W. Emerging roles of lncRNA in cancer and therapeutic opportunities. Am J Cancer Res. 2019;9(7):1354–66.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Wang Y, Tian Y, Liu S, Wang Z, Xing Q. Prognostic value and immunological role of AXL gene in clear cell renal cell carcinoma associated with identifying LncRNA/RBP/AXL mRNA networks. Cancer Cell Int. 2021;21(1):625.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Yu YZ, Lv DJ, Wang C, Song XL, Xie T, Wang T, Li ZM, Guo JD, Fu DJ, Li KJ, et al. Hsa_circ_0003258 promotes prostate cancer metastasis by complexing with IGF2BP3 and sponging miR-653-5p. Mol Cancer. 2022;21(1):12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Contributions

PAC: manuscript writing/editing; SLC: data collection or management; XMC: data analysis; QST, LMM: protocol/project development. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Qingsheng Tang or Limin Ma.

Ethics declarations

Ethics approval and consent to participate

None declared.

Consent for publications

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chang, P., Chen, S., Chang, X. et al. EXTL3 could serve as a potential biomarker of prognosis and immunotherapy for prostate cancer and its potential mechanisms. Eur J Med Res 27, 115 (2022). https://doi.org/10.1186/s40001-022-00740-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-022-00740-w

Keywords