Citations

Abstract Background Breast cancer is etiologically heterogeneous, but which risk factors differ in their associations across tumor subtypes remains unclear. We conducted a large, pooled analysis to evaluate independent, dose-response associations between breast cancer risk factors and quantitative tumor features. Methods Analyses of 15,731 invasive breast cancers from 24 studies evaluated associations (p-trend) between reproductive and hormonal factors, body mass index (BMI), alcohol, smoking, and family history in relation to quantitative immunohistochemistry measures on tissue microarrays (ER, PR, HER2, KI67, TP53) and tumor grade. Analyses in a subset of 10 population-based studies estimated subtype-specific odds ratios (ORs) comparing cases to controls. A Bayesian False Discovery Probability (BFDP) <0.2 was used to identify associations with strong statistical evidence. Results Nulliparity and later age at menopause were associated with higher ER-positivity (p-trend=0.021 and 0.001, respectively), with corresponding OR[ER+] (95% CI) = 1.49 (1.16-1.90) for nulliparous vs. parous and 1.07 (1.03-1.11) per 5 years. Current combined menopausal hormone therapy (MHT) use was associated with lower grade (p-trend<0.001), with OR [grade1] = 3.37 (2.69-4.21) for current vs. never users. Higher BMI was associated with lower ER-positivity and higher grade in premenopausal women (p-trend<0.001 and <0.001), with OR[ER+] = 0.80 (0.74-0.87) and OR[grade1] = 0.75 (0.63-0.88) per 5 units, and with higher PR-positivity and higher grade in postmenopausal women (p-trend<0.001 and <0.001), with OR[PR+] = 1.08 (1.03-1.14) and OR[grade 3] = 1.10 (1.04-1.17) per 5 units. Conclusion This pooled analysis of 15,731 cases showed that nulliparity, age at menopause, MHT, and BMI have independent, dose-response associations with ER, PR, and grade, clarifying patterns of etiologic heterogeneity. Associations with HER2, KI67 and TP53, or other risk factors did not meet our threshold for strong evidence.

Abstract Polygenic risk scores (PRSs), which quantify inherited susceptibility to complex traits and diseases, have emerged as valuable tools for risk stratification and precision medicine. Despite their promise, PRS developed on European cohorts often demonstrate substantially reduced predictive accuracy in non-European populations, due to differences in genetic architecture. The disproportionate representation of European ancestry cohorts in genome-wide association studies (GWAS) leads to inequitable deployment of PRS technologies across diverse populations. Here, we introduce PRANA (Polygenic Risk Adaptation via Neural-network Architecture), a deep learning framework that adapts an existing PRS developed on one population to other ancestries. Unlike methods that require large-scale GWAS in the target population, PRANA leverages pre-trained PRS models derived from European cohorts and adapts them using modestly sized cohorts from the target population. We evaluated PRANA on seven complex traits in South Asian, East Asian and Ashkenazi Jewish populations, as well as in selected smaller East Asian subpopulations where the scarcity of training data poses a particular challenge. PRANA mostly improved predictive performance of the baseline PRS models by 5%-20% in terms of effect size (β) and Nagelkerke’s R², and, in most cases, outperformed existing cross-ancestry multi-PRS approaches. These results highlight PRANA as a scalable and practical strategy to reduce disparities in genomic risk prediction and advance the equitable application of PRS in diverse populations.

Abstract Background Breast cancer genome-wide association studies (GWAS) have identified more than 200 susceptibility loci, but most studies are dominated by European and East Asian populations. Methods We analyzed breast cancer GWAS summary statistics from African (AFR), East Asian (EAS), European (EUR), and Hispanic/Latina (H/L) samples (159,297 cases and 212,102 controls). We estimated logit-scale SNP-based heritability, polygenicity, and cross-ancestry genetic correlation, partitioned heritability across functional annotations, and integrated GWAS results with the Tabula Sapiens single-cell atlas using scDRS+. Results The logit-scale heritability of breast cancer ranged from ℎ 2 =0.47 (SE = 0.07) in EAS to AFR ℎ 2 =0.61 (SE = 0.10), with no significant differences across ancestries (p=0.63). The estimated number of susceptibility markers in a sparse normal-mixture effects model also varied from 4,446 (SE = 3,100) in EAS to 8,308 (SE = 2,751) in AFR, but differences were not significant across ancestries (p=0.55). Cross-sample genetic correlations varied, with the strongest correlation between EUR and EAS (𝜌 = 0.79, SE = 0.08) and weakest between AFR and H/L (𝜌 = 0.26, SE = 0.24). Regulatory annotations were enriched for breast cancer heritability across samples. Integration with single-cell expression profiles implicated ancestry-shared associations with innate immune, secretory epithelial, and stromal cell types. Conclusion These results indicate substantial cross-ancestry sharing of breast cancer polygenic architecture, highlight a consistent contribution of regulatory variation, and identify convergent cellular contexts that motivate functional follow-up and inform expectations for the transferability and attainable performance of common-variant risk prediction across populations.

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 Introduction: Breast cancer genome-wide association studies (GWAS) have identified over 200 susceptibility loci, many replicated in diverse populations. However, cross-ancestry evaluation of breast cancer genetic architecture remains limited. We examined breast cancer genetic architecture using GWAS summary results from European (EUR; cases (ca) = 188,474, controls (co) = 96,201), East Asian (EAS; ca = 20,393, co = 86,329), African American (AA; ca = 9,235, co = 10,184), and US Hispanic/Latina and Latin American (H/L; ca = 2,396, co =7,468) studies. Methods: GWAS results were derived from the Breast Cancer Association Consortium, the African Ancestry Breast Cancer Genetics Consortium (AABCG), Biobank Japan, and a meta-analysis of five studies of H/L women with breast cancer. Linkage disequilibrium (LD) scores were generated for EUR, EAS, and AMR using the 1000 Genomes Project, while AA LD scores were derived from a subset of AABCG controls. Heritability was estimated for each ancestry group using LD score regression. Genetic correlations between populations were assessed via Popcorn. Polygenicity was analyzed using GENESIS; however, due to the limited sample size in the H/L studies, these estimates were restricted to AA, EAS, and EUR. Enrichment of heritability across regulatory elements was evaluated using stratified LD score regression across all populations. Results: The logit-scale heritability and corresponding standard error (SE) were 0.466 (0.066) for EAS, 0.501 (0.050) for EUR, 0.588 (0.360) for H/L, and 0.614 (0.095) for AA. The estimated number of independent susceptibility loci was 4,446 for EAS, 5,235 for EUR, and 8,308 for AA. Using a clumping and thresholding approach, an optimal set of common variants were projected to explain 38.6% (EUR), 39.4% (EAS), and 26.2% (AA) of genetic variance for samples of 100,000 cases and 100,000 controls with AUCs from corresponding estimated polygenic risk scores (PRSs) of 0.621 (EUR), 0.634 (EAS), and 0.611 (AA). Genetic correlations were strongest between EUR and EAS (ρ = 0.79, SE = 0.08) and EUR and H/L (ρ = 0.68, SE = 0.21), and weakest between AA and EAS (ρ = 0.42, SE = 0.14) and AA and H/L (ρ = 0.26, SE = 0.24). Among 73 genomic features, we found significant (P < 0.05/73) enrichment in heritability in EUR for ‘ancient promoter’ regions, transcription factor binding sites, H3K4me3, ‘super-enhancers’, H3K4me1, H3K27ac; in EAS for ‘super-enhancers’; and in AA for H3K27ac. No genomic features were significantly enriched in H/L, likely due to limited power. Conclusion: These findings suggest a shared breast cancer genetic architecture across diverse populations, as well as the potential for a similar level of breast cancer risk stratification by PRS in these populations. Expanding GWAS in underrepresented populations is essential to improve genetic risk predictions and foster equitable cancer prevention. Citation Format: James L. Li, Maria Zanti, Jacob Williams, Om Jahagirdar, Guochong Jia, Qiang Hu, Jean-Tristan Brandenburg, Li Yan, Weang-Kee Ho, Jingmei Li, José P. Miranda, Devika Godbole, Julie-Alexia Dias, Leila Dorling, Wenlong C. Chen, Nicholas Boddicker, Ying Wang, Alicia Martin, Martin J. Zhang, Yan Zhang, Joe Dennis, Esther M. John, Gabriela Torres-Mejia, Larry Kushi, Jeffrey Weitzel, Susan L. Neuhausen, Luis Carvajal-Carmona, Christopher Haiman, Elad Ziv, Laura Fejerman, Wei Zheng, Dezheng Huo, Douglas Easton, Nilanjan Chatterjee, Peter Kraft, Montserrat Garcia-Closas, Wendy Wong, Kyriaki Michailidou, Qianqian Zhu, Diptavo Dutta, Thomas U. Ahearn, Haoyu Zhang. Genetic architecture of breast cancer across diverse populations: Assessing heritability, genetic correlation, and polygenicity [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 2282.

Abstract Studies have identified genetic and epidemiologic factors associated with mammographic density (MD) phenotypes. However, MD-associated genetic variants only account for a small proportion of the total estimated heritability. Interrogating interactions between genetic and epidemiologic factors could potentially identify additional MD-associated loci, expand our understanding of the genetic basis of MD phenotypes, and clarify how epidemiologic factors modulate relationships between genetic variants and MD. We conducted six separate genome-wide, gene–environment (GxE) interaction analyses, applying 2 degrees of freedom (df) and 1df interaction tests, for each of three MD phenotypes (percent density, dense area (DA), and nondense area (NDA)). The six epidemiologic factors considered were height, ever parous, parity, ever menopausal hormone therapy, ever breastfeeding, and months of breastfeeding. We included European ancestry participants from multiple studies within the Markers of Density consortium and the Breast Cancer Association Consortium (n = 4895-16 218 depending on specific analyses). We identified 11 loci with genome-wide significant (P < 5 × 10−8) interaction tests including two novel common genetic signals interacting with parity (8p21.2) and ever breastfeeding (19p13.2) for NDA. Our results suggest that epidemiologic risk factors might influence relationships between common genetic variants and MD phenotypes at particular genomic loci.

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.