Robust BRCA1-like classification of copy number profiles of samples repeated across different datasets and platforms

Long-term outcomes of young, node-negative, chemotherapy-naïve, triple-negative breast cancer patients according to BRCA1 status

BACKGROUND

Due to the abundant usage of chemotherapy in young triple-negative breast cancer (TNBC) patients, the unbiased prognostic value of BRCA1-related biomarkers in this population remains unclear. In addition, whether BRCA1-related biomarkers modify the well-established prognostic value of stromal tumor-infiltrating lymphocytes (sTILs) is unknown. This study aimed to compare the outcomes of young, node-negative, chemotherapy-naïve TNBC patients according to BRCA1 status, taking sTILs into account.

METHODS

We included 485 Dutch women diagnosed with node-negative TNBC under age 40 between 1989 and 2000. During this period, these women were considered low-risk and did not receive chemotherapy. BRCA1 status, including pathogenic germline BRCA1 mutation (gBRCA1m), somatic BRCA1 mutation (sBRCA1m), and tumor BRCA1 promoter methylation (BRCA1-PM), was assessed using DNA from formalin-fixed paraffin-embedded tissue. sTILs were assessed according to the international guideline. Patients' outcomes were compared using Cox regression and competing risk models.

RESULTS

Among the 399 patients with BRCA1 status, 26.3% had a gBRCA1m, 5.3% had a sBRCA1m, 36.6% had tumor BRCA1-PM, and 31.8% had BRCA1-non-altered tumors. Compared to BRCA1-non-alteration, gBRCA1m was associated with worse overall survival (OS) from the fourth year after diagnosis (adjusted HR, 2.11; 95% CI, 1.18-3.75), and this association attenuated after adjustment for second primary tumors. Every 10% sTIL increment was associated with 16% higher OS (adjusted HR, 0.84; 95% CI, 0.78-0.90) in gBRCA1m, sBRCA1m, or BRCA1-non-altered patients and 31% higher OS in tumor BRCA1-PM patients. Among the 66 patients with tumor BRCA1-PM and ≥ 50% sTILs, we observed excellent 15-year OS (97.0%; 95% CI, 92.9-100%). Conversely, among the 61 patients with gBRCA1m and < 50% sTILs, we observed poor 15-year OS (50.8%; 95% CI, 39.7-65.0%). Furthermore, gBRCA1m was associated with higher (adjusted subdistribution HR, 4.04; 95% CI, 2.29-7.13) and tumor BRCA1-PM with lower (adjusted subdistribution HR, 0.42; 95% CI, 0.19-0.95) incidence of second primary tumors, compared to BRCA1-non-alteration.

CONCLUSIONS

Although both gBRCA1m and tumor BRCA1-PM alter BRCA1 gene transcription, they are associated with different outcomes in young, node-negative, chemotherapy-naïve TNBC patients. By combining sTILs and BRCA1 status for risk classification, we were able to identify potential subgroups in this population to intensify and optimize adjuvant treatment.

Identifying the BRCA1 c.-107A > T variant in Dutch patients with a tumor BRCA1 promoter hypermethylation

An inherited single nucleotide variant (SNV) in the 5'UTR of the BRCA1 gene c.-107A > T was identified to be related to BRCA1 promoter hypermethylation and a hereditary breast and ovarian cancer phenotype in two UK families. We investigated whether this BRCA1 variant was also present in a Dutch cohort of breast and ovarian cancer patients with tumor BRCA1 promoter hypermethylation. We selected all breast and ovarian cancer cases that tested positive for tumor BRCA1 promoter hypermethylation at the Netherlands Cancer Institute and Sanger sequenced the specific mutation in the tumor DNA. In total, we identified 193 tumors with BRCA1 promoter hypermethylation in 178 unique patients. The wild-type allele was identified in 100% (193/193) of sequenced tumor samples. In a large cohort of 178 patients, none had tumors harboring the previously identified c.-107A > T SNV in BRCA1. We therefore can conclude that the germline SNV is not pervasive in patients with tumor BRCA1 promoter hypermethylation.

Chromosomal Instability, Selection and Competition: Factors That Shape the Level of Karyotype Intra-Tumor Heterogeneity

Simple Summary

Each cancer consists of billions of cells. These cells are far from identical; hence, the population of cells that constitute a tumor is heterogeneous. A salient property that varies between cells in a tumor is their karyotype, the number and configuration of the chromosomes. The level of karyotype heterogeneity can be used to predict the survival of a patient. In this review, we describe the processes that shape the level of karyotype heterogeneity in a cancer.

Abstract

Intra-tumor heterogeneity (ITH) is a pan-cancer predictor of survival, with high ITH being correlated to a dismal prognosis. The level of ITH is, hence, a clinically relevant characteristic of a malignancy. ITH of karyotypes is driven by chromosomal instability (CIN). However, not all new karyotypes generated by CIN are viable or competitive, which limits the amount of ITH. Here, we review the cellular processes and ecological properties that determine karyotype ITH. We propose a framework to understand karyotype ITH, in which cells with new karyotypes emerge through CIN, are selected by cell intrinsic and cell extrinsic selective pressures, and propagate through a cancer in competition with other malignant cells. We further discuss how CIN modulates the cell phenotype and immune microenvironment, and the implications this has for the subsequent selection of karyotypes. Together, we aim to provide a comprehensive overview of the biological processes that shape the level of karyotype heterogeneity.

Functional RECAP (REpair CAPacity) assay identifies homologous recombination deficiency undetected by DNA-based BRCAness tests

Germline BRCA1/2 mutation status is predictive for response to Poly-[ADP-Ribose]-Polymerase (PARP) inhibitors in breast cancer (BC) patients. However, non-germline BRCA1/2 mutated and homologous recombination repair deficient (HRD) tumors are likely also PARP-inhibitor sensitive. Clinical validity and utility of various HRD biomarkers are under investigation. The REpair CAPacity (RECAP) test is a functional method to select HRD tumors based on their inability to form RAD51 foci. We investigated whether this functional test defines a similar group of HRD tumors as DNA-based tests. An HRD enriched cohort (n = 71; 52 primary and 19 metastatic BCs) selected based on the RECAP test (26 RECAP-HRD; 37%), was subjected to DNA-based HRD tests (i.e., Classifier of HOmologous Recombination Deficiency (CHORD) and BRCA1/2-like classifier). Whole genome sequencing (WGS) was carried out for 38 primary and 19 metastatic BCs. The RECAP test identified all bi-allelic BRCA deficient samples (n = 15) in this cohort. RECAP status partially correlated with DNA-based HRD test outcomes (70% concordance for both RECAP-CHORD and RECAP-BRCA1/2-like classifier). RECAP selected additional samples unable to form RAD51 foci, suggesting that this functional assay identified deficiencies in other DNA repair genes, which could also result in PARP-inhibitor sensitivity. Direct comparison of these HRD tests in clinical trials will be required to evaluate the optimal predictive test for clinical decision making.

Effect of HIPEC according to HRD/BRCAwt genomic profile in stage III ovarian cancer - results from the phase III OVHIPEC trial

The addition of hyperthermic intraperitoneal chemotherapy (HIPEC) with cisplatin to interval cytoreductive surgery improves recurrence-free (RFS) and overall survival (OS) in patients with stage III ovarian cancer. Homologous recombination deficient (HRD) ovarian tumors are usually more platinum sensitive. Since hyperthermia impairs BRCA1/2 protein function, we hypothesized that HRD tumors respond best to treatment with HIPEC. We analyzed the effect of HIPEC in patients in the OVHIPEC trial, stratified by HRD status and BRCAm status. Clinical data and tissue samples were collected from patients included in the randomized, phase III OVHIPEC-1 trial. DNA copy number variation (CNV) profiles, HRD-related pathogenic mutations, and BRCA1 promotor hypermethylation were determined. CNV-profiles were categorized as HRD or non-HRD, based on a previously validated algorithm-based BRCA1-like classifier. Hazard ratios (HR) and corresponding 99% confidence intervals (CI) for the effect of RFS and OS of HIPEC in the BRCAm, the HRD/BRCAwt and the non-HRD group were estimated using Cox proportional hazard models. DNA was available from 200/245 (82%) patients. Seventeen (9%) tumors carried a pathogenic mutation in BRCA1 and 14 (7%) in BRCA2. Ninety-one (46%) tumors classified as BRCA1-like. The effect of HIPEC on RFS and OS was absent in BRCAm tumors (HR 1.25; 99%CI 0.48-3.29), and most present in HRD/BRCAwt (HR 0.44; 99%CI 0.21-0.91), and non-HRD/BRCAwt tumors (HR 0.82; 99%CI 0.48-1.42), interaction p-value: 0.024. Patients with HRD tumors without pathogenic BRCA1/2 mutation appear to benefit most from treatment with HIPEC, while benefit in patients with BRCA1/2 pathogenic mutations and patients without HRD seems less evident.

Adjuvant capecitabine-containing chemotherapy benefit and homologous recombination deficiency in early-stage triple-negative breast cancer patients

Background

The addition of adjuvant capecitabine to standard chemotherapy of early-stage triple-negative breast cancer (TNBC) patients has improved survival in a few randomised trials and in meta-analyses. However, many patients did not benefit. We evaluated the BRCA1-like DNA copy number signature, indicative of homologous recombination deficiency, as a predictive biomarker for capecitabine benefit in the TNBC subgroup of the FinXX trial.

Methods

Early-stage TNBC patients were randomised between adjuvant capecitabine-containing (TX + CEX: capecitabine-docetaxel, followed by cyclophosphamide-epirubicin-capecitabine) and conventional chemotherapy (T + CEF: docetaxel, followed by cyclophosphamide-epirubicin-fluorouracil). Tumour BRCA1-like status was determined on low-coverage, whole genome next-generation sequencing data using an established DNA comparative genomic hybridisation algorithm.

Results

For 129/202 (63.9%) patients the BRCA1-like status could be determined, mostly due to lack of tissue. During a median follow-up of 10.7 years, 35 recurrences and 32 deaths occurred. Addition of capecitabine appears to improve recurrence-free survival more among 61 (47.3%) patients with non-BRCA1-like tumours (HR 0.23, 95% CI 0.08–0.70) compared to 68 (52.7%) patients with BRCA1-like tumours (HR 0.66, 95% CI 0.24–1.81) (P-interaction = 0.17).

Conclusion

Based on our data, patients with non-BRCA1-like TNBC appear to benefit from the addition of capecitabine to adjuvant chemotherapy. Patients with BRCA1-like TNBC may also benefit. Additional research is needed to define the subgroup within BRCA1-like TNBC patients who may not benefit from adjuvant capecitabine.

PARP Inhibitors: A Major Therapeutic Option in Endocrine-Receptor Positive Breast Cancers

Simple Summary

OlympiAD and EMBRACA trials demonstrated the efficacy of PARPi, compared to chemotherapy, in patients with HER2-negative metastatic breast cancers (mBC) carrying a germline BRCA mutation. Patients with ER+/HER2-BRCA-mutated mBC seemed to have a higher risk of early disease progression while on CDK4/6 inhibitors and benefit from PARPi, especially when prescribed before chemotherapy. Importantly, the frequency of BRCA pathogenic variant (PV) carriers among ER+/HER2- breast cancer patients has been underestimated, and 50% of all BRCA1/2 mutated breast cancers are actually of ER+/HER2- subtype. Recent studies also highlight the benefit of PARPi in BRCA wild type mBC with HRD representing up to 20% of ER+/HER2- breast cancers. The OLYMPIA trial also demonstrated PARPi utility in patients with ER+/HER2- early breast cancers with BRCA PV at high risk of relapse. Consequently, implementation of early genotyping and new strategies for identifying patients with high-risk ER+/HER2- HRD breast cancers likely to benefit from PARPi is of high importance.

Abstract

Recently, OlympiAD and EMBRACA trials demonstrated the favorable efficacy/toxicity ratio of PARPi, compared to chemotherapy, in patients with HER2-negative metastatic breast cancers (mBC) carrying a germline BRCA mutation. PARPi have been largely adopted in triple-negative metastatic breast cancer, but their place has been less clearly defined in endocrine-receptor positive, HER2 negative (ER+/ HER2-) mBC. The present narrative review aims at addressing this question by identifying the patients that are more likely benefit from PARPi. Frequencies of BRCA pathogenic variant (PV) carriers among ER+/HER2- breast cancer patients have been underestimated, and many experts assume than 50% of all BRCA1/2 mutated breast cancers are of ER+/HER2- subtype. Patients with ER+/HER2- BRCA-mutated mBC seemed to have a higher risk of early disease progression while on CDK4/6 inhibitors and PARPi are effective especially when prescribed before exposure to chemotherapy. The OLYMPIA trial also highlighted the utility of PARPi in patients with early breast cancers at high risk of relapse and carrying PV of BRCA. PARPi might also be effective in patients with HRD diseases, representing up to 20% of ER+/HER2- breast cancers. Consequently, the future implementation of early genotyping strategies for identifying the patients with high-risk ER+/HER2- HRD breast cancers likely to benefit from PARPi is of high importance.

Ovarian Cancer-Specific BRCA-like Copy-Number Aberration Classifiers Detect Mutations Associated with Homologous Recombination Deficiency in the AGO-TR1 Trial

PURPOSE

Previously, we developed breast cancer BRCA1-like and BRCA2-like copy-number profile shrunken centroid classifiers predictive for mutation status and response to therapy, targeting homologous recombination deficiency (HRD). Therefore, we investigated BRCA1- and BRCA2-like classification in ovarian cancer, aiming to acquire classifiers with similar properties as those in breast cancer.Experimental Design: We analyzed DNA copy-number profiles of germline BRCA1- and BRCA2-mutant ovarian cancers and control tumors and observed that existing breast cancer classifiers did not sufficiently predict mutation status. Hence, we trained new shrunken centroid classifiers on this set and validated them in the independent The Cancer Genome Atlas dataset. Subsequently, we assessed BRCA1/2-like classification and obtained germline and tumor mutation and methylation status of cancer predisposition genes, among them several involved in HR repair, of 300 ovarian cancer samples derived from the consecutive cohort trial AGO-TR1 (NCT02222883).

RESULTS

The detection rate of the BRCA1-like classifier for BRCA1 mutations and promoter hypermethylation was 95.6%. The BRCA2-like classifier performed less accurately, likely due to a smaller training set. Furthermore, three quarters of the BRCA1/2-like tumors could be explained by (epi)genetic alterations in BRCA1/2, germline RAD51C mutations and alterations in other genes involved in HR. Around half of the non-BRCA-mutated ovarian cancer cases displayed a BRCA-like phenotype.

CONCLUSIONS

The newly trained classifiers detected most BRCA-mutated and methylated cancers and all tumors harboring a RAD51C germline mutations. Beyond that, we found an additional substantial proportion of ovarian cancers to be BRCA-like.

A survey of cancer genome signatures identifies genes connected to distinct chromosomal instability phenotypes

Certain breast and ovarian cancers are characterised by high levels of chromosomal instability. We established a suite of eleven SNP array-based signatures of various forms of chromosomal instability. To understand what biological mechanisms might underpin these signatures, we generated and assembled genetic-feature data including allele-specific expression, fusion genes, gene expression, methylation, somatic coding mutations and protein expression. For each signature, we extracted a compendium of significantly associated genetic features using machine learning. We established an association between subchromosomal allelic imbalance-based measures and DNA repair genes. Numerical chromosomal instability and chromothripsis were found to have distinct genetic associations but shared a relationship to mitotic processes, while whole-genome doubling was characterised by TP53 mutation, and high AURKA and GINS1 expression. Furthermore, we identified two genetically distinct forms of tandem duplicator phenotypes. These findings identify potentially novel genomic targets for validation and drug development for specific subsets of breast and ovarian cancer.

Assessment of structural chromosomal instability phenotypes as biomarkers of carboplatin response in triple negative breast cancer: the TNT trial

BACKGROUND

In the TNT trial of triple negative breast cancer (NCT00532727), germline BRCA1/2 mutations were present in 28% of carboplatin responders. We assessed quantitative measures of structural chromosomal instability (CIN) to identify a wider patient subgroup within TNT with preferential benefit from carboplatin over docetaxel.

PATIENTS AND METHODS

Copy number aberrations (CNAs) were established from 135 formalin-fixed paraffin-embedded primary carcinomas using Illumina OmniExpress SNP-arrays. Seven published [allelic imbalanced CNA (AiCNA); allelic balanced CNA (AbCNA); copy number neutral loss of heterozygosity (CnLOH); number of telomeric allelic imbalances (NtAI); BRCA1-like status; percentage of genome altered (PGA); homologous recombination deficiency (HRD) scores] and two novel [Shannon diversity index (SI); high-level amplifications (HLAMP)] CIN-measurements were derived. HLAMP was defined based on the presence of at least one of the top 5% amplified cytobands located on 1q, 8q and 10p. Continuous CIN-measurements were divided into tertiles. All nine CIN-measurements were used to analyse objective response rate (ORR) and progression-free survival (PFS).

RESULTS

Patients with tumours without HLAMP had a numerically higher ORR and significantly longer PFS in the carboplatin (C) than in the docetaxel (D) arm [56% (C) versus 29% (D), PHLAMP,quiet = 0.085; PFS 6.1 months (C) versus 4.1 months (D), Pinteraction/HLAMP = 0.047]. In the carboplatin arm, patients with tumours showing intermediate telomeric NtAI and AiCNA had higher ORR [54% (C) versus 20% (D), PNtAI,intermediate = 0.03; 62% (C) versus 33% (D), PAiCNA,intermediate = 0.076]. Patients with high AiCNA and PGA had shorter PFS in the carboplatin arm [3.4 months (high) versus 5.7 months (low/intermediate); and 3.8 months (high) versus 5.6 months (low/intermediate), respectively; Pinteraction/AiCNA = 0.027, Padj.interaction/AiCNA = 0.125 and Pinteraction/PGA = 0.053, Padj.interaction/PGA = 0.176], whilst no difference was observed in the docetaxel arm.

CONCLUSIONS

Patients with tumours lacking HLAMP and demonstrating intermediate CIN-measurements formed a subgroup benefitting from carboplatin relative to docetaxel treatment within the TNT trial. This suggests a complex and paradoxical relationship between the extent of genomic instability in primary tumours and treatment response in the metastatic setting.

An immune-centric exploration of BRCA1 and BRCA2 germline mutation related breast and ovarian cancers

Background

BRCA1/2 germline mutation related cancers are candidates for new immune therapeutic interventions. This study was a hypothesis generating exploration of genomic data collected at diagnosis for 19 patients. The prominent tumor mutation burden (TMB) in hereditary breast and ovarian cancers in this cohort was not correlated with high global immune activity in their microenvironments. More information is needed about the relationship between genomic instability, phenotypes and immune microenvironments of these hereditary tumors in order to find appropriate markers of immune activity and the most effective anticancer immune strategies.

Methods

Mining and statistical analyses of the original DNA and RNA sequencing data and The Cancer Genome Atlas data were performed. To interpret the data, we have used published literature and web available resources such as Gene Ontology, The Cancer immunome Atlas and the Cancer Research Institute iAtlas.

Results

We found that BRCA1/2 germline related breast and ovarian cancers do not represent a unique phenotypic identity, but they express a range of phenotypes similar to sporadic cancers. All breast and ovarian BRCA1/2 related tumors are characterized by high homologous recombination deficiency (HRD) and low aneuploidy. Interestingly, all sporadic high grade serous ovarian cancers (HGSOC) and most of the subtypes of triple negative breast cancers (TNBC) also express a high degree of HRD.

Conclusions

TMB is not associated with the magnitude of the immune response in hereditary BRCA1/2 related breast and ovarian cancers or in sporadic TNBC and sporadic HGSOC. Hereditary tumors express phenotypes as heterogenous as sporadic tumors with various degree of “BRCAness” and various characteristics of the immune microenvironments. The subtyping criteria developed for sporadic tumors can be applied for the classification of hereditary tumors and possibly also characterization of their immune microenvironment. A high HRD score may be a good candidate biomarker for response to platinum, and potentially PARP-inhibition.

Trial registration

Phase I Study of the Oral PI3kinase Inhibitor BKM120 or BYL719 and the Oral PARP Inhibitor Olaparib in Patients With Recurrent TNBC or HGSOC (NCT01623349), first posted on June 20, 2012. The design and the outcome of the clinical trial is not in the scope of this study.

EZH2 Is Overexpressed in BRCA1-like Breast Tumors and Predictive for Sensitivity to High-Dose Platinum-Based Chemotherapy

PURPOSE

BRCA1-deficient breast cancers carry a specific DNA copy-number signature ("BRCA1-like") and are hypersensitive to DNA double-strand break (DSB) inducing compounds. Here, we explored whether (i) EZH2 is overexpressed in human BRCA1-deficient breast tumors and might predict sensitivity to DSB-inducing drugs; (ii) EZH2 inhibition potentiates cisplatin efficacy in Brca1-deficient murine mammary tumors.

EXPERIMENTAL DESIGN

EZH2 expression was analyzed in 497 breast cancers using IHC or RNA sequencing. We classified 370 tumors by copy-number profiles as BRCA1-like or non-BRCA1-like and examined its association with EZH2 expression. Additionally, we assessed BRCA1 loss through mutation or promoter methylation status and investigated the predictive value of EZH2 expression in a study population of breast cancer patients treated with adjuvant high-dose platinum-based chemotherapy compared with standard anthracycline-based chemotherapy. To explore whether EZH2 inhibition by GSK126 enhances sensitivity to platinum drugs in EZH2-overexpressing breast cancers we used a Brca1-deficient mouse model.

RESULTS

The highest EZH2 expression was found in BRCA1-associated tumors harboring a BRCA1 mutation, BRCA1-promoter methylation or were classified as BRCA1 like. We observed a greater benefit from high-dose platinum-based chemotherapy in BRCA1-like and non-BRCA1-like patients with high EZH2 expression. Combined treatment with the EZH2 inhibitor GSK126 and cisplatin decreased cell proliferation and improved survival in Brca1-deficient mice in comparison with single agents.

CONCLUSIONS

Our findings demonstrate that EZH2 is expressed at significantly higher levels in BRCA1-deficient breast cancers. EZH2 overexpression can identify patients with breast cancer who benefit significantly from intensified DSB-inducing platinum-based chemotherapy independent of BRCA1-like status. EZH2 inhibition improves the antitumor effect of platinum drugs in Brca1-deficient breast tumors in vivo.

Immune induction strategies in metastatic triple-negative breast cancer to enhance the sensitivity to PD-1 blockade: the TONIC trial

The efficacy of programmed cell death protein 1 (PD-1) blockade in metastatic triple-negative breast cancer (TNBC) is low^1-5, highlighting a need for strategies that render the tumor microenvironment more sensitive to PD-1 blockade. Preclinical research has suggested immunomodulatory properties for chemotherapy and irradiation^6-13. In the first stage of this adaptive, non-comparative phase 2 trial, 67 patients with metastatic TNBC were randomized to nivolumab (1) without induction or with 2-week low-dose induction, or with (2) irradiation (3 × 8 Gy), (3) cyclophosphamide, (4) cisplatin or (5) doxorubicin, all followed by nivolumab. In the overall cohort, the objective response rate (ORR; iRECIST¹⁴) was 20%. The majority of responses were observed in the cisplatin (ORR 23%) and doxorubicin (ORR 35%) cohorts. After doxorubicin and cisplatin induction, we detected an upregulation of immune-related genes involved in PD-1-PD-L1 (programmed death ligand 1) and T cell cytotoxicity pathways. This was further supported by enrichment among upregulated genes related to inflammation, JAK-STAT and TNF-α signaling after doxorubicin. Together, the clinical and translational data of this study indicate that short-term doxorubicin and cisplatin may induce a more favorable tumor microenvironment and increase the likelihood of response to PD-1 blockade in TNBC. These data warrant confirmation in TNBC and exploration of induction treatments prior to PD-1 blockade in other cancer types.

Therapeutic targeting of macrophages enhances chemotherapy efficacy by unleashing type I interferon response

Recent studies have revealed a role for macrophages and neutrophils in limiting chemotherapy efficacy; however, the mechanisms underlying the therapeutic benefit of myeloid-targeting agents in combination with chemotherapy are incompletely understood. Here, we show that targeting tumour-associated macrophages by colony-stimulating factor-1 receptor (CSF-1R) blockade in the K14cre;Cdh1^F/F;Trp53^F/F transgenic mouse model for breast cancer stimulates intratumoural type I interferon (IFN) signalling, which enhances the anticancer efficacy of platinum-based chemotherapeutics. Notably, anti-CSF-1R treatment also increased intratumoural expression of type I IFN-stimulated genes in patients with cancer, confirming that CSF-1R blockade is a powerful strategy to trigger an intratumoural type I IFN response. By inducing an inflamed, type I IFN-enriched tumour microenvironment and by further targeting immunosuppressive neutrophils during cisplatin therapy, antitumour immunity was activated in this poorly immunogenic breast cancer mouse model. These data illustrate the importance of breaching multiple layers of immunosuppression during cytotoxic therapy to successfully engage antitumour immunity in breast cancer.

Molecular and epigenetic profiles of BRCA1-like hormone-receptor-positive breast tumors identified with development and application of a copy-number-based classifier

Background

BRCA1-mutated cancers exhibit deficient homologous recombination (HR) DNA repair, resulting in extensive copy number alterations and genome instability. HR deficiency can also arise in tumors without a BRCA1 mutation. Compared with other breast tumors, HR-deficient, BRCA1-like tumors exhibit worse prognosis but selective chemotherapeutic sensitivity. Presently, patients with triple negative breast cancer (TNBC) who do not respond to hormone endocrine-targeting therapy are given cytotoxic chemotherapy. However, more recent evidence showed a similar genomic profile between BRCA1-deficient TNBCs and hormone-receptor-positive tumors. Characterization of the somatic alterations of BRCA1-like hormone-receptor-positive breast tumors as a group, which is currently lacking, can potentially help develop biomarkers for identifying additional patients who might respond to chemotherapy.

Methods

We retrained and validated a copy-number-based support vector machine (SVM) classifier to identify HR-deficient, BRCA1-like breast tumors. We applied this classifier to The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast tumors. We assessed mutational profiles and proliferative capacity by covariate-adjusted linear models and identified differentially methylated regions using DMRcate in BRCA1-like hormone-receptor-positive tumors.

Results

Of the breast tumors in TCGA and METABRIC, 22% (651/2925) were BRCA1-like. Stratifying on hormone-receptor status, 13% (302/2405) receptor-positive and 69% (288/417) triple-negative tumors were BRCA1-like. Among the hormone-receptor-positive subgroup, BRCA1-like tumors showed significantly increased mutational burden and proliferative capacity (both P < 0.05). Genome-scale DNA methylation analysis of BRCA1-like tumors identified 202 differentially methylated gene regions, including hypermethylated BRCA1. Individually significant CpGs were enriched for enhancer regions (P < 0.05). The hypermethylated gene sets were enriched for DNA and chromatin conformation (all Bonferroni P < 0.05).

Conclusions

To provide insights into alternative classification and potential therapeutic targeting strategies of BRCA1-like hormone-receptor-positive tumors we developed and applied a novel copy number classifier to identify BRCA1-like hormone-receptor-positive tumors and their characteristic somatic alteration profiles.

Electronic supplementary material

The online version of this article (10.1186/s13058-018-1090-z) contains supplementary material, which is available to authorized users.

Targeting homologous repair deficiency in breast and ovarian cancers: Biological pathways, preclinical and clinical data

Mutation or epigenetic silencing of homologous recombination (HR) repair genes is characteristic of a growing proportion of triple-negative breast cancers (TNBCs) and high-grade serous ovarian carcinomas. Defects in HR lead to genome instability, allowing cells to acquire the multiple genetic alterations essential for cancer development. However, this deficiency can also be exploited by using DNA damaging agents or by targeting compensatory repair pathways. A noteworthy example is treatment of TNBC and epithelial ovarian cancer harboring BRCA1/2 germline mutations using platinum salts and/or PARP inhibitors. Dramatic responses to PARP inhibitors may support a wider use in the HR-deficient population beyond those with mutated germline BRCA1 and 2. In this review, we discuss HR deficiency hallmarks as predictive biomarkers for platinum salt and PARP inhibitor sensitivity for selecting patients affected by TNBC or epithelial ovarian cancer who could benefit from these therapeutic options.

Computational Investigation of Homologous Recombination DNA Repair Deficiency in Sporadic Breast Cancer

BRCAness has important implications in the management and treatment of patients with breast and ovarian cancer. In this study, we propose a computational framework to measure the BRCAness of breast and ovarian tumor samples based on their gene expression profiles. We define a characteristic profile for BRCAness by comparing gene expression differences between BRCA1/2 mutant familial tumors and sporadic breast cancer tumors while adjusting for relevant clinical factors. With this BRCAness profile, our framework calculates sample-specific BRCA scores, which indicates homologous recombination (HR)-mediated DNA repair pathway activity of samples. We found that in sporadic breast cancer high BRCAness score is associated with aberrant copy number of HR genes rather than somatic mutation and other genomic features. Moreover, we observed significant correlations of BRCA score with genome instability and neoadjuvant chemotherapy. More importantly, BRCA score provides significant prognostic value in both breast and ovarian cancers after considering established clinical variables. In summary, the inferred BRCAness from our framework can be used as a robust biomarker for the prediction of prognosis and treatment response in breast and ovarian cancers.

BRCA1-like profile is not significantly associated with survival benefit of non-myeloablative intensified chemotherapy in the GAIN randomized controlled trial

PURPOSE

The BRCA1-like profile identifies tumors with a defect in homologous recombination due to inactivation of BRCA1. This profile has been shown to predict which stage III breast cancer patients benefit from myeloablative, DNA double-strand-break-inducing chemotherapy. We tested the predictive potential of the BRCA1-like profile for adjuvant non-myeloablative, intensified dose-dense chemotherapy in the GAIN trial.

METHODS

Lymph node positive breast cancer patients were randomized to 3 × 3 dose-dense cycles of intensified epirubicin, paclitaxel, and cyclophosphamide (ETC) or 4 cycles concurrent epirubicin and cyclophosphamide followed by 10 cycles of weekly paclitaxel combined with 4 cycles capecitabine (EC-TX). Only triple negative breast cancer patients (TNBC) for whom tissue was available were included in these planned analyses. BRCA1-like or non-BRCA1-like copy number profiles were derived from low coverage sequencing data.

RESULTS

119 out of 163 TNBC patients (73%) had a BRCA1-like profile. After median follow-up of 83 months, disease free survival (DFS) was not significantly different between BRCA1-like and non-BRCA1-like patients [adjusted hazard ratio (adj.HR) 1.02; 95% confidence interval (CI) 0.55-1.86], neither was overall survival (OS; adj.HR 1.26; 95% CI 0.58-2.71). When split by BRCA1-like status, DFS and OS were not significantly different between treatments. However, EC-TX seemed to result in a trend to an improvement in DFS in patients with a BRCA1-like tumor, while the reverse accounted for ETC treatment in patients with a non-BRCA1-like tumor (p for interaction = 0.094).

CONCLUSIONS

The BRCA1-like profile is not associated with survival benefit for a non-myeloablative, intensified regimen in this study population. Considering the limited cohort size, capecitabine might have additional benefit for TNBC patients.

Relationship of immunohistochemistry, copy number aberrations and epigenetic disorders with BRCAness pattern in hereditary and sporadic breast cancer

The study aims to identify the relevance of immunohistochemistry (IHC), copy number aberrations (CNA) and epigenetic disorders in BRCAness breast cancers (BCs). We studied 95 paraffin included BCs, of which 41 carried BRCA1/BRCA2 germline mutations and 54 were non hereditary (BRCAX/Sporadic). Samples were assessed for BRCA1ness and CNAs by Multiplex Ligation-dependent Probe Amplification (MLPA); promoter methylation (PM) was assessed by methylation-specific-MLPA and the expression of miR-4417, miR-423-3p, miR-590-5p and miR-187-3p by quantitative RT-PCR. IHC markers Ki67, ER, PR, HER2, CK5/6, EGFR and CK18 were detected with specific primary antibodies (DAKO, Denmark). BRCAness association with covariates was performed using multivariate binary logistic regression (stepwise backwards Wald option). BRCA1/2 mutational status (p = 0.027), large tumor size (p = 0.041) and advanced histological grade (p = 0.017) among clinic-pathological variables; ER (p < 0.001) among IHC markers; MYC (p < 0.001) among CNA; APC (p = 0.065), ATM (p = 0.014) and RASSF1 (p = 0.044) among PM; and miR-590-5p (p = 0.001), miR-4417 (p = 0.019) and miR-423 (p = 0.013) among microRNA expression, were the selected parameters significantly related with the BRCAness status. The logistic regression performed with all these parameters selected ER+ as linked with the lack of BRCAness (p = 0.001) and MYC CNA, APC PM and miR-590-5p expression with BRCAness (p = 0.014, 0.045 and 0.007, respectively). In conclusion, the parameters ER expression, APC PM, MYC CNA and miR-590-5p expression, allowed detection of most BRCAness BCs. The identification of BRCAness can help establish a personalized medicine addressed to predict the response to specific treatments.

Identification of BRCA1-like triple-negative breast cancers by quantitative multiplex-ligation-dependent probe amplification (MLPA) analysis of BRCA1-associated chromosomal regions: a validation study

Background

Triple-negative breast cancer (TNBC) with a BRCA1-like molecular signature has been demonstrated to remarkably respond to platinum-based chemotherapy and might be suited for a future treatment with poly(ADP-ribose)polymerase (PARP) inhibitors. In order to rapidly assess this signature we have previously developed a multiplex-ligation-dependent probe amplification (MLPA)-based assay. Here we present an independent validation of this assay to confirm its important clinical impact.

Methods

One-hundred-forty-four TNBC tumor specimens were analysed by the MLPA-based “BRCA1-like” test. Classification into BRCA1-like vs. non-BRCA1-like samples was performed by our formerly established nearest shrunken centroids classifier. Data were subsequently compared with the BRCA1-mutation/methylation status of the samples. T-lymphocyte infiltration and expression of the main target of PARP inhibitors, PARP1, were assessed on a subset of samples by immunohistochemistry. Data acquisition and interpretation was performed in a blinded manner.

Results

In the studied TNBC cohort, 63 out of 144 (44 %) tumors were classified into the BRCA1-like category. Among these, the MLPA test correctly predicted 15 out of 18 (83 %) samples with a pathogenic BRCA1-mutation and 20 of 22 (91 %) samples exhibiting BRCA1-promoter methylation. Five false-negative samples were observed. We identified high lymphocyte infiltration as one possible basis for misclassification. However, two falsely classified BRCA1-mutated tumors were also characterized by rather non-BRCA1-associated histopathological features such as borderline ER expression. The BRCA1-like vs. non-BRCA1-like signature was specifically enriched in high-grade (G3) cancers (90 % vs. 58 %, p = 0.0004) and was also frequent in tumors with strong (3+) nuclear PARP1 expression (37 % vs. 16 %; p = 0.087).

Conclusions

This validation study confirmed the good performance of the initial MLPA assay which might thus serve as a valuable tool to select patients for platinum-based chemotherapy regimens. Moreover, frequent PARP1 upregulation in BRCA1-like tumors may also point to susceptibility to treatment with PARP inhibitors. Limitations are the requirement of high tumor content and high-quality DNA.

Smoking, second-hand smoke exposure and smoking cessation in relation to leukocyte telomere length and mortality

OBJECTIVES

To investigate the link between smoking exposure, telomere length and mortality, with emphasis on second-hand smoke (SHS) exposure and the duration of smoking cessation.

RESULTS

A total of 1,018 participants died during follow-up (mean: 10.3 years). A 50 base-pair decrease in LTL was shown among cotinine-confirmed current versus never smokers. The 90th quantile of LTL decreased with increasing cotinine among never smokers, indicating a role of SHS. Longer telomeres with smoking cessation were indicated but limited to a 3-16 year period of abstaining smoking. When assessing mortality, we observed a lower risk of all-cause death for the second quintile compared to the first among never smokers (HR: 0.67, 95% CI: 0.52-0.87), and a higher risk was found among current smokers (HR: 1.89, 1.19-2.92).

MATERIALS AND METHODS

We studied 6,456 nationally representative U.S. respondents with mortality follow-up through to 31 December 2011. Smoking status was assessed by interviews and cotinine levels. Relative leukocyte telomere length (LTL) was quantified by polymerase chain reaction (PCR). Multivariable linear regression was performed to examine LTL by smoking exposure, adjusted for age, sex, race/ethnicity, socioeconomic status, education, body mass index, alcohol consumption, and physical activity. We further estimated the association of LTL with cotinine levels using quantile regression, and with smoking cessation dynamics. Cox regression was used to estimate mortality by smoking status and LTL.

CONCLUSION

Our findings indicated a complex association between smoking, telomere length, and mortality. LTL alterations with SHS and smoking cessation warrant further investigation for translation to public health measures.

RUBIC identifies driver genes by detecting recurrent DNA copy number breaks

The frequent recurrence of copy number aberrations across tumour samples is a reliable hallmark of certain cancer driver genes. However, state-of-the-art algorithms for detecting recurrent aberrations fail to detect several known drivers. In this study, we propose RUBIC, an approach that detects recurrent copy number breaks, rather than recurrently amplified or deleted regions. This change of perspective allows for a simplified approach as recursive peak splitting procedures and repeated re-estimation of the background model are avoided. Furthermore, we control the false discovery rate on the level of called regions, rather than at the probe level, as in competing algorithms. We benchmark RUBIC against GISTIC2 (a state-of-the-art approach) and RAIG (a recently proposed approach) on simulated copy number data and on three SNP6 and NGS copy number data sets from TCGA. We show that RUBIC calls more focal recurrent regions and identifies a much larger fraction of known cancer genes.

Analysis of cancer genome sequencing data has been used to predict genes associated with the pathogenesis of cancer. Here, the authors propose a new algorithm entitled RUBIC that predicts breaks in DNA as opposed to previously published methods that predict amplifications and deletions of DNA.

Comprehensive Screening of Gene Copy Number Aberrations in Formalin-Fixed, Paraffin-Embedded Solid Tumors Using Molecular Inversion Probe-Based Single-Nucleotide Polymorphism Array

Gene copy number aberrations (CNAs) represent a major class of cancer-related genomic alterations that drive solid tumors. Comprehensive and sensitive detection of CNAs is challenging because of often low quality and quantity of DNA isolated from the formalin-fixed, paraffin-embedded (FFPE) solid tumor samples. Here, in a clinical molecular diagnostic laboratory, we tested the utility and validated a molecular inversion probe-based (MIP) array to routinely screen for CNAs in solid tumors. Using low-input FFPE DNA, the array detects genome-wide CNAs with a special focus on 900 cancer-related genes. A cohort of 76 solid tumors of various types and tumor cellularity (20% to 100%), and four cancer cell lines were used. These harbored CNAs in clinically important genes (ERBB2, EGFR, FGFR1, KRAS, MYC) as detected by orthogonal techniques like next-generation sequencing or fluorescence in situ hybridization. Results of the MIP array were concordant with results from orthogonal techniques, and also provided additional information regarding the allelic nature of the CNAs. Limit-of-detection and assay reproducibility studies showed a high degree of sensitivity and reproducibility of detection, respectively. FFPE compatibility, ability to detect CNAs with high sensitivity, accuracy, and provide valuable information such as loss of heterozygosity along with relatively short turnaround times makes the MIP array a desirable clinical platform for routine screening of solid tumors in a clinical laboratory.

BRCA1-like profile predicts benefit of tandem high dose epirubicin-cyclophospamide-thiotepa in high risk breast cancer patients randomized in the WSG-AM01 trial

BRCA1 is an important protein in the repair of DNA double strand breaks (DSBs), which are induced by alkylating chemotherapy. A BRCA1-like DNA copy number signature derived from tumors with a BRCA1 mutation is indicative for impaired BRCA1 function and associated with good outcome after high dose (HD) and tandem HD DSB inducing chemotherapy. We investigated whether BRCA1-like status was a predictive biomarker in the WSG AM 01 trial. WSG AM 01 randomized high-risk breast cancer patients to induction (2× epirubicin-cyclophosphamide) followed by tandem HD chemotherapy with epirubicin, cyclophosphamide and thiotepa versus dose dense chemotherapy (4× epirubicin-cyclophospamide followed by 3× cyclophosphamide-methotrexate-5-fluorouracil). We generated copy number profiles for 143 tumors and classified them as being BRCA1-like or non-BRCA1-like. Twenty-six out of 143 patients were BRCA1-like. BRCA1-like status was associated with high grade and triple negative tumors. With regard to event-free-survival, the primary endpoint of the trial, patients with a BRCA1-like tumor had a hazard rate of 0.2, 95% confidence interval (CI): 0.07-0.63, p = 0.006. In the interaction analysis, the combination of BRCA1-like status and HD chemotherapy had a hazard rate of 0.19, 95% CI: 0.067-0.54, p = 0.003. Similar results were observed for overall survival. These findings suggest that BRCA1-like status is a predictor for benefit of tandem HD chemotherapy with epirubicin-thiotepa-cyclophosphamide.

Robust BRCA1-like classification of copy number profiles of samples repeated across different datasets and platforms

Breast cancers with BRCA1 germline mutation have a characteristic DNA copy number (CN) pattern. We developed a test that assigns CN profiles to be 'BRCA1-like' or 'non-BRCA1-like', which refers to resembling a BRCA1-mutated tumor or resembling a tumor without a BRCA1 mutation, respectively. Approximately one third of the BRCA1-like breast cancers have a BRCA1 mutation, one third has hypermethylation of the BRCA1 promoter and one third has an unknown reason for being BRCA1-like. This classification is indicative of patients' response to high dose alkylating and platinum containing chemotherapy regimens, which targets the inability of BRCA1 deficient cells to repair DNA double strand breaks. We investigated whether this classification can be reliably obtained with next generation sequencing and copy number platforms other than the bacterial artificial chromosome (BAC) array Comparative Genomic Hybridization (aCGH) on which it was originally developed. We investigated samples from 230 breast cancer patients for which a CN profile had been generated on two to five platforms, comprising low coverage CN sequencing, CN extraction from targeted sequencing panels (CopywriteR), Affymetrix SNP6.0, 135K/720K oligonucleotide aCGH, Affymetrix Oncoscan FFPE (MIP) technology, 3K BAC and 32K BAC aCGH. Pairwise comparison of genomic position-mapped profiles from the original aCGH platform and other platforms revealed concordance. For most cases, biological differences between samples exceeded the differences between platforms within one sample. We observed the same classification across different platforms in over 80% of the patients and kappa values of at least 0.36. Differential classification could be attributed to CN profiles that were not strongly associated to one class. In conclusion, we have shown that the genomic regions that define our BRCA1-like classifier are robustly measured by different CN profiling technologies, providing the possibility to retro- and prospectively investigate BRCA1-like classification across a wide range of CN platforms.