Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

This corrects the article DOI: 10.1038/ncomms5999.

PALB2, CHEK2 and ATM rare variants and cancer risk: data from COGS

BACKGROUND

The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study.

METHODS

We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant.

RESULTS

For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10-5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10-8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants.

CONCLUSIONS

This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

A polymorphism in the base excision repair gene PARP2 is associated with differential prognosis by chemotherapy among postmenopausal breast cancer patients

BACKGROUND

Personalized therapy considering clinical and genetic patient characteristics will further improve breast cancer survival. Two widely used treatments, chemotherapy and radiotherapy, can induce oxidative DNA damage and, if not repaired, cell death. Since base excision repair (BER) activity is specific for oxidative DNA damage, we hypothesized that germline genetic variation in this pathway will affect breast cancer-specific survival depending on treatment.

METHODS

We assessed in 1,408 postmenopausal breast cancer patients from the German MARIE study whether cancer specific survival after adjuvant chemotherapy, anthracycline chemotherapy, and radiotherapy is modulated by 127 Single Nucleotide Polymorphisms (SNPs) in 21 BER genes. For SNPs with interaction terms showing p<0.1 (likelihood ratio test) using multivariable Cox proportional hazard analyses, replication in 6,392 patients from nine studies of the Breast Cancer Association Consortium (BCAC) was performed.

RESULTS

rs878156 in PARP2 showed a differential effect by chemotherapy (p=0.093) and was replicated in BCAC studies (p=0.009; combined analysis p=0.002). Compared to non-carriers, carriers of the variant G allele (minor allele frequency=0.07) showed better survival after chemotherapy (combined allelic hazard ratio (HR)=0.75, 95% 0.53-1.07) and poorer survival when not treated with chemotherapy (HR=1.42, 95% 1.08-1.85). A similar effect modification by rs878156 was observed for anthracycline-based chemotherapy in both MARIE and BCAC, with improved survival in carriers (combined allelic HR=0.73, 95% CI 0.40-1.32). None of the SNPs showed significant differential effects by radiotherapy.

CONCLUSIONS

Our data suggest for the first time that a SNP in PARP2, rs878156, may together with other genetic variants modulate cancer specific survival in breast cancer patients depending on chemotherapy. These germline SNPs could contribute towards the design of predictive tests for breast cancer patients.

Fine-mapping identifies two additional breast cancer susceptibility loci at 9q31.2

We recently identified a novel susceptibility variant, rs865686, for estrogen-receptor positive breast cancer at 9q31.2. Here, we report a fine-mapping analysis of the 9q31.2 susceptibility locus using 43 160 cases and 42 600 controls of European ancestry ascertained from 52 studies and a further 5795 cases and 6624 controls of Asian ancestry from nine studies. Single nucleotide polymorphism (SNP) rs676256 was most strongly associated with risk in Europeans (odds ratios [OR] = 0.90 [0.88-0.92]; P-value = 1.58 × 10(-25)). This SNP is one of a cluster of highly correlated variants, including rs865686, that spans ∼14.5 kb. We identified two additional independent association signals demarcated by SNPs rs10816625 (OR = 1.12 [1.08-1.17]; P-value = 7.89 × 10(-09)) and rs13294895 (OR = 1.09 [1.06-1.12]; P-value = 2.97 × 10(-11)). SNP rs10816625, but not rs13294895, was also associated with risk of breast cancer in Asian individuals (OR = 1.12 [1.06-1.18]; P-value = 2.77 × 10(-05)). Functional genomic annotation using data derived from breast cancer cell-line models indicates that these SNPs localise to putative enhancer elements that bind known drivers of hormone-dependent breast cancer, including ER-α, FOXA1 and GATA-3. In vitro analyses indicate that rs10816625 and rs13294895 have allele-specific effects on enhancer activity and suggest chromatin interactions with the KLF4 gene locus. These results demonstrate the power of dense genotyping in large studies to identify independent susceptibility variants. Analysis of associations using subjects with different ancestry, combined with bioinformatic and genomic characterisation, can provide strong evidence for the likely causative alleles and their functional basis.

Prediction of breast cancer risk based on profiling with common genetic variants

BACKGROUND

Data for multiple common susceptibility alleles for breast cancer may be combined to identify women at different levels of breast cancer risk. Such stratification could guide preventive and screening strategies. However, empirical evidence for genetic risk stratification is lacking.

METHODS

We investigated the value of using 77 breast cancer-associated single nucleotide polymorphisms (SNPs) for risk stratification, in a study of 33 673 breast cancer cases and 33 381 control women of European origin. We tested all possible pair-wise multiplicative interactions and constructed a 77-SNP polygenic risk score (PRS) for breast cancer overall and by estrogen receptor (ER) status. Absolute risks of breast cancer by PRS were derived from relative risk estimates and UK incidence and mortality rates.

RESULTS

There was no strong evidence for departure from a multiplicative model for any SNP pair. Women in the highest 1% of the PRS had a three-fold increased risk of developing breast cancer compared with women in the middle quintile (odds ratio [OR] = 3.36, 95% confidence interval [CI] = 2.95 to 3.83). The ORs for ER-positive and ER-negative disease were 3.73 (95% CI = 3.24 to 4.30) and 2.80 (95% CI = 2.26 to 3.46), respectively. Lifetime risk of breast cancer for women in the lowest and highest quintiles of the PRS were 5.2% and 16.6% for a woman without family history, and 8.6% and 24.4% for a woman with a first-degree family history of breast cancer.

CONCLUSIONS

The PRS stratifies breast cancer risk in women both with and without a family history of breast cancer. The observed level of risk discrimination could inform targeted screening and prevention strategies. Further discrimination may be achievable through combining the PRS with lifestyle/environmental factors, although these were not considered in this report.

Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

GWAS have identified a breast cancer susceptibility locus on 2q35. Here we report the fine mapping of this locus using data from 101,943 subjects from 50 case-control studies. We genotype 276 SNPs using the 'iCOGS' genotyping array and impute genotypes for a further 1,284 using 1000 Genomes Project data. All but two, strongly correlated SNPs (rs4442975 G/T and rs6721996 G/A) are excluded as candidate causal variants at odds against >100:1. The best functional candidate, rs4442975, is associated with oestrogen receptor positive (ER+) disease with an odds ratio (OR) in Europeans of 0.85 (95% confidence interval=0.84-0.87; P=1.7 × 10(-43)) per t-allele. This SNP flanks a transcriptional enhancer that physically interacts with the promoter of IGFBP5 (encoding insulin-like growth factor-binding protein 5) and displays allele-specific gene expression, FOXA1 binding and chromatin looping. Evidence suggests that the g-allele confers increased breast cancer susceptibility through relative downregulation of IGFBP5, a gene with known roles in breast cell biology.

Genetic variation at CYP3A is associated with age at menarche and breast cancer risk: a case-control study

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.

A large-scale assessment of two-way SNP interactions in breast cancer susceptibility using 46,450 cases and 42,461 controls from the breast cancer association consortium

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.

Fine-scale mapping of the FGFR2 breast cancer risk locus: putative functional variants differentially bind FOXA1 and E2F1

The 10q26 locus in the second intron of FGFR2 is the locus most strongly associated with estrogen-receptor-positive breast cancer in genome-wide association studies. We conducted fine-scale mapping in case-control studies genotyped with a custom chip (iCOGS), comprising 41 studies (n = 89,050) of European ancestry, 9 Asian ancestry studies (n = 13,983), and 2 African ancestry studies (n = 2,028) from the Breast Cancer Association Consortium. We identified three statistically independent risk signals within the locus. Within risk signals 1 and 3, genetic analysis identified five and two variants, respectively, highly correlated with the most strongly associated SNPs. By using a combination of genetic fine mapping, data on DNase hypersensitivity, and electrophoretic mobility shift assays to study protein-DNA binding, we identified rs35054928, rs2981578, and rs45631563 as putative functional SNPs. Chromatin immunoprecipitation showed that FOXA1 preferentially bound to the risk-associated allele (C) of rs2981578 and was able to recruit ERα to this site in an allele-specific manner, whereas E2F1 preferentially bound the risk variant of rs35054928. The risk alleles were preferentially found in open chromatin and bound by Ser5 phosphorylated RNA polymerase II, suggesting that the risk alleles are associated with changes in transcription. Chromatin conformation capture demonstrated that the risk region was able to interact with the promoter of FGFR2, the likely target gene of this risk region. A role for FOXA1 in mediating breast cancer susceptibility at this locus is consistent with the finding that the FGFR2 risk locus primarily predisposes to estrogen-receptor-positive disease.

The 5p12 breast cancer susceptibility locus affects MRPS30 expression in estrogen-receptor positive tumors

Genome-wide association studies have identified numerous loci linked to breast cancer susceptibility, but the mechanism by which variations at these loci influence susceptibility is usually unknown. Some variants are only associated with particular clinical subtypes of breast cancer. Understanding how and why these variants influence subtype-specific cancer risk contributes to our understanding of cancer etiology. We conducted a genome-wide expression Quantitative Trait Locus (eQTL) study in a discovery set of 287 breast tumors and 97 normal mammary tissue samples and a replication set of 235 breast tumors. We found that the risk-associated allele of rs7716600 in the 5p12 estrogen receptor-positive (ER-positive) susceptibility locus was associated with elevated expression of the nearby gene MRPS30 exclusively in ER-positive tumors. We replicated this finding in 235 independent tumors. Further, we showed the rs7716600 risk genotype was associated with decreased MRPS30 promoter methylation exclusively in ER-positive breast tumors. In vitro studies in MCF-7 cells carrying the protective genotype showed that estrogen stimulation decreased MRPS30 promoter chromatin availability and mRNA levels. In contrast, in 600MPE cells carrying the risk genotype, estrogen increased MRPS30 expression and did not affect promoter availability. Our data suggest the 5p12 risk allele affects MRPS30 expression in estrogen-responsive tumor cells after tumor initiation by a mechanism affecting chromatin availability. These studies emphasize that the genetic architecture of breast cancer is context-specific, and integrated analysis of gene expression and chromatin remodeling in normal and tumor tissues will be required to explain the mechanisms of risk alleles.

Abstract P5-05-03: The 5p12 breast cancer susceptibility locus is associated with MRPS30 expression in estrogen receptor - positive tumors

Genome-wide association studies (GWAS) have identified numerous loci linked to breast cancer susceptibility, but few of these loci have been convincingly linked to causal variants. We identified genes whose mRNA expression is linked to germline variation in normal mammary tissue and breast tumors by performing a genome-wide expression Quantitative Trait Locus (eQTL) analysis in 97 samples of normal mammary tissue and 286 breast adenocarcinomas. In the tumors we found 164 cis-acting loci, which affect expression of nearby genes. Twenty nine of these loci are previously unreported in other eQTL studies, including a cis-acting locus at 5p12 affecting expression of the mitochondrial ribosomal 28S subunit protein MRPS30 in estrogen receptor-positive but not estrogen receptor-negative tumors. This locus has been associated with susceptibility to estrogen receptor-positive breast cancer in several GWAS, with MRPS30 being the suggested candidate gene. One of these studies found that the rs771660 variant was most significantly associated with breast cancer in the 5p12 region, and we found that rs771660 was most strongly associated with MRPS30 gene expression. We provide the first evidence that this breast cancer susceptibility locus affects MRPS30 expression specifically in estrogen receptor-positive tumors and suggest a mechanism for estrogen activation of MRPS30 via proteins in the AP-1 complex.

Citation Information: Cancer Res 2012;72(24 Suppl):Abstract nr P5-05-03.

CHEK2*1100delC heterozygosity in women with breast cancer associated with early death, breast cancer-specific death, and increased risk of a second breast cancer

PURPOSE

We tested the hypotheses that CHEK2*1100delC heterozygosity is associated with increased risk of early death, breast cancer-specific death, and risk of a second breast cancer in women with a first breast cancer.

PATIENTS AND METHODS

From 22 studies participating in the Breast Cancer Association Consortium, 25,571 white women with invasive breast cancer were genotyped for CHEK2*1100delC and observed for up to 20 years (median, 6.6 years). We examined risk of early death and breast cancer-specific death by estrogen receptor status and risk of a second breast cancer after a first breast cancer in prospective studies.

RESULTS

CHEK2*1100delC heterozygosity was found in 459 patients (1.8%). In women with estrogen receptor-positive breast cancer, multifactorially adjusted hazard ratios for heterozygotes versus noncarriers were 1.43 (95% CI, 1.12 to 1.82; log-rank P = .004) for early death and 1.63 (95% CI, 1.24 to 2.15; log-rank P < .001) for breast cancer-specific death. In all women, hazard ratio for a second breast cancer was 2.77 (95% CI, 2.00 to 3.83; log-rank P < .001) increasing to 3.52 (95% CI, 2.35 to 5.27; log-rank P < .001) in women with estrogen receptor-positive first breast cancer only.

CONCLUSION

Among women with estrogen receptor-positive breast cancer, CHEK2*1100delC heterozygosity was associated with a 1.4-fold risk of early death, a 1.6-fold risk of breast cancer-specific death, and a 3.5-fold risk of a second breast cancer. This is one of the few examples of a genetic factor that influences long-term prognosis being documented in an extensive series of women with breast cancer.

Comparison of 6q25 breast cancer hits from Asian and European Genome Wide Association Studies in the Breast Cancer Association Consortium (BCAC)

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER-) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26-1.48), p = 7.6 × 10(-14) in Asians and 1.09 (95% CI 1.07-1.11), p = 6.8 × 10(-18) in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19-1.41), p = 1.2 × 10(-9) in Asians and 1.12 (95% CI 1.08-1.17), p = 3.8 × 10(-9) in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER- than ER+ tumours in Europeans [OR (ER-) = 1.20 (95% CI 1.15-1.25), p = 1.8 × 10(-17) versus OR (ER+) = 1.07 (95% CI 1.04-1.1), p = 1.3 × 10(-7), p(heterogeneity) = 5.1 × 10(-6)]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER- tumours.

19p13.1 is a triple-negative-specific breast cancer susceptibility locus

The 19p13.1 breast cancer susceptibility locus is a modifier of breast cancer risk in BRCA1 mutation carriers and is also associated with the risk of ovarian cancer. Here, we investigated 19p13.1 variation and risk of breast cancer subtypes, defined by estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor-2 (HER2) status, using 48,869 breast cancer cases and 49,787 controls from the Breast Cancer Association Consortium (BCAC). Variants from 19p13.1 were not associated with breast cancer overall or with ER-positive breast cancer but were significantly associated with ER-negative breast cancer risk [rs8170 OR, 1.10; 95% confidence interval (CI), 1.05-1.15; P = 3.49 × 10(-5)] and triple-negative (ER-, PR-, and HER2-negative) breast cancer (rs8170: OR, 1.22; 95% CI, 1.13-1.31; P = 2.22 × 10(-7)). However, rs8170 was no longer associated with ER-negative breast cancer risk when triple-negative cases were excluded (OR, 0.98; 95% CI, 0.89-1.07; P = 0.62). In addition, a combined analysis of triple-negative cases from BCAC and the Triple Negative Breast Cancer Consortium (TNBCC; N = 3,566) identified a genome-wide significant association between rs8170 and triple-negative breast cancer risk (OR, 1.25; 95% CI, 1.18-1.33; P = 3.31 × 10(-13)]. Thus, 19p13.1 is the first triple-negative-specific breast cancer risk locus and the first locus specific to a histologic subtype defined by ER, PR, and HER2 to be identified. These findings provide convincing evidence that genetic susceptibility to breast cancer varies by tumor subtype and that triple-negative tumors and other subtypes likely arise through distinct etiologic pathways.

Confirmation of 5p12 as a susceptibility locus for progesterone-receptor-positive, lower grade breast cancer

BACKGROUND

The single-nucleotide polymorphism (SNP) 5p12-rs10941679 has been found to be associated with risk of breast cancer, particularly estrogen receptor (ER)-positive disease. We aimed to further explore this association overall, and by tumor histopathology, in the Breast Cancer Association Consortium.

METHODS

Data were combined from 37 studies, including 40,972 invasive cases, 1,398 cases of ductal carcinoma in situ (DCIS), and 46,334 controls, all of white European ancestry, as well as 3,007 invasive cases and 2,337 controls of Asian ancestry. Associations overall and by tumor invasiveness and histopathology were assessed using logistic regression.

RESULTS

For white Europeans, the per-allele OR associated with 5p12-rs10941679 was 1.11 (95% CI = 1.08-1.14, P = 7 × 10(-18)) for invasive breast cancer and 1.10 (95% CI = 1.01-1.21, P = 0.03) for DCIS. For Asian women, the estimated OR for invasive disease was similar (OR = 1.07, 95%CI = 0.99-1.15, P = 0.09). Further analyses suggested that the association in white Europeans was largely limited to progesterone receptor (PR)-positive disease (per-allele OR = 1.16, 95% CI = 1.12-1.20, P = 1 × 10(-18) vs. OR = 1.03, 95% CI = 0.99-1.07, P = 0.2 for PR-negative disease; P(heterogeneity) = 2 × 10(-7)); heterogeneity by ER status was not observed (P = 0.2) once PR status was accounted for. The association was also stronger for lower grade tumors [per-allele OR (95% CI) = 1.20 (1.14-1.25), 1.13 (1.09-1.16), and 1.04 (0.99-1.08) for grade 1, 2, and 3/4, respectively; P(trend) = 5 × 10(-7)].

CONCLUSION

5p12 is a breast cancer susceptibility locus for PR-positive, lower grade breast cancer.

IMPACT

Multicenter fine-mapping studies of this region are needed as a first step to identifying the causal variant or variants.

Low penetrance breast cancer susceptibility loci are associated with specific breast tumor subtypes: findings from the Breast Cancer Association Consortium

Breast cancers demonstrate substantial biological, clinical and etiological heterogeneity. We investigated breast cancer risk associations of eight susceptibility loci identified in GWAS and two putative susceptibility loci in candidate genes in relation to specific breast tumor subtypes. Subtypes were defined by five markers (ER, PR, HER2, CK5/6, EGFR) and other pathological and clinical features. Analyses included up to 30 040 invasive breast cancer cases and 53 692 controls from 31 studies within the Breast Cancer Association Consortium. We confirmed previous reports of stronger associations with ER+ than ER- tumors for six of the eight loci identified in GWAS: rs2981582 (10q26) (P-heterogeneity = 6.1 × 10(-18)), rs3803662 (16q12) (P = 3.7 × 10(-5)), rs13281615 (8q24) (P = 0.002), rs13387042 (2q35) (P = 0.006), rs4973768 (3p24) (P = 0.003) and rs6504950 (17q23) (P = 0.002). The two candidate loci, CASP8 (rs1045485, rs17468277) and TGFB1 (rs1982073), were most strongly related with the risk of PR negative tumors (P = 5.1 × 10(-6) and P = 4.1 × 10(-4), respectively), as previously suggested. Four of the eight loci identified in GWAS were associated with triple negative tumors (P ≤ 0.016): rs3803662 (16q12), rs889312 (5q11), rs3817198 (11p15) and rs13387042 (2q35); however, only two of them (16q12 and 2q35) were associated with tumors with the core basal phenotype (P ≤ 0.002). These analyses are consistent with different biological origins of breast cancers, and indicate that tumor stratification might help in the identification and characterization of novel risk factors for breast cancer subtypes. This may eventually result in further improvements in prevention, early detection and treatment.