Troester, Melissa A.

Background: Discovering risk biomarkers in small benign breast disease (BBD) biopsies is constrained by scarce tissue and microanatomic heterogeneity of terminal duct lobular units (TDLUs). We tested whether tissue-sparing tissue microarray (TMA)–based Digital Spatial Profiling (DSP) can deliver reproducible, biologically coherent protein measurements across benign lobules and breast cancers (BCs), and how well DSP aligns with standard immunoassays. Methods: We performed a pilot using tissues from the Mayo Clinic BBD cohort using TMAs representing four contexts: terminal duct lobular units (TDLUs) from BBD biopsies preceding BC and matched BBD-controls, subsequent BCs, and BC-associated TDLUs. We profiled 79 proteins by DSP (37 retained after QC) and benchmarked against chromogenic IHC and OPAL immunofluorescence. Reproducibility was evaluated using intraclass correlation coefficients (ICCs), cross-platform agreement (weighted kappa), marker correlations, and mixed-effects models with false-discovery-rate (FDR) control. Results: We analyzed 368 BBD-TDLU cores (88 cases; 88 controls), 204 BC cores and 110 BC-associated TDLU cores. ICCs were highest in BC tissues, and lower in BC-associated TDLUs and BBD-TDLUs. Agreement was slight–to-fair in TDLUs but moderate (ER/PR) to substantial (BCL2) in BC. DSP recapitulated expected immunologic correlations (CD45 with T-cell, B-cell, and macrophage markers) and tissue-type gradients (BC > BC-associated TDLUs > BBD-TDLUs). Exploratory case–control differences in BBD-TDLUs did not persist after FDR control. Conclusions: TMA-based DSP is feasible in archival breast tissues and yields biologically coherent, cross-platform-benchmarked profiles that are particularly robust in BC, while conserving scarce TDLUS and clarifying current limits of single-marker risk stratification from benign lobules. These data provide a foundation for refined sampling and expanded panels in future TDLU-focused studies.

ImportanceEtiologic heterogeneity in breast carcinogenesis needs to be well characterized for targeted prevention. Associations between menopausal hormonal therapy (MHT) and oral contraceptive (OC) use and breast cancer intrinsic-like subtypes are not well understood.ObjectiveTo examine whether exogenous hormone use is differentially associated with breast cancer subtypes and to evaluate heterogeneity by intrinsic-like subtypes.Design, Setting, and ParticipantsThis study pooled data from 31 nested and population-based case-control studies involved in the Breast Cancer Association Consortium. The study population included individuals with breast cancer and control participants from 13 case-control studies nested in prospective cohorts (recruited between 1982 and 2011) and 18 population-based case-control studies (recruited between 1990 and 2013). Data analysis was performed in June 2024.ExposureMHT use (estrogen-progestin therapy [EPT] or estrogen-only therapy [ET]) in postmenopausal women and OC use in premenopausal women (never, past use, or current use).Main Outcomes and MeasuresBreast cancer intrinsic-like subtypes (luminal A–like, luminal B–like, luminal B–ERBB2 [formerly HER2 or HER2/neu]-like, ERBB2 enriched–like, or triple-negative) were determined by immunohistochemistry of tumor sections. Polytomous logistic regression was performed to estimate the association between exogenous hormones and risk of breast cancer by intrinsic-like subtypes. Analyses by subtypes were stratified by body mass index (BMI [calculated as weight in kilograms divided by height in meters squared]; healthy weight, 18.5-<25; overweight, 25-<30; or obesity, ≥30).ResultsThis study included 42 269 individuals with breast cancer (11 901 [28.2%] premenopausal and 30 368 [71.8%] postmenopausal; 23 353 [55.2%] had a known intrinsic-like subtype) and 71 072 control participants. The mean (SD) age of all participants was 57.9 (10.9) years. In postmenopausal women, associations between current MHT use (EPT or ET) and breast cancer differed by subtype. Current EPT users with healthy weight were more likely to be diagnosed with luminal A–like (odds ratio [OR], 2.51 [95% CI, 2.26-2.80]) or luminal B–ERBB2-like (OR, 1.95 [95% CI, 1.61-2.37]) subtypes. These associations were attenuated but remained for individuals with overweight (OR, 1.40 [95% CI, 1.02-1.92]) or obesity (OR, 1.68 [95% CI, 1.01-2.78]). EPT use increased the odds of being diagnosed with luminal B–like tumors solely in women with healthy weight (OR, 1.47 [95% CI, 1.17-1.86]). Current ET use was positively associated with luminal A–like disease in women with healthy weight only (OR, 1.16 [95% CI, 1.01-1.32]), showing inverse associations with higher BMI (obesity: OR, 0.65 [95% CI, 0.50-0.85]). In premenopausal women, recent OC use was associated with luminal B–ERBB2-like (OR, 1.50 [95% CI, 1.09-2.08]), ERBB2 enriched–like (OR, 2.33 [95% CI, 1.55-3.51]), and triple-negative (OR, 1.75 [95% CI, 1.33-2.29]; P < .04 for heterogeneity) tumors.Conclusions and RelevanceIn this study, clear differences were observed in associations between current EPT use and luminal-like breast cancer subtypes and other subtypes. EPT users with healthy weight were more likely to be diagnosed with luminal-like breast cancer compared with nonusers. Subtype heterogeneity was less apparent in associations of OC and ET use. Future studies on contemporary formulations, patterns of use, and routes of administration of exogenous hormone usage are warranted.

Abstract Background Most breast biopsies are diagnosed as benign breast disease, with 1.5- to 4-fold increased breast cancer risk. Apart from pathologic diagnoses of atypical hyperplasia, few factors aid in breast cancer risk assessment of these patients. We assessed whether a 313–single nucleotide variation (formerly single-nucleotide polymorphism) polygenic risk score stratifies risk of benign breast disease patients. Methods We pooled data from 5 Breast Cancer Association Consortium case-control studies (mean age = 59.9 years), including 6706 participants with breast cancer and 8488 participants without breast cancer. Using logistic regression, we estimated breast cancer risk associations by self-reported benign breast disease history and strata of polygenic risk score, with median polygenic risk score category among women without benign breast disease as the referent. We assessed interactions and mediation of benign breast disease and polygenic risk score with breast cancer risk. Results Benign breast disease history was associated with increased breast cancer risk (odds ratio [OR] = 1.48, 95% confidence interval [CI] = 1.37 to 1.60; P < .001). Polygenic risk score increased breast cancer risk, irrespective of benign breast disease history (Pinteraction = .48), with minimal evidence of mediation of either factor by the other. Women with benign breast disease and polygenic risk score in the highest tertile had more than twofold increased odds of breast cancer (OR = 2.73, 95% CI = 2.41 to 3.09), and those with benign breast disease and polygenic risk score in the lowest tertile experienced reduced breast cancer risk (OR = 0.79, 95% CI = 0.70 to 0.91) compared with the referent group. Women with benign breast disease and polygenic risk score in the highest decile had a 3.7-fold increase (95% CI = 3.00 to 4.61) compared with those with median polygenic risk score without benign breast disease. Conclusion Breast cancer risks are elevated among women with benign breast disease and increase progressively with polygenic risk score, suggesting that optimal combinations of these factors may improve risk stratification.

A large number of variants identified through clinical genetic testing in disease susceptibility genes are of uncertain significance (VUS). Following the recommendations of the American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP), the frequency in case-control datasets (PS4 criterion) can inform their interpretation. We present a novel case-control likelihood ratio-based method that incorporates gene-specific age-related penetrance. We demonstrate the utility of this method in the analysis of simulated and real datasets. In the analysis of simulated data, the likelihood ratio method was more powerful compared to other methods. Likelihood ratios were calculated for a case-control dataset of BRCA1 and BRCA2 variants from the Breast Cancer Association Consortium (BCAC) and compared with logistic regression results. A larger number of variants reached evidence in favor of pathogenicity, and a substantial number of variants had evidence against pathogenicity—findings that would not have been reached using other case-control analysis methods. Our novel method provides greater power to classify rare variants compared with classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and preformatted Excel calculators for implementation of the method for rare variants in BRCA1, BRCA2, and other high-risk genes with known penetrance.

Background Polygenic risk score (PRS), calculated based on genome-wide association studies (GWASs), can improve breast cancer (BC) risk assessment. To date, most BC GWASs have been performed in individuals of European (EUR) ancestry, and the generalisation of EUR-based PRS to other populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women. Methods We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel. Results In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28). Conclusions Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women.

Abstract Background Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC‐specific survival (BCSS) compared to non‐carriers. Aim To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Methods Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow‐up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi‐state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. Results There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55–0.78)]. No association was observed with radiotherapy. Results from the multi‐state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non‐carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09–1.56)]. Conclusion Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.

Abstract Over one million biopsies are performed annually in the U.S. to evaluate suspicious breast lesions. Most biopsies lead to a diagnosis of benign breast disease (BBD) and these women have a 1.5- to fourfold increase in subsequent breast cancer (BC) risk compared with the general population. A polygenic risk score (PRS) including 313 common variants was developed to predict BC risk in the general population and has improved performance when added to existing BC prediction models. However, the value of PRS in risk prediction models among women with BBD has not been evaluated. To evaluate the role of the established PRS, BBD, and BC risk, we pooled data from 15,194 white women (6,706 BC cases and 8,488 controls) participating in five Breast Cancer Association Consortium case-control studies. BC risk associations were evaluated by self-reported BBD history and combinations of PRS expressed in tertiles using odds ratios (OR), considering the lowest category of PRS among women without BBD as the referent. Associations were also assessed in an independent case-cohort study (831 BC cases and 1,845 controls) using hazard ratios (HR). Among the case-control studies, the PRS was strongly associated with BC risk (OR=2.50 per standard deviation increase), and associations were similar in women with (OR=2.60) and without a history of BBD (OR=2.44, p-interaction=0.23). Women with BBD have a higher risk of BC than those without BBD within each level of PRS. Compared to women without BBD and in the lowest PRS category, those with BBD and PRS in the highest tertile had 5-fold increased odds of BC (OR = 5.38; 95% CI 4.73-6.13). Associations were similar in an independent case-cohort study. These findings indicate that BBD history and PRS are independent contributors to BC risk, and consideration of both simultaneously could improve breast cancer risk prediction. Future studies are needed to integrate PRS with other risk measures to refine absolute BC risk prediction among BBD patients. Associations of BC risk with combinations of PRS (tertiles) and history of BBD in five case-control studies and a case-cohort study from the Breast Cancer Association Consortium. Case-control studies Case-cohort study PRS score range History of BBD n BC cases n controls OR (95% CI)1 n BC cases n controls HR (95% CI)2 Tertile 1 (≤ -0.57) No 1009 2692 1.00 (Ref) 117 584 1.00 (Ref) Yes 508 902 1.54 (1.35-1.76) 69 163 2.05 (1.45-2.89) Tertile 2 (-0.58, -0.03) No 1469 2001 1.97 (1.78-2.17) 181 475 1.93 (1.48-2.51) Yes 760 798 2.65 (2.34-3.01) 108 170 2.99 (2.18-4.09) Tertile 3 (≥ -0.02) No 1910 1551 3.32 (3.01-3.67) 214 336 3.28 (2.52-4.26) Yes 1050 544 5.38 (4.73-6.13) 142 117 4.83 4.25-8.02) 1.Logistic regression analysis, adjusted for study site/genomic platform and age at enrollment. OR, odds ratio; CI, confidence interval.2.Cox proportional hazards regression analysis, accounting for the case-cohort study design and adjusted for age at enrollment. HR, hazard ratio; CI, confidence interval. Citation Format: Stacey J. Winham, Mark Sherman, Robert Vierkant, Bryan McCauley, Christopher Scott, Mia Gaudet, Melissa Troester, Sandhya Pruthi, Derek Radisky, Amy Degnim, Kristan Aronson, Rachel Murphy, Pascal Guenel, Therese Truong, Jenny Chang-Claude, Heiko Becher, Montserrat Garcia-Closas, Stephen Chanock, Thomas Ahearn, Xiaohong Yang, Doug Easton, Manjeet Bolla, Celine Vachon, Breast Cancer Association Consortium. Association of polygenic risk scores with breast cancer risk among women with benign breast disease. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5241.

Abstract Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. We aimed to assess the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Analyses were based on 82,701 women diagnosed with invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations of treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR(95%CI): 0.66 (0.55–0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR(95%CI) :1.30 (1.09–1.56)]. In conclusion, systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk. (Main MS: 3201 words)

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 Background Mammographic breast density (MBD) decline post-tamoxifen initiation is a favorable prognostic factor in estrogen receptor (ER)–positive breast cancer (BC) and has potential utility as a biomarker of tamoxifen response. However, the prognostic value of MBD decline may vary by molecular characteristics among ER–positive patients. Methods We investigated associations between MBD decline (≥10% vs <10%) and breast cancer–specific mortality (BCSM) among ER–positive breast cancer patients aged 36-87 years at diagnosis treated with tamoxifen at Kaiser Permanente Northwest (1990-2008). Patients who died of BC (case patients; n = 62) were compared with those who did not (control patients; n = 215) overall and by tumor molecular characteristics (immunohistochemistry [IHC]–based subtype [luminal A–like: ER–positive/progesterone receptor [PR]–positive/HER2–negative/low Ki67; luminal B–like: ER–positive and 1 or more of PR–negative, HER2–positive, high Ki67] and modified IHC [mIHC]–based recurrence score of ER/PR/Ki67). Percent MBD was measured in the unaffected breast at baseline mammogram (mean = 6 months before tamoxifen initiation) and follow-up (mean = 12 months post-tamoxifen initiation). Adjusted odds ratios (ORs) and 95% confidence intervals (CIs) were computed from logistic regression models. All statistical tests were 2-sided. Results MBD decline was statistically significantly associated with reduced risk of BCSM overall (OR = 0.38, 95% CI = 0.15 to 0.92). This association was, however, stronger among women with aggressive tumor characteristics including luminal B–like (OR = 0.17, 95% CI = 0.04 to 0.73) vs A–like (OR = 0.74, 95% CI = 0.19 to 2.92); large (OR = 0.26, 95% CI = 0.08 to 0.78) vs small (OR = 0.41, 95% CI = 0.04 to 3.79) tumors; PR–negative (OR = 0.02, 95% CI = 0.001 to 0.37) vs PR–positive (OR = 0.50, 95% CI = 0.18 to 1.40) disease; and high (OR = 0.25, 95% CI = 0.07 to 0.93) vs low (OR = 0.44, 95% CI = 0.10 to 2.09) mIHC3 score. Conclusion The findings support MBD decline as a prognostic marker of tamoxifen response among patients with aggressive ER–positive BC phenotypes, for whom understanding treatment effectiveness is critical.