Association of Age with Non-muscle-invasive Bladder Cancer: Unearthing a Biological Basis for Epidemiological Disparities?

BACKGROUND: Age disparity in patients with non-muscle-invasive bladder cancer (NMIBC) exists. Whether this is due to differences in adequate cancer care or tumour biology is unclear. OBJECTIVE: To investigate age disparities in NMIBC using the Surveillance, Epidemiology, and End Results (SEER)-Medicare and UROMOL datasets. DESIGN, SETTING, AND PARTICIPANTS: The SEER-Medicare data were used to identify patients with clinical stage Ta, Tis, and T1 NMIBC between 2005 and 2017 (n = 32 225). Using the UROMOL cohort (n = 834), age disparities across transcriptomic, genomic, and spatial proteomic domains were assessed. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: For the SEER-Medicare data, multivariable competing-risk regression was used to examine the association between age and recurrence, progression, and bladder cancer-specific mortality (BCSM). For the UROMOL cohort, multivariable general linear model and multinomial logistic regression were performed to evaluate the association between age and tumour biology. RESULTS AND LIMITATIONS: An analysis of the SEER-Medicare cohort revealed 5-yr recurrence rates of 55.2%, 57.4%, and 58.9%; 5-yr progression rates of 25.6%, 29.2%, and 36.9%; and 5-yr BCSM rates of 3.9%, 5.8%, and 11.8% in patients aged 66-70, 71-80, and ≥81 yr, respectively. After multivariable adjustment, age ≥81 yr was associated with a higher risk of recurrence (hazard ratio [HR] 1.07, 95% confidence interval [CI] 1.03-1.12; p = 0.001), progression (HR 1.32, p < 0.001), and BCSM (HR 2.58, p < 0.001). UROMOL2021 transcriptomic class 2a was most frequently observed in patients with advanced age (34.0% in ≥76 yr vs 21.6% in ≤65 yr; p = 0.004), a finding confirmed on multivariable analysis (risk ratio [RR] 3.86, p = 0.002). UROMOL2021 genomic class 3 was observed more frequently in patients aged ≥76 yr (4.9% vs 24.2%; p = 0.001). Limitations include the definitions used for recurrence and progression, which may lead to under- or overestimation of true rates. CONCLUSIONS: Among SEER-Medicare patients with NMIBC, advanced age is associated with inferior oncological outcomes. These results reflect age-related molecular biological differences observed across transcriptomic and genomic domains, providing further evidence that innate tumour biology contributes to observed disparities in NMIBC outcomes. PATIENT SUMMARY: Older patients with non-muscle-invasive bladder cancer have worse oncological outcomes than younger patients. Some of this age disparity may be due to differences in tumour biology.

Improved protocol for single-nucleus RNA-sequencing of frozen human bladder tumor biopsies.

This paper provides a laboratory workflow for single-nucleus RNA-sequencing (snRNA-seq) including a protocol for gentle nuclei isolation from fresh frozen tumor biopsies, making it possible to analyze biobanked material. To develop this protocol, we used non-frozen and frozen human bladder tumors and cell lines. We tested different lysis buffers (IgePal and Nuclei EZ) and incubation times in combination with different approaches for tissue and cell dissection: sectioning, semi-automated dissociation, manual dissociation with pestles, and semi-automated dissociation combined with manual dissociation with pestles. Our results showed that a combination of IgePal lysis buffer, tissue dissection by sectioning, and short incubation time was the best conditions for gentle nuclei isolation applicable for snRNA-seq, and we found limited confounding transcriptomic changes based on the isolation procedure. This protocol makes it possible to analyze biobanked material from patients with well-described clinical and histopathological information and known clinical outcomes with snRNA-seq.

Single-nucleus and Spatially Resolved Intratumor Subtype Heterogeneity in Bladder Cancer.

Background: Current bulk transcriptomic classification systems for bladder cancer do not consider the level of intratumor subtype heterogeneity. Objective: To investigate the extent and possible clinical impact of intratumor subtype heterogeneity across early and more advanced stages of bladder cancer. Design setting and participants: We performed single-nucleus RNA sequencing (RNA-seq) of 48 bladder tumors and additional spatial transcriptomics for four of these tumors. Total bulk RNA-seq and spatial proteomics data were available from the same tumors for comparison, along with detailed clinical follow-up of the patients. Outcome measurements and statistical analysis: The primary outcome was progression-free survival for non-muscle-invasive bladder cancer. Cox regression analysis, log-rank tests, Wilcoxon rank-sum tests, Spearman correlation, and Pearson correlation were used for statistical analysis. Results and limitations: We found that the tumors exhibited varying levels of intratumor subtype heterogeneity and that the level of subtype heterogeneity can be estimated from both single-nucleus and bulk RNA-seq data, with high concordance between the two. We found that a higher class 2a weight estimated from bulk RNA-seq data is associated with worse outcome for patients with molecular high-risk class 2a tumors. The sparsity of the data generated using the DroNc-seq sequencing protocol is a limitation. Conclusions: Our results indicate that discrete subtype assignments from bulk RNA-seq data may lack biological granularity and that continuous class scores may improve clinical risk stratification of patients with bladder cancer. Patient summary: We found that several molecular subtypes can exist within a single bladder tumor and that continuous subtype scores can be used to identify a subgroup of patients with poor outcomes. Use of these subtype scores may improve risk stratification for patients with bladder cancer, which can help in making decisions on treatment.

Cell-Free Urine and Plasma DNA Mutational Analysis Predicts Neoadjuvant Chemotherapy Response and Outcome in Patients with Muscle-Invasive Bladder Cancer.

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.

Circulating Tumor DNA Analysis in Advanced Urothelial Carcinoma: Insights from Biological Analysis and Extended Clinical Follow-up.

PURPOSE: To investigate whether circulating tumor DNA (ctDNA) assessment in patients with muscle-invasive bladder cancer predicts treatment response and provides early detection of metastatic disease. EXPERIMENTAL DESIGN: We present full follow-up results (median follow-up: 68 months) from a previously described cohort of 68 neoadjuvant chemotherapy (NAC)-treated patients who underwent longitudinal ctDNA testing (712 plasma samples). In addition, we performed ctDNA evaluation of 153 plasma samples collected before and after radical cystectomy (RC) in a separate cohort of 102 NAC-naïve patients (median follow-up: 72 months). Total RNA sequencing of tumors was performed to investigate biological characteristics of ctDNA shedding tumors. RESULTS: Assessment of ctDNA after RC identified metastatic relapse with a sensitivity of 94% and specificity of 98% using the expanded follow-up data for the NAC-treated patients. ctDNA dynamics during NAC was independently associated with patient outcomes when adjusted for pathologic downstaging (HR = 4.7; P = 0.029). For the NAC-naïve patients, ctDNA was a prognostic predictor before (HR = 3.4; P = 0.0005) and after RC (HR = 17.8; P = 0.0002). No statistically significant difference in recurrence-free survival for patients without detectable ctDNA at diagnosis was observed between the cohorts. Baseline ctDNA positivity was associated with the Basal/Squamous (Ba/Sq) subtype and enrichment of epithelial-to-mesenchymal transition and cell cycle-associated gene sets. CONCLUSIONS: ctDNA is prognostic in NAC-treated and NAC-naïve patients with more than 5 years follow-up and outperforms pathologic downstaging in predicting treatment efficacy. Patients without detectable ctDNA at diagnosis may benefit significantly less from NAC, but additional studies are needed.

Stage-stratified molecular profiling of non-muscle-invasive bladder cancer enhances biological, clinical, and therapeutic insight.

Understanding the molecular determinants that underpin the clinical heterogeneity of non-muscle-invasive bladder cancer (NMIBC) is essential for prognostication and therapy development. Stage T1 disease in particular presents a high risk of progression and requires improved understanding. We present a detailed multi-omics study containing gene expression, copy number, and mutational profiles that show relationships to immune infiltration, disease recurrence, and progression to muscle invasion. We compare expression and genomic subtypes derived from all NMIBCs with those derived from the individual disease stages Ta and T1. We show that sufficient molecular heterogeneity exists within the separate stages to allow subclassification and that this is more clinically meaningful for stage T1 disease than that derived from all NMIBCs. This provides improved biological understanding and identifies subtypes of T1 tumors that may benefit from chemo- or immunotherapy.