Multiomics profiling of urothelial carcinoma in situ reveals CIS-specific gene signature and immune characteristics.

Urothelial carcinoma in situ (CIS) is an aggressive phenotype of non-muscle-invasive bladder cancer. Molecular features unique to CIS compared to high-grade papillary tumors are underexplored. RNA sequencing of CIS, papillary tumors, and normal urothelium showed lower immune marker expression in CIS compared to papillary tumors. We identified a 46-gene expression signature in CIS samples including selectively upregulated known druggable targets MTOR, TYK2, AXIN1, CPT1B, GAK, and PIEZO1 and selectively downregulated BRD2 and NDUFB2. High expression of selected genes was significantly associated with CIS in an independent dataset. Mutation analysis of matched CIS and papillary tumors revealed shared mutations between samples across time points and mutational heterogeneity. CCDC138 was the most frequently mutated gene in CIS. The immunological landscape showed higher levels of PD-1-positive cells in CIS lesions compared to papillary tumors. We identified CIS lesions to have distinct characteristics compared to papillary tumors potentially contributing to the aggressive phenotype.

Field Cancerization Is Associated with Tumor Development, T-cell Exhaustion, and Clinical Outcomes in 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.

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.

Elevated T-cell Exhaustion and Urinary Tumor DNA Levels Are Associated with Bacillus Calmette-Guérin Failure in Patients with Non-muscle-invasive Bladder Cancer.

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.

Refinement of an Established Procedure and Its Application for Identification of Hypoxia in Prostate Cancer Xenografts.

BACKGROUND: This pre-clinical study was designed to refine a dissection method for validating the use of a 15-gene hypoxia classifier, which was previously established for head and neck squamous cell carcinoma (HNSCC) patients, to identify hypoxia in prostate cancer. METHODS: PC3 and DU-145 adenocarcinoma cells, in vitro, were gassed with various oxygen concentrations (0-21%) for 24 h, followed by real-time PCR. Xenografts were established in vivo, and the mice were injected with the hypoxic markers [18F]-FAZA and pimonidazole. Subsequently, tumors were excised, frozen, cryo-sectioned, and analyzed using autoradiography ([18F]-FAZA) and immunohistochemistry (pimonidazole); the autoradiograms used as templates for laser capture microdissection of hypoxic and non-hypoxic areas, which were lysed, and real-time PCR was performed. RESULTS: In vitro, all 15 genes were increasingly up-regulated as oxygen concentrations decreased. With the xenografts, all 15 genes were up-regulated in the hypoxic compared to non-hypoxic areas for both cell lines, although this effect was greater in the DU-145. CONCLUSIONS: We have developed a combined autoradiographic/laser-guided microdissection method with broad applicability. Using this approach on fresh frozen tumor material, thereby minimizing the degree of RNA degradation, we showed that the 15-gene hypoxia gene classifier developed in HNSCC may be applicable for adenocarcinomas such as prostate cancer.

An integrated multi-omics analysis identifies prognostic molecular subtypes of non-muscle-invasive bladder cancer.

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.

Maternally Contributed Folate Receptor 1 Is Expressed in Ovarian Follicles and Contributes to Preimplantation Development.

Folates have been shown to play a crucial role for proper development of the embryo as folate deficiency has been associated with reduced developmental capacity such as increased risk of fetal neural tube defects and spontanous abortion. Transcripts encoding the reduced folate carrier RFC1 (SLC19A1 protein) and the high-affinity folate receptor FOLR1 are expressed in oocytes and preimplantation embryos, respectively. In this study, we observed maternally contributed FOLR1 protein during mouse and human ovarian follicle development, and 2-cell mouse embryos. In mice, FOLR1 was highly enriched in oocytes from primary, secondary and tertiary follicles, and in the surrounding granulosa cells. Interestingly, during human follicle development, we noted a high and specific presence of FOLR1 in oocytes from primary and intermediate follicles, but not in the granulosa cells. The distribution of FOLR1 in follicles was noted as membrane-enriched but also seen in the cytoplasm in oocytes and granulosa cells. In 2-cell embryos, FOLR1-eGFP fusion protein was detected as cytoplasmic and membrane-associated dense structures, resembling the distribution pattern observed in ovarian follicle development. Knock-down of Folr1 mRNA function was accomplished by microinjection of short interference (si)RNA targeting Folr1, into mouse pronuclear zygotes. This revealed a reduced capacity of Folr1 siRNA-treated embryos to develop to blastocyst compared to the siRNA-scrambled control group, indicating that maternally contributed protein and zygotic transcripts sustain embryonic development combined. In summary, maternally contributed FOLR1 protein appears to maintain ovarian functions, and contribute to preimplantation development combined with embryonically synthesized FOLR1.