Spatial molecular profiling of mixed invasive ductal and lobular breast cancers reveals heterogeneity in intrinsic molecular subtypes, oncogenic signatures, and mutations
Mixed invasive ductal and lobular carcinoma (MDLC) is a rare histologic subtype of breast cancer displaying both E-cadherin positive ductal and E-cadherin negative lobular morphologies within the same tumor, posing challenges with regard to anticipated clinical management. It remains unclear whether these distinct morphologies also have distinct biology and risk of recurrence. Our spatially resolved transcriptomic, genomic, and single-cell profiling revealed clinically significant differences between ductal and lobular tumor regions including distinct intrinsic subtype heterogeneity – e.g., MDLC with triple-negative breast cancer (TNBC) or basal ductal and estrogen receptor positive (ER+) luminal lobular regions, distinct enrichment of cell cycle arrest/senescence and oncogenic (ER and MYC ) signatures, genetic and epigenetic CDH1 inactivation in lobular but not ductal regions, and single-cell ductal and lobular subpopulations with unique oncogenic signatures further highlighting intraregional heterogeneity. Altogether, we demonstrated that the intratumoral morphological/histological heterogeneity within MDLC is underpinned by intrinsic subtype and oncogenic heterogeneity which may result in prognostic uncertainty and therapeutic dilemma.
A Multiparameter Molecular Classifier to Predict Response to Neoadjuvant Lapatinib plus Trastuzumab without Chemotherapy in HER2+ Breast Cancer
Abstract Purpose: Clinical trials reported 25% to 30% pathologic complete response (pCR) rates in HER2+ patients with breast cancer treated with anti-HER2 therapies without chemotherapy. We hypothesize that a multiparameter classifier can identify patients with HER2-“addicted” tumors who may benefit from a chemotherapy-sparing strategy. Experimental Design: Baseline HER2+ breast cancer specimens from the TBCRC023 and PAMELA trials, which included neoadjuvant treatment with lapatinib and trastuzumab, were used. In the case of estrogen receptor–positive (ER+) tumors, endocrine therapy was also administered. HER2 protein and gene amplification (ratio), HER2-enriched (HER2-E), and PIK3CA mutation status were assessed by dual gene protein assay (GPA), research-based PAM50, and targeted DNA-sequencing. GPA cutoffs and classifier of response were constructed in TBCRC023 using a decision tree algorithm, then validated in PAMELA. Results: In TBCRC023, 72 breast cancer specimens had GPA, PAM50, and sequencing data, of which 15 had pCR. Recursive partitioning identified cutoffs of HER2 ratio ≥ 4.6 and %3+ IHC staining ≥ 97.5%. With PAM50 and sequencing data, the model added HER2-E and PIK3CA wild-type (WT). For clinical implementation, the classifier was locked as HER2 ratio ≥ 4.5, %3+ IHC staining ≥ 90%, and PIK3CA-WT and HER2-E, yielding 55% and 94% positive (PPV) and negative (NPV) predictive values, respectively. Independent validation using 44 PAMELA cases with all three biomarkers yielded 47% PPV and 82% NPV. Importantly, our classifier's high NPV signifies its strength in accurately identifying patients who may not be good candidates for treatment deescalation. Conclusions: Our multiparameter classifier differentially identifies patients who may benefit from HER2-targeted therapy alone from those who need chemotherapy and predicts pCR to anti-HER2 therapy alone comparable with chemotherapy plus dual anti-HER2 therapy in unselected patients.
Molecular Subclasses of Clear Cell Ovarian Carcinoma and Their Impact on Disease Behavior and Outcomes
Abstract Purpose: To identify molecular subclasses of clear cell ovarian carcinoma (CCOC) and assess their impact on clinical presentation and outcomes. Experimental Design: We profiled 421 primary CCOCs that passed quality control using a targeted deep sequencing panel of 163 putative CCOC driver genes and whole transcriptome sequencing of 211 of these tumors. Molecularly defined subgroups were identified and tested for association with clinical characteristics and overall survival. Results: We detected a putative somatic driver mutation in at least one candidate gene in 95% (401/421) of CCOC tumors including ARID1A (in 49% of tumors), PIK3CA (49%), TERT (20%), and TP53 (16%). Clustering of cancer driver mutations and RNA expression converged upon two distinct subclasses of CCOC. The first was dominated by ARID1A-mutated tumors with enriched expression of canonical CCOC genes and markers of platinum resistance; the second was largely comprised of tumors with TP53 mutations and enriched for the expression of genes involved in extracellular matrix organization and mesenchymal differentiation. Compared with the ARID1A-mutated group, women with TP53-mutated tumors were more likely to have advanced-stage disease, no antecedent history of endometriosis, and poorer survival, driven by their advanced stage at presentation. In women with ARID1A-mutated tumors, there was a trend toward a lower rate of response to first-line platinum-based therapy. Conclusions: Our study suggests that CCOC consists of two distinct molecular subclasses with distinct clinical presentation and outcomes, with potential relevance to both traditional and experimental therapy responsiveness. See related commentary by Lheureux, p. 4838
A Distinctive Adnexal (Usually Paratubal) Neoplasm Often Associated With Peutz-Jeghers Syndrome and Characterized by STK11 Alterations (STK11 Adnexal Tumor)
We describe 22 examples of a novel, usually paratubal, adnexal tumor associated with Peutz-Jeghers syndrome in nearly 50% of cases that harbored STK11 alterations in all tested (n=21). The patients ranged from 17 to 66 years (median=39 y) and the tumors from 4.5 to 25.5 cm (median=11 cm). Most (n=18) were paratubal, with metastases noted in 11/22 (50%) and recurrences in 12/15 (80%). Morphologically, they were characterized by interanastomosing cords and trabeculae of predominantly epithelioid cells, set in a variably prominent myxoid to focally edematous stroma, that often merged to form tubular, cystic, cribriform, and microacinar formations, reminiscent of salivary gland-type tumors. The tumor cells were uniformly atypical, often with prominent nucleoli and a variable mitotic index (median=9/10 HPFs). The tumors were usually positive to a variable extent for epithelial (CAM5.2, AE1/AE3, cytokeratin 7), sex cord (calretinin, inhibin, WT1), and mesothelial (calretinin, D2-40) markers, as well as hormone receptors. PAX8, SF1, and GATA-3 were rarely positive, while claudin-4, FOXL2, and TTF-1 were consistently negative. All sequenced tumors (n=21) harbored alterations in STK11, often with a loss of heterozygosity event. There were no other recurrently mutated genes. Recurrent copy number alterations included loss of 1p and 11q, and gain of 1q, 15q, and 15p. Despite an extensive morphologic, immunohistochemical, and molecular evaluation, we are unable to determine with certainty the histogenesis of this unique tumor. Wolffian, sex cord stromal, epithelial, and mesothelial origins were considered. We propose the term STK11 adnexal tumor to describe this novel entity and emphasize the importance of genetic counseling in these patients as a significant number of neoplasms occur in association with Peutz-Jeghers syndrome.
A multiparameter classifier to predict response to lapatinib plus trastuzumab (LT) without chemotherapy in HER2+ breast cancer (BC).
1011 Background: Several trials have shown 25-30% pathologic complete response (pCR) rates in patients with HER2+ BC treated with LT therapy (+ endocrine therapy if ER+), but no chemotherapy (CTX). We hypothesize that a multiparameter classifier, comprised of HER2 gene and protein levels, intratumor heterogeneity (ITH), HER2-enriched (E) subtype, and PIK3CA mutation status can identify patients whose tumors are “addicted” to HER2 signaling and are likely to achieve pCR from a CTX-sparing de-escalation strategy. Methods: Baseline specimens from 2 trials (TBCRC023 [NCT00999804] and PAMELA [NCT01973660]) of neoadjuvant CTX-sparing LT (+ endocrine therapy if ER+) in HER2+ BC were used. HER2 protein and ITH (scored for % 3+ by IHC), and gene amplification (HER2:CEP17 ratio and copy number (CN) by CISH) were measured on the same slide by the dual gene protein assay (GPA). HER2-E and PIK3CA mutation status were assessed by research-based PAM50 and MSK-IMPACT platforms, respectively. A decision tree algorithm was used to determine the GPA cutoffs and to construct the classifier of response (by pCR) in TBCRC023, which was then validated in PAMELA. Results: Of the evaluable patients from TBCRC023 (N = 130) and PAMELA (N = 151), GPA data were available for 121 and 94 cases, respectively. Both cohorts exhibited similar distributions for HER2 ratio, CN, and % 3+, and a strong correlation between HER2 ratio and CN (R > 0.92). In TBCRC023, 73 cases had data from GPA, PAM50, and IMPACT, of which 15 had pCR. Recursive partitioning identified cutoffs of HER2 ratio > 4.6 and % 3+ > 97.5% in both the GPA data cohort (N = 121) and complete data cohort (N = 73). With PAM50 and IMPACT data, the model added HER2-E and PIK3CA wild-type (wt). For practical reasons, the classifier was locked as HER2 ratio ≥ 4.5 AND % 3+ ≥ 90% AND PIK3CA-wt AND HER2-E, which yielded a PPV of 55% and NPV of 94%. Validation in PAMELA using 45 cases with data for all 3 assays yielded PPV of 44% and NPV of 82%. 29 TBCRC023 cases without IMPACT data could be predicted to be non-pCR, of which 26 were correct (NPV = 89%). In PAMELA, 66 additional cases could be predicted to be non-pCR, of which 54 were correct (NPV = 81%). Conclusions: We have constructed a multiparameter classifier that can predict pCR with targeted therapy alone that compare to pCR rates of CTX + dual anti-HER2 in unselected patients. Prospective validation in a clinical trial is warranted. [Table: see text]
Abstract PD1-15: The landscape of somatic genetic alterations in breast cancers from ATM germline mutation carriers
Abstract Introduction:Pathogenic and/or founder germline variants in the ataxia-telangiectasia mutated (ATM) gene confer an increased breast cancer (BC) risk. The protein kinase ATM plays a central role inDNA double-strand break-repair and in the activation of downstream targets such as p53 and BRCA1. We sought to define the repertoire of somatic genetic alterations of BCs from patients with pathogenic germline ATM mutations and whether somatic loss of heterozygosity (LOH) of ATM would be present in these cancers. Methods: 21 BCs from ATM germline mutation carriers were microdissected. Tumor and normal DNA samples were subjected to whole-exome sequencing (WES, n=12) or massively parallel sequencing targeting all coding regions and selected intronic and regulatory regions of 410 key cancer genes (n=9). Somatic mutations, copy number alterations, cancer cell fractions, large-scale state transitions (LSTs) and mutational signatures were defined using state-of-the-art bioinformatics algorithms. ABSOLUTE and FACETS were employed to assess LOH of the wild-type allele of ATM. Results: Of the patients included in this study, 71%, 24% and 5% of cases harbored ATM missense (all but one p.V2424G), frame-shift and nonsense germline mutations, respectively. All tumors were ER-positive and four (19%) were HER2-positive. The median age of the patients was 46 years (32–79 years). Our analyses revealed biallelic inactivation of ATM through LOH of the wild-type allele in 16 of 21 cases (76%), and second somatic ATM mutations were not found. The median number of non-synonymous somatic mutations was 38 (range 15-113) and 2 (range 0-8)in tumors subjected to WES and targeted sequencing, respectively. The repertoire of somatic genetic alterations of ATM-associated BCs was found to be heterogeneous, including clonal PIK3CA mutations (24%), GATA3 mutations (19%), FANCI amplifications (19%) and CCND1 amplifications (14%). Importantly, however, no somatic mutations affecting TP53 were found. Analysis of the WES data revealed that 5 (42%) ATM-associated BCs displayed high LST scores, all of which harbored bi-allelic ATM inactivation. In contrast to BRCA1- and BRCA2-associated BCs, which frequently display the mutational signature 3 associated with defective homologous recombination DNA repair, the ATM-associated BCs studied displayed the ageing mutational signature (i.e. signature 1). Comparison of the mutational profiles of the ATM--associated BCs subjected to WES (n=12) with those of BRCA1- (n=11) and BRCA2-associated (n=10) BCs from The Cancer Genome Atlas revealed that TP53 was more frequently mutated in BCs from BRCA1 germline mutation carriers (0% vs 72%, P<0.001), while no differences with BRCA2-associated BCs were found. Conclusion: ATM-associated BCs frequently display bi-allelic ATM inactivation through LOH of the wild-type allele and a subset of these cases displayed high levels of LSTs. These findings suggest that at least in a subset of ATM-associated BCs, biallelic inactivation of ATM rather than a dominant negative effect of the germline mutation may be the mechanism of inactivation of this tumor suppressor gene. The repertoire of somatic genetic alterations of ATM-associated BCs is heterogeneous, with a noticeable lack of TP53 somatic mutations. Citation Format: Weigelt B, Bi R, Kumar R, James PA, Thorne H, Couch FJ, Eccles DM, Blows F, Geyer FC, Li A, Selenica P, Lim RS, Blecua P, Shen R, Wen H, Robson ME, Reis-Filho JS, Chenevix-Trench G. The landscape of somatic genetic alterations in breast cancers from ATM germline mutation carriers [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr PD1-15.
The Landscape of Somatic Genetic Alterations in Breast Cancers From ATM Germline Mutation Carriers
Abstract S4-03: A functional assay for homologous recombination (HR) DNA repair and whole exome sequencing reveal that HR-defective sporadic breast cancers are enriched for genetic alterations in DNA repair genes
Abstract Background: Germline mutations in the BRCA1 and BRCA2 genes lead to hereditary breast cancers that are defective in homologous recombination (HR) repair and sensitive to DNA damaging agents. HR deficiency (HRD) also occurs in sporadic breast cancers, but its incidence and etiology are unclear. Genomic signatures of HRD recently were employed as biomarkers in clinical trials with modest success. We posited that sporadic breast cancers displaying functional HR deficiency would harbor genetic alterations affecting HR DNA repair genes. To test this hypothesis, we applied a functional assay to define lack of competent HR DNA repair and sequenced the exomes of consecutive sporadic breast cancers. Methods: We developed an assay to assess the ability of cancer cells to localize RAD51 into sub nuclear foci in response to ex-vivo irradiation (IR) in fresh sporadic breast cancer tissue specimens from 60 patients. RAD51 focus formation was compared between mock and IR conditions to determine relative fold induction. Twenty-nine tumors with sufficient DNA underwent whole-exome sequencing. Structural genomic signatures of HRD (i.e. large state transitions (LST), telomeric imbalance (NtAI), loss-of-heterozygositiy (LOH)) and a previously reported mutational signature related to BRCA1/2 hereditary breast cancers were assessed. HR deficient tumors were defined as those with both RAD51 foci defects and a genomic signature of HR deficiency (LST>15 or presence of a BRCA mutational signature). Somatic, germ-line, and copy number changes in HR genes were investigated. BRCA1 methylation status was determined. Results: Seventeen of 60 (28%) tumors displayed defective RAD51 recruitment following ex-vivo IR (RAD51-DEF). RAD51-DEF was seen in all breast cancer subtypes , including 7 of 33 (21%) ER+/HER2-, 4 of 14 (29%) HER2+, and 6 of 13 (48%) triple-negative cases. Of the 29 sequenced tumors, 13 (45%) were RAD51-DEF and 16 (55%) were competent for inducing RAD51 foci. LST was elevated in 10 tumors (LST >15) and associated with RAD51-DEF (p=0.02), whereas NtAI (p=0.10) and LOH (p=0.052) did not show a significant association with RAD51-DEF. The BRCA1/2 mutational signature was evident in 4 tumors, all were RAD51-DEF (p=0.03) and 2 were BRCA2 mutated. Nine of 29 (31%) sequenced tumors were determined to have HRD by the RAD51 assay and presence of a genomic scar. Eight of these 9 (88%) cases with HRD had a genetic alteration of both alleles of a bona fide HR gene due to a pathogenic mutation (somatic or germline) coupled with loss of heterozygosity or a homozygous deletion compared to 1 (5%) tumor without HRD (p<0.001). Both alleles of a gene were affected for BRCA2 (n=4), FAAP100 (n=2), CHEK2 (n=1), TP53BP1 (n=1) and BRCA1 (n=1). BRCA1 gene promoter methylation was found not to be significantly associated with HRD. Conclusion: Combined functional and genomic analyses of breast tumors demonstrated that genetic loss of an HR gene may underpin HRD in sporadic breast cancers. Our findings warrant further comprehensive genetic assessment (somatic, germline, and copy number) of HR genes as potential biomarker for HR-directed therapies. Citation Format: Powell SN, Riaz N, Mutter RW, Ng CKY, Delsite R, Piscuoglio S, King TA, Martelotto L, Sakr R, Brogi E, Edelweiss M, Lim R, Higginson D, Weigelt B, Lee W, Reis-Filho JS. A functional assay for homologous recombination (HR) DNA repair and whole exome sequencing reveal that HR-defective sporadic breast cancers are enriched for genetic alterations in DNA repair genes. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr S4-03.