Love, Michael I.

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: Breast cancers in recently postpartum women may have worse outcomes, but studies examining tumor molecular features by pregnancy recency have shown conflicting results. Methods: This analysis used Carolina Breast Cancer Study data to examine clinical and molecular tumor features among women less than 50 years of age who were recently (≤10 years prior) or remotely (>10 years prior) postpartum, or nulliparous. Prevalence odds ratios (POR) and 95% confidence intervals (CI) were estimated using multivariable models. Results: Recently postpartum women (N = 618) were more frequently lymph node–positive [POR (95% CI): 1.66 (1.26–2.19)], estrogen receptor (ER)-negative [1.37 (1.02–1.83)], and IHC-based triple negative [1.57 (1.00–2.47)] compared with nulliparous (N = 360) women. Some differences were identified between recent versus remotely postpartum; smaller tumor size [0.67 (0.52–0.86)], p53 wildtype [0.53 (0.36–0.77)], and non–basal-like phenotype [0.53 (0.33–0.84)] were more common among recently postpartum. Recently postpartum (vs. nulliparous) had significant enrichment for adaptive immunity, T cells, B cells, CD8 T cells, activated CD8 T cells/natural killer (NK) cells, and T follicular helper (Tfh) cells and higher overall immune cell scores. These differences were attenuated in remotely (compared with recently) postpartum women. Conclusions: These results suggest a dominant effect of parity (vs. nulliparity) and a lesser effect of pregnancy recency on tumor molecular features, although tumor immune microenvironments were altered in association with pregnancy recency. Impact: Our study is unique in examining tumor immune microenvironment and RNA-based markers according to time since last childbirth. Future studies should examine the interplay between tumor features, postdiagnostic treatment, and outcomes among recently postpartum women. See related commentary by McDonald et al., p. 518