Functional annotation of the 2q35 breast cancer risk locus implicates a structural variant in influencing activity of a long-range enhancer element.

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

CYP3A7*1C allele: linking premenopausal oestrone and progesterone levels with risk of hormone receptor-positive breast cancers.

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.

Evidence that the 5p12 Variant rs10941679 Confers Susceptibility to Estrogen-Receptor-Positive Breast Cancer through FGF10 and MRPS30 Regulation.

Genome-wide association studies (GWASs) have revealed increased breast cancer risk associated with multiple genetic variants at 5p12. Here, we report the fine mapping of this locus using data from 104,660 subjects from 50 case-control studies in the Breast Cancer Association Consortium (BCAC). With data for 3,365 genotyped and imputed SNPs across a 1 Mb region (positions 44,394,495-45,364,167; NCBI build 37), we found evidence for at least three independent signals: the strongest signal, consisting of a single SNP rs10941679, was associated with risk of estrogen-receptor-positive (ER+) breast cancer (per-g allele OR ER+ = 1.15; 95% CI 1.13-1.18; p = 8.35 × 10-30). After adjustment for rs10941679, we detected signal 2, consisting of 38 SNPs more strongly associated with ER-negative (ER-) breast cancer (lead SNP rs6864776: per-a allele OR ER- = 1.10; 95% CI 1.05-1.14; p conditional = 1.44 × 10-12), and a single signal 3 SNP (rs200229088: per-t allele OR ER+ = 1.12; 95% CI 1.09-1.15; p conditional = 1.12 × 10-05). Expression quantitative trait locus analysis in normal breast tissues and breast tumors showed that the g (risk) allele of rs10941679 was associated with increased expression of FGF10 and MRPS30. Functional assays demonstrated that SNP rs10941679 maps to an enhancer element that physically interacts with the FGF10 and MRPS30 promoter regions in breast cancer cell lines. FGF10 is an oncogene that binds to FGFR2 and is overexpressed in ∼10% of human breast cancers, whereas MRPS30 plays a key role in apoptosis. These data suggest that the strongest signal of association at 5p12 is mediated through coordinated activation of FGF10 and MRPS30, two candidate genes for breast cancer pathogenesis.

Combined genetic and splicing analysis of BRCA1 c.[594-2A>C; 641A>G] highlights the relevance of naturally occurring in-frame transcripts for developing disease gene variant classification algorithms.

A recent analysis using family history weighting and co-observation classification modeling indicated that BRCA1 c.594-2A > C (IVS9-2A > C), previously described to cause exon 10 skipping (a truncating alteration), displays characteristics inconsistent with those of a high risk pathogenic BRCA1 variant. We used large-scale genetic and clinical resources from the ENIGMA, CIMBA and BCAC consortia to assess pathogenicity of c.594-2A > C. The combined odds for causality considering case-control, segregation and breast tumor pathology information was 3.23 × 10-8 Our data indicate that c.594-2A > C is always in cis with c.641A > G. The spliceogenic effect of c.[594-2A > C;641A > G] was characterized using RNA analysis of human samples and splicing minigenes. As expected, c.[594-2A > C; 641A > G] caused exon 10 skipping, albeit not due to c.594-2A > C impairing the acceptor site but rather by c.641A > G modifying exon 10 splicing regulatory element(s). Multiple blood-based RNA assays indicated that the variant allele did not produce detectable levels of full-length transcripts, with a per allele BRCA1 expression profile composed of ≈70-80% truncating transcripts, and ≈20-30% of in-frame Δ9,10 transcripts predicted to encode a BRCA1 protein with tumor suppression function.We confirm that BRCA1c.[594-2A > C;641A > G] should not be considered a high-risk pathogenic variant. Importantly, results from our detailed mRNA analysis suggest that BRCA-associated cancer risk is likely not markedly increased for individuals who carry a truncating variant in BRCA1 exons 9 or 10, or any other BRCA1 allele that permits 20-30% of tumor suppressor function. More generally, our findings highlight the importance of assessing naturally occurring alternative splicing for clinical evaluation of variants in disease-causing genes.

Fine-scale mapping of the 5q11.2 breast cancer locus reveals at least three independent risk variants regulating MAP3K1.

Genome-wide association studies (GWASs) have revealed SNP rs889312 on 5q11.2 to be associated with breast cancer risk in women of European ancestry. In an attempt to identify the biologically relevant variants, we analyzed 909 genetic variants across 5q11.2 in 103,991 breast cancer individuals and control individuals from 52 studies in the Breast Cancer Association Consortium. Multiple logistic regression analyses identified three independent risk signals: the strongest associations were with 15 correlated variants (iCHAV1), where the minor allele of the best candidate, rs62355902, associated with significantly increased risks of both estrogen-receptor-positive (ER(+): odds ratio [OR] = 1.24, 95% confidence interval [CI] = 1.21-1.27, ptrend = 5.7 × 10(-44)) and estrogen-receptor-negative (ER(-): OR = 1.10, 95% CI = 1.05-1.15, ptrend = 3.0 × 10(-4)) tumors. After adjustment for rs62355902, we found evidence of association of a further 173 variants (iCHAV2) containing three subsets with a range of effects (the strongest was rs113317823 [pcond = 1.61 × 10(-5)]) and five variants composing iCHAV3 (lead rs11949391; ER(+): OR = 0.90, 95% CI = 0.87-0.93, pcond = 1.4 × 10(-4)). Twenty-six percent of the prioritized candidate variants coincided with four putative regulatory elements that interact with the MAP3K1 promoter through chromatin looping and affect MAP3K1 promoter activity. Functional analysis indicated that the cancer risk alleles of four candidates (rs74345699 and rs62355900 [iCHAV1], rs16886397 [iCHAV2a], and rs17432750 [iCHAV3]) increased MAP3K1 transcriptional activity. Chromatin immunoprecipitation analysis revealed diminished GATA3 binding to the minor (cancer-protective) allele of rs17432750, indicating a mechanism for its action. We propose that the cancer risk alleles act to increase MAP3K1 expression in vivo and might promote breast cancer cell survival.

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.

Identification and characterization of novel associations in the CASP8/ALS2CR12 region on chromosome 2 with breast cancer risk.

Previous studies have suggested that polymorphisms in CASP8 on chromosome 2 are associated with breast cancer risk. To clarify the role of CASP8 in breast cancer susceptibility, we carried out dense genotyping of this region in the Breast Cancer Association Consortium (BCAC). Single-nucleotide polymorphisms (SNPs) spanning a 1 Mb region around CASP8 were genotyped in 46 450 breast cancer cases and 42 600 controls of European origin from 41 studies participating in the BCAC as part of a custom genotyping array experiment (iCOGS). Missing genotypes and SNPs were imputed and, after quality exclusions, 501 typed and 1232 imputed SNPs were included in logistic regression models adjusting for study and ancestry principal components. The SNPs retained in the final model were investigated further in data from nine genome-wide association studies (GWAS) comprising in total 10 052 case and 12 575 control subjects. The most significant association signal observed in European subjects was for the imputed intronic SNP rs1830298 in ALS2CR12 (telomeric to CASP8), with per allele odds ratio and 95% confidence interval [OR (95% confidence interval, CI)] for the minor allele of 1.05 (1.03-1.07), P = 1 × 10(-5). Three additional independent signals from intronic SNPs were identified, in CASP8 (rs36043647), ALS2CR11 (rs59278883) and CFLAR (rs7558475). The association with rs1830298 was replicated in the imputed results from the combined GWAS (P = 3 × 10(-6)), yielding a combined OR (95% CI) of 1.06 (1.04-1.08), P = 1 × 10(-9). Analyses of gene expression associations in peripheral blood and normal breast tissue indicate that CASP8 might be the target gene, suggesting a mechanism involving apoptosis.

CYP3A7*1C allele is associated with reduced levels of 2-hydroxylation pathway oestrogen metabolites.

BACKGROUND: Endogenous sex hormones are well-established risk factors for breast cancer; the contribution of specific oestrogen metabolites (EMs) and/or ratios of specific EMs is less clear. We have previously identified a CYP3A7*1C allele that is associated with lower urinary oestrone (E1) levels in premenopausal women. The purpose of this analysis was to determine whether this allele was associated with specific pathway EMs. METHODS: We measured successfully 12 EMs in mid-follicular phase urine samples from 30 CYP3A7*1C carriers and 30 non-carriers using HPLC-MS/MS. RESULTS: In addition to having lower urinary E1 levels, CYP3A7*1C carriers had significantly lower levels of four of the 2-hydroxylation pathway EMs that we measured (2-hydroxyestrone, P=1.1 × 10-12; 2-hydroxyestradiol, P=2.7 × 10-7; 2-methoxyestrone, P=1.9 × 10-12; and 2-methoxyestradiol, P=0.0009). By contrast, 16α-hydroxylation pathway EMs were slightly higher in carriers and significantly so for 17-epiestriol (P=0.002). CONCLUSIONS: The CYP3A7*1C allele is associated with a lower urinary E1 levels, a more pronounced reduction in 2-hydroxylation pathway EMs and a lower ratio of 2-hydroxylation:16α-hydroxylation EMs in premenopausal women. To further characterise the association between parent oestrogens, EMs and subsequent risk of breast cancer, characterisation of additional genetic variants that influence oestrogen metabolism and large prospective studies of a broad spectrum of EMs will be required.

Unbiased analysis of potential targets of breast cancer susceptibility loci by Capture Hi-C.

Genome-wide association studies have identified more than 70 common variants that are associated with breast cancer risk. Most of these variants map to non-protein-coding regions and several map to gene deserts, regions of several hundred kilobases lacking protein-coding genes. We hypothesized that gene deserts harbor long-range regulatory elements that can physically interact with target genes to influence their expression. To test this, we developed Capture Hi-C (CHi-C), which, by incorporating a sequence capture step into a Hi-C protocol, allows high-resolution analysis of targeted regions of the genome. We used CHi-C to investigate long-range interactions at three breast cancer gene deserts mapping to 2q35, 8q24.21, and 9q31.2. We identified interaction peaks between putative regulatory elements ("bait fragments") within the captured regions and "targets" that included both protein-coding genes and long noncoding (lnc) RNAs over distances of 6.6 kb to 2.6 Mb. Target protein-coding genes were IGFBP5, KLF4, NSMCE2, and MYC; and target lncRNAs included DIRC3, PVT1, and CCDC26. For one gene desert, we were able to define two SNPs (rs12613955 and rs4442975) that were highly correlated with the published risk variant and that mapped within the bait end of an interaction peak. In vivo ChIP-qPCR data show that one of these, rs4442975, affects the binding of FOXA1 and implicate this SNP as a putative functional variant.

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.

Genetic predisposition to in situ and invasive lobular carcinoma of the breast.

Invasive lobular breast cancer (ILC) accounts for 10-15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09-1.18), P = 6.0 × 10(-10); P-het for ILC vs IDC ER+ tumors = 1.8 × 10(-4)). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P<0.05. Two SNPs showed significantly stronger associations for ILC than LCIS (rs2981579/10q26/FGFR2, P-het = 0.04 and rs889312/5q11/MAP3K1, P-het = 0.03); and two showed stronger associations for LCIS than ILC (rs6678914/1q32/LGR6, P-het = 0.001 and rs1752911/6q14, P-het = 0.04). In addition, seven of the 75 known loci showed significant differences between ER+ tumors with IDC and ILC histology, three of these showing stronger associations for ILC (rs11249433/1p11, rs2981579/10q26/FGFR2 and rs10995190/10q21/ZNF365) and four associated only with IDC (5p12/rs10941679; rs2588809/14q24/RAD51L1, rs6472903/8q21 and rs1550623/2q31/CDCA7). In conclusion, we have identified one novel lobular breast cancer specific predisposition polymorphism at 7q34, and shown for the first time that common breast cancer polymorphisms predispose to LCIS. We have shown that many of the ER+ breast cancer predisposition loci also predispose to ILC, although there is some heterogeneity between ER+ lobular and ER+ IDC tumors. These data provide evidence for overlapping, but distinct etiological pathways within ER+ breast cancer between morphological subtypes.

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.

A common coding variant in CASP8 is associated with breast cancer risk.

The Breast Cancer Association Consortium (BCAC) has been established to conduct combined case-control analyses with augmented statistical power to try to confirm putative genetic associations with breast cancer. We genotyped nine SNPs for which there was some prior evidence of an association with breast cancer: CASP8 D302H (rs1045485), IGFBP3 -202 C --> A (rs2854744), SOD2 V16A (rs1799725), TGFB1 L10P (rs1982073), ATM S49C (rs1800054), ADH1B 3' UTR A --> G (rs1042026), CDKN1A S31R (rs1801270), ICAM5 V301I (rs1056538) and NUMA1 A794G (rs3750913). We included data from 9-15 studies, comprising 11,391-18,290 cases and 14,753-22,670 controls. We found evidence of an association with breast cancer for CASP8 D302H (with odds ratios (OR) of 0.89 (95% confidence interval (c.i.): 0.85-0.94) and 0.74 (95% c.i.: 0.62-0.87) for heterozygotes and rare homozygotes, respectively, compared with common homozygotes; P(trend) = 1.1 x 10(-7)) and weaker evidence for TGFB1 L10P (OR = 1.07 (95% c.i.: 1.02-1.13) and 1.16 (95% c.i.: 1.08-1.25), respectively; P(trend) = 2.8 x 10(-5)). These results demonstrate that common breast cancer susceptibility alleles with small effects on risk can be identified, given sufficiently powerful studies.

Inherited variants in the inner centromere protein (INCENP) gene of the chromosomal passenger complex contribute to the susceptibility of ER-negative breast cancer.

The chromosomal passenger complex (CPC) plays a pivotal role in the regulation of cell division. Therefore, inherited CPC variability could influence tumor development. The present candidate gene approach investigates the relationship between single nucleotide polymorphisms (SNPs) in genes encoding key CPC components and breast cancer risk. Fifteen SNPs in four CPC genes (INCENP, AURKB, BIRC5 and CDCA8) were genotyped in 88 911 European women from 39 case-control studies of the Breast Cancer Association Consortium. Possible associations were investigated in fixed-effects meta-analyses. The synonymous SNP rs1675126 in exon 7 of INCENP was associated with overall breast cancer risk [per A allele odds ratio (OR) 0.95, 95% confidence interval (CI) 0.92-0.98, P = 0.007] and particularly with estrogen receptor (ER)-negative breast tumors (per A allele OR 0.89, 95% CI 0.83-0.95, P = 0.0005). SNPs not directly genotyped were imputed based on 1000 Genomes. The SNPs rs1047739 in the 3' untranslated region and rs144045115 downstream of INCENP showed the strongest association signals for overall (per T allele OR 1.03, 95% CI 1.00-1.06, P = 0.0009) and ER-negative breast cancer risk (per A allele OR 1.06, 95% CI 1.02-1.10, P = 0.0002). Two genotyped SNPs in BIRC5 were associated with familial breast cancer risk (top SNP rs2071214: per G allele OR 1.12, 95% CI 1.04-1.21, P = 0.002). The data suggest that INCENP in the CPC pathway contributes to ER-negative breast cancer susceptibility in the European population. In spite of a modest contribution of CPC-inherited variants to the total burden of sporadic and familial breast cancer, their potential as novel targets for breast cancer treatment should be further investigated.

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.

Associations of common variants at 1p11.2 and 14q24.1 (RAD51L1) with breast cancer risk and heterogeneity by tumor subtype: findings from the Breast Cancer Association Consortium.

A genome-wide association study (GWAS) identified single-nucleotide polymorphisms (SNPs) at 1p11.2 and 14q24.1 (RAD51L1) as breast cancer susceptibility loci. The initial GWAS suggested stronger effects for both loci for estrogen receptor (ER)-positive tumors. Using data from the Breast Cancer Association Consortium (BCAC), we sought to determine whether risks differ by ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), grade, node status, tumor size, and ductal or lobular morphology. We genotyped rs11249433 at 1p.11.2, and two highly correlated SNPs rs999737 and rs10483813 (r(2)= 0.98) at 14q24.1 (RAD51L1), for up to 46 036 invasive breast cancer cases and 46 930 controls from 39 studies. Analyses by tumor characteristics focused on subjects reporting to be white women of European ancestry and were based on 25 458 cases, of which 87% had ER data. The SNP at 1p11.2 showed significantly stronger associations with ER-positive tumors [per-allele odds ratio (OR) for ER-positive tumors was 1.13, 95% CI = 1.10-1.16 and, for ER-negative tumors, OR was 1.03, 95% CI = 0.98-1.07, case-only P-heterogeneity = 7.6 × 10(-5)]. The association with ER-positive tumors was stronger for tumors of lower grade (case-only P= 6.7 × 10(-3)) and lobular histology (case-only P= 0.01). SNPs at 14q24.1 were associated with risk for most tumor subtypes evaluated, including triple-negative breast cancers, which has not been described previously. Our results underscore the need for large pooling efforts with tumor pathology data to help refine risk estimates for SNP associations with susceptibility to different subtypes of breast cancer.

A genome-wide association study to identify genetic susceptibility loci that modify ductal and lobular postmenopausal breast cancer risk associated with menopausal hormone therapy use: a two-stage design with replication.

Menopausal hormone therapy (MHT) is associated with an elevated risk of breast cancer in postmenopausal women. To identify genetic loci that modify breast cancer risk related to MHT use in postmenopausal women, we conducted a two-stage genome-wide association study (GWAS) with replication. In stage I, we performed a case-only GWAS in 731 invasive breast cancer cases from the German case-control study Mammary Carcinoma Risk Factor Investigation (MARIE). The 1,200 single nucleotide polymorphisms (SNPs) showing the lowest P values for interaction with current MHT use (within 6 months prior to breast cancer diagnosis), were carried forward to stage II, involving pooled case-control analyses including additional MARIE subjects (1,375 cases, 1,974 controls) as well as 795 cases and 764 controls of a Swedish case-control study. A joint P value was calculated for a combined analysis of stages I and II. Replication of the most significant interaction of the combined stage I and II was performed using 5,795 cases and 5,390 controls from nine studies of the Breast Cancer Association Consortium (BCAC). The combined stage I and II yielded five SNPs on chromosomes 2, 7, and 18 with joint P values <6 × 10(-6) for effect modification of current MHT use. The most significant interaction was observed for rs6707272 (P = 3 × 10(-7)) on chromosome 2 but was not replicated in the BCAC studies (P = 0.21). The potentially modifying SNPs are in strong linkage disequilibrium with SNPs in TRIP12 and DNER on chromosome 2 and SETBP1 on chromosome 18, previously linked to carcinogenesis. However, none of the interaction effects reached genome-wide significance. The inability to replicate the top SNP × MHT interaction may be due to limited power of the replication phase. Our study, however, suggests that there are unlikely to be SNPs that interact strongly enough with MHT use to be clinically significant in European women.

Genome-wide association study identifies novel breast cancer susceptibility loci.

Breast cancer exhibits familial aggregation, consistent with variation in genetic susceptibility to the disease. Known susceptibility genes account for less than 25% of the familial risk of breast cancer, and the residual genetic variance is likely to be due to variants conferring more moderate risks. To identify further susceptibility alleles, we conducted a two-stage genome-wide association study in 4,398 breast cancer cases and 4,316 controls, followed by a third stage in which 30 single nucleotide polymorphisms (SNPs) were tested for confirmation in 21,860 cases and 22,578 controls from 22 studies. We used 227,876 SNPs that were estimated to correlate with 77% of known common SNPs in Europeans at r2 > 0.5. SNPs in five novel independent loci exhibited strong and consistent evidence of association with breast cancer (P < 10(-7)). Four of these contain plausible causative genes (FGFR2, TNRC9, MAP3K1 and LSP1). At the second stage, 1,792 SNPs were significant at the P < 0.05 level compared with an estimated 1,343 that would be expected by chance, indicating that many additional common susceptibility alleles may be identifiable by this approach.

Gene-body hypermethylation of ATM in peripheral blood DNA of bilateral breast cancer patients.

Bilaterality of breast cancer is an indicator of constitutional cancer susceptibility; however, the molecular causes underlying this predisposition in the majority of cases is not known. We hypothesize that epigenetic misregulation of cancer-related genes could partially account for this predisposition. We have performed methylation microarray analysis of peripheral blood DNA from 14 women with bilateral breast cancer compared with 14 unaffected matched controls throughout 17 candidate breast cancer susceptibility genes including BRCA1, BRCA2, CHEK2, ATM, ESR1, SFN, CDKN2A, TP53, GSTP1, CDH1, CDH13, HIC1, PGR, SFRP1, MLH1, RARB and HSD17B4. We show that the majority of methylation variability is associated with intragenic repetitive elements. Detailed validation of the tiled region around ATM was performed by bisulphite modification and pyrosequencing of the same samples and in a second set of peripheral blood DNA from 190 bilateral breast cancer patients compared with 190 controls. We show significant hypermethylation of one intragenic repetitive element in breast cancer cases compared with controls (P = 0.0017), with the highest quartile of methylation associated with a 3-fold increased risk of breast cancer (OR 3.20, 95% CI 1.78-5.86, P = 0.000083). Increased methylation of this locus is associated with lower steady-state ATM mRNA level and correlates with age of cancer patients but not controls, suggesting a combined age-phenotype-related association. This research demonstrates the potential for gene-body epigenetic misregulation of ATM and other cancer-related genes in peripheral blood DNA that may be useful as a novel marker to estimate breast cancer risk. ACCESSION NUMBERS: The microarray data and associated .BED and .WIG files can be accessed through Gene Expression Omnibus accession number: GSE14603.

Association of ESR1 gene tagging SNPs with breast cancer risk.

We have conducted a three-stage, comprehensive single nucleotide polymorphism (SNP)-tagging association study of ESR1 gene variants (SNPs) in more than 55,000 breast cancer cases and controls from studies within the Breast Cancer Association Consortium (BCAC). No large risks or highly significant associations were revealed. SNP rs3020314, tagging a region of ESR1 intron 4, is associated with an increase in breast cancer susceptibility with a dominant mode of action in European populations. Carriers of the c-allele have an odds ratio (OR) of 1.05 [95% Confidence Intervals (CI) 1.02-1.09] relative to t-allele homozygotes, P = 0.004. There is significant heterogeneity between studies, P = 0.002. The increased risk appears largely confined to oestrogen receptor-positive tumour risk. The region tagged by SNP rs3020314 contains sequence that is more highly conserved across mammalian species than the rest of intron 4, and it may subtly alter the ratio of two mRNA splice forms.