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.

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.

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.

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 < 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.

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.

Abstract Background: It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype. Methods: We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer–specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype. Results: There was no evidence of heterogeneous associations between risk factors and mortality by subtype (Padj &gt; 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5–25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06–1.34)]; current versus never smoking [1.37 (1.27–1.47)], high versus low physical activity [0.43 (0.21–0.86)], age ≥30 years versus &lt;20 years at first pregnancy [0.79 (0.72–0.86)]; &gt;0–&lt;5 years versus ≥10 years since last full-term birth [1.31 (1.11–1.55)]; ever versus never use of oral contraceptives [0.91 (0.87–0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen–progestin therapy [0.61 (0.54–0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02–1.21) for current versus never smoking. Conclusions: We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype. Impact: Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.

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.

AbstractPrevious transcriptome‐wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome‐wide association studies (GWAS), but analyses of breast cancer subtype‐specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta‐analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER−). We further compared associations with ER+ and ER− subtypes, using a case‐only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER– breast cancer. We further identified 30 TWAS‐significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast‐cancer gene in three of six regions harboring multiple TWAS‐significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER− breast cancer.