Skip to main content

Prognostic factors of non-muscle invasive bladder cancer: a study based on next-generation sequencing



To investigate the genetic prognostic factors for the recurrence of non-muscle invasive bladder cancer.

Materials and methods

The patients underwent transurethral resection of bladder tumor and received bacillus Calmette–Guérin (BCG) or epirubicin. Next-generation sequencing was performed and alterations of genes, pathways, and tumor mutation burden were recorded. Associations between these clinicopathological and genetic variants were estimated, and prognostic factor identified.


A total of 58 cases were included in our study, and 46 patients underwent treatment with BCG. FGFR3 was the most frequently altered gene (48%), and more commonly detected in intermediate-risk patients. Univariate Cox analysis demonstrated that 10 genes were significantly correlated with BCG failure, while NEB, FGFR1 and SDHC were independent recurrence predictors. Besides, epigenetic-related gene pathway mutations were negatively correlated with recurrence (hazard ratio: 0.198, P = 0.023). DNA damage response and repair gene alterations were positively correlated with tumor burden, while altered TP53 was most frequent among these genes and significant correlated with high tumor burden.


BCG instillation significantly reduced the rate of recurrence compared with epirubicin in this population. Potential biomarkers and therapeutic targets were found with the help of next-generation sequencing; correlations between DDR genes alterations and high tumor mutation burden were also demonstrated.


Approximately 75% of patients with bladder cancer have non-muscle invasive bladder cancer (NMIBC; Ta or T1 disease) at the time of primary diagnosis [1]. In one study, more than 70% of patients suffered recurrences in the course of treatment and surveillance, and disease progression occurred in 33% of them [2]. For these reasons, many efforts have focused on exploring prognostic factors and optimized treatments for NMIBC. Prior treatment, prior recurrence, tumor numbers, tumor size, T stage, nuclear grade and presence of carcinoma in situ (CIS) were shown to be the most important prognostic factors for NMIBC [3, 4]. The European Association of Urology (EAU) stratified these patients into three risk groups: low, intermediate, and high-risk, according to these results [5].

For intermediate- and high-risk patients, TUR-BT and subsequent intravesical drug instillation are the standard treatments. It has been confirmed that treatment with Bacillus Calmette–Guérin (BCG) following TUR-BT is superior to TUR-BT alone or combined with chemotherapy in Caucasian patients [6, 7]. However, despite the administration of maintenance therapy with BCG, 32.6% and 13.4% of patients experience recurrence and tumor progression, respectively [4].

Many attempts on bladder tumor markers have been made and markers were developed to detect tumor associated antigens, growth factors, extracellular matrix proteins [8]. Even some of them were supposed to replace urinary cytology or cystoscopy, however, no prospective data supported their impact on prognosis, especially for prediction of BCG failure [8, 9]. The revolutionary technology called next generation sequencing (NGS) enables fast and cost-effective generation of genetic sequence data with high accuracy [10]. NGS can acquire substantial information on genetic alterations, such that potential tumor markers, such as biomarkers related to failure of intravesical therapy, are more likely to be identified. Pietzak et al. used targeted NGS with cancer-associated gene panel and demonstrated that NMIBC tumors had at least one potentially actionable mutation, and DNA damage repair (DDR) gene alterations were associated with increased mutational load [11]. In an NGS study of urothelial cancer, DDR alterations were demonstrated to be independently associated with effect of immune therapy [12]. In this study, we analyzed the prognostic factors of recurrence based on a randomized controlled trial of BCG and epirubicin (EPI), examined the predictive value of prognostic models, and investigated the genetic alterations of this cohort to define the potential clinical implications of these variants.

Patients and methods

Patients and treatment

Patients were collected from a prospective randomized study of intravesical therapy after TUR-BT: those diagnosed as intermediate- or high-risk NMIBC were randomized to receive BCG or EPI for one-year regime by a ratio of 4:1 (phase IV clinical trial: CTR20150840 for Center for Drug Evaluation of National Medical Products Administration, China). The inclusion criteria were: age 20–75 years; TUR-BT with completely tumor resection. BT was completely resected with pathologically proven intermediate or high-risk non-muscle invasive bladder urothelium carcinoma according to the EAU guidelines for risk stratification of NMIBC (low-risk tumors were defined as primary, solitary, TaG1 (papillary urothelial neoplasm of low malignant potential and low-grade carcinomas), < 3 cm, no CIS; intermediate risk tumors were those between the category of low and high risk; characteristics of high-risk tumors were any of the following: T1, G3 tumor, CIS, multiple, recurrent and large (> 3 cm) tumors), and adjacent detrusor muscle was examined separately or together with BT [5]. Exclusion criteria were as follows: (a) Eastern Cooperative Oncology Group performance status score > 1; (b) active tuberculosis or treatment for anti-tuberculous; (c) immune deficiency or immunosuppressive therapy; (d) severe complication (e.g., serious cardiovascular and cerebrovascular disease), or presence of other types of cancer; e. previous diagnosis with muscle-invasive bladder cancer; (e) history of treatment (e.g., chemotherapy, radiotherapy or immunotherapy) during the previous 4 weeks that may influence the research results; (f) serious intraoperative and postoperative complications (e.g., bladder perforation, serious postoperative hematuria, bladder irritation, etc.); and (g) contraindications to treatment or inability to participate in our trial due to pregnancy, severe disability, or serious psychological problems.

Follow-up and outcomes

For intermediate-risk patients, cystoscopy was performed 3 and 6 months after surgery, then biannually for 5 years. For high-risk patients, cystoscopy was performed trimonthly until 2 years after surgery, then biannually for 5 years. The initiation time was the day TUR-BT was performed. The primary end point was recurrence during follow-up, while the secondary end point was progression to muscle-invasive bladder cancer or therapy-related discontinuation of treatment. Recurrence-free survival (RFS) and progression-free survival (PFS) times were recorded. BCG failure and intolerance were recorded according to the definitions of the EAU [5].

Data collection

Clinicopathological variables included gender, age, tumor size, site, number, tumor histopathological type, T stage (8th American Joint Committee on Cancer TNM classification system) [13], tumor grade (2004 World Health Organization grading system), and CIS status. All pathological reports were provided by the Pathology Department of West China Hospital, Sichuan University.

Paraffin-embedded tissues obtained from the BCG treatment group were collected for NGS. Patients provided written informed consent. Sequencing was performed on a Nextseq 500 sequencer (Illumina, Inc., San Diego, CA, USA) for 520 cancer-related genes, including the whole exons of 312 genes, as well as critical exons, introns, and promoter regions of the remaining 208 genes. Resultant sequences were analyzed for gene alternations. Nonsynonymous mutations, deletions, nonsense mutations, transcoding mutations, and splicing mutations were considered significant. In addition, the tumor mutation burden (TMB) was also recorded. Genomic pathways were analyzed in qualitative and quantitative dimensions (i.e., whether mutations are present and the number of mutations genes in each pathway).

Statistical methods

One-way analysis of variance and the Chi-squared test were used for comparisons between groups. The Chi-squared test was used to analyze categorical variants. The Pearson Chi-squared test was used when n ≥ 40 and T ≥ 5; the continuously correcting Chi-squared test was used when n ≥ 40 and 1 ≤ T<5; if n < 40 or T < 0, Fisher’s exact test was used. Kruskal–Wallis and Wilcoxon tests were used to analyze continuous variants. RFS and PFS were estimated using the Kaplan–Meier method. Differences in survival between two or more subgroups were evaluated using log-rank tests. Univariate and multivariate Cox regression analyses were performed to determine the clinicopathological parameters associated with the recurrence of patients with NMIBC. Patients were stratified using the European Organisation for Research and Treatment of Cancer (EORTC) risk tables [3] and the Spanish Urological Club for Oncological Treatment (CUETO) scoring model [4].

Statistical analyses were performed using the SPSS Statistics version 25 (IBM, Armonk, NY, USA). P-value < 0.05 denoted statistical significance.


Patient baseline data and prognosis

The clinicopathological data from 58 patients are shown in Table 1. The mean age was 61.91 years (standard deviation, SD: 7.67 years), and the median follow-up time was 32.86 months (interquartile range, IQR: 21.00–44.49 months). There were no significant differences between BCG and EPI subgroups in terms of baseline clinicopathological data except for T stage. For the BCG subgroup, recurrence occurred in nine patients, and six of these recurred during course of intravesical BCG therapy. Two-year RFS rate for the BCG group was 82.2%; while for EPI subgroup, seven patients experienced recurrences, five of whom did so during the EPI treatment course. The two-year RFS rate for EPI was 58.3%. Univariate analysis revealed significantly better RFS in the BCG-treated subgroup (P = 0.020) and in those with incipient tumor (P = 0.001 for all 58 patients and P = 0.014 for 46 BCG-treated patients) (Additional file 1: Table S1).

Table 1 Clinicopathological data of patients with NMIBC

Genomic landscape of patients with NMIBC and relationship with prognosis

NGS was performed on the total 58 patients, the most frequently altered genes were FGFR3 (48%), KDM6A (47%), KMT2D (43%), KMT2C (34%), and STAG2 (31%). Oncoprints for these genes are presented in Additional file 2, Figures S1–S7). Significant mutations of FGFR3 were observed in 73.1% of intermediate-risk patients (19/26) and 15.6% of high-risk patients (5/32). For STAG2, these numbers were 46.2% (12/26) and 18.8% (6/32), respectively. Mutations of FGFR3, STAG2, and PRKDC were significantly associated with tumor risk stage (P < 0.05, Additional file 1: Table S2). Mutations of 17 genes (such as KDM6A and ARID1A) significant correlated with higher EORTC scores, while mutations of other 17 genes (such as TP53) significantly correlated with higher CUETO scores (Additional file 1: Table S3).

Cox univariate analysis was performed to identify correlations of genomic mutations with recurrence. In all 58 cases, 18 genes were significantly correlated with recurrence (Table 2 and Additional file 1: Table S4); all these mutations were positively associated with poor outcome. However, among these genes, only mutations of NEB, MLH1, GATA3, FGFR1 and RAF1 occurred in > 5% of patients. Similar results were observed in the BCG subgroup (46 cases), with ten genes significantly associated with recurrence (P < 0.05, Table 2 and Additional file 1: Table S5). The multivariate Cox analysis included all genes significantly associated with recurrence, the results demonstrated that NEB, MLH1, FGF12, and FGFR1 independently predicted poor outcome in all 58 patients (P = 0.001, 0.001, 0.001, and 0.007, respectively), while NEB, FGFR1 and SDHC were independent prognostic predictor for recurrence in patients treated with BCG (P = 0.001, 0.004 and 0.017 respectively) (Additional file 1: Table S6).

Table 2 Univariate analysis of genetic mutations and recurrence survival rate

Genomic pathways and recurrence

Canonical genomic pathways were analyzed to define whether mutations were associated with prognosis (Fig. 1). Mutations of the receptor tyrosine kinase/RAS/phosphatidylinositol 3-kinase (RTK/RAS/PI3K) pathway were found in 87.9% of patients, and were associated with higher risk stage (P = 0.013, Additional file 1: Table S7). Epigenetic-related gene mutations were also frequently detected in patients with NMIBC (49/58, 84.4%), with markedly lower rates of mutations in the TP53/cell cycle pathway, the switch/sucrose nonfermentable pathway, the DNA damage pathway, and the alternative splicing pathway (50.0%, 31.0%, 32.8%, and 34.5%, respectively). Correlations between the number of pathway mutations and prediction models were also analyzed: larger numbers of mutated genes in the epigenetic pathway, histone modification pathway, and the switch/sucrose nonfermentable pathway significantly correlated with a poorer EORTC score, whereas there was no significant correlation with the CUETO score (Additional file 1: Table S7).

Fig. 1

Oncoprint of genetic pathways. BCG: Bacillus Calmette–Guérin treated group; EPI: epirubicin treated group. CUETO model: Spanish Urological Club for Oncological Treatment scoring model; EORTC: European Organization for Research and Treatment of Cancer; CN: copy number; RTK: receptor tyrosine kinases; PI3K: phosphatidylinositol 3-kinase; SWI/SNF: Switch/sucrose nonfermentable

Correlations of pathway mutations with recurrence were analyzed using Cox regression in qualitative and quantitative dimensions (Table 3). In qualitative analysis, mutations in epigenetic-related genes, the RTK-PI3K pathway, the SWI/SNF pathway and the alternative splicing pathway were protective factors (hazard ratio [HR] < 1), while TP53/cell cycle pathway mutations were risk factor, in all 58 patients and in the BCG subgroup. Mutations in epigenetic-related genes were significantly associated with recurrence in qualitative analysis (P = 0.020 and 0.023 for total patients and for the BCG subgroup, respectively). Nevertheless, when quantitative analysis was performed, we found that no pathway mutations significantly correlated with recurrence (P > 0.05, Table 3).

Table 3 Univariate analysis of mutated genetic pathways, DDR genes, and recurrence survival rate

DDR genes and TMB

DDR genes have been associated with increased mutation load and with the effect of immunotherapy against urothelial carcinoma [12]. We also analyzed the correlation between DDR, TMB, and prognosis (the list of DDR genes is provided in Additional file 1: Table S8 and the alteration status is shown in Additional file 2: Figure S8). The percentages of DDR mutations in 58 cases with NMIBC and BCG subgroups were 86.2% and 84.8%, respectively (numbers of mutations were 165 and 141, respectively). TP53 was the most frequently mutated gene among the DDR genes (19/165 in the total cases and 17/141 in the BCG subgroups). Of note, missense variant was the most common alteration in TP53 (16/19 and 14/17 in the total and BCG group, respectively) (Additional file 3: Figure S8). In our previous analysis, TP53 mutation correlated with poorer CUETO score; however, there were no correlations with tumor risk group stratification, EORTC score, or RFS (Additional file 1: Tables S2–5). When considering mutations in all DDR genes, there were no correlations demonstrated with the predicted models or RFS (P > 0.05) (Table 3 and Additional file 1: Table S7).

In the analysis of TMB, significant correlations were found with the EORTC model (P-values for correlation with risk group stratification, CUETO, and EORTC scores were 0.485, 0.706, and 0.037, respectively). There was no significant correlation between TMB and RFS (HR: 0.972, P = 0.316 for 58 NMIBC cases and HR: 0.976, P = 0.497 for 46 BCG subgroup cases). However, significant correlations were demonstrated between alterations in DDR and TMB, and between DDR mutated numbers and TMB (Figs. 2a, b). In addition, TP53 mutations were also significantly associated with TMB (Fig. 2c).

Fig. 2

Correlation analysis between DNA damage response and repair genes and TMB. a Correlation of the TMB and DDR alteration status. b Correlation of the TMB and DDR mutational numbers. DDR: DNA damage response and repair. c Correlation of the TMB and TP53 alteration status. TMB: tumor mutational burden


Recurrence of NMIBC is affected by numerous factors, including tumor size, number of tumors, prior recurrence rate, T stage, grade, and presence of CIS [3]. Since first reported, the efficacy of BCG irrigation has been widely demonstrated and became the standard therapy for intermediate- and high-risk NMIBC [14]. However, recurrence continues to occur in many patients treated with BCG; 32.6% as reported by Jesus et al. [4]. Even after a full course of maintenance therapy with BCG, 15.7% and 26.3% of patients continue to experience early and late recurrence, respectively [15].

Pretreatment biomarkers can help to predict the therapeutic effect of BCG against NMIBC [2]. In fact, these biomarkers can be used as predictors and carry therapeutic value. For example, the protein encoded by the MTOR gene belongs to the family of phosphatidylinositol kinases, which mediate DNA damage and nutrient deficiency [16]. Hansel et al. reported that abnormal activation of phosphorylated MTOR and its downstream target was associated with increased pathological stage and poorer survival [17]. In our study, we also found a correlation between MTOR mutations and tumor recurrence. Therefore, we believe that MTOR can serve as a potential biomarker for recurrence of NMIBC. In addition, it was reported that inhibition of MTOR expression induced apoptosis of bladder cancer cells and inhibited cell cycle progression, cell growth, angiogenesis, and endothelial cell proliferation, indicating its therapeutic value [18]. In fact, MTOR protein is a key downstream protein kinase in the PI3K/AKT signaling pathway, with mutations occurring in approximately 40% of bladder cancer cases [19]. However, there was no correlation found between mutations of the RTK-PI3K pathway and recurrence in our study, and the underlying mechanism warrants further research.

Epigenetic regulations, which include DNA methylation, histone modification and chromatin remodeling, control gene expression and tumor biology [20]. Alterations of histone modification occurred in nearly 89% of invasive bladder cancer cases [21], while DNA hypermethylation correlated with grade and stage [22, 23]. In vitro, tumor suppressor genes can be restored after treatment with DNA methyltransferases and histone deacetylase [24, 25], suggesting the therapeutic significance of the epigenetic pathway in bladder cancer. Moreover, chromatin remodeling genes, such as ARID1A and KDM6A were frequently mutated in bladder cancer [21, 26]. Pietzak et al., found that mutated ARID1A correlated with recurrence and could be used as a biomarker [27]. However, in our study, a negative correlation was found between mutations of epigenetic-related genes and recurrence, and there was no prognostic significance demonstrated for ARID1A. In fact, epigenetic dysregulations occur in different stages of various tumors; this can be either a cause or an effect of tumor recurrence and progression. Further studies of the underlying mechanism are warranted to understand the significance of epigenetic regulation in NMIBC.

In our study, mutated FGFR1 and NEB were shown to be independent predictors of recurrence, despite their low frequency of alterations observed in the present cohort (5.2%). Rehab et al. found that overexpression of FGFR1 was common in malignant bladder cancer, and the sensitivity and specificity of combined FGFR1 and CK20 for the differentiation between MIBC and NMIBC can be as high as 97.5% and 92.5%, respectively [28]. These results suggest that FGFR1 can serve as a biomarker of tumor invasiveness. Although we believe this is also the case for NEB, further investigations are warranted. However, except for the aforementioned independent predictors, we did not consider all the statistically significant gene mutations with low mutation rates as clinically significant.

The clinical importance of some genes could not be ignored even in the absence of prognostic significance. The telomerase reverse transcriptase gene (TERT) is involved in DNA protection and cellular ageing processes [29], and mutations of this gene occur in > 70% of bladder cancer cases [27, 30]. Although there was no correlation found between TERT mutation and tumor stage, grade, and prognosis in our and previous studies, this gene was also considered a potential biomarker for screening and surveillance owing to its high mutation rate in patients with bladder cancer [30, 31]. FGFR3 is also clinically significant, whose mutations have been identified as a key to tumor initiation [32]. FGFR3 can serve as a potential biomarker, based on the negative correlation of mutated FGFR3 and tumor malignant degree. FGFR3 mutations were more commonly found in NMIBC than MIBC in a previous study [11]. In our study, these mutations were also more common in the intermediate-risk than the high-risk NMIBC group. FGFR3 is also a potential treatment target as it can activate many downstream kinase pathways, including the RTK/RAS pathway [32]. The antitumor activity of BGJ398, a selective inhibitor of FGFR, has been demonstrated in patients with FGFR3 mutations [33]. There is an ongoing clinical trial of BGJ398 in patients with NMIBC, and results are expected in the future.

DDR-related genes were frequently altered in high-grade NMIBC [27], and these alterations were associated with increased mutation burden and better efficacy of chemotherapy and immunotherapy [12, 34]. Min et al. studied DDR alterations in 60 patients with advanced urothelial cancer treated with adjuvant PD-1/PD-L1 blockade therapy, they found that a higher response rate to immune therapy in the DDR-altered group (67.9% vs. 18.8%), as well as longer PFS and overall survival [12]. Therefore, we believe that DDR alterations may influence the efficacy of BCG instillation, also based on the activation of the autoimmune system. Nevertheless, we found that DDR alterations were only significantly correlated with TMB and there was no correlation with risk stage or prognosis.

TP53 was the most frequently mutated DDR gene in our study. DNA damage can trigger the activation of p53, which maintains genetic stability and controls the transcription of DNA repair-related genes [35]. A recent meta-analysis suggested that p53 overexpression in NMIBC might be associated with poor RFS in an Asian population [36]. However, there was no prognostic significance of TP53 mutation or the mutated TP53/cell cycle pathway observed in our study. Interestingly, mutated TP53 gene was associated with larger TMB, while these two variants correlated with higher CUETO and EORTC scores. This evidence suggests that prognostic models may serve as bridges for understanding the correlation between DDR alterations and NMIBC recurrence.

There were some limitations in our study. Although it was based on a randomized controlled trial with limited heterogeneity among the subjects, the sample size was relatively small for the analysis of the association between prognosis and various clinicopathological and genomic variants. Furthermore, the follow-up period was limited, with a small number of patients reaching the follow-up endpoint; therefore, PFS and overall survival could not be calculated.


Maintenance treatment with BCG provided better control of recurrence in Chinese patients with NMIBC after TUR-BT versus intravesical therapy with EPI. NGS of these patients showed that FGFR1 and NEB are biomarkers of recurrence, while mutations in the MTOR gene and epigenetic pathway are potential treatment targets for patients who fail treatment with BCG. DDR alterations were positively correlated with higher TMB, TP53 may be the key to understanding the relationship between DDR and recurrence after treatment with BCG (Additional file 4).

Availability of data and materials

The dataset supporting the conclusions of this article is included within the article and its supplementary materials.


  1. 1.

    Burger M, Catto JW, Dalbagni G, et al. Epidemiology and risk factors of urothelial bladder cancer. Eur Urol. 2013;63(2):234–41.

    PubMed  Article  Google Scholar 

  2. 2.

    Chamie K, Litwin MS, Bassett JC, et al. Recurrence of high-risk bladder cancer: a population-based analysis. Cancer. 2013;119(17):3219–27.

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Sylvester RJ, van der Meijden AP, Oosterlinck W, et al. Predicting recurrence and progression in individual patients with stage Ta T1 bladder cancer using EORTC risk tables: a combined analysis of 2596 patients from seven EORTC trials. Eur Urol. 2006;49(3):466 (discussion 475–467).

    PubMed  Article  Google Scholar 

  4. 4.

    Fernandez-Gomez J, Madero R, Solsona E, et al. Predicting nonmuscle invasive bladder cancer recurrence and progression in patients treated with bacillus Calmette-Guerin: the CUETO scoring model. J Urol. 2009;182(5):2195–203.

    PubMed  Article  Google Scholar 

  5. 5.

    Babjuk M, Burger M, Comperat EM, et al. European Association of Urology Guidelines on Non-muscle-invasive Bladder Cancer (TaT1 and Carcinoma In Situ)-2019 Update. Eur Urol. 2019;76(5):639–57.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Sylvester RJ, Brausi MA, Kirkels WJ, et al. Long-term efficacy results of EORTC genito-urinary group randomized phase 3 study 30911 comparing intravesical instillations of epirubicin, bacillus Calmette-Guerin, and bacillus Calmette-Guerin plus isoniazid in patients with intermediate- and high-risk stage Ta T1 urothelial carcinoma of the bladder. Eur Urol. 2010;57(5):766–73.

    PubMed  Article  Google Scholar 

  7. 7.

    Pan J, Liu M, Zhou X. Can intravesical bacillus Calmette-Guerin reduce recurrence in patients with non-muscle invasive bladder cancer? An update and cumulative meta-analysis. Front Med. 2014;8(2):241–9.

    PubMed  Article  Google Scholar 

  8. 8.

    Lokeshwar VB, Habuchi T, Grossman HB, et al. Bladder tumor markers beyond cytology: International Consensus Panel on bladder tumor markers. Urology. 2005;66(6 Suppl 1):35–63.

    PubMed  Article  Google Scholar 

  9. 9.

    Liou LS. Urothelial cancer biomarkers for detection and surveillance. Urology. 2006;67(3 Suppl 1):25–33 (discussion 33–24).

    PubMed  Article  Google Scholar 

  10. 10.

    Xuan J, Yu Y, Qing T, Guo L, Shi L. Next-generation sequencing in the clinic: promises and challenges. Cancer Lett. 2013;340(2):284–95.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Audenet F, Attalla K, Sfakianos JP. The evolution of bladder cancer genomics: what have we learned and how can we use it? Urol Oncol. 2018;36(7):313–20.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Teo MY, Seier K, Ostrovnaya I, et al. Alterations in DNA damage response and repair genes as potential marker of clinical benefit from PD-1/PD-L1 blockade in advanced urothelial cancers. J Clin Oncol. 2018;36(17):1685–94.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Deh BHB, Efstathiou JA, KonetyB, Lee CT, Jam es M. M cKiernan, Plimack ER, Reuter VE, Srídhar S, Vikram R, and Stadler WM. AJCC Cancer Staging Manual 8th Edition-Urinary Bladder. 2017.

  14. 14.

    Chang SS, Boorjian SA, Chou R, et al. Diagnosis and treatment of non-muscle invasive bladder cancer: AUA/SUO Guideline. J Urol. 2016;196(4):1021–9.

    PubMed  Article  Google Scholar 

  15. 15.

    Cambier S, Sylvester RJ, Collette L, et al. EORTC nomograms and risk groups for predicting recurrence, progression, and disease-specific and overall survival in non-muscle-invasive Stage Ta-T1 urothelial bladder cancer patients treated with 1–3 years of maintenance Bacillus Calmette-Guerin. Eur Urol. 2016;69(1):60–9.

    PubMed  Article  Google Scholar 

  16. 16.

    Sathe A, Nawroth R. Targeting the PI3K/AKT/mTOR pathway in bladder cancer. Methods Mol Biol. 2018;1655:335–50.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Hansel DE, Platt E, Orloff M, et al. Mammalian target of rapamycin (mTOR) regulates cellular proliferation and tumor growth in urothelial carcinoma. Am J Pathol. 2010;176(6):3062–72.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Pinto-Leite R, Arantes-Rodrigues R, Sousa N, Oliveira PA, Santos L. mTOR inhibitors in urinary bladder cancer. Tumour Biol. 2016;37(9):11541–51.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Carneiro BA, Meeks JJ, Kuzel TM, Scaranti M, Abdulkadir SA, Giles FJ. Emerging therapeutic targets in bladder cancer. Cancer Treat Rev. 2015;41(2):170–8.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Esteller M. Epigenetics in cancer. N Engl J Med. 2008;358(11):1148–59.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Robertson AG, Kim J, Al-Ahmadie H, et al. Comprehensive molecular characterization of muscle-invasive bladder cancer. Cell. 2017;171(3):540–56.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Sacristan R, Gonzalez C, Fernandez-Gomez JM, Fresno F, Escaf S, Sanchez-Carbayo M. Molecular classification of non-muscle-invasive bladder cancer (pTa low-grade, pT1 low-grade, and pT1 high-grade subgroups) using methylation of tumor-suppressor genes. J Mol Diagn. 2014;16(5):564–72.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Olkhov-Mitsel E, Savio AJ, Kron KJ, et al. Epigenome-wide DNA methylation profiling identifies differential methylation biomarkers in high-grade bladder cancer. Transl Oncol. 2017;10(2):168–77.

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Faleiro I, Leao R, Binnie A, de Mello RA, Maia AT, Castelo-Branco P. Epigenetic therapy in urologic cancers: an update on clinical trials. Oncotarget. 2017;8(7):12484–500.

    PubMed  Article  Google Scholar 

  25. 25.

    Shang D, Liu Y, Matsui Y, et al. Demethylating agent 5-aza-2′-deoxycytidine enhances susceptibility of bladder transitional cell carcinoma to Cisplatin. Urology. 2008;71(6):1220–5.

    PubMed  Article  Google Scholar 

  26. 26.

    Gui Y, Guo G, Huang Y, et al. Frequent mutations of chromatin remodeling genes in transitional cell carcinoma of the bladder. Nat Genet. 2011;43(9):875–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Pietzak EJ, Bagrodia A, Cha EK, et al. Next-generation sequencing of nonmuscle invasive bladder cancer reveals potential biomarkers and rational therapeutic targets. Eur Urol. 2017;72(6):952–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Abdul-Maksoud RS, Shalaby SM, Elsayed WS, Elkady S. Fibroblast growth factor receptor 1 and cytokeratin 20 expressions and their relation to prognostic variables in bladder cancer. Gene. 2016;591(2):320–6.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Zhang X, Zhang Y. Bladder cancer and genetic mutations. Cell Biochem Biophys. 2015;73(1):65–9.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Allory Y, Beukers W, Sagrera A, et al. Telomerase reverse transcriptase promoter mutations in bladder cancer: high frequency across stages, detection in urine, and lack of association with outcome. Eur Urol. 2014;65(2):360–6.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Hurst CD, Platt FM, Knowles MA. Comprehensive mutation analysis of the TERT promoter in bladder cancer and detection of mutations in voided urine. Eur Urol. 2014;65(2):367–9.

    CAS  PubMed  Article  Google Scholar 

  32. 32.

    Al Hussain TO, Akhtar M. Molecular basis of urinary bladder cancer. Adv Anat Pathol. 2013;20(1):53–60.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Nogova L, Sequist LV, Perez Garcia JM, et al. Evaluation of BGJ398, a fibroblast growth factor receptor 1–3 kinase inhibitor, in patients with advanced solid tumors harboring genetic alterations in fibroblast growth factor receptors: results of a global phase I, dose-escalation and dose-expansion study. J Clin Oncol. 2017;35(2):157–65.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Plimack ER, Dunbrack RL, Brennan TA, et al. Defects in DNA repair genes predict response to neoadjuvant cisplatin-based chemotherapy in muscle-invasive bladder cancer. Eur Urol. 2015;68(6):959–67.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Ciccarese C, Massari F, Blanca A, et al. Tp53 and its potential therapeutic role as a target in bladder cancer. Expert Opin Ther Targets. 2017;21(4):401–14.

    CAS  PubMed  Article  Google Scholar 

  36. 36.

    Zhou X, Zhang G, Tian Y. p53 Status correlates with the risk of recurrence in non-muscle invasive bladder cancers treated with Bacillus Calmette-Guerin: a meta-analysis. PLoS ONE. 2015;10(3):e0119476.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

Download references


We would like to thanks Dr Ranran Liu and Dr Xuelin Zhang for their help of data collecting and patients’ follow-up; we also thanks Dr Tian Qin and Dr Hao Liu of Burning Rock Biotech (Guangzhou, China) for their helping on genetic sequencing work.


This study was supported by the National Natural Science Foundation of China (Reference Number: 81672552) and the 1.3.5 Project for Disciplines of Excellence, West China Hospital, Sichuan University (Reference number: ZY2016104).

Author information




XL designed the study, guided experiments, carried out analysis and reviewed the manuscript. YS, XH and ZY contributed to the study design, collected the data, performed analysis, and wrote the manuscript. TL, WY, KW and SR contributed to the data collection and analysis and reviewed the manuscript. SX, WD, SF, YW, YL and KW helped with data collation and statistical analysis and reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiang Li.

Ethics declarations

Ethics approval and consent to participate

This study was approved by Center for Drug Evaluation of National Medical Products Administration (CDE-NMPA, China) and Ethics Committee of our institute, with the whole process supervised. The registration number in CDE-NMPA was CTR20150840 (website of CDE: All the research registration and ethical approval documents were added below. Patients and their authorized family members had been fully informed before follow-up work was performed, with informed consent signed.

Consent for publication

Consent for publication was obtained from all participants.

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.

Supplementary information

Additional file 1:

Tables S1 to S8.

Additional file 2:

Figures S1 to S7.

Additional file 3:

Figure S8.

Additional file 4:

Original dataset.

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

Shao, Y., Hu, X., Yang, Z. et al. Prognostic factors of non-muscle invasive bladder cancer: a study based on next-generation sequencing. Cancer Cell Int 21, 23 (2021).

Download citation


  • Bladder cancer
  • Bacillus Calmette–Guérin
  • Predictive model
  • Next-generation sequencing
  • DNA damage response and repair