Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

The fate of BRCA1-related germline mutations in triple-negative breast tumors

The preservation of pathogenic BRCA1/2 germline mutations in tumor tissues is usually not questioned, while it remains unknown whether these interact with somatic genotypes for patient outcome. Herein we compared germline and tumor genotypes in operable triple-negative breast cancer (TNBC) and evaluated their combined effects on prognosis. We analyzed baseline germline and primary tumor genotype data obtained by Sanger and Next Generation Sequencing in 194 TNBC patients. We also performed multiple tests interrogating the preservation of germline mutations in matched tumors and breast tissue from carriers with available material. Patients had been treated within clinical trials with adjuvant anthracyclines-taxanes based chemotherapy. We identified 50 (26%) germline mutation carriers (78% in BRCA1) and 136 (71%) tumors with somatic mutations (83% in TP53). Tumor mutation patterns differed between carriers and non-carriers (P<0.001); PIK3CA mutations were exclusively present in non-carriers (P=0.007). Germline BRCA1/2 mutations were not detected in matched tumors and breast tissues from 14 out of 33 (42%) evaluable carriers. Microsatellite markers revealed tumor loss of the germline mutant allele in one case only. Tumors that had lost the germline mutation demonstrated a higher incidence of somatic TP53 mutations as compared to tumors with preserved germline mutations (P=0.036). Germline mutation status significantly interacted with tumor TP53 mutations for patient disease-free survival (interaction P=0.026): In non-carriers, tumor TP53 mutations did not affect outcome; In carriers, those with mutated TP53 tumors experienced more relapses compared to those with wild-type TP53 tumors (36% vs. 9% relapse rate, respectively). In conclusion, we show that loss of germline BRCA1/2 mutations is not a rare event in TNBC. This finding, the observed differences in tumor genotypes with respect to germline status and the prognostic interaction between germline BRCA1-related and tumor TP53 mutation status prompt for combined germline and tumor genotyping for the classification of TNBC, particularly in the context of clinical trials evaluating synthetic lethality drugs.

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.

rs2735383, located at a microRNA binding site in the 3'UTR of NBS1, is not associated with breast cancer risk

NBS1, also known as NBN, plays an important role in maintaining genomic stability. Interestingly, rs2735383 G > C, located in a microRNA binding site in the 3'-untranslated region (UTR) of NBS1, was shown to be associated with increased susceptibility to lung and colorectal cancer. However, the relation between rs2735383 and susceptibility to breast cancer is not yet clear. Therefore, we genotyped rs2735383 in 1,170 familial non-BRCA1/2 breast cancer cases and 1,077 controls using PCR-based restriction fragment length polymorphism (RFLP-PCR) analysis, but found no association between rs2735383CC and breast cancer risk (OR = 1.214, 95% CI = 0.936-1.574, P = 0.144). Because we could not exclude a small effect size due to a limited sample size, we further analyzed imputed rs2735383 genotypes (r2 > 0.999) of 47,640 breast cancer cases and 46,656 controls from the Breast Cancer Association Consortium (BCAC). However, rs2735383CC was not associated with overall breast cancer risk in European (OR = 1.014, 95% CI = 0.969-1.060, P = 0.556) nor in Asian women (OR = 0.998, 95% CI = 0.905-1.100, P = 0.961). Subgroup analyses by age, age at menarche, age at menopause, menopausal status, number of pregnancies, breast feeding, family history and receptor status also did not reveal a significant association. This study therefore does not support the involvement of the genotype at NBS1 rs2735383 in breast cancer susceptibility.

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Genome-Wide Meta-Analyses of Breast, Ovarian, and Prostate Cancer Association Studies Identify Multiple New Susceptibility Loci Shared by at Least Two Cancer Types

Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10(-8) seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type-specific expression quantitative trait locus and enhancer-gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10(-5) in the three-cancer meta-analysis.

SIGNIFICANCE

We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052-67. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 932.

An intergenic risk locus containing an enhancer deletion in 2q35 modulates breast cancer risk by deregulating IGFBP5 expression

Breast cancer is the most diagnosed malignancy and the second leading cause of cancer mortality in females. Previous association studies have identified variants on 2q35 associated with the risk of breast cancer. To identify functional susceptibility loci for breast cancer, we interrogated the 2q35 gene desert for chromatin architecture and functional variation correlated with gene expression. We report a novel intergenic breast cancer risk locus containing an enhancer copy number variation (enCNV; deletion) located approximately 400Kb upstream to IGFBP5, which overlaps an intergenic ERα-bound enhancer that loops to the IGFBP5 promoter. The enCNV is correlated with modified ERα binding and monoallelic-repression of IGFBP5 following oestrogen treatment. We investigated the association of enCNV genotype with breast cancer in 1,182 cases and 1,362 controls, and replicate our findings in an independent set of 62,533 cases and 60,966 controls from 41 case control studies and 11 GWAS. We report a dose-dependent inverse association of 2q35 enCNV genotype (percopy OR = 0.68 95%CI 0.55-0.83, P = 0.0002; replication OR = 0.77 95% CI 0.73-0.82, P = 2.1 × 10-19) and identify 13 additional linked variants (r2 > 0.8) in the 20Kb linkage block containing the enCNV (P = 3.2 × 10-15 - 5.6 × 10-17). These associations were independent of previously reported 2q35 variants, rs13387042/rs4442975 and rs16857609, and were stronger for ER-positive than ER-negative disease. Together, these results suggest that 2q35 breast cancer risk loci may be mediating their effect through IGFBP5.

Fine-Mapping of the 1p11.2 Breast Cancer Susceptibility Locus

The Cancer Genetic Markers of Susceptibility genome-wide association study (GWAS) originally identified a single nucleotide polymorphism (SNP) rs11249433 at 1p11.2 associated with breast cancer risk. To fine-map this locus, we genotyped 92 SNPs in a 900kb region (120,505,799-121,481,132) flanking rs11249433 in 45,276 breast cancer cases and 48,998 controls of European, Asian and African ancestry from 50 studies in the Breast Cancer Association Consortium. Genotyping was done using iCOGS, a custom-built array. Due to the complicated nature of the region on chr1p11.2: 120,300,000-120,505,798, that lies near the centromere and contains seven duplicated genomic segments, we restricted analyses to 429 SNPs excluding the duplicated regions (42 genotyped and 387 imputed). Per-allelic associations with breast cancer risk were estimated using logistic regression models adjusting for study and ancestry-specific principal components. The strongest association observed was with the original identified index SNP rs11249433 (minor allele frequency (MAF) 0.402; per-allele odds ratio (OR) = 1.10, 95% confidence interval (CI) 1.08-1.13, P = 1.49 x 10-21). The association for rs11249433 was limited to ER-positive breast cancers (test for heterogeneity P≤8.41 x 10-5). Additional analyses by other tumor characteristics showed stronger associations with moderately/well differentiated tumors and tumors of lobular histology. Although no significant eQTL associations were observed, in silico analyses showed that rs11249433 was located in a region that is likely a weak enhancer/promoter. Fine-mapping analysis of the 1p11.2 breast cancer susceptibility locus confirms this region to be limited to risk to cancers that are ER-positive.

Haplotype analysis reveals that the recurrent BRCA1 deletion of exons 23 and 24 is a Greek founder mutation

A recurrent large genomic rearrangement (LGR) encompassing exons 23 and 24 of the BRCA1 gene has been identified in breast-ovarian cancer families of Greek origin. Its breakpoints have been determined as c.5406 + 664_*8273del11052 (RefSeq: NM_007294.3) and a diagnostic polymerase chain reaction (PCR) has been set up for rapid screening. In a series of 2,092 high-risk families completely screened for BRCA1 and BRCA2 germline mutations, we have found the deletion in 35 families (1.68%), representing 7.83% of the mutations identified in both genes and 10.3% of the total BRCA1 mutations. In order to characterize this deletion as a founder mutation, haplotype analysis was conducted in 60 carriers from 35 families, using three BRCA1 intragenic microsatellite markers and four markers surrounding the BRCA1 locus. Our results demonstrate a common shared core disease-associated haplotype of 2.89Mb. Our calculations estimate that the deletion has originated from a common ancestor 1450 years ago, which most probably inhabited the Asia Minor area. The particular (LGR) is the third mutation of such type that is proven to have a Greek founder effect in the Greek population, illustrating the necessity for LGRs testing in individuals of Greek descent.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND

Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk.

METHOD

We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation.

RESULTS

Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05.

CONCLUSION

This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk.

RAD51B in Familial Breast Cancer

Common variation on 14q24.1, close to RAD51B, has been associated with breast cancer: rs999737 and rs2588809 with the risk of female breast cancer and rs1314913 with the risk of male breast cancer. The aim of this study was to investigate the role of RAD51B variants in breast cancer predisposition, particularly in the context of familial breast cancer in Finland. We sequenced the coding region of RAD51B in 168 Finnish breast cancer patients from the Helsinki region for identification of possible recurrent founder mutations. In addition, we studied the known rs999737, rs2588809, and rs1314913 SNPs and RAD51B haplotypes in 44,791 breast cancer cases and 43,583 controls from 40 studies participating in the Breast Cancer Association Consortium (BCAC) that were genotyped on a custom chip (iCOGS). We identified one putatively pathogenic missense mutation c.541C>T among the Finnish cancer patients and subsequently genotyped the mutation in additional breast cancer cases (n = 5259) and population controls (n = 3586) from Finland and Belarus. No significant association with breast cancer risk was seen in the meta-analysis of the Finnish datasets or in the large BCAC dataset. The association with previously identified risk variants rs999737, rs2588809, and rs1314913 was replicated among all breast cancer cases and also among familial cases in the BCAC dataset. The most significant association was observed for the haplotype carrying the risk-alleles of all the three SNPs both among all cases (odds ratio (OR): 1.15, 95% confidence interval (CI): 1.11-1.19, P = 8.88 x 10-16) and among familial cases (OR: 1.24, 95% CI: 1.16-1.32, P = 6.19 x 10-11), compared to the haplotype with the respective protective alleles. Our results suggest that loss-of-function mutations in RAD51B are rare, but common variation at the RAD51B region is significantly associated with familial breast cancer risk.

Identification of four novel susceptibility loci for oestrogen receptor negative breast cancer

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

Genes associated with histopathologic features of triple negative breast tumors predict molecular subtypes

Distinct subtypes of triple negative (TN) breast cancer have been identified by tumor expression profiling. However, little is known about the relationship between histopathologic features of TN tumors, which reflect aspects of both tumor behavior and tumor microenvironment, and molecular TN subtypes. The histopathologic features of TN tumors were assessed by central review and 593 TN tumors were subjected to whole genome expression profiling using the Illumina Whole Genome DASL array. TN molecular subtypes were defined based on gene expression data associated with histopathologic features of TN tumors. Gene expression analysis yielded signatures for four TN subtypes (basal-like, androgen receptor positive, immune, and stromal) consistent with previous studies. Expression analysis also identified genes significantly associated with the 12 histological features of TN tumors. Development of signatures using these markers of histopathological features resulted in six distinct TN subtype signatures, including an additional basal-like and stromal signature. The additional basal-like subtype was distinguished by elevated expression of cell motility and glucose metabolism genes and reduced expression of immune signaling genes, whereas the additional stromal subtype was distinguished by elevated expression of immunomodulatory pathway genes. Histopathologic features that reflect heterogeneity in tumor architecture, cell structure, and tumor microenvironment are related to TN subtype. Accounting for histopathologic features in the development of gene expression signatures, six major subtypes of TN breast cancer were identified.

Height and Breast Cancer Risk: Evidence From Prospective Studies and Mendelian Randomization

BACKGROUND

Epidemiological studies have linked adult height with breast cancer risk in women. However, the magnitude of the association, particularly by subtypes of breast cancer, has not been established. Furthermore, the mechanisms of the association remain unclear.

METHODS

We performed a meta-analysis to investigate associations between height and breast cancer risk using data from 159 prospective cohorts totaling 5216302 women, including 113178 events. In a consortium with individual-level data from 46325 case patients and 42482 control patients, we conducted a Mendelian randomization analysis using a genetic score that comprised 168 height-associated variants as an instrument. This association was further evaluated in a second consortium using summary statistics data from 16003 case patients and 41335 control patients.

RESULTS

The pooled relative risk of breast cancer was 1.17 (95% confidence interval [CI] = 1.15 to 1.19) per 10cm increase in height in the meta-analysis of prospective studies. In Mendelian randomization analysis, the odds ratio of breast cancer per 10cm increase in genetically predicted height was 1.22 (95% CI = 1.13 to 1.32) in the first consortium and 1.21 (95% CI = 1.05 to 1.39) in the second consortium. The association was found in both premenopausal and postmenopausal women but restricted to hormone receptor-positive breast cancer. Analyses of height-associated variants identified eight new loci associated with breast cancer risk after adjusting for multiple comparisons, including three loci at 1q21.2, DNAJC27, and CCDC91 at genome-wide significance level P < 5×10(-8).

CONCLUSIONS

Our study provides strong evidence that adult height is a risk factor for breast cancer in women and certain genetic factors and biological pathways affecting adult height have an important role in the etiology of breast cancer.

Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Abstract 1924: Genes associated with histopathologic features of triple negative breast tumors predict molecular subtypes

Distinct subtypes of triple negative (TN) breast cancer have currently been identified by tumor expression profiling and are thought to be influenced by the tumor microenvironment. However, little is known about the relationship between distinguishing histopathologic features of TN tumors, which reflect aspects of both tumor behavior and tumor microenvironment, and molecular TN subtypes. We analyzed whole genome expression profiles of 593 TN tumors from the Triple Negative Breast Cancer Consortium (TNBCC) using the Illumina Whole Genome DASL array. Histopathologic review of mitotic index, grade, trabecular pattern, growth pattern, involution, presence of ductal carcinoma in situ, and degree of lymphoid infiltrate, necrosis, fibrosis, and lobulitis was completed for 310 of these tumors. Clusters were identified both agnostically and using genes associated with these histopathologic features, and subtypes were evaluated for enriched pathways and associations with histopathologic features. Agnostic clustering identified four stable clusters which represented a basal-like (BL) signature, a luminal signature with overexpression of androgen receptor and targets, an immune signature, and a stromal (STR) signature. Overall, the TN tumors tended to have little or no necrosis (91%), high grade (88%), a high mitotic index (80%), fibrosis (77%), an epicentric growth pattern (62%), incomplete lobular involution (55%), trabecular structures (37%), lymphocytic infiltration (25%), adjacent ductal carcinoma in situ (24%), and variable amounts of lobulitis. Using a combined total of 2,776 unique probes individually associated with the ten histopathologic features measured, we identified six distinct TN subtypes, with an additional BL and STR signature compared to the four agnostic subtypes. Our findings suggest that histopathologic features that reflect heterogeneity in tumor cell structure and behavior as well as the tumor microenvironment, particularly the presence of lymphocytic infiltrate, are critical for TN subtype identification. Further work is needed to validate these signatures and to better understand the immunologic pathways associated with each subtype.

Citation Format: Kristen Purrington, Daniel Visscher, Drakoulis Yannoukakos, Jane Carpenter, Heli Nevanlinna, Arto Mannerma, Xianshu Wang, Graham Giles, Wei Zheng, Angela Cox, Hiltrud Brauch, Ute Hamann, Diana Eccles, Celine M. Vachon, Fergus J. Couch. Genes associated with histopathologic features of triple negative breast tumors predict molecular subtypes. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1924. doi:10.1158/1538-7445.AM2015-1924

Polymorphisms in a Putative Enhancer at the 10q21.2 Breast Cancer Risk Locus Regulate NRBF2 Expression

Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82-0.88]) and ER-negative (OR = 0.87 [0.82-0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91-0.95] and OR = 1.06 [1.03-1.09]) and ER-negative (OR = 0.95 [0.91-0.98] and OR = 1.08 [1.04-1.13]) disease. There was weaker evidence for iCHAV4, located 5' of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90-0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1-4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.

Epidemiological and clinicopathological characteristics of BRCA-positive and BRCA-negative breast cancer patients in Greece

PURPOSE

BRCA mutation carriers can benefit from targeted clinical interventions. On the other hand, families with evident aggregation of breast cancer (BC) cases and a BRCA-negative genetic test can still be considered as of elevated risk, since the underlying genetic factor remains unidentified. In the present study, we compared clinical and demographic characteristics between BRCA1 mutation carriers (BRCA1mut) and non-carriers (non-BRCA1) in a Greek group of BC patients (n=321).

METHODS

Data were collected and analyzed from 321 women with BC, with 131 patients screened for pathogenic mutations in the high-penetrant genes BRCA1 and BRCA2. Collected data included demographics, pedigrees, tumor histopathology and immunohistochemistry findings.

RESULTS

In BRCA1mut patients, their mothers and grand- mothers were diagnosed at a younger age compared to non-BRCA1-carriers. Additionally, BRCA1mut patients were diagnosed with mainly estrogen receptor (ER) negative (p<0.001), Her-2 negative (p<0.05) and triple negative (p<0.01) tumors. The youngest generation was diagnosed with familial breast cancer (FBC) 9.7 years earlier than their mothers (p<0.001). Age at BC diagnosis negatively correlated with the nuclear grade of breast tumors (r=-0.3, p<0.05). Among parous individuals, the number of full-term pregnancies significantly correlated with the age at BC onset (r=0.19, p<0.05).

CONCLUSION

Despite their similarities, FBC cases with identified BRCA1 mutations exhibit a clearly distinct profile. We have identified an anticipation effect in FBC patients, with significantly reduced age at diagnosis in younger generations. Increased parity seems to prevent early BC onset. This is the first study comparing clinical and demographic characteristics of FBC BRCA1mut and non-carriers in a Greek cohort.

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.

Abstract P4-12-03: Triple-negative breast cancer: Frequency of inherited mutations in breast cancer susceptibility genes

Background: Guidelines recommend germline mutation testing of breast cancer predisposition genes in triple negative (TN) breast cancer cases with a family history of breast or ovarian cancer or when diagnosed under age 60. However, the prevalence of mutations in these genes among TN cases unselected for family history of breast or ovarian cancer is not known.

Methods: To assess the frequency of mutations in 16 predisposition genes in TN cases we screened a large cohort of TN patients (n=1824) unselected for family history of breast or ovarian cancer from 12 centers and 824 study matched unaffected controls for mutations using a panel-based sequencing approach.

Results: Deleterious mutations were identified in 15% of TN patients: 8.5% had BRCA1, 2.7% had BRCA2, and 3.6% had mutations in 12 other genes. Mutations in non-BRCA1/2 genes encoding proteins implicated in homologous recombination repair of DNA double strand breaks were detected at the same frequency as in breast cancer families. TN cases with mutations had high-grade tumors and were diagnosed at an earlier age than non-mutated cases. However, 10% of TN cases diagnosed at ≥60 years and 5% with no family history of cancer were also found to carry mutations. Inactivating mutations in non-BRCA1/2 predisposition genes were associated with moderate to high risks of TN breast cancer.

Conclusions: National Comprehensive Cancer Network (NCCN) guidelines support clinical genetic testing of breast cancer predisposition genes in 95% of TN breast cancer patients carrying mutations in susceptibility genes. In contrast, National Institute of Health and Care Excellence (NICE) guidelines in the U.K. do not support genetic testing of a substantial proportion of TN patients with predisposing alleles. Frequency tables for inherited mutations in known predisposition genes based on age of diagnosis and family history of cancer will allow for selection of TN patients most likely to carry mutations in the predisposition genes.

Citation Format: Fergus J Couch, Steven N Hart, Priyanka Sharma, Amanda Ewart Toland, Penelope Miron, Janet E Olson, Andrew Godwin, Vernon S Pankratz, Curtis Olswold, Seth Slettedahl, Lucia Guidugli, Matthias W Beckmann, Brigitte Rack, Arif B Ekici, Irene Konstantopoulou, Florentia Fostira, George Fountzilas, Liisa M Pelttari, Song Yao, Judy Garber, Angela Cox, Hiltrud Brauch, Christine Ambrosone, Heli Nevanlinna, Drakoulis Yannoukakos, Susan L Slager, Celine M Vachon, Diana M Eccles, Peter A Fasching. Triple-negative breast cancer: Frequency of inherited mutations in breast cancer susceptibility genes [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P4-12-03.

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.

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

INTRODUCTION

Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers.

METHODS

We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals.

RESULTS

We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk.

CONCLUSIONS

This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

Association of type and location of BRCA1 and BRCA2 mutations with risk of breast and ovarian cancer

IMPORTANCE

Limited information about the relationship between specific mutations in BRCA1 or BRCA2 (BRCA1/2) and cancer risk exists.

OBJECTIVE

To identify mutation-specific cancer risks for carriers of BRCA1/2.

DESIGN, SETTING, AND PARTICIPANTS

Observational study of women who were ascertained between 1937 and 2011 (median, 1999) and found to carry disease-associated BRCA1 or BRCA2 mutations. The international sample comprised 19,581 carriers of BRCA1 mutations and 11,900 carriers of BRCA2 mutations from 55 centers in 33 countries on 6 continents. We estimated hazard ratios for breast and ovarian cancer based on mutation type, function, and nucleotide position. We also estimated RHR, the ratio of breast vs ovarian cancer hazard ratios. A value of RHR greater than 1 indicated elevated breast cancer risk; a value of RHR less than 1 indicated elevated ovarian cancer risk.

EXPOSURES

Mutations of BRCA1 or BRCA2.

MAIN OUTCOMES AND MEASURES

Breast and ovarian cancer risks.

RESULTS

Among BRCA1 mutation carriers, 9052 women (46%) were diagnosed with breast cancer, 2317 (12%) with ovarian cancer, 1041 (5%) with breast and ovarian cancer, and 7171 (37%) without cancer. Among BRCA2 mutation carriers, 6180 women (52%) were diagnosed with breast cancer, 682 (6%) with ovarian cancer, 272 (2%) with breast and ovarian cancer, and 4766 (40%) without cancer. In BRCA1, we identified 3 breast cancer cluster regions (BCCRs) located at c.179 to c.505 (BCCR1; RHR = 1.46; 95% CI, 1.22-1.74; P = 2 × 10(-6)), c.4328 to c.4945 (BCCR2; RHR = 1.34; 95% CI, 1.01-1.78; P = .04), and c. 5261 to c.5563 (BCCR2', RHR = 1.38; 95% CI, 1.22-1.55; P = 6 × 10(-9)). We also identified an ovarian cancer cluster region (OCCR) from c.1380 to c.4062 (approximately exon 11) with RHR = 0.62 (95% CI, 0.56-0.70; P = 9 × 10(-17)). In BRCA2, we observed multiple BCCRs spanning c.1 to c.596 (BCCR1; RHR = 1.71; 95% CI, 1.06-2.78; P = .03), c.772 to c.1806 (BCCR1'; RHR = 1.63; 95% CI, 1.10-2.40; P = .01), and c.7394 to c.8904 (BCCR2; RHR = 2.31; 95% CI, 1.69-3.16; P = .00002). We also identified 3 OCCRs: the first (OCCR1) spanned c.3249 to c.5681 that was adjacent to c.5946delT (6174delT; RHR = 0.51; 95% CI, 0.44-0.60; P = 6 × 10(-17)). The second OCCR spanned c.6645 to c.7471 (OCCR2; RHR = 0.57; 95% CI, 0.41-0.80; P = .001). Mutations conferring nonsense-mediated decay were associated with differential breast or ovarian cancer risks and an earlier age of breast cancer diagnosis for both BRCA1 and BRCA2 mutation carriers.

CONCLUSIONS AND RELEVANCE

Breast and ovarian cancer risks varied by type and location of BRCA1/2 mutations. With appropriate validation, these data may have implications for risk assessment and cancer prevention decision making for carriers of BRCA1 and BRCA2 mutations.