Non-muscle-invasive bladder cancer (NMIBC) is a heterogeneous disease with widely different outcomes. We performed a comprehensive transcriptional analysis of 460 early-stage urothelial carcinomas and showed that NMIBC can be subgrouped into three major classes with basal- and luminal-like characteristics and different clinical outcomes. Large differences in biological processes such as the cell cycle, epithelial-mesenchymal transition, and differentiation were observed. Analysis of transcript variants revealed frequent mutations in genes encoding proteins involved in chromatin organization and cytoskeletal functions. Furthermore, mutations in well-known cancer driver genes (e.g., TP53 and ERBB2) were primarily found in high-risk tumors, together with APOBEC-related mutational signatures. The identification of subclasses in NMIBC may offer better prognostication and treatment selection based on subclass assignment.

Different molecular subtypes can be identified in muscle-invasive bladder cancer. Although cisplatin-based neoadjuvant chemotherapy improves patient outcomes, we identified that the benefit is highest in patients with basal tumors. Our newly discovered classifier can identify these molecular subtypes in a single patient and could be integrated into routine clinical practice after further validation.

The molecular landscape in non-muscle-invasive bladder cancer (NMIBC) is characterized by large biological heterogeneity with variable clinical outcomes. Here, we perform an integrative multi-omics analysis of patients diagnosed with NMIBC (n = 834). Transcriptomic analysis identifies four classes (1, 2a, 2b and 3) reflecting tumor biology and disease aggressiveness. Both transcriptome-based subtyping and the level of chromosomal instability provide independent prognostic value beyond established prognostic clinicopathological parameters. High chromosomal instability, p53-pathway disruption and APOBEC-related mutations are significantly associated with transcriptomic class 2a and poor outcome. RNA-derived immune cell infiltration is associated with chromosomally unstable tumors and enriched in class 2b. Spatial proteomics analysis confirms the higher infiltration of class 2b tumors and demonstrates an association between higher immune cell infiltration and lower recurrence rates. Finally, the independent prognostic value of the transcriptomic classes is documented in 1228 validation samples using a single sample classification tool. The classifier provides a framework for biomarker discovery and for optimizing treatment and surveillance in next-generation clinical trials.

Analyzing hematuria patients for the risk of bladder cancer using novel molecular markers may lead to a reduction in diagnostic cystoscopy. Combining methylation analysis (OTX1, ONECUT2 and TWIST1) with mutation analysis (FGFR3, TERT and HRAS) and patient age resulted in a validated accurate prediction model.