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.

No abstract

Discrete bladder cancer molecular subtypes exhibit differential clinical aggressiveness and therapeutic response, which may have significant implications for identifying novel treatments for this common malignancy. However, research is hindered by the lack of suitable models to study each subtype. To address this limitation, we classified bladder cancer cell lines into molecular subtypes using publically available data in the Cancer Cell Line Encyclopedia (CCLE), guided by genomic characterization of bladder cancer by The Cancer Genome Atlas (TCGA). This identified a panel of bladder cancer cell lines which exhibit genetic alterations and gene expression patterns consistent with luminal and basal molecular subtypes of human disease. A subset of bladder cancer cell lines exhibit in vivo histomorphologic patterns consistent with luminal and basal subtypes, including papillary architecture and squamous differentiation. Using the molecular subtype assignments, and our own RNA-seq analysis, we found overexpression of GATA3 and FOXA1 cooperate with PPARɣ activation to drive transdifferentiation of a basal bladder cancer cells to a luminial phenotype. In summary, our analysis identified a set of human cell lines suitable for the study of molecular subtypes in bladder cancer, and furthermore indicates a cooperative regulatory network consisting of GATA3, FOXA1, and PPARɣ drive luminal cell fate.