Perou, Charles M.

AbstractBackgroundReproductive factors have been shown to be differentially associated with risk of estrogen receptor (ER)-positive and ER-negative breast cancer. However, their associations with intrinsic-like subtypes are less clear.MethodsAnalyses included up to 23 353 cases and 71 072 controls pooled from 31 population-based case-control or cohort studies in the Breast Cancer Association Consortium across 16 countries on 4 continents. Polytomous logistic regression was used to estimate the association between reproductive factors and risk of breast cancer by intrinsic-like subtypes (luminal A-like, luminal B-like, luminal B-HER2–like, HER2-enriched–like, and triple-negative breast cancer) and by invasiveness. All statistical tests were 2-sided.ResultsCompared with nulliparous women, parous women had a lower risk of luminal A-like, luminal B-like, luminal B-HER2–like, and HER2-enriched–like disease. This association was apparent only after approximately 10 years since last birth and became stronger with increasing time (odds ratio [OR] = 0.59, 95% confidence interval [CI] = 0.49 to 0.71; and OR = 0.36, 95% CI = 0.28 to 0.46 for multiparous women with luminal A-like tumors 20 to less than 25 years after last birth and 45 to less than 50 years after last birth, respectively). In contrast, parous women had a higher risk of triple-negative breast cancer right after their last birth (for multiparous women: OR = 3.12, 95% CI = 2.02 to 4.83) that was attenuated with time but persisted for decades (OR = 1.03, 95% CI = 0.79 to 1.34, for multiparous women 25 to less than 30 years after last birth). Older age at first birth (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) and breastfeeding (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) were associated with lower risk of triple-negative breast cancer but not with other disease subtypes. Younger age at menarche was associated with higher risk of all subtypes; older age at menopause was associated with higher risk of luminal A-like but not triple-negative breast cancer. Associations for in situ tumors were similar to luminal A-like.ConclusionsThis large and comprehensive study demonstrates a distinct reproductive risk factor profile for triple-negative breast cancer compared with other subtypes, with implications for the understanding of disease etiology and risk prediction.

Abstract Introduction: The ALTERNATE trial randomized postmenopausal women with ER Allred 6-8 HER2- breast cancer to 6 months of NET with anastrozole (A), fulvestrant (F) or the combination (A+F). Biopsies were taken preNET and after 4-weeks(wks). Patients with Ki67 values >10% at 4-wks were offered triage to neoadjuvant chemotherapy. Patients with on-treatment Ki67 ≤ 10% who completed NET underwent surgery and Ki67 was reassessed. The primary endpoint was endocrine-sensitive disease rate (ESDR). ESD is defined as pCR or PEPI-0 residual disease (pT1-2, pN0, Ki67 ≤ 2.7%). We previously reported that the ESDR difference between the F-containing arms and the A arm was not >10% (ASCO 2020) and that baseline RNA-seq-based intrinsic subtypes predicted outcomes overall (SABCS 2021). Herein we describe relationships between PAM50 intrinsic subtype and Ki67 values by treatment arm because comparative drug effectiveness in adjuvant endocrine therapy studies in ER+ HER2- breast cancer can be predicted by the degree of Ki67 suppression (PMC3518447). Methods: 743 of the 1297 eligible patients (A: 264; F: 231; A+F: 248) had RNA extracted from preNET frozen tumor biopsies with >50% tumor content and subjected to RNA seq. Intrinsic subtypes were then assigned as LumA, LumB, and NonLum (Basal or HER2-E) using open-source PAM50-based informatics. Differences in the proportion with wk4 Ki67 > 10%, % change in wk4 ki67, and surgical CCCA (Ki67 ≤ 2.7%) rate (sxCCCA) between treatments and by intrinsic subtype was assessed using stratified logistic regression, Wilcoxon rank sum test, and Fisher’s exact test, respectively. Analysis of sxCCCA excluded those who failed to complete NET for reasons other than disease progression or early Ki67 >10%. Results: Amongst the 358 LumA cases there were no significant differences in Ki67-based endpoints between treatments. Among the 292 LumB cases, the wk4 ki67 > 10% rate was lower with A+F (19.4%) than A (43%) (P=0.0002) and was somewhat lower in F (31%) versus A (P=0.076). The % change in wk4 Ki67 in LumB cases, adjusted for baseline Ki67, showed markedly superior suppression for A+F versus A (-90% vs. -77%; P=<0.0001) and versus F (-90% vs. -80%; P=0.0026). Furthermore sxCCCA rates were significantly higher with A+F than A (53% vs. 25% P = <0.0001) and somewhat higher for F (37%) than A (p=0.068), indicating that superior antiproliferative effects for A+F persist after 6 months on therapy. Lack of Ki67 suppression in response to treatment was observed in the majority of 43 NonLum samples regardless of treatment. Conclusion: The combination of A+F was significantly more effective than either drug alone for the control of LumB breast cancer cell proliferation. This suggests that A+F may be a more effective adjuvant endocrine therapy than A alone in LumB disease. The lower Ki67 suppression with A alone also suggests that poorer outcome in some LumB tumors may be due to insufficient ER targeting rather than ER-independent tumor growth Support: U10CA180821, U10CA180882, U24CA196171, UG1CA189856, U10CA180868 (NRG), NCI BIQSFP, BCRF, Genentech, AstraZeneca. https://acknowledgments.alliancefound.org. (MJE) CPRIT RR140033, P50-CA186784, P50-CA58223, U01-CA214125, U24-CA210954, Gift from Ralph and Lisa Eads, McNair Scholarship. ClinicalTrials.gov Identifier: NCT01953588 Citation Format: Matthew J. Ellis, Meenakshi Anurag, Jeremy Hoog, Aranzazu Fernandez-Martinez, Cheng Fan, Richard Gibbs, Souzan Sanati, Kiran Vij, Mark Watson, Travis Dockter, Olwen Hahn, Joseph Guenther, Abigail Caudle, Erica Crouch, Amy Tiersten, Monica Mita, Wajeeha Razaq, Tina J. Hieken, Yang Wang, A. Marilyn Leitch, Gary W. Unzeitig, Eric Winer, Anna Weiss, Kelly Hunt, Ann H. Partridge, Charles M. Perou, Vera Suman, Cynthia X. Ma, Lisa A. Carey. The effect of intrinsic subtype on inhibition of tumor growth by anastrozole vs. fulvestrant vs. the combination: Results from the Alliance neoadjuvant endocrine therapy (NET) ALTERNATE trial [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr CT026.

Abstract Background: Neoadjuvant ET (NET) offers an opportunity to assess ET sensitivity for ER+ HER2- BC and potentially to tailor therapy. Ki67 >10% on biopsy after 2-4 weeks (wks) of NET identifies patients (pts) with intrinsic ET resistance; while pathologic complete response (pCR) and modified preoperative endocrine prognostic index of 0 (mPEPI 0: pT1-2N0, Ki67 ≤2.7%) at surgery indicates sensitivity to ET. However, on-NET biopsy is not always acceptable or feasible and delays the ET sensitivity determination. PAM50 ROR score and intrinsic subtypes by tumor RNA profiling are prognostic in pts with early stage ER+ HER2- BC, and predict pCR rates to neoadjuvant chemotherapy (NCT) (PMC2667820). We therefore hypothesized that PAM50 analysis on pre-NET biopsies could predict the likelihood of a) a high on-NET Ki67, b) mPEPI-0 or pCR at surgery and, c) pCR for pts triaged to NCT. Methods: The ALTERNATE trial is a phase III study that randomized postmenopausal pts with clinical stage II/III ER+ (Allred score 6-8) HER2- BC to receive neoadjuvant anastrozole, fulvestrant, or both for 6 months before surgery. Research biopsy was required at pre-NET and wk 4, then optional at wk 12. Pts with Ki67 >10% on biopsy at wk 4 or 12 discontinued NET and were offered NCT. PAM50 intrinsic subtype and ROR-P values were generated from mRNA sequencing (RNASeq) analysis on pre-NET biopsies using open-source informatics (PMC7723687) and evaluated for prediction of on-NET Ki67 >10% at wk 4 or 12, pCR or mPEPI-0 post NET, and pCR post NCT. Results: 749 of 1,297 eligible trial pts were included in the analyses, after excluding 548 pts due to insufficient pre-NET tumor for RNASeq (n=511) or PAM50 normal subtype (n=37). Similar to the entire ALTERNATE population, the rate of Ki67 >10% at wk 4 or 12 was 24.4% (95% CI: 21.4-27.7%) and the rate of mPEPI-0/pCR post NET was 19.8% (95% CI: 17.0-22.8%). There were 393 (52.5%) Lum A, 302 (40.3%) Lum B, and 54 (7.2%) non-Lum (9 Basal, 45 HER2-E) BCs. These included 196 (26.2%) ROR-P low, 354 (47.3%) ROR-P medium and 199 (26.6%) ROR-P high BCs. Both the rates of Ki67 >10% at wk 4 or 12 and mPEPI-0/pCR differed significantly with respect to PAM50 subtype or ROR-P category, such that Lum A or ROR-P low BCs were least likely to have a Ki67 >10% at wk 4 or 12 and most likely to achieve mPEPI-0/pCR (Table). 93 of 168 (55.4%) pts triaged to NCT had RNA-seq results, yielding 26 Lum A, 49 Lum B, 4 Basal and 14 HER2-E, with the pCR rates of 0%, 6.1%, 0%, and 21.4%, respectively. There were 10 ROR-P low, 39 medium, and 44 high tumors, with a pCR rate of 0%, 5.1% and 9.1%, respectively. Conclusion: These data indicate that both baseline ROR-P and intrinsic subtype are predictive of early on-NET Ki67 > 10% and mPEPI 0/pCR at surgery after NET. For pts triaged to NCT based on an early on-NET Ki67 >10%, the HER2-E group had the highest pCR rate (20%) and no pCRs were observed in Lum A. These data may be useful for directing neoadjuvant therapy in postmenopausal pts with ER+ HER2- BC. Support: U10CA180821, U10CA180882, U24CA196171, UG1CA189856, U10CA180868 (NRG), NCI BIQSFP, BCRF, Genentech, AstraZeneca. https://acknowledgments.alliancefound.org. (MJE) CPRIT RR140033, P50CA186784, P50-CA58223, U01 CA214125, U24CA210954, Gift from Ralph and Lisa Eads, McNair Scholarship. Trials.gov Identifier: NCT01953588. Table 1.Rates of Ki67 >10% and mPEPI-0/pCR post NET by PAM50 subtype and ROR-P categoryKi67 >10% at wk 4 or 12mPEPI 0/pCR post NETPAM50 SubtypenYes, n (%)PnNo, n (%)PLum A37251 (13.7%) 95% CI: 10.4-17.6%<0.0001393104 (26.5%) 95%CI: 22.2-31.1%<0.0001Lum B29394 (32.1%) 95% CI: 26.8-37.8%30243 (14.2%) 95%CI: 10.5-18.7%Non-luminal (Basal and HER2-E)5338 (71.7%) 95%CI: 57.6-83.2%541 (1.9%) 95%CI: 0.05-9.9%ROR-P CategorynYes, n (%)PnNo, n (%)PLow18018 (10.0%) 95%CI: 6.0-15.3%<0.000119660 (30.6%) 95%CI: 24.2-37.6%<0.0001Intermediate34474 (21.5%) 95%CI: 17.3-26.2%35471 (20.1%) 95%CI: 16.0-24.6%High19491 (46.9%) 95%CI: 39.7-54.2%19917 (8.5%) 95%CI: 5.1-13.3% Citation Format: Cynthia X Ma, Meenakshi Anurag, Travis Dockter, Jeremy Hoog, Aranzazu Fernandez-Martinez, Cheng Fan, Richard Gibbs, Souzan Sanati, Kiran Vij, Mark Watson, Olwen Hahn, Joseph Guenther, Abigail Caudle, Erika Crouch, Amy Tiersten, Monica Mita, Wajeeha Razaq, Tina J Hieken, Yang Wang, A. Marilyn Leitch, Gary W Unzeitig, Anna Weiss, Eric P Winer, Kelly Hunt, Ann H Partridge, Lisa A Carey, Charles M Perou, Matthew J Ellis, Vera Suman. Pam50 intrinsic subtype and risk of recurrence score (ROR) for the prediction of endocrine (ET) sensitivity and pathologic response to chemotherapy in postmenopausal women with clinical stage II/III estrogen receptor positive (ER+) and HER2 negative (HER2-) breast cancer (BC) in the alternate trial (Alliance A011106) [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr PD9-03.

PURPOSE CALGB 40603 ( NCT00861705 ), a 2 × 2 randomized phase II trial, demonstrated that adding carboplatin or bevacizumab to weekly paclitaxel (wP) followed by doxorubicin and cyclophosphamide significantly increased the pathologic complete response (pCR) rate in stage II-III triple-negative breast cancer. We now report long-term outcomes (LTOs) and correlative science end points. PATIENTS AND METHODS The Kaplan-Meier method was used to estimate LTOs in 443 patients who initiated study treatment. Log-rank tests and Cox proportional hazards models evaluated the impact of clinical characteristics, pathologic response, calculated residual cancer burden (RCB) in patients with residual disease (RD), treatment assignment, and dose delivery during wP on LTOs, including event-free survival (EFS). Genomic predictors of treatment response and outcomes were assessed on pretreatment tumor samples by mRNA sequencing. RESULTS Among baseline characteristics, only the clinical stage was associated with LTOs. At a median follow-up of 7.9 years, LTOs were not significantly improved with either carboplatin or bevacizumab, overall or in patients with basal-like subtype cancers by genomic analysis. Patients with pCR (n = 205, 46.3%) had significantly higher 5-year EFS (85.5% v 56.6%, log-rank P < .0001) and overall survival (87.9% v 63.4%, P < .0001) rates compared with patients with RD, even those with RCB class I. Among clinical and genomic features, evidence of immune activation, including tumor-infiltrating lymphocytes and low B-cell receptor evenness, was associated with pCR and improved EFS. CONCLUSION Despite higher pCR rates, neither carboplatin nor bevacizumab appeared to improve LTOs although the study was not powered to assess these secondary end points. pCR was associated with superior LTOs even when compared with minimal RD. Markers of immune activation in pretreatment tumor biopsies were independently associated with higher pCR rates and improved survival.

Abstract Patients with metastatic breast cancer (MBC) typically have short survival times and their successful treatment represents one of most challenging aspects of patient care. This poor prognostic behavior is in part due to molecular features including increased tumor cell clonal heterogeneity, multiple drug resistance mechanisms, and alterations of the tumor microenvironment. The AURORA US Metastasis Project was established with the goal to identify molecular features specifically associated with metastasis. We therefore collected and molecularly characterized specimens from 55 metastatic breast cancer (BC) patients representing 51 primary cancers and 102 metastases. The 153 unique tumors were assayed using RNAseq, tumor/germline DNA exomes and low pass whole genome sequencing, and global DNA methylation microarrays. We found intrinsic molecular subtype differences between primary tumors and their matched metastases to be rare in triple negative breast cancer (TNBC)/Basal-like subtype tumors. Conversely, tumor subtype changes were relatively frequent in estrogen receptor positive (ER+) cancers where ~30% of Luminal A cases switched to Luminal B or HER2-enriched (HER2E) subtypes. Clonal evolution studies identified changes in expression subtype coincident with DNA clonality shifts, especially involving HER2 amplification and/or the HER2E expression subtype. We further found evidence for ER-mediated downregulation of genes involved in cell-cell adhesion in metastases. Microenvironment differences varied according to tumor subtype where ER+/Luminal metastases had lower fibroblast and endothelial cell content, while TNBC/Basal-like metastases showed a dramatic decrease in B cells and T cells. In 17% of metastatic tumors, we identified DNA hypermethylation and/or focal DNA deletions near HLA-A that were associated with its’ significantly reduced expression, and with lower immune cell infiltrates. We also identified low immune cell features in brain and liver metastases when compared to other metastatic sites, even within the same patient. These findings have direct implications for the treatment of metastatic breast cancer patients with immune- and HER2-targeting therapies and suggest potential novel therapeutic avenues for the improvement of outcomes for some patients with MBC.

Abstract Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values<5x10-8 as genome-wide significant, and p-values<1x10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values<1x105. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

Abstract BackgroundCopy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data.ResultsGene burden tests detected the strongest association for deletions in BRCA1 (P= 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P= 0.0008), ATM (P= 0.002) and BRCA2 (P= 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci.ConclusionsThis is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Abstract Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values<5x10-8 as genome-wide significant, and p-values<1x10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values<1x105. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

PURPOSE CALGB 40601 assessed whether dual versus single human epidermal growth factor receptor 2 (HER2) –targeting drugs added to neoadjuvant chemotherapy increased pathologic complete response (pCR). Here, we report relapse-free survival (RFS), overall survival (OS), and gene expression signatures that predict pCR and survival. PATIENTS AND METHODS Three hundred five women with untreated stage II and III HER2-positive breast cancer were randomly assigned to receive weekly paclitaxel combined with trastuzumab plus lapatinib (THL), trastuzumab (TH), or lapatinib (TL). The primary end point was pCR, and secondary end points included RFS, OS, and gene expression analyses. mRNA sequencing was performed on 264 pretreatment samples. RESULTS One hundred eighteen patients were randomly allocated to THL, 120 to TH, and 67 to TL. At more than 7 years of follow-up, THL had significantly better RFS and OS than did TH (RFS hazard ratio, 0.32; 95% CI, 0.14 to 0.71; P = .005; OS hazard ratio, 0.34; 95% CI, 0.12 to 0.94; P = .037), with no difference between TH and TL. Of 688 previously described gene expression signatures, significant associations were found in 215 with pCR, 45 with RFS, and only 22 with both pCR and RFS (3.2%). Specifically, eight immune signatures were significantly correlated with a higher pCR rate and better RFS. Among patients with residual disease, the immunoglobulin G signature was an independent, good prognostic factor, whereas the HER2-enriched signature, which was associated with a higher pCR rate, showed a significantly shorter RFS. CONCLUSION In CALGB 40601, dual HER2-targeting resulted in significant RFS and OS benefits. Integration of intrinsic subtype and immune signatures allowed for the prediction of pCR and RFS, both overall and within the residual disease group. These approaches may provide means for rational escalation and de-escalation treatment strategies in HER2-positive breast cancer.

Abstract We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast cancer.