Skip to main content

Identification of NTRK3 as a potential prognostic biomarker associated with tumor mutation burden and immune infiltration in bladder cancer



Bladder cancer (BLCA) is a common malignant tumor of urinary system with high morbidity and mortality. In recent years, immunotherapy has played a significant role in the treatment of BLCA. Tumor mutation burden (TMB) has been reported to be a powerful biomarker for predicting tumor prognosis and efficacy of immunotherapy. Our study aimed to explore the relationship between TMB, prognosis and immune infiltration to identify the key genes in BLCA.


Clinical information, somatic mutation and gene expression data of BLCA patients were downloaded from The Cancer Genome Atlas (TCGA) database. Patients were divided into high and low TMB groups according to their calculated TMB scores. Gene Set Enrichment Analysis (GSEA) was performed to screen for significantly enriched pathways. Differentially expressed genes (DEGs) between the two groups were identified. Univariate Cox analysis and Kaplan-Meier survival analysis were applied for screening key genes. Immune infiltration was performed for TMB groups and NTRK3.


Higher TMB scores were related with poor survival in BLCA. After filtering, 36 DEGs were identified. NTRK3 had the highest hazard ratio and significant prognostic value. Co-expressed genes of NTRK3 were mainly involved in several pathways, including DNA replication, basal transcription factors, complement and coagulation cascades, and ribosome biogenesis in eukaryotes. There was a significant correlation among TMB scores, NTRK3 expression and immune infiltration.


Our results suggest that NTRK3 is a TMB-related prognostic biomarker, which lays the foundation for further research on the immunomodulatory effect of NTRK3 in BLCA.

Peer Review reports


Bladder cancer (BLCA) is one of the most common urinary malignancies and its incidence is gradually increasing, accounting for nearly 550,000 new cases and 200,000 deaths annually worldwide [1]. 75% of newly diagnosed patients are nonmuscle-invasive bladder cancers (NMIBCs), and although the 5-year overall survival (OS) rate is as high as 90% after treatment, NMIBC has a high recurrence rate, presenting a large social and economic burden [2, 3]. Muscle-invasive bladder cancer (MIBC) represents approximately 20% of newly diagnosed cases, and approximately 15 to 20% of NMIBCs can progress to MIBCs [4]. Despite radical cystectomy and neoadjuvant chemotherapy are available, the MIBC patients have a poor prognosis with a 5-year OS rate of less than 50%. In addition, about half of the patients will eventually progress to distant metastases [5, 6].

In the past 30 years, the systemic treatment for BLCA has remained unchanged and the curative effect has not made a breakthrough, especially in advanced disease [7]. Nowadays, with the rapid evolution of immunotherapy, the use of immune checkpoint inhibitors (ICIs) is becoming emerging as a new treatment strategy for advanced BLCA [8]. Treatment with ICIs such as anti-programmed cell death protein 1(PD-1), anti-programmed death-ligand 1 (PD-L1), and anti-cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) can significantly increase OS and result in durable remission for many advanced cancers, including melanoma, non-small cell lung cancer, renal cancer and urothelial cancer [9,10,11,12]. However, these effects only exist in a small subset of patients who can respond to ICIs. Many studies have shown that the tumor mutation burden (TMB) is a powerful biomarker for predicting the efficacy of ICIs and can be used to identify patients who will benefit from immunotherapy [13]. Tumors can be caused by the accumulation of somatic mutations in cells as a result of carcinogens. Some mutant cells revert to normal via DNA self-modification, some die, and a minority of them express neoantigens on surface. Tumor cells can prevent these mutant cells from being recognized by the immune system through abnormal expressions of antigens or regulating the tumor microenvironment (TME), thereby achieving immune escape [14]. A higher TMB represents more neoantigens and increases the probability of activating the body’s anti-tumor immune response, which can improve the patients’ response to ICIs [15, 16]. Thus, exploring the relationship between TMB and immunoregulation and discovering the related genes or biological mechanisms will help to better understand the role of immunotherapy in cancer treatment.

In our study, mutation information and expression profiles were acquired from The Cancer Genome Atlas (TCGA) database. Bioinformatic analysis was applied to calculate TMB scores and evaluate their relationships with prognosis and immune infiltration. Finally, we identified one key potential biomarker with prognostic value, which might help to filter BLCA patients suitable for immunotherapy.


Data source and data processing

Transcriptome data and clinical information of BLCA patients were downloaded from the TCGA data portal ( Validation cohort GSE48276 was downloaded from the GEO database ( Then we downloaded the simple nucleotide variation data from the “Masked Somatic Mutation” category processed with VarScan software in TCGA. The “maftools” R package was used to analyze the Mutation Annotation Format (MAF) file and visualize the somatic mutation data [17].

Calculation of TMB scores and prognostic evaluation

TMB is also defined as mutation frequency. We calculated TMB scores with the number of all nonsynonymous mutations, divided by exon length (38 Mb). We used the X-title software to calculate the best cut-off value of TMB scores and divided the patients into high and low TMB groups [18]. Survival analysis was conducted while clinicopathological characteristics were compared between high and low TMB groups. The Wilcoxon rank-sum test was for comparing two groups and the Kruskal–Wallis test was for more groups.

Gene set enrichment analysis (GSEA) between high and low TMB groups

GSEA version was 4.0.0 based on JAVA platform and the reference gene set (c2.cp.kegg.v7.0.symbols.gmt) was downloaded from the Molecular Signatures Database (MSigDB) ( [19]. The pathways were considered to be statistically enriched with P < 0.05 and FDR (false discovery rate) < 0.25.

Screening of key genes associated with TMB and survival

We identified the differentially expressed genes (DEGs) between high and low TMB groups. The “edgeR,” “limma,” and “DEseq2” packages in R software were used for this analysis. Univariate Cox and Kaplan-Meier (KM) methods were performed for survival analysis on every DEG to obtain the target genes. |log2FC| > 1 and P value < 0.05 were set as the thresholds.

Coexpression analysis of NTRK3

LinkedOmics ( is a publicly available portal that includes multi-omics data from all 32 TCGA cancer types and provides a web-based platform for analysis [20]. Coexpression of NTRK3 was analyzed statistically on the website using the Pearson correlation coefficient and the results were visualized by volcano plot and heat map. The LinkInterpreter module performed GSEA to obtain the related Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The rank criterion was P < 0.05 and 500 simulations were tested.

Immune infiltration analysis

The “estimate” R package was used to calculate scores for tumor purity, the level of stromal cells presents, and the infiltration level of immune cells in tumor tissues based on expression data. The CIBERSORT script was used to analyze 22 types of immune cell fractions in each patient [21]. Samples with a P value ≥0.05 were excluded from the next step. Wilcoxon rank-sum tests were conducted to analyze the different abundances of 22 immune cell types between high and low TMB groups. The TIMER database ( was used to analyze the correlation between NTRK3 expression and immune cell infiltration level [22]. The TISIDB web portal ( was applied to explore the association between NTRK3 and multiple immune regulatory factors [23].


Overview of mutation information in BLCA patients

We analyzed and visualized the somatic mutation information of BLCA patients using the “maftools” package. Among all the mutation types, missense mutation is the most common, with a frequency far higher than that of other types (Fig. 1a). Single nucleotide polymorphism (SNP) represented the largest fraction in the variant type than insertion or deletion (Fig. 1b). The most common single nucleotide variant (SNV) was C > T, followed by C > G and C > A (Fig. 1c). The waterfall plot showed the top 15 most frequently mutated genes (Fig. 1d). Figure 1e shows the cooccurrence and exclusive associations between mutated genes, in which TP53 was tightly associated with FGFR3 and RB1.

Fig. 1
figure 1

Overview of mutation information in BLCA patients. a-c Variant classifications, variant types and SNV classes in BLCA samples. SNP: single nucleotide polymorphism, INS: insertion, DEL: deletion, SNV: single nucleotide variants. d Waterfall plot showing the mutation information for top 15 genes. e The coincident and exclusive associations across the top 30 mutated genes

Clinical relevance and pathway analysis for TMB

The TMB score of each BLCA patient was calculated and the best cut-off value of TMB was 6.4 computed by X-title software (Fig. 2a). Then we divided the patients into a high TMB group with 164 cases and a low TMB group with 244 cases based on this cut-off. The result of KM survival analysis showed that the overall survival of the high TMB group was significantly better than the low TMB group (Fig. 2b, log-rank P < 0.0001). Clinical correlation analysis revealed that high-grade advanced tumors had higher TMB scores than low-grade tumors (Fig. 2e, P < 0.01). However, no significant associations were found between TMB scores and clinical stage, T stage, or smoking level (Figs. 2c, d and f).

Fig. 2
figure 2

Clinical relevance and pathways analysis for TMB. a The best cut-off value of TMB was 6.4 computed by X-title. b KM curves of overall survival of high and low TMB groups. c-f The correlation between TMB scores and clinical stage, T stage, pathology grade and smoking level. g The top 4 results of GSEA between high and low TMB patients

In addition, GSEA was performed between high and low TMB groups to explore the biological pathways that may be affected by TMB. The results showed that BLCA samples in the low TMB group were significantly enriched for 41 biological processes, and the following four biological processes with high normalized enrichment scores and large sizes of gene counts were selected, including ECM (extracellular matrix) receptor interaction, focal adhesion, leukocyte transendothelial migration and MAPK signaling pathway (Fig. 2g). But for the high TMB group, no significant enrichment was observed.

Screening and identification of key DEGs between high and low TMB groups

Three methods were applied to screen for DEGs between high and low TMB groups. The intersection among the results of the three methods is presented in Venn diagram and Upset plot (Fig. 3a). Thirty-six genes were selected and the expression level of each gene in high and low TMB groups was shown in the form of heatmap (Fig. 3b). Univariate Cox analysis for each of the 36 genes showed that 10 genes were significantly related to overall survival (P < 0.05, Table 1). Then we performed KM survival analysis for these 10 genes and only four genes were with log-rank P < 0.05, including NCAM1, SPON1, NTRK3 and PRND (Table 1). The Forest plot showed the result of univariate Cox analysis for these genes (Fig. 3c). Figure 3d-g shows the KM survival curves for them. Among them, NTRK3 had the biggest hazard ratio (HR) value, indicating that the high expression of NTRK3 may be a risk factor. In addition, we validated the prognostic value of NTRK3 in another independent cohort (GSE48276, Fig. 3h). Therefore, we selected NTRK3 for further analysis.

Fig. 3
figure 3

Screening of DEGs between high and low TMB groups associated with overall survival. a Upset plot and Venn plot of the results of three DEG-analysis methods. b The heatmap showed expression levels of 36 DEGs between high and low TMB groups. c The forest plot showed the results of univariate Cox analysis for 4 genes selected. d-g KM survival curves of 4 genes selected. (H) Validation of prognostic value for NTRK3 in GSE48276

Table 1 The 10 significant genes of univariate Cox analysis and their KM survival analysis results. (HR: Hazard Ratio; CI: Confidence Interval)

NTRK3 coexpression in BLCA

To explore the biological significance of NTRK3 in BLCA, the coexpression profiles of NTRK3 in the BLCA cohort were examined by LinkedOmics. As shown in Fig. 4C, 1,640 genes were significantly negatively correlated with NTRK3, and 5027 genes were significantly positively correlated with NTRK3 (FDR < 0.01). The top 50 significant genes positively and negatively correlated with NTRK3 are shown in heatmaps (Figs. 4a and b). GSEA enrichment results based on NTRK3 were dealt with “Weighted set cover” in LinkedOmics. The significant KEGG pathways identified included DNA replication, basal transcription factors, complement and coagulation cascades, ribosome biogenesis in eukaryotes (Fig. 4d).

Fig. 4
figure 4

NTRK3 coexpression analysis in HCC (LinkedOmics). a, b Heat maps showing top 50 genes positively and negatively correlated with NTRK3 in BLCA. Red and blue represent positive and negative correlations, respectively. c The total NTRK3 highly correlated genes identified by Pearson test in BLCA cohort. d The significant KEGG pathways of NTRK3 in BLCA cohort

Relationship between TMB scores, NTRK3 expression level and ESTIMATE results

In the low TMB group, though the P > 0.05, the expression of NTRK3 was obviously higher than that of the high TMB group (Fig. 5a). Then, we calculated immune and stromal scores for the patients by ESTIMATE algorithm. In Fig. 5b, patients were divided into two groups according to the median value of NTRK3 expression. We found that the high expression group had significantly higher immune and stromal scores than the low expression group. Next, we measured the immune and stromal scores between high and low TMB groups. Compared with the high TMB group, the low TMB group had higher immune and stromal scores (Figs. 5c and d). All these results indicated that there was an evident correlation between TMB, NTRK3 and the infiltration level of immune cells and stromal cells. There may be a potential mechanism of interaction and mutual influence among them.

Fig. 5
figure 5

Relationship between TMB scores, NTRK3 expression level and ESTIMATE results. a Comparison of NTRK3 expression level between high and low TMB groups. b Comparison of immune and stromal scores between high and low groups based on NTRK3 expression. c, d Comparison of immune and stromal scores between high and low TMB groups

Comparison of immune cell infiltration between high and low TMB groups and analysis of NTRK3 immune correlation

Immune cell fractions for each sample were estimated by the CIBERSORT algorithm based on transcriptome data. After excluding samples with P values ≥0.05, 134 low-TMB samples and 102 high-TMB samples were selected for comparing fractions of 22 leukocyte subtypes between high and low TMB groups. The results showed that high-TMB samples had higher fractions of CD 8+ T cells, memory activated CD4+ T cells and T follicular helper cells (Tfh) while low-TMB samples had higher fractions of memory resting CD4+ T cells and resting mast cells (Fig. 6a). TIMER was used to assess the correlation between NTRK3 expression and common immune infiltrating cells. The expression level of NTRK3 was positively correlated with infiltrating levels of immune cells (B cell, CD8+ T cell, CD4+ T cell, macrophage, neutrophil, and dendritic cell) (Fig. 6b). Besides, we also evaluated the relationship between NTRK3 and immunomodulators through the TISIDB portal (Fig. 7a-d). The results showed that NTRK3 in BLCA has a significant positive correlation with chemokine CCL14, immunoinhibitor ADORA2A, and immunostimulator CXCL12 (Figs. 6c, d and e). Figure 7e-g shows the relationship between the three immunomodulators and overall survival in BLCA. Thus it can be seen that NTRK3 interacts with immune regulation in TME and may become a potential biomarker which has an important impact on the development of tumors and prognosis of patients.

Fig. 6
figure 6

Immune infiltration analysis in BLCA. a Comparison of the infiltration of 22 leukocyte types between high and low TMB groups. b The correlation between NTRK3 and immune infiltrating cells. c-e The correlation between the expression of NTRK3 and the expression of CCL14, ADORA2A, and CXCL12 respectively

Fig. 7
figure 7

Relationship between NTRK3 expression and immunoregulators. a-d Spearman correlations between expression of NTRK3 and immunoregulators (Y axis) across human cancers (X axis). e-g The KM survival curves of CCL14, ADORA2A and CXCL12 in BLCA

Validation for NTRK3 by Oncomine database and immunohistochemistry

We searched the Oncomine database ( to compare the mRNA expression levels of NTRK3 between tumor and normal tissues. The result suggested that there was no significant difference between tumor and normal tissues in most studies (Supplementary Fig. S1). The immunohistochemistry data of NTRK3 was acquired from the Human Protein Atlas database ( As shown in Fig. 8, NTRK3 staining was higher in tumor samples than in normal tissue, which was consistent with the result of survival analysis, indicating that high expression of NTRK3 is a risk factor in bladder cancer.

Fig. 8
figure 8

Immunohistochemistry of NTRK3 based on the Human Protein Atlas. a IHC staining of NTRK3 in bladder normal tissue. b IHC staining of NTRK3 in bladder tumor tissue


Prior to the advent of ICIs, surgery and platinum-based systemic chemotherapy were the standards of care in BLCA. However, only 20% of MIBC patients are fit to receive chemotherapy and almost half of them may have sequelae with poor prognosis [24]. About 30–50% of metastatic MIBCs cannot receive cisplatin due to comorbidities, and until recently, platinum-based regimens were still the only treatment option, with no inspiring effect on clinical outcomes [25]. Nowadays, immunotherapy represented by ICIs has gradually revolutionized the treatment paradigm of BLCA. Indeed, the intravesical use of bacillus Calmette–Guérin (BCG) for treating NMIBCs has proven the immune response characteristic of BLCA since the 1970s [26]. Tumorigenesis is a process of mutation accumulation. BLCA has the third-highest level of somatic mutations and is highly antigenic which may facilitate the immunological recognition [27]. Consequently, it is crucial to fully understand and explore the role of TMB in BLCA immunotherapy.

In the present study, we found that BLCA patients with higher TMB had more survival benefits. This is consistent with the results of previous studies in other cancers such as cutaneous melanoma and ovarian carcinoma [28]. For patients treated with ICIs, high TMB has been confirmed to be associated with good response and better OS [29]. However, Cao et al. performed a pooled analysis of 103,078 patients to evaluate the predictive efficiency of TMB and found that the prognostic effect of TMB varied in patients with or without ICIs treatment, which might be related to tumor types [30]. Another research also indicated that the patients without ICIs treatment in some tumors may suffer a worse prognosis despite having high TMB [31]. Therefore, for the real application of TMB, it is still necessary to determine and assess its exact role in different tumors by more clinical trials. Moreover, meaningful GSEA results were all enriched in the low TMB group and most were immune-related, indicating that low TMB might enhance the tumor heterogeneity in BLCA. Of the GSEA results, ECM receptor interaction plays a crucial role in regulating invasiveness and progression of tumors. Cancer cells in the ECM of TME can release signals to mislead immune cells to avoid being attacked [32]. Besides ECM, the overactivation of focal adhesion pathway can result in dysfunction of cell migration and influence immune cell chemotaxis, which leads to tumor metastasis [33].

In terms of immune cell infiltration, the high proportion of CD4+ T cells, CD8+ T cells and Tfhs may be an important factor leading to better OS in the high TMB group. Tumor-infiltrating immune cells are a hallmark of immune surveillance and an integral part of complex microenvironment regulating tumor progression [34]. CD8+ T cells have strong cytotoxic activity for killing cancer cells, being considered as main drivers of anti-tumor immunity. The depletion in numbers and dysfunction of CD8+ T cells in TME create a favorable condition for cancer cell proliferation and metastasis [35]. CD4+ T cells and Tfhs also play a prominent role in anti-tumor immunity. CD4+ T cells can directly eliminate tumor cells through cytolysis or indirectly regulate the TME to target tumor cells [36]. And Tfh may conduce to the formation of tertiary lymphoid structures in primary site and thereby promote intertumoral immune response of CD8+ T cells and B cells [35]. ICIs can help restore T lymphocyte activity and break through the physical barrier of TME to promote T cell homing, thus activating anti-tumor immunity and improving the effect of immunotherapy [37].

After strict screening, NTRK3 was finally selected as a potential biomarker due to its excellent prognostic value and immune infiltration correlation. NTRK3 encodes the TRKC protein, a member of neurotrophic tropomyosin receptor kinase (TRK) family which regulates many aspects of neuronal development and function. After binding to the ligand neurotrophin-3 (NT-3), TRKC autophosphorylates and motivates various signaling pathways such as MAPK and PI3K/AKT pathways which can regulate cellular growth and differentiation [38]. In recent years, many studies have reported that TRK pathway aberrations such as single nucleotide variation, gene fusion and gene overexpression are involved in the pathogenesis of many cancers, among which NTRK3 gene fusion is the most fully verified carcinogenic event [39]. Except for the functional relevance of TRKC in the nervous system, the overexpression of TRKC is observed in many types of tumors, including neuroblastoma, breast cancer, hepatocellular carcinoma and metastatic melanoma. TRKC plays an important role in regulating angiogenesis, inducing tumor growth, preventing apoptosis and promoting metastasis. Abnormal activation of NTRK3 and its fusion proteins may regulate the epithelial-mesenchymal transition (EMT) process, tumor growth rate and tumorigenicity by activating several signaling pathways [40].

In our results, the expression level of NTRK3 was significantly positively correlated with the patients’ immune and stromal scores, suggesting that NTRK3 may be inherently regulated by the TME. Moreover, NTRK3 also showed good correlation with a variety of immune lymphocytes. Among the related immunomodulators, we were particularly interested in CCL14 and CXCL12, both of which were strongly positively coexpressed with NTRK3 and had prognostic value. CCL14 represents a C-C type chemokine with a high concentration in human plasma, mainly involved in the transport of lymphocytes and inflammatory cells [41]. The role of CCL14 in tumor progression is unclear and underreported, especially in bladder cancer. Yu et al. found that high expression of CCL14 played a protective role which can promote apoptosis of cancer cells and improve survival time in hepatic carcinoma [42]. However, Li et al. found that inhibiting the expression of CCL14 could effectively suppress the metastatic potential and angiogenesis of breast cancer [43]. In addition, there is an interesting theory that tumors can actively release chemokines to regulate the microenvironment and alter the host immune response from immunogenic to tolerogenic, thereby achieving immune escape and promoting tumorigenesis and metastasis. Shields et al. expounded this view in their research on CCL21 [44]. Thus, the underlying mechanism of the interaction between CCL14 and NTRK3 and how they jointly affect the progression of BLCA requires further investigation. As for CXCL12, it has been extensively studied in the field of cancers. CXCL12, also known as stromal cell derived factor-1 (SDF-1), was first characterized as a pre-B cell growth factor and combines its receptor CXCR4 as a signaling axis which mainly participated in many physiological processes such as hematopoiesis, immune responses and vascular formation [45]. Studies have shown that cancer cells in TME can induce overexpression of CXCL12 by autocrine or paracrine, and then activate downstream pathways to affect immune status and promote cancer cell proliferation and distant metastasis. Blocking the CXCL12/CXCR4 signaling axis may inhibit tumor growth and provide new ideas for immunotherapy [46]. Many biology processes in cancers are regulated by adenosine which exerts its control on both tumor and immune cells, modulating the main characteristics such as proliferation, metastasis and immune escape. A2A adenosine subtype controlled by ADORA2A is mainly present on the cell membrane of lymphocytes and able to reduce the activation of the machinery triggered by T cell receptor in immune cells, thus inducing a series of immunosuppressive events and promoting progression in tumors [47]. Cekic et al. found that A2A receptor (A2AR) regulated CD8+ T cells in TME and the application of its inhibitors could enhance the effect of immunotherapy in melanoma and bladder cancer cells [48]. A2AR antagonist showed monotherapy activity in a clinical trial in which patients were almost resistant or refractory to anti-PD (L)1 antibodies [49].

The current study has several limitations. First, our study is retrospective based on public databases and problems such as insufficient and limited data are inevitable. Second, the relationship between TMB and the response to ICIs cannot be evaluated in BLCA due to a lack of immunotherapy information of patients. Third, our results need to be verified in other independent patient cohorts with mutation information and by laboratory or clinical experiments.


In conclusion, the present study systematically explored the relationship between TMB and prognosis as well as immune infiltration in BLCA. With further exploration, our data revealed that the expression of NTRK3 was closely correlated with survival and immune response. Accordingly, NTRK3 may prove to be a novel TMB-related biomarker and contribute to the development of immunotherapy for bladder cancer.

Availability of data and materials

All the data in our study can be accessed from the online public database and all the hyperlinks to publicly archived datasets were listed in the “Methods” part.



Bladder cancer


Tumor mutation burden


The cancer genome atlas


Gene set enrichment analysis


Differentially expressed genes


Nonmuscle-invasive bladder cancer


Muscle-invasive bladder cancer


Immune checkpoint inhibitor


Tumor microenvironment


False discovery rate


Kyoto encyclopedia of genes and genomes


Single nucleotide polymorphism


Single nucleotide variant


Overall survival


Extracellular matrix


NT-3 growth factor receptor


  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  Google Scholar 

  2. Magers MJ, Lopez-Beltran A, Montironi R, Williamson SR, Kaimakliotis HZ, Cheng L. Staging of bladder cancer. Histopathology. 2019;74(1):112–34.

    Article  PubMed  Google Scholar 

  3. Lotan Y, Kamat AM, Porter MP, Robinson VL, Shore N, Jewett M, et al. Key concerns about the current state of bladder cancer: a position paper from the bladder Cancer think tank, the bladder Cancer advocacy network, and the society of urologic oncology. Cancer. 2009;115(18):4096–103.

    Article  PubMed  Google Scholar 

  4. Knowles MA, Hurst CD. Molecular biology of bladder cancer: new insights into pathogenesis and clinical diversity. Nat Rev Cancer. 2015;15(1):25–41.

    Article  CAS  PubMed  Google Scholar 

  5. Kamat AM, Hahn NM, Efstathiou JA, Lerner SP, Malmström P-U, Choi W, et al. Bladder cancer. Lancet. 2016;388(10061):2796–810.

    Article  PubMed  Google Scholar 

  6. Cao R, Yuan L, Ma B, Wang G, Tian Y. Tumour microenvironment (TME) characterization identified prognosis and immunotherapy response in muscle-invasive bladder cancer (MIBC). Cancer Immunol Immunother. 2021;70(1):1–18.

  7. Scarpato KR, Morgans AK, Moses KA. Optimal management of muscle-invasive bladder cancer - a review. Res Rep Urol. 2015;7:143–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Patel VG, Oh WK, Galsky MD. Treatment of muscle-invasive and advanced bladder cancer in 2020. CA Cancer J Clin. 2020;70(5):404-23.

  9. Buchbinder EI, McDermott DF. Cytotoxic T-lymphocyte antigen-4 blockade in melanoma. Clin Ther. 2015;37(4):755–63.

    Article  CAS  PubMed  Google Scholar 

  10. Bansal P, Osman D, Gan GN, Simon GR, Boumber Y. Recent advances in immunotherapy in metastatic NSCLC. Front Oncol. 2016;6:239.

    PubMed  PubMed Central  Google Scholar 

  11. Zibelman M, Plimack ER. Integrating immunotherapy into the Management of Renal Cell Carcinoma. J Natl Compr Cancer Netw. 2017;15(6):841–7.

    Article  CAS  Google Scholar 

  12. Rosenberg JE, Hoffman-Censits J, Powles T, van der Heijden MS, Balar AV, Necchi A, et al. Atezolizumab in patients with locally advanced and metastatic urothelial carcinoma who have progressed following treatment with platinum-based chemotherapy: a single-arm, multicentre, phase 2 trial. Lancet. 2016;387(10031):1909–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Arora S, Velichinskii R, Lesh RW, Ali U, Kubiak M, Bansal P, et al. Existing and emerging biomarkers for immune checkpoint immunotherapy in solid tumors. Adv Ther. 2019;36(10):2638–78.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Riaz N, Morris L, Havel JJ, Makarov V, Desrichard A, Chan TA. The role of neoantigens in response to immune checkpoint blockade. Int Immunol. 2016;28(8):411–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Chan TA, Yarchoan M, Jaffee E, Swanton C, Quezada SA, Stenzinger A, et al. Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic. Ann Oncol. 2019;30(1):44–56.

    Article  CAS  PubMed  Google Scholar 

  16. Killock D. TMB - a histology-agnostic predictor of the efficacy of ICIs? Nat Rev Clin Oncol. 2020;17(12):718.

  17. Mayakonda A, Lin DC, Assenov Y, Plass C, Koeffler HP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 2018;28(11):1747–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;10(21):7252–9.

    Article  CAS  PubMed  Google Scholar 

  19. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63.

    Article  CAS  PubMed  Google Scholar 

  21. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Ru B, Wong CN, Tong Y, Zhong JY, Zhong SSW, Wu WC, et al. TISIDB: an integrated repository portal for tumor-immune system interactions. Bioinformatics. 2019;35(20):4200–2.

    Article  CAS  PubMed  Google Scholar 

  24. Galsky MD, Hahn NM, Rosenberg J, Sonpavde G, Hutson T, Oh WK, et al. Treatment of patients with metastatic urothelial cancer "unfit" for Cisplatin-based chemotherapy. J Clin Oncol. 2011;29(17):2432–8.

    Article  PubMed  Google Scholar 

  25. Stenehjem DD, Tran D, Nkrumah MA, Gupta S. PD1/PDL1 inhibitors for the treatment of advanced urothelial bladder cancer. Onco Targets Ther. 2018;11:5973–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Morales A, Eidinger D, Bruce AW. Intracavitary Bacillus Calmette-guerin in the treatment of superficial bladder tumors. J Urol. 1976;116(2):180–2.

    Article  CAS  PubMed  Google Scholar 

  27. Alexandrov LB, Nik-Zainal S, Wedge DC, Aparicio SA, Behjati S, Biankin AV, et al. Signatures of mutational processes in human cancer. Nature. 2013;500(7463):415–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Klebanov N, Artomov M, Goggins WB, Daly E, Daly MJ, Tsao H. Burden of unique and low prevalence somatic mutations correlates with cancer survival. Sci Rep. 2019;9(1):4848.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Samstein RM, Lee CH, Shoushtari AN, Hellmann MD, Shen R, Janjigian YY, et al. Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat Genet. 2019;51(2):202–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Cao D, Xu H, Xu X, Guo T, Ge W. High tumor mutation burden predicts better efficacy of immunotherapy: a pooled analysis of 103078 cancer patients. Oncoimmunology. 2019;8(9):e1629258.

    Article  PubMed  PubMed Central  Google Scholar 

  31. McNamara MG, Jacobs T, Lamarca A, Hubner RA, Valle JW, Amir E. Impact of high tumor mutational burden in solid tumors and challenges for biomarker application. Cancer Treat Rev. 2020;89:102084.

    Article  CAS  PubMed  Google Scholar 

  32. Mushtaq MU, Papadas A, Pagenkopf A, Flietner E, Morrow Z, Chaudhary SG, et al. Tumor matrix remodeling and novel immunotherapies: the promise of matrix-derived immune biomarkers. J Immunother Cancer. 2018;6(1):65.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Luo M, Guan JL. Focal adhesion kinase: a prominent determinant in breast cancer initiation, progression and metastasis. Cancer Lett. 2010;289(2):127–39.

    Article  CAS  PubMed  Google Scholar 

  34. Domingues P, Gonzalez-Tablas M, Otero A, Pascual D, Miranda D, Ruiz L, et al. Tumor infiltrating immune cells in gliomas and meningiomas. Brain Behav Immun. 2016;53:1–15.

    Article  CAS  PubMed  Google Scholar 

  35. van der Leun AM, Thommen DS, Schumacher TN. CD8(+) T cell states in human cancer: insights from single-cell analysis. Nat Rev Cancer. 2020;20(4):218–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Borst J, Ahrends T, Babala N, Melief CJM, Kastenmuller W. CD4(+) T cell help in cancer immunology and immunotherapy. Nat Rev Immunol. 2018;18(10):635–47.

    Article  CAS  PubMed  Google Scholar 

  37. Zhang J, Endres S, Kobold S. Enhancing tumor T cell infiltration to enable cancer immunotherapy. Immunotherapy. 2019;11(3):201–13.

    Article  CAS  PubMed  Google Scholar 

  38. Okamura R, Boichard A, Kato S, Sicklick JK, Bazhenova L, Kurzrock R. Analysis of NTRK alterations in pan-Cancer adult and pediatric malignancies: implications for NTRK-targeted therapeutics. JCO Precis Oncol. 2018;2018.

  39. Khotskaya YB, Holla VR, Farago AF, Mills Shaw KR, Meric-Bernstam F, Hong DS. Targeting TRK family proteins in cancer. Pharmacol Ther. 2017;173:58–66.

    Article  CAS  PubMed  Google Scholar 

  40. Jin W. Roles of TrkC signaling in the regulation of tumorigenicity and metastasis of cancer. Cancers (Basel). 2020;12(1):147.

  41. Forssmann U, Mägert HJ, Adermann K, Escher SE, Forssmann WG. Hemofiltrate CC chemokines with unique biochemical properties: HCC-1/CCL14a and HCC-2/CCL15. J Leukoc Biol. 2001;70(3):357–66.

    CAS  PubMed  Google Scholar 

  42. Gu Y, Li X, Bi Y, Zheng Y, Wang J, Li X, et al. CCL14 is a prognostic biomarker and correlates with immune infiltrates in hepatocellular carcinoma. Aging. 2020;12(1):784–807.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Li Q, Shi L, Gui B, Yu W, Wang J, Zhang D, et al. Binding of the JmjC demethylase JARID1B to LSD1/NuRD suppresses angiogenesis and metastasis in breast cancer cells by repressing chemokine CCL14. Cancer Res. 2011;71(21):6899–908.

    Article  CAS  PubMed  Google Scholar 

  44. Shields JD, Kourtis IC, Tomei AA, Roberts JM, Swartz MA. Induction of lymphoidlike stroma and immune escape by tumors that express the chemokine CCL21. Science. 2010;328(5979):749–52.

    Article  CAS  PubMed  Google Scholar 

  45. Nazari A, Khorramdelazad H, Hassanshahi G. Biological/pathological functions of the CXCL12/CXCR4/CXCR7 axes in the pathogenesis of bladder cancer. Int J Clin Oncol. 2017;22(6):991–1000.

    Article  CAS  PubMed  Google Scholar 

  46. Meng W, Xue S, Chen Y. The role of CXCL12 in tumor microenvironment. Gene. 2018;641:105–10.

    Article  CAS  PubMed  Google Scholar 

  47. Merighi S, Battistello E, Giacomelli L, Varani K, Vincenzi F, Borea PA, et al. Targeting A3 and A2A adenosine receptors in the fight against cancer. Expert Opin Ther Targets. 2019;23(8):669–78.

    Article  CAS  PubMed  Google Scholar 

  48. Cekic C, Linden J. Adenosine A2A receptors intrinsically regulate CD8+ T cells in the tumor microenvironment. Cancer Res. 2014;74(24):7239–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Fong L, Hotson A, Powderly JD, Sznol M, Heist RS, Choueiri TK, et al. Adenosine 2A receptor blockade as an immunotherapy for treatment-refractory renal cell Cancer. Cancer Discov. 2020;10(1):40–53.

    Article  CAS  PubMed  Google Scholar 

Download references


Not applicable.


This work was supported by the National Natural Science Foundation of China (grant no. 82071750, 81772713, 81472411), Taishan Scholar Program of Shandong Province (grant no. tsqn20161077), Major Science and technology innovation project of Shandong Province (grant no. 2019JZZY021002), Key projects of Qingdao Science and Technology Program (grant no. 18–6–1-64-nsh), Research and Development Program of Shandong Province (grant no. 2018GSF118197).

Author information

Authors and Affiliations



ZZ and HTN contributed to conception and design. YBY and PFZ contributed to collection and assembly of data. ZZ, GFM and MXZ contributed to data analysis and interpretation. YL and WJ contributed to manuscript writing. HTN contributed to the manuscript revision and finalization. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Wei Jiao or Haitao Niu.

Ethics declarations

Ethics approval and consent to participate

Any repository data used in this study are open access and do not require any permissions. Ethics approval and consent to participate are not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Fig. S1

. Comparison of NTRK3 mRNA expression between tumor and normal tissues across multi-studies by Oncomine.

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 The Creative Commons Public Domain Dedication waiver ( 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

Verify currency and authenticity via CrossMark

Cite this article

Zhang, Z., Yu, Y., Zhang, P. et al. Identification of NTRK3 as a potential prognostic biomarker associated with tumor mutation burden and immune infiltration in bladder cancer. BMC Cancer 21, 458 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Tumor mutation burden
  • Immune infiltration
  • NTRK3
  • Prognosis
  • TCGA