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.

Part of the substantial unexplained familial aggregation of breast cancer may be due to interactions between common variants, but few studies have had adequate statistical power to detect interactions of realistic magnitude. We aimed to assess all two-way interactions in breast cancer susceptibility between 70,917 single nucleotide polymorphisms (SNPs) selected primarily based on prior evidence of a marginal effect. Thirty-eight international studies contributed data for 46,450 breast cancer cases and 42,461 controls of European origin as part of a multi-consortium project (COGS). First, SNPs were preselected based on evidence (P < 0.01) of a per-allele main effect, and all two-way combinations of those were evaluated by a per-allele (1 d.f.) test for interaction using logistic regression. Second, all 2.5 billion possible two-SNP combinations were evaluated using Boolean operation-based screening and testing, and SNP pairs with the strongest evidence of interaction (P < 10(-4)) were selected for more careful assessment by logistic regression. Under the first approach, 3277 SNPs were preselected, but an evaluation of all possible two-SNP combinations (1 d.f.) identified no interactions at P < 10(-8). Results from the second analytic approach were consistent with those from the first (P > 10(-10)). In summary, we observed little evidence of two-way SNP interactions in breast cancer susceptibility, despite the large number of SNPs with potential marginal effects considered and the very large sample size. This finding may have important implications for risk prediction, simplifying the modelling required. Further comprehensive, large-scale genome-wide interaction studies may identify novel interacting loci if the inherent logistic and computational challenges can be overcome.

Part of the substantial unexplained familial aggregation of breast cancer may be due to interactions between common variants, but few studies have had adequate statistical power to detect interactions of realistic magnitude. We aimed to assess all two-way interactions in breast cancer susceptibility between 70,917 single nucleotide polymorphisms (SNPs) selected primarily based on prior evidence of a marginal effect. Thirty-eight international studies contributed data for 46,450 breast cancer cases and 42,461 controls of European origin as part of a multi-consortium project (COGS). First, SNPs were preselected based on evidence (P < 0.01) of a per-allele main effect, and all two-way combinations of those were evaluated by a per-allele (1 d.f.) test for interaction using logistic regression. Second, all 2.5 billion possible two-SNP combinations were evaluated using Boolean operation-based screening and testing, and SNP pairs with the strongest evidence of interaction (P < 10(-4)) were selected for more careful assessment by logistic regression. Under the first approach, 3277 SNPs were preselected, but an evaluation of all possible two-SNP combinations (1 d.f.) identified no interactions at P < 10(-8). Results from the second analytic approach were consistent with those from the first (P > 10(-10)). In summary, we observed little evidence of two-way SNP interactions in breast cancer susceptibility, despite the large number of SNPs with potential marginal effects considered and the very large sample size. This finding may have important implications for risk prediction, simplifying the modelling required. Further comprehensive, large-scale genome-wide interaction studies may identify novel interacting loci if the inherent logistic and computational challenges can be overcome.

Background

Despite a modest association between tobacco smoking and breast cancer risk reported by recent epidemiological studies, it is still equivocal whether smoking is causally related to breast cancer risk.

Methods

We applied Mendelian randomisation (MR) to evaluate a potential causal effect of cigarette smoking on breast cancer risk. Both individual-level data as well as summary statistics for 164 single-nucleotide polymorphisms (SNPs) reported in genome-wide association studies of lifetime smoking index (LSI) or cigarette per day (CPD) were used to obtain MR effect estimates. Data from 108,420 invasive breast cancer cases and 87,681 controls were used for the LSI analysis and for the CPD analysis conducted among ever-smokers from 26,147 cancer cases and 26,072 controls. Sensitivity analyses were conducted to address pleiotropy.

Results

Genetically predicted LSI was associated with increased breast cancer risk (OR 1.18 per SD, 95% CI: 1.07-1.30, P = 0.11 × 10^-2), but there was no evidence of association for genetically predicted CPD (OR 1.02, 95% CI: 0.78-1.19, P = 0.85). The sensitivity analyses yielded similar results and showed no strong evidence of pleiotropic effect.

Conclusion

Our MR study provides supportive evidence for a potential causal association with breast cancer risk for lifetime smoking exposure but not cigarettes per day among smokers.

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10^-31).

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.

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 (P _adj > 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5-25 kg/m² [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 <20 years at first pregnancy [0.79 (0.72-0.86)]; >0-<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.

Background

Our recent genome-wide association study identified a novel breast cancer susceptibility locus at 9q31.2 (rs865686).

Methods

To further investigate the rs865686-breast cancer association, we conducted a replication study within the Breast Cancer Association Consortium, which comprises 37 case-control studies (48,394 cases, 50,836 controls).

Results

This replication study provides additional strong evidence of an inverse association between rs865686 and breast cancer risk [study-adjusted per G-allele OR, 0.90; 95% confidence interval (CI), 0.88; 0.91, P = 2.01 × 10(-29)] among women of European ancestry. There were ethnic differences in the estimated minor (G)-allele frequency among controls [0.09, 0.30, and 0.38 among, respectively, Asians, Eastern Europeans, and other Europeans; P for heterogeneity (P(het)) = 1.3 × 10(-143)], but no evidence of ethnic differences in per allele OR (P(het) = 0.43). rs865686 was associated with estrogen receptor-positive (ER(+)) disease (per G-allele OR, 0.89; 95% CI, 0.86-0.91; P = 3.13 × 10(-22)) but less strongly, if at all, with ER-negative (ER(-)) disease (OR, 0.98; 95% CI, 0.94-1.02; P = 0.26; P(het) = 1.16 × 10(-6)), with no evidence of independent heterogeneity by progesterone receptor or HER2 status. The strength of the breast cancer association decreased with increasing age at diagnosis, with case-only analysis showing a trend in the number of copies of the G allele with increasing age at diagnosis (P for linear trend = 0.0095), but only among women with ER(+) tumors.

Conclusions

This study is the first to show that rs865686 is a susceptibility marker for ER(+) breast cancer.

Impact

The findings further support the view that genetic susceptibility varies according to tumor subtype.

Background

Given the high heterogeneity among breast tumors, associations between common germline genetic variants and survival that may exist within specific subgroups could go undetected in an unstratified set of breast cancer patients.

Methods

We performed genome-wide association analyses within 15 subgroups of breast cancer patients based on prognostic factors, including hormone receptors, tumor grade, age, and type of systemic treatment. Analyses were based on 91,686 female patients of European ancestry from the Breast Cancer Association Consortium, including 7531 breast cancer-specific deaths over a median follow-up of 8.1 years. Cox regression was used to assess associations of common germline variants with 15-year and 5-year breast cancer-specific survival. We assessed the probability of these associations being true positives via the Bayesian false discovery probability (BFDP < 0.15).

Results

Evidence of associations with breast cancer-specific survival was observed in three patient subgroups, with variant rs5934618 in patients with grade 3 tumors (15-year-hazard ratio (HR) [95% confidence interval (CI)] 1.32 [1.20, 1.45], P = 1.4E-08, BFDP = 0.01, per G allele); variant rs4679741 in patients with ER-positive tumors treated with endocrine therapy (15-year-HR [95% CI] 1.18 [1.11, 1.26], P = 1.6E-07, BFDP = 0.09, per G allele); variants rs1106333 (15-year-HR [95% CI] 1.68 [1.39,2.03], P = 5.6E-08, BFDP = 0.12, per A allele) and rs78754389 (5-year-HR [95% CI] 1.79 [1.46,2.20], P = 1.7E-08, BFDP = 0.07, per A allele), in patients with ER-negative tumors treated with chemotherapy.

Conclusions

We found evidence of four loci associated with breast cancer-specific survival within three patient subgroups. There was limited evidence for the existence of associations in other patient subgroups. However, the power for many subgroups is limited due to the low number of events. Even so, our results suggest that the impact of common germline genetic variants on breast cancer-specific survival might be limited.

Background

Previous studies have suggested that breast cancer risk factors are associated with estrogen receptor (ER) and progesterone receptor (PR) expression status of the tumors.

Methods

We pooled tumor marker and epidemiological risk factor data from 35,568 invasive breast cancer case patients from 34 studies participating in the Breast Cancer Association Consortium. Logistic regression models were used in case-case analyses to estimate associations between epidemiological risk factors and tumor subtypes, and case-control analyses to estimate associations between epidemiological risk factors and the risk of developing specific tumor subtypes in 12 population-based studies. All statistical tests were two-sided.

Results

In case-case analyses, of the epidemiological risk factors examined, early age at menarche (≤12 years) was less frequent in case patients with PR(-) than PR(+) tumors (P = .001). Nulliparity (P = 3 × 10(-6)) and increasing age at first birth (P = 2 × 10(-9)) were less frequent in ER(-) than in ER(+) tumors. Obesity (body mass index [BMI] ≥ 30 kg/m(2)) in younger women (≤50 years) was more frequent in ER(-)/PR(-) than in ER(+)/PR(+) tumors (P = 1 × 10(-7)), whereas obesity in older women (>50 years) was less frequent in PR(-) than in PR(+) tumors (P = 6 × 10(-4)). The triple-negative (ER(-)/PR(-)/HER2(-)) or core basal phenotype (CBP; triple-negative and cytokeratins [CK]5/6(+) and/or epidermal growth factor receptor [EGFR](+)) accounted for much of the heterogeneity in parity-related variables and BMI in younger women. Case-control analyses showed that nulliparity, increasing age at first birth, and obesity in younger women showed the expected associations with the risk of ER(+) or PR(+) tumors but not triple-negative (nulliparity vs parity, odds ratio [OR] = 0.94, 95% confidence interval [CI] = 0.75 to 1.19, P = .61; 5-year increase in age at first full-term birth, OR = 0.95, 95% CI = 0.86 to 1.05, P = .34; obesity in younger women, OR = 1.36, 95% CI = 0.95 to 1.94, P = .09) or CBP tumors.

Conclusions

This study shows that reproductive factors and BMI are most clearly associated with hormone receptor-positive tumors and suggest that triple-negative or CBP tumors may have distinct etiology.