Stone, Jennifer

Background: Breast cancer polygenic risk scores (PRS) and traditional risk models (e.g., the Gail model [Gail]) are known to contribute largely independent information, but it is unclear how the overlap varies by ancestry, age, disease type (invasive breast cancer, DCIS), and risk threshold. Methods: In a retrospective case–control study, we evaluated risk prediction performance in 180,398 women (161,849 of European ancestry; 18,549 of Asian ancestry). Odds ratios (ORs) from logistic regression models and the area under the receiver operating characteristic curve (AUC) were estimated. Results: PRS for invasive disease showed a stronger association in younger (<50 years) women (OR = 2.51, AUC = 0.622) than in women ≥ 50 years (OR = 2.06, AUC = 0.653) of European ancestry. PRS performance in Asians was lower (OR range = 1.62–1.64, AUC = 0.551–0.600). Gail performance was modest across groups and poor in younger Asian women (OR = 0.94–0.99, AUC = 0.523–0.533). Age interactions were observed for both PRS (p < 0.001) and Gail (p < 0.001) in Europeans, whereas in Asians, age interaction was observed only for Gail (invasive: p < 0.001; DCIS: p = 0.002). PRS identified more high-risk individuals than Gail in Asian populations, especially ≥50 years, while Gail identified more in Europeans. Overlap between PRS, Gail, and family history was limited at higher thresholds. Calibration analysis, comparing empirical and model-based ROC curves, showed divergence for both PRS and Gail (p < 0.001), which indicates miscalibration. In Europeans, family history and prior biopsies drove Gail discrimination. In younger Asians, age at first live birth was influential. Conclusions: PRS adds value to risk stratification beyond traditional tools, especially in younger women and Asian ancestry 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-&amp;lt;25; overweight, 25-&amp;lt;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 &amp;lt; .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 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 &lt; 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.

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.

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 &lt; .001 for triple-negative compared with luminal A-like breast cancer) and breastfeeding (Pheterogeneity &lt; .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.

Genome-wide association studies (GWAS) have identified more than 200 susceptibility loci for breast cancer, but these variants explain less than a fifth of the disease risk. Although gene–environment interactions have been proposed to account for some of the remaining heritability, few studies have empirically assessed this. We obtained genotype and risk factor data from 46,060 cases and 47,929 controls of European ancestry from population-based studies within the Breast Cancer Association Consortium (BCAC). We built gene expression prediction models for 4,864 genes with a significant (P &lt; 0.01) heritable component using the transcriptome and genotype data from the Genotype-Tissue Expression (GTEx) project. We leveraged predicted gene expression information to investigate the interactions between gene-centric genetic variation and 14 established risk factors in association with breast cancer risk, using a mixed-effects score test. After adjusting for number of tests using Bonferroni correction, no interaction remained statistically significant. The strongest interaction observed was between the predicted expression of the C13orf45 gene and age at first full-term pregnancy (PGXE = 4.44 × 10−6). In this transcriptome-informed genome-wide gene–environment interaction study of breast cancer, we found no strong support for the role of gene expression in modifying the associations between established risk factors and breast cancer risk. Our study suggests a limited role of gene–environment interactions in breast cancer risk.

PURPOSEDiagnosis of AIDS lymphoma in low-resource settings, like South Africa, is often delayed, leaving patients with limited treatment options. In tuberculosis (TB) endemic regions, overlapping signs and symptoms often lead to diagnostic delays. Assessment of plasma cell-free DNA (cfDNA) by next-generation sequencing (NGS) may expedite the diagnosis of lymphoma but requires high-quality cfDNA.METHODSPeople living with HIV with newly diagnosed aggressive B-cell lymphoma and those with newly diagnosed TB seeking care at Chris Hani Baragwanath Academic Hospital and its surrounding clinics, in Soweto, South Africa, were enrolled in this study. Each participant provided a whole blood specimen collected in cell-stabilizing tubes. Quantity and quality of plasma cfDNA were assessed. NGS of the immunoglobulin heavy chain was performed.RESULTSNine HIV+ patients with untreated lymphoma and eight HIV+ patients with TB, but without lymphoma, were enrolled. All cfDNA quantity and quality metrics were similar between the two groups, except that cfDNA accounted for a larger fraction of recovered plasma DNA in patients with lymphoma. The concentration of cfDNA in plasma also trended higher in patients with lymphoma. NGS of immunoglobulin heavy chain showed robust amplification of DNA, with large amplicons (> 250 bp) being more readily detected in patients with lymphoma. Clonal sequences were detected in five of nine patients with lymphoma, and none of the patients with TB.CONCLUSIONThis proof-of-principle study demonstrates that whole blood collected for cfDNA in a low-resource setting is suitable for sophisticated sequencing analyses, including clonal immunoglobulin NGS. The detection of clonal sequences in more than half of patients with lymphoma shows promise as a diagnostic marker that may be explored in future studies.

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.