Trine Strandgaard

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.

PURPOSE: To investigate the use of plasma and urine DNA mutation analysis for predicting neoadjuvant chemotherapy (NAC) response and oncological outcome in patients with muscle-invasive bladder cancer. EXPERIMENTAL DESIGN: Whole-exome sequencing of tumor and germline DNA was performed for 92 patients treated with NAC followed by radical cystectomy (RC). A custom NGS-panel capturing approximately 50 mutations per patient was designed and used to track mutated tumor DNA in plasma and urine. A total of 447 plasma samples, 281 urine supernatants, and 123 urine pellets collected before, during, and after treatment were analyzed. Patients were enrolled from 2013 to 2019, with a median follow-up time of 41.3 months after RC. RESULTS: We identified tumor DNA before NAC in 89% of urine supernatants, 85% of urine pellets, and 43% of plasma samples. Tumor DNA levels were higher in urine supernatants and urine pellets compared with plasma samples (P < 0.001). In plasma, detection of circulating tumor DNA (ctDNA) before NAC was associated with a lower NAC response rate (P < 0.001). Detection of tumor DNA after NAC was associated with lower response rates in plasma, urine supernatant, and urine pellet (P < 0.001, P = 0.03, P = 0.002). Tumor DNA dynamics during NAC was predictive of NAC response and outcome in urine supernatant and plasma (P = 0.006 and P = 0.002). A combined measure from plasma and urine supernatant tumor DNA dynamics stratified patients by outcome (P = 0.003). CONCLUSIONS: Analysis of tumor DNA in plasma and urine samples both separately and combined has a potential to predict treatment response and outcome.

BACKGROUND: The functional status of immune cells in the tumor microenvironment and tumor characteristics may explain bacillus Calmette-Guérin (BCG) failure in high-risk non-muscle-invasive bladder cancer (NMIBC). OBJECTIVE: To characterize molecular correlates of post-BCG high-grade (HG) recurrence using multiomics analysis. DESIGN, SETTING, AND PARTICIPANTS: Patients with BCG-treated NMIBC (n = 156) were included in the study. Metachronous tumors were analyzed using RNA sequencing (n = 170) and whole-exome sequencing (n = 195). Urine samples were analyzed for immuno-oncology-related proteins (n = 190) and tumor-derived DNA (tdDNA; n = 187). OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: The primary endpoint was post-BCG HG recurrence. Cox regression and Wilcoxon rank-sum, t, and Fisher's exact tests were used for analyses. RESULTS AND LIMITATIONS: BCG induced activation of the immune system regardless of clinical response; however, immunoinhibitory proteins were observed in the urine of patients with post-BCG HG recurrence (CD70, PD1, CD5). Post-BCG HG recurrence was associated with post-BCG T-cell exhaustion (p = 0.002). Pre-BCG tumors from patients with post-BCG T-cell exhaustion had high expression of genes related to cell division and immune function. A high predicted post-BCG exhaustion score for pre-BCG tumors was associated with worse post-BCG HG recurrence-free survival (HGRFS; p = 0.002). This was validated in independent cohorts. Pre-BCG class 2a and 2b tumors (UROMOL2021 scheme) were associated with worse post-BCG HGRFS (p = 0.015). Post-BCG exhaustion was observed in patients with high pre-BCG neoantigen load (p = 0.017) and MUC4 mutations (p = 0.002). Finally, the absence of post-BCG tdDNA clearance identified patients at high risk of recurrence (p = 0.018). The retrospective design and partial overlap for analyses are study limitations. CONCLUSIONS: Post-BCG HG recurrence may be caused by T-cell exhaustion. Tumor subtype and pre-BCG tumor characteristics may identify patients at high risk of post-BCG HG recurrence. Urinary measurements have potential for real-time assessment of treatment response. PATIENT SUMMARY: A dysfunctional immune response to bacillus Calmette-Guérin (BCG) therapy may explain high-grade recurrences of bladder cancer.

BACKGROUND: Field cancerization is characterized by areas of normal tissue affected by mutated clones. Bladder field cancerization may explain the development and recurrence of bladder cancer and may be associated with treatment outcomes. OBJECTIVE: To investigate the predictive and prognostic roles of field cancerization in patients with high-risk non-muscle-invasive bladder cancer (NMIBC) treated with bacillus Calmette-Guérin (BCG). DESIGN, SETTING, AND PARTICIPANTS: We conducted comprehensive genomic and proteomic analyses for 751 bladder biopsies and 234 urine samples from 136 patients with NMIBC. The samples were collected at multiple time points during the disease course. Field cancerization in normal-appearing bladder biopsies was measured using deep-targeted sequencing and error correction models. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Endpoints included the rates of recurrence and progression. Cox regression and Wilcoxon rank-sum and Fisher's exact tests were used. RESULTS AND LIMITATIONS: A high level of field cancerization was associated with high tumor mutational burden (p = 0.007), high tumor neoantigen load (p = 0.029), and high tumor-associated CD8 T-cell exhaustion (p = 0.017). In addition, high field cancerization was associated with worse short-term outcomes (p = 0.029). Nonsynonymous mutations in bladder cancer-associated genes such as KDM6A, ARID1A, and TP53 were identified as early disease drivers already found in normal-appearing bladder biopsies. Urinary tumor DNA (utDNA) levels reflected the bladder tumor burden and originated from tumors and field cancerization. High levels of utDNA after BCG were associated with worse clinical outcomes (p = 0.027) and with disease progression (p = 0.003). High field cancerization resulted in high urinary levels of proteins associated with angiogenesis and proliferation. Limitations include variation in the number of biopsies and time points analyzed. CONCLUSIONS: Field cancerization levels are associated with tumor development, immune responses, and clinical outcomes. utDNA measurements can be used to monitor disease status and treatment response. PATIENT SUMMARY: Molecular changes in the tissue lining the bladder result in tumor recurrence. Urinary measurements may be used to monitor bladder cancer status and treatment responses.

Understanding the molecular landscape of nonmuscle-invasive bladder cancer (NMIBC) is essential to improve risk assessment and treatment regimens. We performed a comprehensive genomic analysis of patients with NMIBC using whole-exome sequencing (n = 438), shallow whole-genome sequencing (n = 362) and total RNA sequencing (n = 414). A large genomic variation within NMIBC was observed and correlated with different molecular subtypes. Frequent loss of heterozygosity in FGFR3 and 17p (affecting TP53) was found in tumors with mutations in FGFR3 and TP53, respectively. Whole-genome doubling (WGD) was observed in 15% of the tumors and was associated with worse outcomes. Tumors with WGD were genomically unstable, with alterations in cell-cycle-related genes and an altered immune composition. Finally, integrative clustering of multi-omics data highlighted the important role of genomic instability and immune cell exhaustion in disease aggressiveness. These findings advance our understanding of genomic differences associated with disease aggressiveness in NMIBC and may ultimately improve patient stratification.