A generalizable machine learning framework for classifying DNA repair defects using ctDNA exomes

Specific classes of DNA damage repair (DDR) defect can drive sensitivity to emerging therapies for metastatic prostate cancer. However, biomarker approaches based on DDR gene sequencing do not accurately predict DDR deficiency or treatment benefit. Somatic alteration signatures may identify DDR deficiency but historically require whole-genome sequencing of tumour tissue. We assembled whole-exome sequencing data for 155 high ctDNA fraction plasma cell-free DNA and matched leukocyte DNA samples from patients with metastatic prostate or bladder cancer. Labels for DDR gene alterations were established using deep targeted sequencing. Per sample mutation and copy number features were used to train XGBoost ensemble models. Naive somatic features and trinucleotide signatures were associated with specific DDR gene alterations but insufficient to resolve each class. Conversely, XGBoost-derived models showed strong performance including an area under the curve of 0.99, 0.99 and 1.00 for identifying BRCA2, CDK12, and mismatch repair deficiency in metastatic prostate cancer. Our machine learning approach re-classified several samples exhibiting genomic features inconsistent with original labels, identified a metastatic bladder cancer sample with a homozygous BRCA2 copy loss, and outperformed an existing exome-based classifier for BRCA2 deficiency. We present DARC Sign (DnA Repair Classification SIGNatures); a public machine learning tool leveraging clinically-practical liquid biopsy specimens for simultaneously identifying multiple types of metastatic prostate cancer DDR deficiencies. We posit that it will be useful for understanding differential responses to DDR-directed therapies in ongoing clinical trials and may ultimately enable prospective identification of prostate cancers with phenotypic evidence of DDR deficiency.

Somatic Features of Response and Relapse in Non-muscle-invasive Bladder Cancer Treated with Bacillus Calmette-Guérin Immunotherapy

BACKGROUND

High-risk non-muscle-invasive bladder cancer (NMIBC) is treated with bacillus Calmette-Guérin (BCG), but relapse is common. Improvement of patient outcomes requires better understanding of links between BCG resistance and genomic driver alterations.

OBJECTIVE

To validate the prognostic impact of common genomic alterations in NMIBC pretreatment and define somatic changes present in post-BCG relapses.

DESIGN, SETTING, AND PARTICIPANTS

We retrieved tumour tissues and outcomes for 90 patients with BCG-naive NMIBC initiating BCG monotherapy. Post-BCG tissue was available from 34 patients. All tissues underwent targeted sequencing of tumour and matched normal.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Associations between clinical outcomes and genomics were determined using Cox proportional hazard models.

RESULTS AND LIMITATIONS

Of the patients, 58% were relapse free at data cut-off, 24% had NMIBC recurrence, and 18% experienced muscle-invasive progression. The risk of relapse was associated with ARID1A mutation (hazard ratio [HR] = 2.00; p = 0.04) and CCNE1 amplification (HR = 2.61; p = 0.02). Pre- and post-BCG tumours shared truncal driver alterations, with mutations in TERT and chromatin remodelling genes particularly conserved. However, shifts in somatic profiles were common and clinically relevant alterations in FGFR3, PIK3CA, TSC1, and TP53 were temporally variable, despite apparent clonal prevalence at one time point. Limitations include the difficulty of resolving the relative impact of BCG therapy versus surgery on genomics at relapse and biopsy bias.

CONCLUSIONS

Somatic hypermutation and alterations in CCNE1 and ARID1A should be incorporated into future models predicting NMIBC BCG outcomes. Changes in tumour genomics over time highlight the importance of recent biopsy when considering targeted therapies, and suggest that relapse after BCG is due to persisting and evolving precursor populations.

PATIENT SUMMARY

Changes in key cancer genes can predict bladder cancer relapse after treatment with bacillus Calmette-Guérin. Relapses after treatment can be driven by large-scale genetic changes within the cancer. These genetic changes help us understand how superficial bladder cancer can progress to be treatment resistant.

TP53 Alterations Are Associated With Poor Survival in Patients With Primary Mediastinal Nonseminoma Germ Cell Tumors

BACKGROUND

Primary mediastinal nonseminoma germ cell tumors (PMNSGCT) are a subgroup of nonseminoma germ cell tumors (GCT) with poor prognosis. In this study, PMNSGCT-specific genomic landscape was analyzed and correlated with clinical outcomes.

METHODS

DNA was extracted and sequenced from 28 archival tumor tissue of patients with mediastinal GCT (3 seminoma and 25 nonseminoma). Overall survival (OS) and association with gene alterations were estimated using the Kaplan-Meier and univariate Cox regression methods.

RESULTS

Three patients (11%) had a karyotype XXY, 17/28 (61%) tumor samples presented chromosome 12p amplification. Somatic mutations were detected in 19/28 (68%) samples. The most frequently mutated genes were: TP53 (13/28; 46%), KIT (5/28; 18%), and KRAS (5/28; 18%). Deleterious TP53 alterations were associated with significantly reduced overall survival (HR: 7.16; P = .012).

CONCLUSIONS

TP53 alterations are common in PMNSGCT and are associated with reduced overall survival, potentially underlying the poor sensitivity to chemotherapy observed in these patients.

Deep whole-genome ctDNA chronology of treatment-resistant prostate cancer

Circulating tumour DNA (ctDNA) in blood plasma is an emerging tool for clinical cancer genotyping and longitudinal disease monitoring¹. However, owing to past emphasis on targeted and low-resolution profiling approaches, our understanding of the distinct populations that comprise bulk ctDNA is incomplete^2-12. Here we perform deep whole-genome sequencing of serial plasma and synchronous metastases in patients with aggressive prostate cancer. We comprehensively assess all classes of genomic alterations and show that ctDNA contains multiple dominant populations, the evolutionary histories of which frequently indicate whole-genome doubling and shifts in mutational processes. Although tissue and ctDNA showed concordant clonally expanded cancer driver alterations, most individual metastases contributed only a minor share of total ctDNA. By comparing serial ctDNA before and after clinical progression on potent inhibitors of the androgen receptor (AR) pathway, we reveal population restructuring converging solely on AR augmentation as the dominant genomic driver of acquired treatment resistance. Finally, we leverage nucleosome footprints in ctDNA to infer mRNA expression in synchronously biopsied metastases, including treatment-induced changes in AR transcription factor signalling activity. Our results provide insights into cancer biology and show that liquid biopsy can be used as a tool for comprehensive multi-omic discovery.

Cabazitaxel versus abiraterone or enzalutamide in poor prognosis metastatic castration-resistant prostate cancer: a multicentre, randomised, open-label, phase II trial

BACKGROUND

Treatment of poor prognosis metastatic castration-resistant prostate cancer (mCRPC) includes taxane chemotherapy and androgen receptor pathway inhibitors (ARPI). We sought to determine optimal treatment in this setting.

PATIENTS AND METHODS

This multicentre, randomised, open-label, phase II trial recruited patients with ARPI-naive mCRPC and poor prognosis features (presence of liver metastases, progression to mCRPC after <12 months of androgen deprivation therapy, or ≥4 of 6 clinical criteria). Patients were randomly assigned 1 : 1 to receive cabazitaxel plus prednisone (group A) or physician's choice of enzalutamide or abiraterone plus prednisone (group B) at standard doses. Patients could cross over at progression. The primary endpoint was clinical benefit rate for first-line treatment (defined as prostate-specific antigen response ≥50%, radiographic response, or stable disease ≥12 weeks).

RESULTS

Ninety-five patients were accrued (median follow-up 21.9 months). First-line clinical benefit rate was greater in group A versus group B (80% versus 62%, P = 0.039). Overall survival was not different between groups A and B (median 37.0 versus 15.5 months, hazard ratio (HR) = 0.58, P = 0.073) nor was time to progression (median 5.3 versus 2.8 months, HR = 0.87, P = 0.52). The most common first-line treatment-related grade ≥3 adverse events were neutropenia (cabazitaxel 32% versus ARPI 0%), diarrhoea (9% versus 0%), infection (9% versus 0%), and fatigue (7% versus 5%). Baseline circulating tumour DNA (ctDNA) fraction above the cohort median and on-treatment ctDNA increase were associated with shorter time to progression (HR = 2.38, P < 0.001; HR = 4.03, P < 0.001). Patients with >30% ctDNA fraction at baseline had markedly shorter overall survival than those with undetectable ctDNA (HR = 38.22, P < 0.001).

CONCLUSIONS

Cabazitaxel was associated with a higher clinical benefit rate in patients with ARPI-naive poor prognosis mCRPC. ctDNA abundance was prognostic independent of clinical features, and holds promise as a stratification biomarker.

Abstract A18: Plasma circulating tumor DNA is scarce and confounded by clonal hematopoiesis in metastatic renal cell carcinoma

Despite the recent implementation of novel targeted therapies, the 5-year survival rate for metastatic renal cell carcinoma (mRCC) remains dismal. Prognostic and predictive biomarkers are urgently required to avoid unnecessary patient morbidity and financial toxicity. Exploratory studies leveraging archival primary tumor tissue are considered suboptimal in this setting due to the considerable heterogeneity of metastatic lesions. Circulating tumor DNA (ctDNA) permits the noninvasive characterization of metastatic cancers using a simple blood draw. In this study, we sought to establish the utility of a ctDNA-based assay as a tool to profile the somatic genome of patients with mRCC. We collected whole blood from 55 progressing mRCC patients. All patients were systemic therapy naïve at the time of sample collection. Plasma cell-free DNA (cfDNA) and matched leukocyte DNA were subjected to targeted sequencing across 981 cancer-associated genes. Matched tumor tissue from 14 patients was also analyzed. The median cfDNA sequencing depth for this cohort was 938x. 33% of patients had evidence for RCC-derived ctDNA above 1% of total cfDNA; this was significantly lower than prostate or bladder cancer patients analyzed using the same approach. Among ctDNA-positive patients, ctDNA fraction averaged only 3.9% and showed no association with clinical variables or cfDNA yield. In these patients, the most commonly mutated genes were VHL, BAP1, and PBRM1, and matched tissue concordance was 77%. Evidence of somatic expansions unrelated to RCC, such as clonal hematopoiesis of indeterminate potential (CHIP), was detected in 43% of patients. Pathogenic germline mutations in DNA repair genes were detected in 11% of patients. Patients with ctDNA above 1% had shorter overall survival and progression-free survival on first-line therapy. Patients with evidence of CHIP but not ctDNA had an intermediate prognosis compared to ctDNA-positive and ctDNA-negative patients. CfDNA sequencing enables characterization of the somatic RCC genome in only a minority of metastatic RCC patients. Due to low ctDNA abundance and presence of non-RCC derived somatic clones in circulation, cfDNA sequencing may not be a simple pan-patient alternative to tissue biopsy in metastatic RCC.

Citation Format: Jack V. W. Bacon, Matti Annala, Maryam Soleimani, Jean-Michel Lavoie, Alan I. So, Martin E. Gleave, Ladan Fazli, Kim N. Chi, Christian K. Kollmannsberger, Alexander W. Wyatt, Lucia Nappi. Plasma circulating tumor DNA is scarce and confounded by clonal hematopoiesis in metastatic renal cell carcinoma [abstract]. In: Proceedings of the AACR Special Conference on Advances in Liquid Biopsies; Jan 13-16, 2020; Miami, FL. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(11_Suppl):Abstract nr A18.

Plasma Circulating Tumor DNA and Clonal Hematopoiesis in Metastatic Renal Cell Carcinoma

BACKGROUND

There is a lack of molecularly-informed biomarkers for patients with metastatic renal cell carcinoma (RCC). Plasma cell-free DNA (cfDNA) sequencing is a minimally-invasive alternative to tissue for profiling the genome in other cancers but relevance in metastatic RCC remains unclear.

MATERIALS AND METHODS

Whole blood was collected from 55 patients with metastatic RCC. Plasma cfDNA and leukocyte DNA were subjected to targeted sequencing across 981 cancer genes. Matched tumor tissue from 14 patients was analyzed.

RESULTS

Thirty-three percent of patients had evidence for RCC-derived circulating tumor DNA (ctDNA), significantly lower than patients with metastatic prostate or bladder cancer analyzed using the same approach. Among ctDNA-positive patients, ctDNA fraction averaged only 3.9% and showed no strong association with clinical variables. In these patients, the most commonly mutated genes were VHL, BAP1, and PBRM1, and matched tissue concordance was 77%. Evidence of somatic expansions unrelated to RCC, such as clonal hematopoiesis of indeterminate potential, were detected in 43% of patients. Pathogenic germline mutations in DNA repair genes were detected in 11% of patients. CtDNA-positive patients had shorter overall survival and progression-free survival on first-line therapy. Patients with evidence of clonal hematopoiesis of indeterminate potential had an intermediate prognosis compared with ctDNA-positive and -negative patients.

CONCLUSIONS

CfDNA sequencing enables straightforward characterization of the somatic RCC genome in a minority of patients with metastatic RCC. Owing to low ctDNA abundance, and the presence of non-RCC derived somatic clones in circulation, cfDNA sequencing may not be a simple pan-patient alternative to tissue biopsy in metastatic RCC.

Frequent mutation of the FOXA1 untranslated region in prostate cancer

Prostate cancer has a low somatic mutation rate but non-coding regions remain underexplored. We sequenced the untranslated regions (UTRs) of 72 established driver genes in 428 patients with metastatic prostate cancer and identified FOXA1 3'-UTR mutations in 12% of patients. The mutations were predominantly insertions or deletions, covered the entire UTR without motif enrichment, and were not detected in other cancers. FOXA1 lies in head-on orientation with the androgen-regulated non-coding gene AL121790.1, resulting in strong prostate lineage-specific bidirectional transcription across the FOXA1 3'-UTR. This suggests transcriptional activity as a cause for the localized hypermutation. The indel-dominant pattern of somatic mutation extends into the FOXA1 coding region, where it is shaped by clonal selection to yield a cluster of non-frameshift indels inside the forkhead domain. Somatic FOXA1 3'-UTR mutations may prove useful for diagnostic and screening approaches, given their high frequency and lineage specificity.