Nichola Johnson

Genome-wide association studies have led to the identification of more than 150 genomic regions that are associated with breast cancer risk. Translating these findings into a greater understanding of that risk requires identification of functional variants and target genes. Breast cancer progression and metastasis does not depend solely on cancer cell autonomous defects; the stroma, of which fibroblasts comprise a dominant component, also has a functional role. We generated promoter capture Hi-C data in primary and immortalized mammary fibroblasts and identified 28 interaction peaks involving 116 credible causal breast cancer variants and 26 target genes that were exclusive to fibroblasts. Integrating these data with H3K27ac CUT&Tag peaks identified a potentially functional variant (rs17393059) and target gene (filamin A interacting protein 1 like (FILIP1L)) at the 3q12.1 breast cancer risk locus. Using genome-wide functional data in breast-relevant cell types we demonstrate that perturbation of gene expression in mammary fibroblasts may impact risk of breast cancer by a cell non-autonomous mechanism.

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

Germline copy 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. Gene 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. To the best of our knowledge, this 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.

In this study we aim to examine gene–environment interactions (GxEs) between genes involved with estrogen metabolism and environmental factors related to estrogen exposure. GxE analyses were conducted with 1970 Korean breast cancer cases and 2052 controls in the case-control study, the Seoul Breast Cancer Study (SEBCS). A total of 11,555 SNPs from the 137 candidate genes were included in the GxE analyses with eight established environmental factors. A replication test was conducted by using an independent population from the Breast Cancer Association Consortium (BCAC), with 62,485 Europeans and 9047 Asians. The GxE tests were performed by using two-step methods in GxEScan software. Two interactions were found in the SEBCS. The first interaction was shown between rs13035764 of NCOA1 and age at menarche in the GE|2df model (p-2df = 1.2 × 10⁻³). The age at menarche before 14 years old was associated with the high risk of breast cancer, and the risk was higher when subjects had homozygous minor allele G. The second GxE was shown between rs851998 near ESR1 and height in the GE|2df model (p-2df = 1.1 × 10⁻⁴). Height taller than 160 cm was associated with a high risk of breast cancer, and the risk increased when the minor allele was added. The findings were not replicated in the BCAC. These results would suggest specificity in Koreans for breast cancer risk.

The ability to identify regulatory interactions that mediate gene expression changes through distal elements, such as risk loci, is transforming our understanding of how genomes are spatially organized and regulated. Capture Hi-C (CHi-C) is a powerful tool to delineate such regulatory interactions. However, primary analysis and downstream interpretation of CHi-C profiles remains challenging and relies on disparate tools with ad-hoc input/output formats and specific assumptions for statistical modeling. Here we present a data processing and interaction calling toolkit (CHiCANE), specialized for the analysis and meaningful interpretation of CHi-C assays. In this protocol, we demonstrate applications of CHiCANE to region capture Hi-C (rCHi-C) and promoter capture Hi-C (pCHi-C) libraries, followed by quality assessment of interaction peaks, as well as downstream analysis specific to rCHi-C and pCHi-C to aid functional interpretation. For a typical rCHi-C/pCHi-C dataset this protocol takes up to 3 d for users with a moderate understanding of R programming and statistical concepts, although this is dependent on dataset size and compute power available. CHiCANE is freely available at https://cran.r-project.org/web/packages/chicane.

Background: Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk. Methods: We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry. Results: For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (−49.2%, 95% CI −56.1% to −41.1%, P = 3.1 × 10^–18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (−26.7%, 95% CI −39.4% to −11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82–0.91, P = 6.9 × 10^–8). Conclusions: The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women.

Purpose Type 2 diabetes (T2D) has been reported to be associated with an elevated risk of breast cancer. It is unclear, however, whether this association is due to shared genetic factors. Methods We constructed a genetic risk score (GRS) using risk variants from 33 known independent T2D susceptibility loci and evaluated its relation to breast cancer risk using the data from two consortia, including 62,328 breast cancer patients and 83,817 controls of European ancestry. Unconditional logistic regression models were used to derive adjusted odds ratios (ORs) and 95 % confidence intervals (CIs) to measure the association of breast cancer risk with T2D GRS or T2D-associated genetic risk variants. Meta-analyses were conducted to obtain summary ORs across all studies. Results The T2D GRS was not found to be associated with breast cancer risk, overall, by menopausal status, or for estrogen receptor positive or negative breast cancer. Three T2D associated risk variants were individually associated with breast cancer risk after adjustment for multiple comparisons using the Bonferroni method (at p < 0.001), rs9939609 (FTO) (OR 0.94, 95 % CI = 0.92–0.95, p = 4.13E−13), rs7903146 (TCF7L2) (OR 1.04, 95 % CI = 1.02–1.06, p = 1.26E−05), and rs8042680 (PRC1) (OR 0.97, 95 % CI = 0.95–0.99, p = 8.05E−04). Conclusions We have shown that several genetic risk variants were associated with the risk of both T2D and breast cancer. However, overall genetic susceptibility to T2D may not be related to breast cancer risk.

Introduction: We have previously shown that a tag single nucleotide polymorphism (rs10235235), which maps to the CYP3A locus (7q22.1), was associated with a reduction in premenopausal urinary estrone glucuronide levels and a modest reduction in risk of breast cancer in women age ≤50 years. Methods: We further investigated the association of rs10235235 with breast cancer risk in a large case control study of 47,346 cases and 47,570 controls from 52 studies participating in the Breast Cancer Association Consortium. Genotyping of rs10235235 was conducted using a custom Illumina Infinium array. Stratified analyses were conducted to determine whether this association was modified by age at diagnosis, ethnicity, age at menarche or tumor characteristics. Results: We confirmed the association of rs10235235 with breast cancer risk for women of European ancestry but found no evidence that this association differed with age at diagnosis. Heterozygote and homozygote odds ratios (ORs) were OR = 0.98 (95% CI 0.94, 1.01; P = 0.2) and OR = 0.80 (95% CI 0.69, 0.93; P = 0.004), respectively (P_trend = 0.02). There was no evidence of effect modification by tumor characteristics. rs10235235 was, however, associated with age at menarche in controls (P_trend = 0.005) but not cases (P_trend = 0.97). Consequently the association between rs10235235 and breast cancer risk differed according to age at menarche (P_het = 0.02); the rare allele of rs10235235 was associated with a reduction in breast cancer risk for women who had their menarche age ≥15 years (OR_het = 0.84, 95% CI 0.75, 0.94; OR_hom = 0.81, 95% CI 0.51, 1.30; P_trend = 0.002) but not for those who had their menarche age ≤11 years (OR_het = 1.06, 95% CI 0.95, 1.19, OR_hom = 1.07, 95% CI 0.67, 1.72; P_trend = 0.29). Conclusions: To our knowledge rs10235235 is the first single nucleotide polymorphism to be associated with both breast cancer risk and age at menarche consistent with the well-documented association between later age at menarche and a reduction in breast cancer risk. These associations are likely mediated via an effect on circulating hormone levels.

Introduction: Several common breast cancer genetic susceptibility variants have recently been identified. We aimed to determine how these variants combine with a subset of other known risk factors to influence breast cancer risk in white women of European ancestry using case-control studies participating in the Breast Cancer Association Consortium. Methods: We evaluated two-way interactions between each of age at menarche, ever having had a live birth, number of live births, age at first birth and body mass index (BMI) and each of 12 single nucleotide polymorphisms (SNPs) (10q26-rs2981582 (FGFR2), 8q24-rs13281615, 11p15-rs3817198 (LSP1), 5q11-rs889312 (MAP3K1), 16q12-rs3803662 (TOX3), 2q35-rs13387042, 5p12-rs10941679 (MRPS30), 17q23-rs6504950 (COX11), 3p24-rs4973768 (SLC4A7), CASP8-rs17468277, TGFB1-rs1982073 and ESR1-rs3020314). Interactions were tested for by fitting logistic regression models including per-allele and linear trend main effects for SNPs and risk factors, respectively, and single-parameter interaction terms for linear departure from independent multiplicative effects. Results: These analyses were applied to data for up to 26,349 invasive breast cancer cases and up to 32,208 controls from 21 case-control studies. No statistical evidence of interaction was observed beyond that expected by chance. Analyses were repeated using data from 11 population-based studies, and results were very similar. Conclusions: The relative risks for breast cancer associated with the common susceptibility variants identified to date do not appear to vary across women with different reproductive histories or body mass index (BMI). The assumption of multiplicative combined effects for these established genetic and other risk factors in risk prediction models appears justified.

Introduction: Sex steroids, insulin-like growth factors (IGFs) and prolactin are breast cancer risk factors but whether their effects are mediated through mammographic density, one of the strongest risk factors for breast cancer, is unknown. If such a hormonal basis of mammographic density exists, hormones may underlie ethnic differences in both mammographic density and breast cancer incidence rates. Methods: In a cross-sectional study of 270 postmenopausal Caucasian and Afro-Caribbean women attending a population-based breast screening service in London, UK, we investigated whether plasma biomarkers (oestradiol, oestrone, sex hormone binding globulin (SHBG), testosterone, prolactin, leptin, IGF-I, IGF-II and IGF binding protein 3 (IGFBP3)) were related to and explained ethnic differences in mammographic percent density, dense area and nondense area, measured in Cumulus using the threshold method. Results: Mean levels of oestrogens, leptin and IGF-I:IGFBP3 were higher whereas SHBG and IGF-II:IGFBP3 were lower in Afro-Caribbean women compared with Caucasian women after adjustment for higher mean body mass index (BMI) in the former group (by 3.2 kg/m² (95% confidence interval (CI): 1.8, 4.5)). Age-adjusted percent density was lower in Afro-Caribbean compared with Caucasian women by 5.4% (absolute difference), but was attenuated to 2.5% (95% CI: -0.2, 5.1) upon BMI adjustment. Despite ethnic differences in biomarkers and in percent density, strong ethnic-age-adjusted inverse associations of oestradiol, leptin and testosterone with percent density were completely attenuated upon adjustment for BMI. There were no associations of IGF-I, IGF-II or IGFBP3 with percent density or dense area. We found weak evidence that a twofold increase in prolactin and oestrone levels were associated, respectively, with an increase (by 1.7% (95% CI: -0.3, 3.7)) and a decrease (by 2.0% (95% CI: 0, 4.1)) in density after adjustment for BMI. Conclusions: These findings suggest that sex hormone and IGF levels are not associated with BMI-adjusted percent mammographic density in cross-sectional analyses of postmenopausal women and thus do not explain ethnic differences in density. Mammographic density may still, however, be influenced by much higher premenopausal hormone levels.