Abstract
Although the intravesical instillation of Bacillus Calmette-Guerin (BCG) is widely used as adjuvant treatment for nonmuscle-invasive bladder cancers, the clinical benefit is variable across patients, and the molecular mechanisms underlying the sensitivity to BCG administration and disease progression are poorly understood. To establish the molecular signatures that predict the responsiveness and disease progression of bladder cancers treated with BCG, we performed transcriptome sequencing (RNA-seq) for 13 treatment-naïve and 22 post-treatment specimens obtained from 14 bladder cancer patients. To overcome disease heterogeneity, we used non-negative matrix factorization to identify the latent molecular features associated with drug responsiveness and disease progression. At least 12 molecular features were present, among which the immune-related feature was associated with drug responsiveness, indicating that pre-treatment anti-cancer immunity might dictate BCG responsiveness. We also identified disease progression-associated molecular features indicative of elevated cellular proliferation in post-treatment specimens. The progression-associated molecular features were validated in an extended cohort of BCG-treated bladder cancers. Our study advances understanding of the molecular mechanisms of BCG activity in bladder cancers and provides clinically relevant gene markers for evaluating and monitoring patients.