Pervasive chromosomal instability and karyotype order in tumour evolution.

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes1,2. The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution1,3,4. Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2+ breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.

Beyond basics: Key mutation selection features for successful tumor-informed ctDNA detection.

Tumor-informed mutation-based approaches are frequently used for detection of circulating tumor DNA (ctDNA). Not all mutations make equally effective ctDNA markers. The objective was to explore if prioritizing mutations using mutational features-such as cancer cell fraction (CCF), multiplicity, and error rate-would improve the success rate of tumor-informed ctDNA analysis. Additionally, we aimed to develop a practical and easily implementable analysis pipeline for identifying and prioritizing candidate mutations from whole-exome sequencing (WES) data. We analyzed WES and ctDNA data from three tumor-informed ctDNA studies, one on bladder cancer (Cohort A) and two on colorectal cancer (Cohorts I and N). The studies included 390 patients. For each patient, a unique set of mutations (median mutations/patient: 6, interquartile 13, range: 1-46, total n = 4023) were used as markers of ctDNA. The tool PureCN was used to assess the CCF and multiplicity of each mutation. High-CCF mutations were detected more frequently than low-CCF mutations (Cohort A: odds ratio [OR] 20.6, 95% confidence interval [CI] 5.72-173, p = 1.73e-12; Cohort I: OR 2.24, 95% CI 1.44-3.52, p = 1.66e-04; and Cohort N: OR 1.78, 95% CI 1.14-2.79, p = 7.86e-03). The detection-likelihood was additionally improved by selecting mutations with multiplicity of two or above (Cohort A: OR 1.55, 95% CI 1. 14-2.11, p = 3.85e-03; Cohort I: OR 1.78, 95% CI 1.23-2.56, p = 1.34e-03; and Cohort N: OR 1.94, 95% CI 1.63-2.31, p = 2.83e-14). Furthermore, selecting the mutations for which the ctDNA detection method had the lowest error rates, additionally improved the detection-likelihood, particularly evident when plasma cell-free DNA tumor fractions were below 0.1% (p = 2.1e-07). Selecting mutational markers with high CCF, high multiplicity, and low error rate significantly improve ctDNA detection likelihood. We provide free access to the analysis pipeline enabling others to perform qualified prioritization of mutations for tumor-informed ctDNA analysis.

Circulating tumor DNA for prognosis assessment and postoperative management after curative-intent resection of colorectal liver metastases.

The recurrence rate of colorectal liver metastases (CRLM) patients treated with curative intent is above 50%. Standard of care surveillance includes intensive computed tomographic (CT) imaging as well as carcinoembryonic antigen (CEA) measurements. Nonetheless, relapse detection often happens too late to resume curative treatment. This longitudinal cohort study enrolled 115 patients with plasma samples (N = 439) prospectively collected before surgery, postoperatively at day 30 and every third month for up to 3 years. Droplet digital PCR (ddPCR) was used to monitor serial plasma samples for somatic mutations. Assessment of ctDNA status either immediately after surgery, or serially during surveillance, stratified the patients into groups of high and low recurrence risk (hazard ratio [HR], 7.6; 95% CI, 3.0-19.7; P < .0001; and HR, 4.3; 95% CI, 2.3-8.1; P < .0001, respectively). The positive predictive value (PPV) of ctDNA was 100% in all postoperative analyses. In multivariable analyses, postoperative ctDNA status was the only consistently significant risk marker associated with relapse (P < .0001). Indeterminate CT findings were observed for 30.8% (21/68) of patients. All patients (9/21) that were ctDNA positive at the time of the indeterminate CT scan later relapsed, contrasting 42.6% (5/12) of those ctDNA negative (P = .0046). Recurrence diagnoses in patients with indeterminate CT findings were delayed (median 2.8 months, P < .0001). ctDNA status is strongly associated with detection of minimal residual disease and early detection of relapse. Furthermore, ctDNA status can potentially contribute to clinical decision-making in case of indeterminate CT findings, reducing time-to-intervention.

Epigenetic and transcriptomic consequences of excess X-chromosome material in 47,XXX syndrome-A comparison with Turner syndrome and 46,XX females.

47,XXX (triple X) and Turner syndrome (45,X) are sex chromosomal abnormalities with detrimental effects on health with increased mortality and morbidity. In karyotypical normal females, X-chromosome inactivation balances gene expression between sexes and upregulation of the X chromosome in both sexes maintain stoichiometry with the autosomes. In 47,XXX and Turner syndrome a gene dosage imbalance may ensue from increased or decreased expression from the genes that escape X inactivation, as well as from incomplete X chromosome inactivation in 47,XXX. We aim to study genome-wide DNA-methylation and RNA-expression changes can explain phenotypic traits in 47,XXX syndrome. We compare DNA-methylation and RNA-expression data derived from white blood cells of seven women with 47,XXX syndrome, with data from seven female controls, as well as with seven women with Turner syndrome (45,X). To address these questions, we explored genome-wide DNA-methylation and transcriptome data in blood from seven females with 47,XXX syndrome, seven females with Turner syndrome, and seven karyotypically normal females (46,XX). Based on promoter methylation, we describe a demethylation of six X-chromosomal genes (AMOT, HTR2C, IL1RAPL2, STAG2, TCEANC, ZNF673), increased methylation for GEMIN8, and four differentially methylated autosomal regions related to four genes (SPEG, MUC4, SP6, and ZNF492). We illustrate how these changes seem compensated at the transcriptome level although several genes show differential exon usage. In conclusion, our results suggest an impact of the supernumerary X chromosome in 47,XXX syndrome on the methylation status of selected genes despite an overall comparable expression profile.

Discordant molecular subtype classification in the basal-squamous subtype of bladder tumors and matched lymph-node metastases.

Molecular subtypes of muscle-invasive bladder tumors have emerged as a promising research tool with potential to stratify patients for neoadjuvant treatment. Prior to radical cystectomy, the utility of molecular classification and biomarkers depend on concordance between tissue from transurethrally resected specimens and disseminated disease. We assess the concordance of molecular subtypes and a large number of potential biomarkers in 67 pairs of muscle-invasive bladder tumors and synchronous lymph-node metastases. Tissue cores were stained for 29 immunohistochemistry markers and immunohistochemistry-based molecular subtype classification was performed. Molecular subtype was determined by mRNA profiling for 57 bladder tumors and 28 matched lymph-node metastases. Full section immunohistochemistry was performed to assess intra-tumor subtype heterogeneity in discordant cases, and exome sequencing was performed for 20 sample pairs. Discordant subtype classification between the bladder tumor and lymph-node metastasis was generally rare (12/67, 18%), but most (7/12, 58%) involved the Basal/Squamous-like subtype. Discordant Basal/Squamous-like tumors showed either Urothelial-like or Genomically Unstable, luminal-like phenotype in the lymph-node metastasis. Full section immunohistochemistry revealed intra-tumor subtype heterogeneity for six discordant cases including four involving the Basal/Squamous-like subtype. Subtype concordance for non- Basal/Squamous-like tumors was 91%. RNA-based classification agreed with immunohistochemistry classification but quantitative assessment is necessary to avoid false detection of subtype shifts. Most high confidence cancer mutations were shared between samples (n = 93, 78%), and bladder tumor private mutations (n = 20, 17%) were more frequent than those private to the lymph-node metastasis (n = 7, 6%). We conclude that bladder tumors and lymph-node metastases have overall similar molecular subtype, biomarker expression, and cancer mutations. The main exception was tumors of the Basal/Squamous-like subtype where most cases showed discordant classification, some with evidence of intra-tumor heterogeneity. The data are of relevance for neoadjuvant treatment stratification and raises questions on the dynamics of molecular subtypes during bladder cancer progression.

Analysis of circulating tumour DNA to monitor disease burden following colorectal cancer surgery.

OBJECTIVE: To develop an affordable and robust pipeline for selection of patient-specific somatic structural variants (SSVs) being informative about radicality of the primary resection, response to adjuvant therapy, incipient recurrence and response to treatment performed in relation to diagnosis of recurrence. DESIGN: We have established efficient procedures for identification of SSVs by next-generation sequencing and subsequent quantification of 3-6 SSVs in plasma. The consequence of intratumour heterogeneity on our approach was assessed. The level of circulating tumour DNA (ctDNA) was quantified in 151 serial plasma samples from six relapsing and five non-relapsing colorectal cancer (CRC) patients by droplet digital PCR, and correlated to clinical findings. RESULTS: Up to six personalised assays were designed for each patient. Our approach enabled efficient temporal assessment of disease status, response to surgical and oncological intervention, and early detection of incipient recurrence. Our approach provided 2-15 (mean 10) months' lead time on detection of metastatic recurrence compared to conventional follow-up. The sensitivity and specificity of the SSVs in terms of detecting postsurgery relapse were 100%. CONCLUSIONS: We show that assessment of ctDNA is a non-invasive, exquisitely specific and highly sensitive approach for monitoring disease load, which has the potential to provide clinically relevant lead times compared with conventional methods. Furthermore, we provide a low-coverage protocol optimised for identifying SSVs with excellent correlation between SSVs identified in tumours and matched metastases. Application of ctDNA analysis has the potential to change clinical practice in the management of CRC.

Identification of expressed and conserved human noncoding RNAs.

The past decade has shown mammalian genomes to be pervasively transcribed and identified thousands of noncoding (nc) transcripts. It is currently unclear to what extent these transcripts are of functional importance, as experimental functional evidence exists for only a small fraction. Here, we characterize the expression and evolutionary conservation properties of 12,115 known and novel nc transcripts, including structural RNAs, long nc RNAs (lncRNAs), antisense RNAs, EvoFold predictions, ultraconserved elements, and expressed nc regions. Expression levels are evaluated across 12 human tissues using a custom-designed microarray, supplemented with RNAseq. Conservation levels are evaluated at both the base level and at the syntenic level. We combine these measures with epigenetic mark annotations to identify subsets of novel nc transcripts that show characteristics similar to known functional ncRNAs. Few novel nc transcripts show both high expression and conservation levels. However, overall, we observe a positive correlation between expression and both conservation and epigenetic annotations, suggesting that a subset of the expressed transcripts are under purifying selection and likely functional. The identified subsets of expressed and conserved novel nc transcripts may form the basis for further functional characterization.

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.

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.

Whole-genome Mutational Analysis for Tumor-informed Detection of Circulating Tumor DNA in Patients with Urothelial Carcinoma.

BACKGROUND AND OBJECTIVE: Circulating tumor DNA (ctDNA) can be used for sensitive detection of minimal residual disease (MRD). However, the probability of detecting ctDNA in settings of low tumor burden is limited by the number of mutations analyzed and the plasma volume available. We used a whole-genome sequencing (WGS) approach for ctDNA detection in patients with urothelial carcinoma. METHODS: We used a tumor-informed WGS approach for ctDNA-based detection of MRD and evaluation of treatment responses. We analyzed 916 longitudinally collected plasma samples from 112 patients with localized muscle-invasive bladder cancer who received neoadjuvant chemotherapy (NAC) before radical cystectomy. Recurrence-free survival (primary endpoint), overall survival, and ctDNA dynamics during NAC were assessed. KEY FINDINGS AND LIMITATIONS: We found that WGS-based ctDNA detection is prognostic for patient outcomes with a median lead time of 131 d over radiographic imaging. WGS-based ctDNA assessment after radical cystectomy identified recurrence with sensitivity of 91% and specificity of 92%. In addition, genomic characterization of post-treatment plasma samples with a high ctDNA level revealed acquisition of platinum therapy-associated mutational signatures and copy number variations not present in the primary tumors. The sequencing depth is a limitation for studying tumor evolution. CONCLUSIONS AND CLINICAL IMPLICATIONS: Our results support the use of WGS for ultrasensitive ctDNA detection and highlight the possibility of plasma-based tracking of tumor evolution. WGS-based ctDNA detection represents a promising option for clinical use owing to the low volume of plasma needed and the ease of performing WGS, eliminating the need for personalized assay design. PATIENT SUMMARY: Detection of tumor DNA in blood samples from patients with cancer of the urinary tract is associated with poorer outcomes. Disease recurrence after surgery can be identified by the presence of tumor DNA in blood before it can be detected on radiography scans.

NGS-Based Tumor-Informed Analysis of Circulating Tumor DNA.

Accurate circulating tumor DNA (ctDNA) detection has an immense biomarker potential in all phases of the cancer disease course. Presence of ctDNA in the blood has been shown to have prognostic value in various cancer types as it may reflect the actual tumor burden. There are two main methods to consider, a tumor-informed and a tumor-agnostic analysis of ctDNA. Both techniques exploit the short half-life of circulating cell-free DNA (cfDNA)/ctDNA for disease monitoring and ultimately future clinical treatment intervention. Urothelial carcinoma is characterized by a high mutation spectrum but very few hotspot mutations. This limits tumor agnostic usability of hotspot mutation or fixed sets of genes for ctDNA detection. Here we focus on a tumor-informed analysis for ultrasensitive patient- and tumor-specific ctDNA detection using personalized mutation panels, probes that bind to specific genomic sequences to enrich for the region of interest. In this chapter, we describe methods for purification of high-quality cfDNA and guidelines for designing tumor-informed customized capture panels for sensitive detection of ctDNA. Furthermore, a detailed protocol for library preparation and panel capture utilizing a double enrichment strategy with low amplification is described.

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.

Expression of TIP60 (tat-interactive protein) and MRE11 (meiotic recombination 11 homolog) predict treatment-specific outcome of localised invasive bladder cancer.

UNLABELLED: What's known on the subject? and What does the study add? Several studies have shown that defects in DNA-damage response are associated with good survival after chemotherapy and radiotherapy. Furthermore, loss of cell cycle regulators may be prognostic indicators of poor survival after cystectomy. However, the potential clinical impact of previous findings is hampered by insufficient validation of significant results in suitable cystectomy and radiotherapy cohorts. Here we use a large cohort of patients receiving radiotherapy to successfully validate the importance of MRE11 as a predictive marker of disease-specific survival (DSS). Furthermore, using two independent patient cohorts we show for the first time that TIP60 is a predictive marker of DSS after cystectomy. We show that combined use of TIP60 and MRE11 may hold the potential to guide treatment decisions. OBJECTIVE: • To determine the association between the proteins: tat-interactive protein 60 kDa (TIP60), p16, meiotic recombination 11 homolog (MRE11), phosphorylated ataxia telangiectasia mutated (ATM), retinoblastoma protein (Rb), Ki67, and p53 and clinical outcome in invasive lymph node-negative bladder cancer. PATIENTS AND METHODS: • Protein expression was measured by immunohistochemistry in cancer specimens from two independent cohorts of patients with bladder cancer treated with cystectomy (162 patients and 273) and one cohort of patients receiving radiotherapy (148). • Disease-specific survival (DSS) was used as the outcome measure, and patients with no disease-specific death were followed for a minimum of 36 months. RESULTS: • TIP60 was significantly correlated with DSS in both cystectomy cohorts (hazard ratio [HR] 0.42, 95% confidence interval [CI] 0.26-0.68, P < 0.001 and HR 0.45, 95% CI 0.28-0.72, P = 0.001). • MRE11 was significantly correlated with DSS in the cohort receiving radiotherapy (HR 0.64, 95% CI 0.47-0.86, P = 0.005). • P16 was significantly correlated with DSS in all three cohorts (HR 0.46, 95% CI 0.30-0.75, P = 0.032; HR 0.60, 95% CI 0.37-0.97, P = 0.032; HR 0.52, 95% CI 0.28-0.96, P = 0.001). • Rb was significantly correlated with DSS in one cystectomy cohort (HR 1.71, 95% CI 1.13-2.75, P = 0.017). • Ki67, p53, and pATM were not significantly correlated with DSS in any of the cohorts. CONCLUSIONS: • TIP60 protein expression was a predictive marker for DSS after cystectomy in two independent cohorts. This novel marker was the strongest predictive factor in multivariate analysis in patients receiving cystectomy. • MRE11 was shown to be a predictive marker for DSS after radiotherapy. • We have shown that TIP60 and MRE11 hold the potential to guide patients with invasive bladder cancer to either cystectomy or radiotherapy. This study was based on retrospective material and consequently we suggest that these markers should be validated in a prospective study.

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.

Increased expression of transcription factor TFAP2α correlates with chemosensitivity in advanced bladder cancer.

BACKGROUND: The standard treatment for patients with advanced transitional cell carcinoma of the bladder is platin based chemotherapy. Only approximately 50% of the patients respond to chemotherapy. Therefore, molecular predictive markers for identification of chemotherapy sensitive subgroups of patients are highly needed. We selected the transcription factor TFAP2α from a previously identified gene expression signature for chemotherapy response. METHODS: TFAP2α expression and localization was assessed by immunohistochemistry using a tissue microarray (TMA) containing 282 bladder cancer tumors from patients with locally advanced (pT2-T4(b) and N(1-3)) or metastatic (M(1)) disease. All patients had received cisplatin containing chemotherapy. Furthermore, QPCR analysis of three TFAP2α isoforms was performed on tumor specimens of advanced muscle invasive bladder cancers (T2-4). Using the bladder cell lines T24 and SW780 the relation of TFAP2α and cisplatin and gemcitabine sensitivity as well as cell proliferation was examined using siRNA directed TFAP2α knockdown. RESULTS: TFAP2α protein expression was analyzed on a TMA with cores from 282 advanced bladder cancer tumors from patients treated with cisplatin based combinational chemotherapy. TFAP2α was identified as a strong independent predictive marker for a good response and survival after cisplatin-containing chemotherapy in patients with advanced bladder cancer. Strong TFAP2α nuclear and cytoplasmic staining predicted good response to chemotherapy in patients with lymph node metastasis, whereas weak TFAP2α nuclear staining predicted good response in patients without lymph node metastasis. In vitro studies showed that siRNA mediated knockdown of TFAP2α increased the proliferation of SW780 cells and rendered the cells less sensitive to cisplatin and gemcitabine. In contrast to that T24 bladder cells with mutated p53 showed to be more drug sensitive upon TFAP2α depletion. CONCLUSIONS: High levels of nuclear and cytoplasmic TFAP2α protein were a predictor of increased overall survival and progression free survival in patients with advanced bladder cancer treated with cisplatin based chemotherapy. TFAP2α knockdown increased the proliferation of the SW780 bladder cells and reduced cisplatin and gemcitabine induced cell death. The inverse effect was observed in the TP53 mutated T24 cell line where TFAP2α silencing augmented cisplatin and gemcitabine sensitivity and did not stimulate proliferation.

SPTAN1, APC, and FGFR3 Mutation Status and APOBEC Mutation Signatures are Predictive of Mitomycin C Response in Non-muscle-invasive Bladder Cancer.

BACKGROUND: Currently, no biomarkers of response to mitomycin C have been identified in non-muscle-invasive bladder cancer patients. Predictive biomarkers could improve the treatment outcome and eliminate adverse events from unnecessary treatment. OBJECTIVE: To identify and validate predictive biomarkers of chemoresection with mitomycin C. DESIGN SETTING AND PARTICIPANTS: The intervention group of a randomised controlled trial was identified for analyses. The study was conducted between January 2018 and June 2019 in two major urological departments in Denmark. Patients had a history of Ta low-grade/high-grade disease and were included upon recurrence. The intervention group (58 patients) received chemoresection with mitomycin C. Tumour and reference germline DNA from prior tumours were analysed by whole exome sequencing. Predictive biomarkers were validated in the context of Ta low-grade tumours from the UROMOL study. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Response to chemotherapy (intervention group from the randomised controlled trial) and recurrence-free survival (UROMOL cohort) were measured. Groups were compared using Fisher's exact test and Wilcoxon rank sum test. RESULTS AND LIMITATIONS: Chemoresponse was associated with the mutation status of SPTAN1, APC, and FGFR3, and the level of APOBEC signature contribution (p = 0.035, p = 0.034, p = 0.055, and p = 0.035, respectively). The main limitations include no biopsy for biomarker discovery immediately prior to chemoresection and the unmatched validation cohort. CONCLUSIONS: Mutation status of APC, SPTAN1, and FGFR3 and the level of mutational contribution from APOBEC-related signatures were identified as potential predictive biomarkers for chemoresection with mitomycin C in non-muscle-invasive bladder cancer patients. A prospective validation study is however needed. PATIENT SUMMARY: We investigated DNA from noninvasive bladder tumours in order to predict treatment response to chemotherapy. Four biomarkers showed promising results, which should be tested in future studies.

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.