García‐Closas, Montserrat

Background: Several breast cancer (BC) risk prediction models have been developed to provide personal risk assessments. Though individually validated, their performance has not been systematically evaluated across a wide range of populations or ages. Methods: We harmonized individual-level baseline questionnaire data and incident BC diagnoses from 21 cohorts from North America, Europe, and Australia participating in the Breast Cancer Risk Prediction Project. Five-year absolute risk of invasive BC was estimated for five established risk prediction models using classical risk factors only. Discrimination was evaluated by area under the curve (AUC). Calibration was assessed using average and risk-decile specific expected to observed (E/O) ratios. Performance metrics were meta-analyzed across cohorts and models. Metaregression tested associations between cohort characteristics and performance metrics. Results: This analysis included 1,595,977 women aged 20-75 years, enrolled in studies between 1976-2015, with 19,062 (1.2%) invasive BC cases ascertained within 5 years from exposure assessment. Age-adjusted AUCs were similar across models and cohorts (pooled AUCs by model: 0.57-0.58), while E/O ratios varied substantially (pooled E/O ratios by model: 0.83-1.25). Overestimation was common among predicted high-risk individuals (>3%). No appreciable differences in model performance by cohort age, birth year, race, and variable missingness emerged. Calibration improved after assigning race-specific incidence rates. Conclusion: Existing BC risk prediction models provided similar risk discrimination across multiple cohorts, although there was overestimation of risk for high-risk individuals. Performance variation across cohorts was not driven by specific characteristics, which supports development of a unified risk model for diverse populations that leverages appropriate incidence rates.

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 BRCA-associated homologous recombination deficiency (HRD) is present in ~ 50% of high-grade serous carcinomas (HGSC) and predicts sensitivity to platinum-based therapy. However, there is little understanding of why some patients with BRCA-deficient tumors experience unexpectedly poor outcomes. We profiled 154 tumors, enriched for patients with BRCA-deficient tumors that experienced short overall survival (≤ 3 years, n = 42), using whole-genome, transcriptome, and methylation analyses. All but one BRCA-deficient tumor exceeded an accepted HRD genomic scarring threshold. However, patients with BRCA1-deficient HGSC with a more elevated HRD score survived significantly longer. Patients with BRCA2-deficient HGSC and loss of NF1 survived twice as long as those without NF1 loss, whereas PIK3CA or RAD21 amplification defined BRCA2-deficient HGSC with exceptionally short survival. BRCA1-deficient tumors in short survivors had evidence of immunosuppressive c-kit signaling and EMT. In a large HGSC cohort (n = 1,389) including 282 individuals with pathogenic germline BRCA variants (gBRCApv), the location of the mutation within functional domains stratified clinical outcomes. Notably, residual disease after primary surgery had limited prognostic effect in gBRCApv-carriers compared to non-carriers. Our findings indicate that tumor HR proficiency in the context of therapy response and survival is not a binary property, and highlight genomic and immune modifiers of outcomes in BRCA-deficient HGSC.

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 Breast cancer risk assessment tools are widely used in clinical practice to guide decisions regarding screening timing and modality, lifestyle interventions, genetic testing, preventive therapy, and risk-reducing surgery. Although a number of tools designed for the general population are used in practice, they face a series of challenges including: (i) modest discriminatory ability due to lack of a unified model that flexibly incorporates a comprehensive set of risk factors; (ii) inability to produce subtype-specific risk estimates, relevant for tailoring screening and/or preventive strategies (e.g. prophylactic endocrine therapy that is effective only for hormone receptor positive tumors); (iii) lack of data from diverse populations to build models that work well for all women; and, (iv) scarcity of opportunities for comparative validation of different models. The Breast Cancer Risk Prediction Project (BCRPP), a collaboration currently including 20 prospective cohorts from the National Cancer Institute’s Cohort Consortium, is addressing these challenges by harmonizing, sharing, and analyzing data on over 1.8 million women, over 115, 000 of whom have developed breast cancer. The BCRPP is developing a comprehensive tool that will predict breast cancer risk, overall and by subtypes, across major racial and ethnic groups in the United States. This tool incorporates known and emerging risk factors (e.g. reproductive history, anthropometry, smoking, alcohol intake, medication use, mammography, family history, and measured genotypes). In parallel to model building and validation efforts, BCRPP investigators are developing privacy-preserving web and mobile interfaces to deliver the risk tools to patients and clinicians. Here we present descriptive data on participating cohorts (e.g. dates of enrollment, ages at enrollment, length of follow up, numbers of cases by subtype), summarize missing data patterns and our approach to model development and validation that accounts for heterogeneity in missingness patterns across cohorts, and describe data governance procedures and data sharing infrastructure (including pooled and federated analyses, to account for policy limitations on sharing individual data). Descriptive information and procedures for requesting and accessing data are available via the BCRPP data platform (https://epidataplatforms.cancer.gov/bcrpp/). Although originally assembled for risk model development, these data on a large and diverse sample of women with extensive harmonized risk factors and covariates can be a resource to address other outstanding questions in breast cancer epidemiology across time, the lifecourse, and diverse populations. Citation Format: Peter Kraft, Laura Beane-Freeman, Julie Palmer, Ellen O'Meara, Jeanine Genkinger, James Lacey, Thomas Rohan, Lauren Teras, Marian L. Neuhouser, Kala Visvanathan, Celine Vachon, Roger Milne, Christopher Haiman, Yu Chen, Heather Eliassen, Cari Kitahara, Katie O'Brien, Emily White, Garnet Anderson, I-Min Lee, Archie Campbell, Renee Fortner, Ylva Lagerros, Sven Sandin, Mia Gaudet, Montserrat Garcia-Closas, Nilanjan Chatterjee. The Breast Cancer Risk Prediction Project: a resource of 1.8 million women with diverse backgrounds from 20 prospective cohorts to improve breast cancer risk prediction and advance breast cancer epidemiology [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 3583.

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 Introduction Incorporation of mammographic density to breast cancer risk models could improve risk stratification to tailor screening and prevention strategies according to risk. Robust evaluation of the value of adding mammographic density to models with comprehensive information on questionnaire-based risk factors and polygenic risk score is needed to determine its effectiveness in improving risk stratification of such models. Methods We used the Individualized Coherent Absolute Risk Estimator (iCARE) tool for risk model building and validation to incorporate density to a previously validated literature-based model with questionnaire-based risk factors and a 313-variant polygenic risk score (PRS). The model was evaluated for calibration and discrimination in three prospective cohorts of European-ancestry women (1,468 cases, 19,104 controls): US-based Nurses’ Health Study (NHS I and II) and Mayo Mammography Health Study (MMHS); and Sweden-based Karolinska Mammography Project for Risk Prediction of Breast Cancer (KARMA) study. Analyses were done separately for women younger (NHS II, KARMA) and older than 50 years (NHS I, MMHS, KARMA). Improvements in terms of risk stratification and reclassification proportions were assessed among European-ancestry women aged 50–70 years in US and Sweden. Results For women younger and older than 50 years, the model with questionnaire-based risk factors, PRS and density was generally well calibrated across risk with some evidence of miscalibration at the extremes of the risk distribution. Incorporation of density led to modest improvements risk discrimination beyond the model with questionnaire-based risk factors and PRS: the area under the curve (AUC) among younger women was 67.0% (95% CI: 63.5–70.6%) vs. 65.6% (95% CI: 61.9–69.3%) for models with and without density; and 66.1% (95% CI 64.4–67.8%) vs. 65.5% (95% CI: 63.8–67.2%) among older women. The model with density identified 18.4% of US women 50–70 years old ≥ 3% 5-year predicted risk (threshold used for recommending risk-reducing medication in the US), with 42.4% of future cases expected to occur in this group. At this threshold, 7.9% of US women were reclassified by adding density to the model, resulting in the identification of 2.8% of additional future cases. The model with density identified 10.3% of Swedish women ≥ 3% 5-year predicted risk, with 29.4% of future cases expected to occur in this group. At this threshold, 5.3% of women were reclassified with the addition of density, leading to the identification of an additional 4.4% of future cases. Conclusion Integrating density with questionnaire-based risk factors and PRS could potentially identify more women of European-ancestry with elevated risk of breast cancer in the United States and Sweden. Further investigations of the integrated model in non-European ancestry populations are needed prior to considering clinical applications.