Skip to main content

Optimizing identification of consensus molecular subtypes in muscle-invasive bladder cancer: a comparison of two sequencing methods and gene sets using FFPE specimens

Abstract

Background

Molecular subtypes predict prognosis in muscle-invasive bladder cancer (MIBC) and are explored as predictive markers. To provide a common base for molecular subtyping and facilitate clinical applications, a consensus classification has been developed. However, methods to determine consensus molecular subtypes require validation, particularly when FFPE specimens are used. Here, we aimed to evaluate two gene expression analysis methods on FFPE samples and to compare reduced gene sets to classify tumors into molecular subtypes.

Methods

RNA was isolated from FFPE blocks of 15 MIBC patients. Massive analysis of 3’ cDNA ends (MACE) and the HTG transcriptome panel (HTP) were used to retrieve gene expression. We used normalized, log2-transformed data to call consensus and TCGA subtypes with the consensusMIBC package for R using all available genes, a 68-gene panel (ESSEN1), and a 48-gene panel (ESSEN2).

Results

Fifteen MACE-samples and 14 HTP-samples were available for molecular subtyping. The 14 samples were classified as Ba/Sq in 7 (50%), LumP in 2 (14.3%), LumU in 1 (7.1%), LumNS in 1 (7.1%), stroma-rich in 2 (14.3%) and NE-like in 1 (7.1%) case based on MACE- or HTP-derived transcriptome data. Consensus subtypes were concordant in 71% (10/14) of cases when comparing MACE with HTP data. Four cases with aberrant subtypes had a stroma-rich molecular subtype with either method. The overlap of the molecular consensus subtypes with the reduced ESSEN1 and ESSEN2 panels were 86% and 100%, respectively, with HTP data and 86% with MACE data.

Conclusion

Determination of consensus molecular subtypes of MIBC from FFPE samples is feasible using various RNA sequencing methods. Inconsistent classification mainly involves the stroma-rich molecular subtype, which may be the consequence of sample heterogeneity with (stroma)-cell sampling bias and highlights the limitations of bulk RNA-based subclassification. Classification is still reliable when analysis is reduced to selected genes.

Peer Review reports

Background

Sequencing techniques have advanced, leading to broad genomic analyses of bladder cancer cohorts and enabling molecular subtyping. Subtyping of muscle-invasive urothelial bladder cancer (MIBC) categorizes heterogenous cancers with similar molecular and biological characteristics, which has significantly contributed to our knowledge in the recent years [1]. Several groups have simultaneously worked on molecular subtyping of different bladder cancer datasets coming to a description of two main types (luminal and basal), that can further be subclassified into 3–10 subtypes [2,3,4,5,6,7,8]. Different nomenclatures, definitions, and numbers of molecular subtypes had hindered further prospective validation, and clinical translation until the description of a consensus classification. The molecular consensus classification used pooled mRNA expression profiles of 1750 fresh frozen and formalin-fixed, paraffin-embedded (FFPE) MIBC samples and identified six molecular classes: luminal papillary (LumP), luminal nonspecified (LumNS), luminal unstable (LumU), stroma-rich, basal/squamous (Ba/Sq), and neuroendocrine-like (NE-like) [3].

The aim of classifying cancer in subgroups is to identify tumors that share similar prognosis and response to various therapies. In several reports molecular subtypes have been described as predictors of response to chemotherapy and immunotherapy [5, 9,10,11,12]. However, the results are conflicting, and to date, the evidence is insufficient to use molecular subtyping or other gene expression signatures for the treatment decisions in patients with urothelial carcinoma. To facilitate the implementation of molecular subtyping into daily clinical routine, gene sets have been reduced to allow quantification with quantitative RT-qPCR, NanoString or immunohistochemistry panels [13,14,15,16,17,18,19,20].

Many molecular profiling studies and molecular subtyping in The Cancer Genome Atlas (TCGA) are based on fresh frozen tissue, which allows high quality transcriptomic analyses based on RNA sequencing. However, fresh frozen samples are rarely available in clinical practice and for retrospective research projects. The use FFPE tissue is the gold-standard for pathological analyses and long-term storage in hospitals. The paraffin material is usually archived for 10 and more years, allowing correlation with long-term patient outcome. The disadvantage of FFPE tissue is that the RNA is highly degraded by fixation and storage, leading to sequencing artifacts and limiting detection of transcripts [21, 22]. However, advances in sequencing techniques also enable molecular profiling of FFPE specimen [23,24,25].

In this study, we tested and compared the feasibility of two RNA sequencing methods with FFPE tissues from MIBCs to determine uniform molecular subtyping. In addition, we used two reduced predefined gene sets to determine molecular subtypes and compared results with the comprehensive transcriptome data.

Methods

Patients and samples

Tumor samples were provided by the University Cancer Center Frankfurt (UCT). Written informed consent was obtained from all patients. The study was approved by the institutional review boards of the UCT and the ethical committee at the University Hospital Frankfurt (project-number: SUG-6–2018, UCT-53–2021), and conducted according to local and national regulations and to the Declaration of Helsinki. Fifteen FFPE samples from MIBC patients were obtained from the Dr. Senckenberg Institute of Pathology.

Immunohistochemistry (IHC)

Samples were stained for CK5/6 (Clone: D5/16 B4; Dako /Agilent, Santa Clara, CA, USA) and GATA3 (Clone: L50-823; Cell Marque, Rocklin, CA, USA) as described before [19].

RNA-isolation

For each RNA isolation, a 1-mm punch was taken from FFPE blocks of an annotated, representative tumor area with at least 50% tumor content. RNA was either extracted using the truXTRAC FFPE total NA Kit (Covaris, Woburn, MA, USA) or by GenXPro GmbH.

HTG transcriptome panel (HTP)

The mRNA expression was determined using the HTP (HTG Molecular Diagnostics, Tuscon, AZ, USA) as describes before [19, 26]. Briefly, target capture was done by the HTG EdgeSeq chemistry with nuclease protection probes on a 96-well plate. Processed samples were used to set up PCR reactions with specially designed primers, referred to as “tags”. These tags share common sequences that are complementary to 5’-end and 3’-end “wing” sequences of the probes and common adaptors required for cluster generation on an Illumina sequencing platform. In addition, each tag contains a unique barcode that is used for sample identification and multiplexing. The library was prepared using a PCR with OneTaq (New England Biolabs, Ipswich, MA, USA) and EdgeSeq PCR tag primers (HTG Molecular Diagnostics). Sequencing was performed on the Illumina NextSeq 550 system (Illumina, San Diego, CA, USA) in accordance with manufacturer’s recommendations but also including two HTG custom sequencing primers. The sequencing data on mRNA expression of target genes were imported into HTG EdgeSeq Parser software for alignment of FASTQ files to the to the probe list and quantification of the reads. The HTG EdgeSeq Reveal Application was utilized to quality check and normalize data. Gene counts were normalized using CPM and median normalization and log2-transformed for further analysis.

MACE Seq

Massive Analysis of cDNA (MACE) is a 3’mRNA sequencing method based on the analysis of Illumina reads derived from fragments that originate from 3’ mRNA ends [24, 27]. Samples were prepared by GenXPro GmbH (Frankfurt, Germany) using the MACE-Kit V2 according to the manual of the manufacturer (GenXPro GmbH). RNA was fragmented and polyadenylated mRNA was enriched by poly-A specific reverse transcription, a specific adapter was integrated at the 5’ ends, and the products were amplified by competitive PCR and sequenced on an Illumina NextSeq 500 instrument. Duplicate reads as determined by the implemented unique molecular identifiers (TrueQuant IDs) were removed from the raw dataset. Low quality sequence-bases were removed by the software Cutadapt (https://github.com/marcelm/cutadapt/) and poly(A)-tails were clipped by an in-house Python-Script. The reads were mapped to the human genome (hg38) and transcripts were quantified by HTSeq [28].

Molecular subtyping

Molecular consensus classes of MIBC were assigned using the consensusMIBC package for R for the nearest-centroid transcriptomic classifier (https://github.com/cit-bioinfo/consensusMIBC), TCGA classes were assigned using the BLCAsubtyping package (https://github.com/cit-bioinfo/BLCAsubtyping) as described by Kamoun et al. [3]. The minimal threshold for best Pearson’s correlation was set to 0.2. Normalized and log2-transformed gene expression values were used. Retrieved data include the consensus class, the Pearson’s correlation coefficient between each sample and each consensus class, the p-values associated to the Pearson’s correlation of the samples with the nearest centroid (correlation p-value) and the separation level.

For comparison, we used the included example data set of the TCGA bladder cancer cohort, which was created from fresh tumor specimen [2].

We reduced the gene set input according to two proposed panels for bladder cancer subtyping (Table S1 – S2) [16, 20]. The “ESSEN1”-panel is a 68-gene set covering tumor and stromal signatures [16]. The above panel was further optimized and condensed to a set of 48 genes (called “ESSEN2” in the present study) [20].

The heatmap was constructed with the open-source Morpheus software (https://software.broadinstitute.org/morpheus/) using log2-transformed normalized gene expression values. For hierarchical clustering we used Pearson’s correlation metric with complete linkage.

Results

Obtained sequencing data used for molecular classification

Characteristics of the patients included in the study are shown in Table 1.

Table 1 Characteristics of 15 urothelial muscle-invasive bladder cancer patients included in the paired comparison

With MACE, 22,729 of 29,716 transcripts were mapped to HUGO Gene Nomenclature Committee (HGNC) symbols. With the HTP, 19,399 mRNA transcripts were detected. 16,645 genes had the same HGNC symbol, and 6084 and 2754 genes were detected with either MACE or HTP only, respectively. One HTP sample failed quality control (after sequencing) because of minimal expression variability as determined by median log2 for negative control probes. The median time from sample collection (time of surgery) to analysis was 5 years (range 3–10 years); seven samples were < 5 years and eight samples were ≥ 5 years. There was no significant correlation between the age of the samples and the overall degree of separation, nor was the association of the Pearson correlation coefficients of the samples with the nearest centroid (correlation p-value), for either method (Figure S1).

Comparison of classification values between FFPE-sequencing methods and TCGA-data

We used MACE-, HTP- as well as TCGA-transcriptome data to assign molecular consensus subtypes. We analyzed the data and compared the output information. The association of the Pearson’s correlation coefficients of the samples with the nearest centroid (correlation p-value) was < 0.0001 for all samples with either method. However, we noticed significantly higher median correlation p-values for subtype calls resulting from HTP vs. MACE vs. TCGA data (Figure S2).

Figure 1 shows the correlation plots as radar-charts for each consensus class for the 406 TCGA, 14 HTP and 15 MACE samples. Using the MACE data resulted in higher correlation values of the sample to the centroids of the six consensus classes, which was only statistically significant for the Ba/Sq subtype (p(Wilcoxon) = 0.0092). With the HTP data, NE-like and Ba/Sq subtype had the highest median separation levels of 0.67 and 0.79, like the TCGA data. Using the MACE data, the stroma-rich subtype had the highest median separation level of 0.71. Based on MACE data, the NE-like subtype had a very low separation level of 0.03. However, this was only based on a single sample. With the HTP and MACE data LumNS subtype and with the TCGA data LumNS and LumU had rather low separation levels, which might be due to other luminal groups.

Fig. 1
figure 1

Radar-charts for each consensus class for TCGA, HTP and MACE data. The data are generated from the consensus classifier. Blue lines indicate the Pearson’s correlation coefficient between each sample and each consensus class. Red lines and values indicate the median Pearson’s correlation between samples and each consensus class. The median separation level gives a measure of how a sample is representative of its consensus class and is calculated as follows: (correlation to nearest centroid—correlation to second nearest centroid) / median difference of sample-to-centroid correlation

Molecular and histological subtypes and IHC-status of patients

Samples were classified as Ba/Sq in 7 (50%), LumP in 2 (14.3%), LumU in 1 (7.1%), LumNS in 1 (7.1%), NE-like in 1 (7.1%) and stroma-rich in 2 (14.3%) cases with HTP data. And classified as Ba/Sq in 8 (53.3%), LumP in 2 (13.3%), LumU in 1 (6.7%), LumNS in 1 (6.7%), NE-like in 1 (6.7%) and stroma-rich in 2 (13.3%) cases with MACE data (Table 2). Concordance of the consensus classification was 71% (10/14). Two cases (48 and 46) classified as stroma-rich using MACE data were Ba/Sq (48) and LumNS (46) with HTP data. Two cases (49 and 20) classified as stroma-rich with HTP data were Ba/Sq (49) and LumNS (20) with MACE data. Thus, all four discordant cases had a stroma-rich subtype involved with either method. Those cases showed a usual (NOS) histological subtype and no or only low CK5/6 protein expression. An example of one discordant case is shown in Fig. 2. Importantly, LumU and LumP subtypes were concordant in 100% and displayed the expected matching histological subtypes (NOS) and protein expression (CK5/6 negative, GATA3 positive). One NE-like case showed neuroendocrine histological subtype and was negative for CK5/6 and GATA3 IHC. The Ba/Sq was concordant in 86% (6/7) and displayed squamous histological subtype in four cases (57%), NOS in two cases (29%), and lymphoepithelial histological subtype in one case, which has been described to be associated with Ba/Sq molecular subtype [29].

Table 2 Molecular subtype, histological subtype and immunohistochemical expression of CK5/6 and GATA3 for each patient
Fig. 2
figure 2

Case 46, which was classified as stroma-rich with the MACE data and LumNS with the HTP data, shown in the HE-overview with tumor annotations for cores for the construction of TMAs and RNA-Isolation (A). Further images show HE-staining (B), positive GATA3-IHC (C) and negative CK5/6-IHC (D), 200x

Reduction of gene set for molecular subtyping

We used two previously described gene sets, that were developed for molecular subtyping of MIBC, named ESSEN1 and ESSEN2 in the present study [16, 20].

With HTP data, samples were classified as Ba/Sq in 6 (42.8%) and 7 (50%), LumP in 3 (21.4%) and 2 (14.3%), LumU in 1 (7.1%) and 1 (7.1%), LumNS in 0 and 1 (7.1%), NE-like in 1 (7.1%) and 1 (7.1%), stroma-rich in 3 (21.4%) and 2 (14.3%) cases using ESSEN1 and ESSEN2 gene sets, respectively. Using MACE data, samples were classified as Ba/Sq in 5 (35.7%) and 8 (57.1%), LumP in 2 (14.3%) and 1 (7.1%), LumU in 1 (7.1%) and 2 (14.3%), LumNS in 1 (7.1%) and 2 (14.3%), NE-like in 1 (7.1%) and 0 and stroma-rich in 3 (21.4%) and 2 (14.3%) cases using ESSEN1 and ESSEN2 gene sets, respectively (Figure S3).

Compared to the comprehensive transcriptomic data, the reduction in gene sets resulted in a significant increase of the correlation p-values associated to the correlation of the samples with the nearest centroids (Figure S4). The median correlation p-value was still < 0.0001. One sample (44) could not be classified with the MACE data and was excluded for the further analyses because the correlation p-value was above the set threshold of 0.2. Of note, this sample showed a lymphoepithelial histological phenotype, and the heavy infiltration of immune cells might have hindered correct classification.

The Pearson’s correlations coefficient for samples and each consensus class stratified for the called consensus classes are included in the supplementary material.

For the ESSEN1 panel, consensus molecular subtypes were concordant using the HTP data in 86% (12/14). One sample (46) changed from LumNS to LumP and one sample (32) changed from Ba/Sq to stroma-rich. Using the MACE data, concordance was 86% (12/14). One sample (32) changed from Ba/Sq to LumNS and one sample (8) from Ba/Sq to stroma-rich (Fig. 3).

Fig. 3
figure 3

Heatmap showing the molecular subtypes according to the consensus classification (A) and the TCGA classification (B) for each sample as identified by all genes vs. different gene sets (all, ESSEN1, ESSEN2) and assay methods (HTP, MACE). C: Gene expression heatmap of MACE data and the ESSEN1 panel. D: Gene expression heatmap of the ESSEN2 panel with molecular consensus classes assigned using the HTP data. The heatmap was constructed with the open-source Morpheus software (https://software.broadinstitute.org/morpheus/) using log2-transformed normalized gene expression values. Hierarchical clustering was performed with Pearson’s correlation metric with complete linkage. Heatmaps also indicate the IHC-expression of CK5/6 and GATA3. Colors of genes correlate to the signature names. CIS: carcinoma in situ; ECM: extracellular matrix; EMT: epithelial-to-mesenchymal transition; SM: smooth muscle

For the ESSEN2 panel concordance of the consensus molecular subtypes was 100% (14/14) with the HTP data (Fig. 3). Using the MACE data with the ESSEN2 panel, concordance was 86% (12/14). One sample (7) changed from LumP to LumU. Of note, one sample (26) with NE-like molecular subtype, that also showed neuroendocrine histological subtype and was negative for immunohistochemical CK5/6 and GATA3 expression, was misclassified by the ESSEN2 panel when using the consensus (but not TCGA classification) with the MACE (but not HTP) sequencing data as Ba/Sq. This sample also showed very low separation level when using all genes from the MACE data as input.

The overlap between subtypes according to the TCGA classification was 86% with the ESSEN1 and 93% with the ESSEN2 panel using the HTP data (Fig. 3B). With the MACE data the overlap was 86% (12/14) with the ESSEN2 panel and 64% with the ESSEN1 panel (two samples changed within luminal subtypes and three samples had discordant subtypes between basal/squamous and luminal subtypes).

Discussion

Molecular analyses of bladder cancer specimens are emerging and provide elementary information to address bladder cancer research questions. One of the major goals of bladder cancer research is to identify subtypes with different sensitivity to systemic therapies such as chemotherapy, immune-checkpoint inhibition, and further targeted treatments. Although molecular subtype classification of bladder cancer has not yet been incorporated into therapeutic decision making, robust methods are important to achieve progress in clinical translation and validation and to improve reproducibility. FFPE samples represent snapshots of the histology and biological information of a tumor at the time of collection, while the patient is being treated and outcome data can be generated. It is important to use the information, that is always collected and stored as FFPE and to overcome limitations of degraded RNA.

In this study, we used FFPE-isolated RNA to determine molecular consensus subtypes using two different sequencing techniques. The overall agreement between molecular subtypes was high for both techniques, although different RNA, library preparation and sequencing facilities were used. We validated two reduced gene sets to determine molecular subtypes with high accuracy, that can be used as a panel-based method at lower cost, which is an important step to introduce subtyping into routine practice.

Efforts are being made to methodologically simplify subtype classification by using reduced gene sets to enable its implementation in clinical practice. Stratifying patients will be important to select patients, that respond to chemo- or immunotherapy to reduce unnecessary toxicities and costs, as soon as prospective and validated evidence is provided. In a previous study, a 47-gene panel (BASE47) was proposed for the discrimination between luminal and basal subtypes using the NanoString platform [14]. Recently, the ESSEN1 (n = 68) and ESSEN2 (n = 48) gene panels were developed to discriminate between 5–6 gene expression-based molecular subtypes according to different classification systems (e.g. TCGA, Consensus, Lund etc.) by using the qRT-PCR method on fresh frozen samples and the NanoString technique on FFPE samples, respectively [16, 20]. These reduced gene sets still allowed the classifier to designate molecular consensus molecular subtypes in all except for one case, which had a lymphoepithelial histological subtype. Our results show high accuracy of 85%-100% of both reduced panels to call consensus molecular subtypes, compared to the comprehensive transcriptomic data. Discordant subtypes were observed between the stroma-rich and the Ba/Sq or the LumNS subtype and between the different luminal subtypes. Divergent subtype calls can be the result of either differences in gene expression or the composition of genes used for calculation. According to our data, the reduced ESSEN2 panel showed an even higher overlap with the comprehensive transcriptome data than the ESSEN1 panel. Thus, the selection of genes might be more important for the classification than just the quantity.

So far as no use-case for the application of molecular subtypes exists, a comprehensive gene expression analysis provides additional information on further genes and enables the analysis of additional immune or stroma signatures, which could find application as predictive markers [30,31,32].

Limitations include the small number of samples. With only one sample for some of the subtypes general conclusions on which method is superior cannot be drawn. Furthermore, we did not perform a direct comparison between FFPE and fresh tissue. Our study lacks the comparison of HTP and MACE to further sequencing techniques used in previous molecular subtyping studies, like Illumina RNA seq or Affymetrix. Most importantly, sequencing was not performed with the exact same RNA, but RNA from neighboring tumor areas. Thus, discordant results can be the result of slight differences in tumor and stroma cell content or differences in deeper tissue layers not represented on the slide. However, this issue might reflect the heterogeneity of bladder cancer itself. This becomes particularly evident with the stroma-rich molecular subtype, since it was present in all divergent classified samples and highlights problems affecting bulk-RNA sequencing of bladder cancer samples in general. Molecular subtyping is performed on RNA derived from tissue (cores) of biopsy specimens representing only a fraction of the tumor. The subtype determined from the isolated RNA might not be the only and/or predominant subtype of a tumor, which is known to be heterogenous, especially in bladder cancer [17, 33]. Furthermore, the subtype assigned depends on the cellular components, of which the RNA is extracted and the proportion of stromal content will influence, if a tumor is called as stroma-rich or infiltrated [17, 31, 34, 35].

Conclusion

The consensus molecular subtypes represent a robust classification and can be determined based on comprehensive or selected FFPE transcriptome data. Using the data of matching pairs, the agreement of subtypes was high, but differences were observed when the stroma-rich molecular subtype was involved. Based on our results, it seems important to further unravel the heterogeneity of bladder cancer and the influence of stromal components on molecular subtypes to reduce sampling bias and allow more accurate assignment.

Availability of data and materials

The transcriptomic dataset generated and analyzed during the current study in the are available in the GEO repository: GSE225376. Further data that support the findings of this study are available from the corresponding author upon reasonable request.

Abbreviations

Ba/Sq:

Basal/squamous

FFPE:

Formalin-fixed, paraffin-embedded

HTG:

HTG Molecular Diagnostics, Inc.

HTP :

HTG transcriptome panel

IQR:

Interquartile range

LumNS:

Luminal nonspecified

LumP:

Luminal papillary

MACE:

Massive Analysis of cDNA Ends

MIBC:

Muscle-invasive bladder cancer

NE-like:

Neuroendocrine-like

RNA:

Ribonucleic acid

RT-qPCR:

Real-time quantitative polymerase chain reaction

TCGA:

The Cancer Genome Atlas

UCT:

University Cancer Center Frankfurt

References

  1. Jalanko T, de Jong JJ, Gibb EA, Seiler R, Black PC. Genomic Subtyping in Bladder Cancer. Curr Urol Rep. 2020;21(2):9.

    Article  PubMed  Google Scholar 

  2. Robertson AG, Kim J, Al-Ahmadie H, Bellmunt J, Guo G, Cherniack AD, Hinoue T, Laird PW, Hoadley KA, Akbani R, et al. Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer. Cell. 2017; 171(3):540–556 e525.

  3. Kamoun A, de Reynies A, Allory Y, Sjodahl G, Robertson AG, Seiler R, Hoadley KA, Groeneveld CS, Al-Ahmadie H, Choi W, et al. A Consensus Molecular Classification of Muscle-invasive Bladder Cancer. Eur Urol. 2020;77(4):420–33.

  4. Choi W, Porten S, Kim S, Willis D, Plimack ER, Hoffman-Censits J, Roth B, Cheng T, Tran M, Lee IL, et al. Identification of distinct basal and luminal subtypes of muscle-invasive bladder cancer with different sensitivities to frontline chemotherapy. Cancer Cell. 2014;25(2):152–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Seiler R, Ashab HAD, Erho N, van Rhijn BWG, Winters B, Douglas J, Van Kessel KE, Fransen van de Putte EE, Sommerlad M, Wang NQ, et al. Impact of Molecular Subtypes in Muscle-invasive Bladder Cancer on Predicting Response and Survival after Neoadjuvant Chemotherapy. Eur Urol. 2017;72(4):544–554.

  6. Damrauer JS, Hoadley KA, Chism DD, Fan C, Tiganelli CJ, Wobker SE, Yeh JJ, Milowsky MI, Iyer G, Parker JS, et al. Intrinsic subtypes of high-grade bladder cancer reflect the hallmarks of breast cancer biology. Proc Natl Acad Sci USA. 2014;111(8):3110–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Sjodahl G, Lauss M, Lovgren K, Chebil G, Gudjonsson S, Veerla S, Patschan O, Aine M, Ferno M, Ringner M, et al. A molecular taxonomy for urothelial carcinoma. Clin Cancer Res. 2012;18(12):3377–86.

    Article  PubMed  Google Scholar 

  8. Aine M, Eriksson P, Liedberg F, Sjodahl G, Hoglund M. Biological determinants of bladder cancer gene expression subtypes. Sci Rep. 2015;5:10957.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Kim J, Kwiatkowski D, McConkey DJ, Meeks JJ, Freeman SS, Bellmunt J, Getz G, Lerner SP. The Cancer Genome Atlas Expression Subtypes Stratify Response to Checkpoint Inhibition in Advanced Urothelial Cancer and Identify a Subset of Patients with High Survival Probability. Eur Urol. 2019;75(6):961–4.

    Article  CAS  PubMed  Google Scholar 

  10. Lotan Y, de Jong JJ, Liu VYT, Bismar TA, Boorjian SA, Huang HC, Davicioni E, Mian OY, Wright JL, Necchi A, et al. Patients with Muscle-Invasive Bladder Cancer with Nonluminal Subtype Derive Greatest Benefit from Platinum Based Neoadjuvant Chemotherapy. J Urol. 2022;207(3):541–50.

  11. Powles T, Assaf ZJ, Davarpanah N, Banchereau R, Szabados BE, Yuen KC, Grivas P, Hussain M, Oudard S, Gschwend JE, et al. ctDNA guiding adjuvant immunotherapy in urothelial carcinoma. Nature. 2021;595(7867):432–7.

    Article  CAS  PubMed  Google Scholar 

  12. Sjodahl G, Abrahamsson J, Holmsten K, Bernardo C, Chebil G, Eriksson P, Johansson I, Kollberg P, Lindh C, Lovgren K, et al. Different Responses to Neoadjuvant Chemotherapy in Urothelial Carcinoma Molecular Subtypes. 2022;81(5):523–32.

  13. Guo CC, Bondaruk J, Yao H, Wang Z, Zhang L, Lee S, Lee JG, Cogdell D, Zhang M, Yang G, et al. Assessment of Luminal and Basal Phenotypes in Bladder Cancer. Sci Rep. 2020;10(1):9743.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Kardos J, Rose TL, Manocha U, Wobker SE, Damrauer JS, Bivalaqua TJ, Kates M, Moore KJ, Parker JS, Kim WY. Development and validation of a NanoString BASE47 bladder cancer gene classifier. PLoS ONE. 2020;15(12): e0243935.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kollberg P, Chebil G, Eriksson P, Sjodahl G, Liedberg F. Molecular subtypes applied to a population-based modern cystectomy series do not predict cancer-specific survival. Urol Oncol. 2019;37(10):791–9.

    Article  CAS  PubMed  Google Scholar 

  16. Olah C, Hahnen C, Nagy N, Musial J, Varadi M, Nyiro G, Gyorffy B, Hadaschik B, Rawitzer J, Ting S, et al. A quantitative polymerase chain reaction based method for molecular subtype classification of urinary bladder cancer-Stromal gene expressions show higher prognostic values than intrinsic tumor genes. Int J Cancer J Int Cancer. 2021;150(5):856–67.

  17. Sjodahl G, Eriksson P, Lovgren K, Marzouka NA, Bernardo C, Nordentoft I, Dyrskjot L, Liedberg F, Hoglund M. Discordant molecular subtype classification in the basal-squamous subtype of bladder tumors and matched lymph-node metastases. Mod Pathol. 2018;31(12):1869–81.

    Article  PubMed  Google Scholar 

  18. Weyerer V, Stoehr R, Bertz S, Lange F, Geppert CI, Wach S, Taubert H, Sikic D, Wullich B, Hartmann A, et al. Prognostic impact of molecular muscle-invasive bladder cancer subtyping approaches and correlations with variant histology in a population-based mono-institutional cystectomy cohort. World J Urol. 2021;39(11):4011–9.

  19. Koll FJ, Schwarz A, Kollermann J, Banek S, Kluth L, Wittler C, Bankov K, Doring C, Becker N, Chun FKH, et al. CK5/6 and GATA3 Defined Phenotypes of Muscle-Invasive Bladder Cancer: Impact in Adjuvant Chemotherapy and Molecular Subtyping of Negative Cases. Front Med (Lausanne). 2022;9: 875142.

    Article  PubMed  Google Scholar 

  20. Olah C, Reis H, Hoffmann MJ, Mairinger F, Ting S, Hadaschik B, Krafft U, Grunwald V, Nyirady P, Varadi M, et al. Predictive value of molecular subtypes and APOBEC3G for adjuvant chemotherapy in urothelial bladder cancer. Cancer Med. 2022; 12(5):5222–32.

  21. von Ahlfen S, Missel A, Bendrat K, Schlumpberger M. Determinants of RNA quality from FFPE samples. PLoS ONE. 2007;2(12): e1261.

    Article  Google Scholar 

  22. Groelz D, Sobin L, Branton P, Compton C, Wyrich R, Rainen L. Non-formalin fixative versus formalin-fixed tissue: a comparison of histology and RNA quality. Exp Mol Pathol. 2013;94(1):188–94.

    Article  CAS  PubMed  Google Scholar 

  23. Adiconis X, Borges-Rivera D, Satija R, DeLuca DS, Busby MA, Berlin AM, Sivachenko A, Thompson DA, Wysoker A, Fennell T, et al. Comparative analysis of RNA sequencing methods for degraded or low-input samples. Nat Methods. 2013;10(7):623–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Boneva S, Schlecht A, Bohringer D, Mittelviefhaus H, Reinhard T, Agostini H, Auw-Haedrich C, Schlunck G, Wolf J, Lange C. 3’ MACE RNA-sequencing allows for transcriptome profiling in human tissue samples after long-term storage. Lab Invest. 2020;100(10):1345–55.

  25. Qi Z, Wang L, Desai K, Cogswell J, Stern M, Lawson B, Kerkar SP, Vitazka P. Reliable Gene Expression Profiling from Small and Hematoxylin and Eosin-Stained Clinical Formalin-Fixed, Paraffin-Embedded Specimens Using the HTG EdgeSeq Platform. J Mol Diagn. 2019;21(5):796–807.

    Article  CAS  PubMed  Google Scholar 

  26. Koll FJ, Metzger E, Hamann J, Ramos-Triguero A, Bankov K, Köllermann J, et al. Overexpression of KMT9α Is Associated with Aggressive Basal-like Muscle-Invasive Bladder Cancer. Cells. 2023;12(4):589. https://doi.org/10.3390/cells12040589.

  27. Muller N, Warwick T, Noack K, Malacarne PF, Cooper AJL, Weissmann N, et al. Reactive Oxygen Species Differentially Modulate the Metabolic and Transcriptomic Response of Endothelial Cells. Antioxidants (Basel). 2022;11(2):434.

  28. Anders S, Pyl PT, Huber W. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.

    Article  CAS  PubMed  Google Scholar 

  29. Manocha U, Kardos J, Selitsky S, Zhou M, Johnson SM, Breslauer C, Epstein JI, Kim WY, Wobker SE. RNA Expression Profiling of Lymphoepithelioma-Like Carcinoma of the Bladder Reveals a Basal-Like Molecular Subtype. Am J Pathol. 2020;190(1):134–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Wang L, Saci A, Szabo PM, Chasalow SD, Castillo-Martin M, Domingo-Domenech J, Siefker-Radtke A, Sharma P, Sfakianos JP, Gong Y, et al. EMT- and stroma-related gene expression and resistance to PD-1 blockade in urothelial cancer. Nat Commun. 2018;9(1):3503.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Sjodahl G, Jackson CL, Bartlett JM, Siemens DR, Berman DM. Molecular profiling in muscle-invasive bladder cancer: more than the sum of its parts. J Pathol. 2019;247(5):563–73.

    Article  PubMed  Google Scholar 

  32. Damrauer JS, Beckabir W, Klomp J, Zhou M, Plimack ER, Galsky MD, Grivas P, Hahn NM, O’Donnell PH, Iyer G, et al. Collaborative study from the Bladder Cancer Advocacy Network for the genomic analysis of metastatic urothelial cancer. Nat Commun. 2022;13(1):6658.

  33. Meeks JJ, Al-Ahmadie H, Faltas BM, Taylor JA 3rd, Flaig TW, DeGraff DJ, Christensen E, Woolbright BL, McConkey DJ, Dyrskjot L. Genomic heterogeneity in bladder cancer: challenges and possible solutions to improve outcomes. Nat Rev Urol. 2020;17(5):259–70.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Sjodahl G, Liedberg F, Hoglund M, Eriksson P. When the Molecular Subtype Is Hidden Behind a Veil of Stroma. Eur Urol. 2021;80(2):160–1.

    Article  PubMed  Google Scholar 

  35. Yuan Y. Spatial Heterogeneity in the Tumor Microenvironment. Cold Spring Harb Perspect Med. 2016;6(8):a026583.

Download references

Acknowledgements

The team of the Senckenberg Biobank is thanked for excellent technical assistance.

Funding

Open Access funding enabled and organized by Projekt DEAL. Florestan Koll was funded by the Mildred Scheel Career Center Frankfurt (Deutsche Krebshilfe).

Author information

Authors and Affiliations

Authors

Contributions

FJK: Conceptualization; Data curation; Formal analysis; Funding acquisition; Investigation; Visualization; Methodology; Writing – original draft; Writing – review & editing. CD: Investigation; Data curation; Formal analysis. CO: Writing – review & editing. TS: Writing – review & editing. JK: Investigation; Writing – review & editing. BH: Writing – review & editing. FC: Writing – review & editing; Supervision; Resources. HR: Conceptualization; Investigation; Supervision; Data curation; Writing – original draft; Writing – review & editing. PJW: Conceptualization; Supervision; Project administration; Resources; Writing – review & editing. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Florestan J. Koll.

Ethics declarations

Ethics approval and consent to participate

Tissue/tumor samples and/or patient data used in this study were provided by the University Cancer Center Frankfurt (UCT). Written informed consent was obtained from all patients and the study was approved by the institutional Review Boards of the UCT and the Ethical Committee at the University Hospital Frankfurt (project-number: SUG-6–2018 and UCT-53–2021). The study was conducted according to local and national regulations and according to the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

HTG Transcriptome Panel was performed as part of a collaboration agreement with HTG Molecular Diagnostics. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. PJW has received consulting fees and honoraria for lectures by Bayer, Janssen-Cilag, Novartis, Roche, MSD, Astellas Pharma, Bristol-Myers Squibb, Thermo Fisher Scientific, Molecular Health, Guardant Health, Sophia Genetics, Qiagen, Eli Lilly, Myriad, Hedera Dx, and AstraZeneca. Research Support was provided by AstraZeneca and Roche. HR was on an advisory board of Bristol-Myers Squibb, received honoraria from Roche, Bristol-Myers Squibb, Janssen-Cilag, Novartis, Astra-Zeneca, MCI, CHOP GmbH and Diaceutics, received travel support from Philips, Roche, and Bristol-Myers Squibb, and received grants from Bristol-Myers Squibb. All other authors indicate that there are 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:

Supplementary material. Table S1. ESSEN1 gene set. Table S2. ESSEN2 gene set. Figure S1. Comparison of correlation p-value between for samples older and younger than five years for MACE and HTG. Figure S2. Pearson’s correlation between samples and each consensus class stratified for the different methods and for the called consensus class. Figure S3. Relative distribution of subtype calls with each method and gene panels. Figure S4. A: Comparison of p-values associated to the Pearson’s correlation of the samples with the nearest centroid for each method (HTG and MACE) for the ESSEN1 and ESSEN2 gene panels. Figure S5. Pearson’s correlation between samples and each consensus class stratified for the different gene sets.

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

Koll, F.J., Döring, C., Olah, C. et al. Optimizing identification of consensus molecular subtypes in muscle-invasive bladder cancer: a comparison of two sequencing methods and gene sets using FFPE specimens. BMC Cancer 23, 504 (2023). https://doi.org/10.1186/s12885-023-11016-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11016-9

Keywords