Abstract
Recent efforts on genome wide profiling of muscle invasive bladder cancer (MIBC) have led to its classification into distinct genomic and transcriptomic molecular subtypes that exhibit variability in prognosis. Evolving evidence from recent immunotherapy trials has demonstrated the significance of pre-existing tumour immune profiles that could guide treatment decisions. To identify immune gene expression patterns associated with the molecular subtypes, we performed a comprehensive in silico immune transcriptomic profiling, utilizing transcriptomic data from 347 MIBC cases from The Cancer Genome Atlas (TCGA). To investigate subtype-associated immune gene expression patterns, we assembled 924 immune response genes and specifically those involved in T-cell cytotoxicity and the Type I/II interferon pathways. A set of 157 ranked genes was able to distinguish the four subtypes in an unsupervised analysis in an original training cohort (n=122) and an expanded, validation cohort (n=225). The most common overrepresented pathways distinguishing the four molecular subtypes, included JAK/STAT signaling, Toll-like receptor signaling, interleukin signaling, and T-cell activation. Some of the most enriched biological processes were responses to IFN-γ, antigen processing and presentation, cytokine mediated signaling, hemopoeisis, cell proliferation and cellular defense response in the TCGA cluster IV. Our novel findings provide further insights into the association between genomic subtypes and immune activation in MIBC and may open novel opportunities for their exploitation towards precise treatment with immunotherapy.