Genome-wide association analysis of more than 120,000 individuals identifies 15 new susceptibility loci for breast cancer

Genome-wide association studies (GWAS) and large-scale replication studies have identified common variants in 79 loci associated with breast cancer, explaining ∼14% of the familial risk of the disease. To identify new susceptibility loci, we performed a meta-analysis of 11 GWAS, comprising 15,748 breast cancer cases and 18,084 controls together with 46,785 cases and 42,892 controls from 41 studies genotyped on a 211,155-marker custom array (iCOGS). Analyses were restricted to women of European ancestry. We generated genotypes for more than 11 million SNPs by imputation using the 1000 Genomes Project reference panel, and we identified 15 new loci associated with breast cancer at P < 5 × 10(-8). Combining association analysis with ChIP-seq chromatin binding data in mammary cell lines and ChIA-PET chromatin interaction data from ENCODE, we identified likely target genes in two regions: SETBP1 at 18q12.3 and RNF115 and PDZK1 at 1q21.1. One association appears to be driven by an amino acid substitution encoded in EXO1.

Improved DNA microarray detection sensitivity through immobilization of preformed in solution streptavidin/biotinylated oligonucleotide conjugates

A novel immobilization approach involving binding of preformed streptavidin/biotinylated oligonucleotide conjugates onto surfaces coated with biotinylated bovine serum albumin is presented. Microarrays prepared according to the proposed method were compared, in terms of detection sensitivity and specificity, with other immobilization schemes employing coupling of biotinylated oligonucleotides onto directly adsorbed surface streptavidin, or sequential coupling of streptavidin and biotinylated oligonucleotides onto a layer of adsorbed biotinylated bovine serum albumin. A comparison was performed employing biotinylated oligonucleotides corresponding to wild- and mutant-type sequences of seven single point mutations of the BRCA1 gene. With respect to the other immobilization protocols, the proposed oligonucleotide immobilization approach offered the highest hybridization signals (at least 5 times higher) and permitted more elaborative washings, thus providing considerably higher discrimination between complimentary and non-complementary DNA sequences for all mutations tested. In addition, the hybridization kinetics were significantly enhanced compared to two other immobilization protocols, permitting PCR sample analysis in less than 40 min. Thus, the proposed oligonucleotide immobilization approach offered improved detection sensitivity and discrimination ability along with considerably reduced analysis time, and it is expected to find wide application in DNA mutation detection.

CHEK2 c.1100delC allele is rarely identified in Greek breast cancer cases

The CHEK2 gene encodes a protein kinase that plays a crucial role in maintenance of genomic integrity and the DNA repair mechanism. CHEK2 germline mutations are associated with increased risk of breast cancer and other malignancies. From a clinical perspective, the most significant mutation identified is the c.1100delC mutation, which is associated with an approximately 25% lifetime breast cancer risk. The distribution of this mutation shows wide geographical variation; it is more prevalent in the Northern European countries and less common, or even absent, in Southern Europe. In order to estimate the frequency of the CHEK2 c.1100delC mutation in Greek breast cancer patients, we genotyped 2,449 patients (2,408 females and 41 males), which was the largest series ever tested for c.1100delC. The mean age of female and male breast cancer diagnosis was 49 and 59 years, respectively. All patients had previously tested negative for the Greek BRCA1 founder and recurrent mutations. The CHEK2 c.1100delC mutation was detected in 0.16% (4 of 2,408) of females, all of whom were diagnosed with breast cancer before the age of 50 years. Only one c.1100delC carrier was reported with breast cancer family history. The present study indicates that the CHEK2 c.1100delC mutation does not contribute substantially to hereditary breast cancer in patients of Greek descent.

Inherited mutations in 17 breast cancer susceptibility genes among a large triple-negative breast cancer cohort unselected for family history of breast cancer

PURPOSE

Recent advances in DNA sequencing have led to the development of breast cancer susceptibility gene panels for germline genetic testing of patients. We assessed the frequency of mutations in 17 predisposition genes, including BRCA1 and BRCA2, in a large cohort of patients with triple-negative breast cancer (TNBC) unselected for family history of breast or ovarian cancer to determine the utility of germline genetic testing for those with TNBC.

PATIENTS AND METHODS

Patients with TNBC (N = 1,824) unselected for family history of breast or ovarian cancer were recruited through 12 studies, and germline DNA was sequenced to identify mutations.

RESULTS

Deleterious mutations were identified in 14.6% of all patients. Of these, 11.2% had mutations in the BRCA1 (8.5%) and BRCA2 (2.7%) genes. Deleterious mutations in 15 other predisposition genes were detected in 3.7% of patients, with the majority observed in genes involved in homologous recombination, including PALB2 (1.2%) and BARD1, RAD51D, RAD51C, and BRIP1 (0.3% to 0.5%). Patients with TNBC with mutations were diagnosed at an earlier age (P < .001) and had higher-grade tumors (P = .01) than those without mutations.

CONCLUSION

Deleterious mutations in predisposition genes are present at high frequency in patients with TNBC unselected for family history of cancer. Mutation prevalence estimates suggest that patients with TNBC, regardless of age at diagnosis or family history of cancer, should be considered for germline genetic testing of BRCA1 and BRCA2. Although mutations in other predisposition genes are observed among patients with TNBC, better cancer risk estimates are needed before these mutations are used for clinical risk assessment in relatives.

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.

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.

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.

Common non-synonymous SNPs associated with breast cancer susceptibility: findings from the Breast Cancer Association Consortium

Candidate variant association studies have been largely unsuccessful in identifying common breast cancer susceptibility variants, although most studies have been underpowered to detect associations of a realistic magnitude. We assessed 41 common non-synonymous single-nucleotide polymorphisms (nsSNPs) for which evidence of association with breast cancer risk had been previously reported. Case-control data were combined from 38 studies of white European women (46 450 cases and 42 600 controls) and analyzed using unconditional logistic regression. Strong evidence of association was observed for three nsSNPs: ATXN7-K264R at 3p21 [rs1053338, per allele OR = 1.07, 95% confidence interval (CI) = 1.04-1.10, P = 2.9 × 10(-6)], AKAP9-M463I at 7q21 (rs6964587, OR = 1.05, 95% CI = 1.03-1.07, P = 1.7 × 10(-6)) and NEK10-L513S at 3p24 (rs10510592, OR = 1.10, 95% CI = 1.07-1.12, P = 5.1 × 10(-17)). The first two associations reached genome-wide statistical significance in a combined analysis of available data, including independent data from nine genome-wide association studies (GWASs): for ATXN7-K264R, OR = 1.07 (95% CI = 1.05-1.10, P = 1.0 × 10(-8)); for AKAP9-M463I, OR = 1.05 (95% CI = 1.04-1.07, P = 2.0 × 10(-10)). Further analysis of other common variants in these two regions suggested that intronic SNPs nearby are more strongly associated with disease risk. We have thus identified a novel susceptibility locus at 3p21, and confirmed previous suggestive evidence that rs6964587 at 7q21 is associated with risk. The third locus, rs10510592, is located in an established breast cancer susceptibility region; the association was substantially attenuated after adjustment for the known GWAS hit. Thus, each of the associated nsSNPs is likely to be a marker for another, non-coding, variant causally related to breast cancer risk. Further fine-mapping and functional studies are required to identify the underlying risk-modifying variants and the genes through which they act.

Genetic variation in mitotic regulatory pathway genes is associated with breast tumor grade

Mitotic index is an important component of histologic grade and has an etiologic role in breast tumorigenesis. Several small candidate gene studies have reported associations between variation in mitotic genes and breast cancer risk. We measured associations between 2156 single nucleotide polymorphisms (SNPs) from 194 mitotic genes and breast cancer risk, overall and by histologic grade, in the Breast Cancer Association Consortium (BCAC) iCOGS study (n = 39 067 cases; n = 42 106 controls). SNPs in TACC2 [rs17550038: odds ratio (OR) = 1.24, 95% confidence interval (CI) 1.16-1.33, P = 4.2 × 10(-10)) and EIF3H (rs799890: OR = 1.07, 95% CI 1.04-1.11, P = 8.7 × 10(-6)) were significantly associated with risk of low-grade breast cancer. The TACC2 signal was retained (rs17550038: OR = 1.15, 95% CI 1.07-1.23, P = 7.9 × 10(-5)) after adjustment for breast cancer risk SNPs in the nearby FGFR2 gene, suggesting that TACC2 is a novel, independent genome-wide significant genetic risk locus for low-grade breast cancer. While no SNPs were individually associated with high-grade disease, a pathway-level gene set analysis showed that variation across the 194 mitotic genes was associated with high-grade breast cancer risk (P = 2.1 × 10(-3)). These observations will provide insight into the contribution of mitotic defects to histological grade and the etiology of breast cancer.

MicroRNA related polymorphisms and breast cancer risk

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88-0.96), rs1052532 (OR 0.97; 95% CI: 0.95-0.99), rs10719 (OR 0.97; 95% CI: 0.94-0.99), rs4687554 (OR 0.97; 95% CI: 0.95-0.99, and rs3134615 (OR 1.03; 95% CI: 1.01-1.05) located in the 3' UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.

Abstract 3266: Expression quantitative trait locus analysis of triple negative breast cancer

Recent studies have shown that associations between common genetic variation and gene expression in breast tumors provide insight into the functional mechanisms underlying breast cancer risk loci. However, the relationship between genetic variation, gene expression, and risk of triple negative (TN) breast cancer remains largely unexplored. We performed a genome-wide expression quantitative trait locus (eQTL) analysis of TN breast cancer using 668 formalin-fixed paraffin embedded TN tumors from the Triple Negative Breast Cancer Consortium (TNBCC) to identify novel genes relevant to TN breast cancer risk and to further explore the biology underlying known TN risk loci. Cis-eQTLs were defined as correlations between single nucleotide polymorphisms (SNP) and expression of genes within 1Mb, and trans-eQTLs were defined as the remainder of all SNP-gene correlations. In a genome-wide analysis, we identified 68,012 cis-eQTLs, representing correlations between 42,193 SNPs and 2,092 genes that were significant at a 10% false discovery rate (FDR). We also identified 70 trans-eQTLs across two unique genes at the same FDR threshold. Five cis-eQTLs involved SNPs previously associated with TN breast cancer risk (p&lt;5.0x10-5) in 3,700 TN cases and 4,700 controls from the TNBCC. Three of these, located on chromosomes 5p (cis-eQTL p=4.8x10-34; risk p=1.1x10-7), 7p (cis-eQTL p=6.1x10-10; risk p=1.3x10-7), and 17q (cis-eQTL p=1.8x10-15; risk p=2.9x10-7), have previously been reported in public databases of normal tissue. The remaining two TN breast cancer cis-eQTLs, located on chromosomes 8q (cis-eQTL p=6.3x10-6, risk p=8.3x10-6) and 14q (cis-eQTL p=3.1x10-6, risk p=1.2x10-5), have not been reported as eQTLs in other tissues. In addition, 12 of 25 known TN breast cancer risk loci (PEX14, MDM4, 2q31.1, ESR1, 11q13.1, 11q24.3, 12p13.1, PTHLH, NTN4, 12q24, 19p13.1, MLK1), contained at least one significant cis-eQTL. In summary, we have identified five cis-eQTLs that may identify novel risk loci for TN breast cancer.

Citation Format: Kristen S. Purrington, Drakoulis Yannoukakos, Jane Carpenter, Heli Nevanlinna, Angela Cox, Gianluca Severi, Christine Ambrosone, Amanda Ewart Toland, Andrew K. Godwin, Hiltrud Brauch, Peter A. Fasching, Penelope Miron, Jenny Chang-Claude, Nicholas G. Martin, Grant W. Montgomery, Vessela Kristensen, Hoda Anton-Culver, Paul Goodfellow, Janet E. Olson, Hugues Sicotte, Naresh Prodduturi, Daniel W. Visscher, Jeanette E. Eckel-Passow, S. Keith Anderson, Seth Slettedahl, Curtis Olswold, Xianshu Wang, V. Shane Pankratz, Susan Slager, Wei Zheng, Arto Mannermaa, Ute Hamann, Diana M. Eccles, Celine M. Vachon, Fergus J. Couch. Expression quantitative trait locus analysis of triple negative breast cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3266. doi:10.1158/1538-7445.AM2014-3266

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.

Breast-cancer risk in families with mutations in PALB2

BACKGROUND

Germline loss-of-function mutations in PALB2 are known to confer a predisposition to breast cancer. However, the lifetime risk of breast cancer that is conferred by such mutations remains unknown.

METHODS

We analyzed the risk of breast cancer among 362 members of 154 families who had deleterious truncating, splice, or deletion mutations in PALB2. The age-specific breast-cancer risk for mutation carriers was estimated with the use of a modified segregation-analysis approach that allowed for the effects of PALB2 genotype and residual familial aggregation.

RESULTS

The risk of breast cancer for female PALB2 mutation carriers, as compared with the general population, was eight to nine times as high among those younger than 40 years of age, six to eight times as high among those 40 to 60 years of age, and five times as high among those older than 60 years of age. The estimated cumulative risk of breast cancer among female mutation carriers was 14% (95% confidence interval [CI], 9 to 20) by 50 years of age and 35% (95% CI, 26 to 46) by 70 years of age. Breast-cancer risk was also significantly influenced by birth cohort (P<0.001) and by other familial factors (P=0.04). The absolute breast-cancer risk for PALB2 female mutation carriers by 70 years of age ranged from 33% (95% CI, 25 to 44) for those with no family history of breast cancer to 58% (95% CI, 50 to 66) for those with two or more first-degree relatives with breast cancer at 50 years of age.

CONCLUSIONS

Loss-of-function mutations in PALB2 are an important cause of hereditary breast cancer, with respect both to the frequency of cancer-predisposing mutations and to the risk associated with them. Our data suggest the breast-cancer risk for PALB2 mutation carriers may overlap with that for BRCA2 mutation carriers. (Funded by the European Research Council and others.).

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.

Genome-wide association study identifies 25 known breast cancer susceptibility loci as risk factors for triple-negative breast cancer

Triple-negative (TN) breast cancer is an aggressive subtype of breast cancer associated with a unique set of epidemiologic and genetic risk factors. We conducted a two-stage genome-wide association study of TN breast cancer (stage 1: 1529 TN cases, 3399 controls; stage 2: 2148 cases, 1309 controls) to identify loci that influence TN breast cancer risk. Variants in the 19p13.1 and PTHLH loci showed genome-wide significant associations (P < 5 × 10(-) (8)) in stage 1 and 2 combined. Results also suggested a substantial enrichment of significantly associated variants among the single nucleotide polymorphisms (SNPs) analyzed in stage 2. Variants from 25 of 74 known breast cancer susceptibility loci were also associated with risk of TN breast cancer (P < 0.05). Associations with TN breast cancer were confirmed for 10 loci (LGR6, MDM4, CASP8, 2q35, 2p24.1, TERT-rs10069690, ESR1, TOX3, 19p13.1, RALY), and we identified associations with TN breast cancer for 15 additional breast cancer loci (P < 0.05: PEX14, 2q24.1, 2q31.1, ADAM29, EBF1, TCF7L2, 11q13.1, 11q24.3, 12p13.1, PTHLH, NTN4, 12q24, BRCA2, RAD51L1-rs2588809, MKL1). Further, two SNPs independent of previously reported signals in ESR1 [rs12525163 odds ratio (OR) = 1.15, P = 4.9 × 10(-) (4)] and 19p13.1 (rs1864112 OR = 0.84, P = 1.8 × 10(-) (9)) were associated with TN breast cancer. A polygenic risk score (PRS) for TN breast cancer based on known breast cancer risk variants showed a 4-fold difference in risk between the highest and lowest PRS quintiles (OR = 4.03, 95% confidence interval 3.46-4.70, P = 4.8 × 10(-) (69)). This translates to an absolute risk for TN breast cancer ranging from 0.8% to 3.4%, suggesting that genetic variation may be used for TN breast cancer risk prediction.

Genetic evaluation based on family history and Her2 status correctly identifies TP53 mutations in very early onset breast cancer cases

Currently, hereditary breast cancer is being attributed to more than 20 genes of differing penetrance. Although BRCA1 and BRCA2 are still the genes of reference for breast cancer susceptibility, extreme breast cancer phenotypes may be the result of deleterious alleles of other genes. Here, we report three families with early-onset breast cancer that were initially referred for BRCA1/BRCA2 genetic testing. They were diagnosed with breast cancer at an extraordinarily early age. On the basis of their extensive family history, which included multiple cancer types, and their Her2 status, they were suspected for Li-Fraumeni syndrome. Indeed, all three probands were found to harbor TP53 tumor suppressor gene mutations. These included p.C275X, described here for the first time, as well as p.R213X and p.Y220C, which have been described in the past. Our conclusion is that decisions on genetic analysis for inherited early onset breast cancer should always be based on detailed pedigree information, combined with Her2 status.

DNA glycosylases involved in base excision repair may be associated with cancer risk in BRCA1 and BRCA2 mutation carriers

Single Nucleotide Polymorphisms (SNPs) in genes involved in the DNA Base Excision Repair (BER) pathway could be associated with cancer risk in carriers of mutations in the high-penetrance susceptibility genes BRCA1 and BRCA2, given the relation of synthetic lethality that exists between one of the components of the BER pathway, PARP1 (poly ADP ribose polymerase), and both BRCA1 and BRCA2. In the present study, we have performed a comprehensive analysis of 18 genes involved in BER using a tagging SNP approach in a large series of BRCA1 and BRCA2 mutation carriers. 144 SNPs were analyzed in a two stage study involving 23,463 carriers from the CIMBA consortium (the Consortium of Investigators of Modifiers of BRCA1 and BRCA2). Eleven SNPs showed evidence of association with breast and/or ovarian cancer at p<0.05 in the combined analysis. Four of the five genes for which strongest evidence of association was observed were DNA glycosylases. The strongest evidence was for rs1466785 in the NEIL2 (endonuclease VIII-like 2) gene (HR: 1.09, 95% CI (1.03-1.16), p = 2.7 × 10(-3)) for association with breast cancer risk in BRCA2 mutation carriers, and rs2304277 in the OGG1 (8-guanine DNA glycosylase) gene, with ovarian cancer risk in BRCA1 mutation carriers (HR: 1.12 95%CI: 1.03-1.21, p = 4.8 × 10(-3)). DNA glycosylases involved in the first steps of the BER pathway may be associated with cancer risk in BRCA1/2 mutation carriers and should be more comprehensively studied.

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.

FGF receptor genes and breast cancer susceptibility: results from the Breast Cancer Association Consortium

BACKGROUND

Breast cancer is one of the most common malignancies in women. Genome-wide association studies have identified FGFR2 as a breast cancer susceptibility gene. Common variation in other fibroblast growth factor (FGF) receptors might also modify risk. We tested this hypothesis by studying genotyped single-nucleotide polymorphisms (SNPs) and imputed SNPs in FGFR1, FGFR3, FGFR4 and FGFRL1 in the Breast Cancer Association Consortium.

METHODS

Data were combined from 49 studies, including 53 835 cases and 50 156 controls, of which 89 050 (46 450 cases and 42 600 controls) were of European ancestry, 12 893 (6269 cases and 6624 controls) of Asian and 2048 (1116 cases and 932 controls) of African ancestry. Associations with risk of breast cancer, overall and by disease sub-type, were assessed using unconditional logistic regression.

RESULTS

Little evidence of association with breast cancer risk was observed for SNPs in the FGF receptor genes. The strongest evidence in European women was for rs743682 in FGFR3; the estimated per-allele odds ratio was 1.05 (95% confidence interval=1.02-1.09, P=0.0020), which is substantially lower than that observed for SNPs in FGFR2.

CONCLUSION

Our results suggest that common variants in the other FGF receptors are not associated with risk of breast cancer to the degree observed for FGFR2.

Functional variants at the 11q13 risk locus for breast cancer regulate cyclin D1 expression through long-range enhancers

Analysis of 4,405 variants in 89,050 European subjects from 41 case-control studies identified three independent association signals for estrogen-receptor-positive tumors at 11q13. The strongest signal maps to a transcriptional enhancer element in which the G allele of the best candidate causative variant rs554219 increases risk of breast cancer, reduces both binding of ELK4 transcription factor and luciferase activity in reporter assays, and may be associated with low cyclin D1 protein levels in tumors. Another candidate variant, rs78540526, lies in the same enhancer element. Risk association signal 2, rs75915166, creates a GATA3 binding site within a silencer element. Chromatin conformation studies demonstrate that these enhancer and silencer elements interact with each other and with their likely target gene, CCND1.

Multiple independent variants at the TERT locus are associated with telomere length and risks of breast and ovarian cancer

TERT-locus SNPs and leukocyte telomere measures are reportedly associated with risks of multiple cancers. Using the Illumina custom genotyping array iCOGs, we analyzed ∼480 SNPs at the TERT locus in breast (n = 103,991), ovarian (n = 39,774) and BRCA1 mutation carrier (n = 11,705) cancer cases and controls. Leukocyte telomere measurements were also available for 53,724 participants. Most associations cluster into three independent peaks. The minor allele at the peak 1 SNP rs2736108 associates with longer telomeres (P = 5.8 × 10(-7)), lower risks for estrogen receptor (ER)-negative (P = 1.0 × 10(-8)) and BRCA1 mutation carrier (P = 1.1 × 10(-5)) breast cancers and altered promoter assay signal. The minor allele at the peak 2 SNP rs7705526 associates with longer telomeres (P = 2.3 × 10(-14)), higher risk of low-malignant-potential ovarian cancer (P = 1.3 × 10(-15)) and greater promoter activity. The minor alleles at the peak 3 SNPs rs10069690 and rs2242652 increase ER-negative (P = 1.2 × 10(-12)) and BRCA1 mutation carrier (P = 1.6 × 10(-14)) breast and invasive ovarian (P = 1.3 × 10(-11)) cancer risks but not via altered telomere length. The cancer risk alleles of rs2242652 and rs10069690, respectively, increase silencing and generate a truncated TERT splice variant.

Genome-wide association study in BRCA1 mutation carriers identifies novel loci associated with breast and ovarian cancer risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7 × 10(-8), HR = 1.14, 95% CI: 1.09-1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4 × 10(-8), HR = 1.27, 95% CI: 1.17-1.38) and 4q32.3 (rs4691139, P = 3.4 × 10(-8), HR = 1.20, 95% CI: 1.17-1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10(-4)). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%-50% compared to 81%-100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.