Investigation of gene-environment interactions between 47 newly identified breast cancer susceptibility loci and environmental risk factors

A large genotyping project within the Breast Cancer Association Consortium (BCAC) recently identified 41 associations between single nucleotide polymorphisms (SNPs) and overall breast cancer (BC) risk. We investigated whether the effects of these 41 SNPs, as well as six SNPs associated with estrogen receptor (ER) negative BC risk are modified by 13 environmental risk factors for BC. Data from 22 studies participating in BCAC were pooled, comprising up to 26,633 cases and 30,119 controls. Interactions between SNPs and environmental factors were evaluated using an empirical Bayes-type shrinkage estimator. Six SNPs showed interactions with associated p-values (pint ) <1.1 × 10(-3) . None of the observed interactions was significant after accounting for multiple testing. The Bayesian False Discovery Probability was used to rank the findings, which indicated three interactions as being noteworthy at 1% prior probability of interaction. SNP rs6828523 was associated with increased ER-negative BC risk in women ≥170 cm (OR = 1.22, p = 0.017), but inversely associated with ER-negative BC risk in women <160 cm (OR = 0.83, p = 0.039, pint = 1.9 × 10(-4) ). The inverse association between rs4808801 and overall BC risk was stronger for women who had had four or more pregnancies (OR = 0.85, p = 2.0 × 10(-4) ), and absent in women who had had just one (OR = 0.96, p = 0.19, pint = 6.1 × 10(-4) ). SNP rs11242675 was inversely associated with overall BC risk in never/former smokers (OR = 0.93, p = 2.8 × 10(-5) ), but no association was observed in current smokers (OR = 1.07, p = 0.14, pint = 3.4 × 10(-4) ). In conclusion, recently identified BC susceptibility loci are not strongly modified by established risk factors and the observed potential interactions require confirmation in independent studies.

Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with prognosis of estrogen receptor-negative breast cancer after chemotherapy

INTRODUCTION

Tumor lymphocyte infiltration is associated with clinical response to chemotherapy in estrogen receptor (ER) negative breast cancer. To identify variants in immunosuppressive pathway genes associated with prognosis after adjuvant chemotherapy for ER-negative patients, we studied stage I-III invasive breast cancer patients of European ancestry, including 9,334 ER-positive (3,151 treated with chemotherapy) and 2,334 ER-negative patients (1,499 treated with chemotherapy).

METHODS

We pooled data from sixteen studies from the Breast Cancer Association Consortium (BCAC), and employed two independent studies for replications. Overall 3,610 single nucleotide polymorphisms (SNPs) in 133 genes were genotyped as part of the Collaborative Oncological Gene-environment Study, in which phenotype and clinical data were collected and harmonized. Multivariable Cox proportional hazard regression was used to assess genetic associations with overall survival (OS) and breast cancer-specific survival (BCSS). Heterogeneity according to chemotherapy or ER status was evaluated with the log-likelihood ratio test.

RESULTS

Three independent SNPs in TGFBR2 and IL12B were associated with OS (P <10⁻³) solely in ER-negative patients after chemotherapy (267 events). Poorer OS associated with TGFBR2 rs1367610 (G > C) (per allele hazard ratio (HR) 1.54 (95% confidence interval (CI) 1.22 to 1.95), P = 3.08 × 10⁻⁴) was not found in ER-negative patients without chemotherapy or ER-positive patients with chemotherapy (P for interaction <10-3). Two SNPs in IL12B (r² = 0.20) showed different associations with ER-negative disease after chemotherapy: rs2546892 (G > A) with poorer OS (HR 1.50 (95% CI 1.21 to 1.86), P = 1.81 × 10⁻⁴), and rs2853694 (A > C) with improved OS (HR 0.73 (95% CI 0.61 to 0.87), P = 3.67 × 10⁻⁴). Similar associations were observed with BCSS. Association with TGFBR2 rs1367610 but not IL12B variants replicated using BCAC Asian samples and the independent Prospective Study of Outcomes in Sporadic versus Hereditary Breast Cancer Study and yielded a combined HR of 1.57 ((95% CI 1.28 to 1.94), P = 2.05 × 10⁻⁵) without study heterogeneity.

CONCLUSIONS

TGFBR2 variants may have prognostic and predictive value in ER-negative breast cancer patients treated with adjuvant chemotherapy. Our findings provide further insights into the development of immunotherapeutic targets for ER-negative breast cancer.

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.

Refined histopathological predictors of BRCA1 and BRCA2 mutation status: a large-scale analysis of breast cancer characteristics from the BCAC, CIMBA, and ENIGMA consortia

INTRODUCTION

The distribution of histopathological features of invasive breast tumors in BRCA1 or BRCA2 germline mutation carriers differs from that of individuals with no known mutation. Histopathological features thus have utility for mutation prediction, including statistical modeling to assess pathogenicity of BRCA1 or BRCA2 variants of uncertain clinical significance. We analyzed large pathology datasets accrued by the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and the Breast Cancer Association Consortium (BCAC) to reassess histopathological predictors of BRCA1 and BRCA2 mutation status, and provide robust likelihood ratio (LR) estimates for statistical modeling.

METHODS

Selection criteria for study/center inclusion were estrogen receptor (ER) status or grade data available for invasive breast cancer diagnosed younger than 70 years. The dataset included 4,477 BRCA1 mutation carriers, 2,565 BRCA2 mutation carriers, and 47,565 BCAC breast cancer cases. Country-stratified estimates of the likelihood of mutation status by histopathological markers were derived using a Mantel-Haenszel approach.

RESULTS

ER-positive phenotype negatively predicted BRCA1 mutation status, irrespective of grade (LRs from 0.08 to 0.90). ER-negative grade 3 histopathology was more predictive of positive BRCA1 mutation status in women 50 years or older (LR = 4.13 (3.70 to 4.62)) versus younger than 50 years (LR = 3.16 (2.96 to 3.37)). For BRCA2, ER-positive grade 3 phenotype modestly predicted positive mutation status irrespective of age (LR = 1.7-fold), whereas ER-negative grade 3 features modestly predicted positive mutation status at 50 years or older (LR = 1.54 (1.27 to 1.88)). Triple-negative tumor status was highly predictive of BRCA1 mutation status for women younger than 50 years (LR = 3.73 (3.43 to 4.05)) and 50 years or older (LR = 4.41 (3.86 to 5.04)), and modestly predictive of positive BRCA2 mutation status in women 50 years or older (LR = 1.79 (1.42 to 2.24)).

CONCLUSIONS

These results refine likelihood-ratio estimates for predicting BRCA1 and BRCA2 mutation status by using commonly measured histopathological features. Age at diagnosis is an important variable for most analyses, and grade is more informative than ER status for BRCA2 mutation carrier prediction. The estimates will improve BRCA1 and BRCA2 variant classification and inform patient mutation testing and clinical management.

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.

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 2188: Genetic variation in angiogenesis-related genes is associated with colorectal cancer risk and prognosis

Introduction

Angiogenesis, the generation of new blood vessels, is crucial in tumor growth, progression, and metastasis. The process involves a variety of factors including signaling, adhesion and chemotactic molecules, ECM (extracellular matrix) proteins, proteinases, transcription factors, growth factors and receptors. Variation in genes encoding these factors potentially influence angiogenic processes. We thus evaluated the association between variants in angiogenesis-related genes and risk as well as prognosis of colorectal cancer (CRC).

Approach

In a candidate pathway approach, we investigated 437 variants in 36 angiogenesis-related genes for association with CRC risk and prognosis in ∼1800 patients and ∼1800 controls of the German DACHS/IMPACT study. Patients were aged 30 years or older and diagnosed between 2003 and 2007. CRC risk was estimated using conditional logistic regression based on the co-dominant inheritance model. In addition, we investigated the association between polymorphisms and overall survival using multivariable Cox regression. Correction for multiple testing was performed using false discovery rates (FDR).

Results

Risk

Several of the investigated variants in ANGTP1, ETS1, FLT4, MMP2, NOTCH4, PDGFRB, TGFB2 were associated with risk of CRC, and one variant in DLL1 (delta-like 1) remained significant after accounting for multiple testing. DLL1 rs9348307 was associated with decreased risk of CRC (ORGC 0.81, 95% CI 0.68-0.96; ORCC 0.62, 95% CI 0.36-1.06, p=0.02).

Survival

After a median follow up time of 5 years, variants in twelve genes (ANGPT1, EFNB2, ETS1, FLT4, JAG1, KDR, MMP2, MMP9, NRP1, NRP2, PDGFRB, TGFB2) were associated with overall survival in CRC patients. When accounting for multiple testing, variants in EFNB2, JAG1 and MMP2 remained significant.

Most of the EFNB2 variants were associated with poorer survival, while one variant in EFNB2 was associated with better survival (rs2391333: HRCT 0.88, 95% CI 0.74-1.05; HRTT 0.60, 95% CI 0.45-0.78, p&lt;0.001). Furthermore, two variants in JAG1 were associated with worse survival in CRC patients (rs3748480: HRTG 1.27, 95% CI 1.05-1.53; HRGG 1.79, 95% CI 1.02-3.13, p=0.02; rs6040062: HRGA 1.26, 95% CI 1.04-1.52; HRAA 1.78, 95% CI 1.02-3.12, p=0.02). Finally, three variants in MMP2 (rs11639960, rs1561219, rs17301608) were associated with poorer survival, while one MMP2 variant was associated with improved survival in CRC patients (rs243847: HRTC 0.87, 95% CI 0.73-1.04, HRCC 0.60, 95% CI 0.46-0.78, p&lt;0.001).

Conclusion

Variants in angiogenesis-related genes were associated with CRC risk and overall survival in CRC patients. Two of the identified genes (DLL1 and JAG1) are involved in NOTCH signaling, which is crucial to cell differentiation, proliferation, apoptosis and angiogenesis, underlining the importance of this oncogenic pathway in colorectal carcinogenesis.

Citation Format: Dominique Scherer, Yesilda Balavarca, Nina Habermann, Katharina Buck, Petra Seibold, Lisanne Kap, Katja Butterbach, Katrin Pfütze, Axel Benner, Michael Hoffmeister, Hermann Brenner, Barbara Burwinkel, Jenny Chang-Claude, Cornelia M. Ulrich. Genetic variation in angiogenesis-related genes is associated with colorectal cancer risk and prognosis. [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 2188. doi:10.1158/1538-7445.AM2014-2188

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.

2q36.3 is associated with prognosis for oestrogen receptor-negative breast cancer patients treated with chemotherapy

Large population-based registry studies have shown that breast cancer prognosis is inherited. Here we analyse single-nucleotide polymorphisms (SNPs) of genes implicated in human immunology and inflammation as candidates for prognostic markers of breast cancer survival involving 1,804 oestrogen receptor (ER)-negative patients treated with chemotherapy (279 events) from 14 European studies in a prior large-scale genotyping experiment, which is part of the Collaborative Oncological Gene-environment Study (COGS) initiative. We carry out replication using Asian COGS samples (n=522, 53 events) and the Prospective Study of Outcomes in Sporadic versus Hereditary breast cancer (POSH) study (n=315, 108 events). Rs4458204_A near CCL20 (2q36.3) is found to be associated with breast cancer-specific death at a genome-wide significant level (n=2,641, 440 events, combined allelic hazard ratio (HR)=1.81 (1.49-2.19); P for trend=1.90 × 10(-9)). Such survival-associated variants can represent ideal targets for tailored therapeutics, and may also enhance our current prognostic prediction capabilities.

Alcohol consumption and survival after a breast cancer diagnosis: a literature-based meta-analysis and collaborative analysis of data for 29,239 cases

BACKGROUND

Evidence for an association of alcohol consumption with prognosis after a diagnosis of breast cancer has been inconsistent. We have reviewed and summarized the published evidence and evaluated the association using individual patient data from multiple case cohorts.

METHODS

A MEDLINE search to identify studies published up to January 2013 was performed. We combined published estimates of survival time for "moderate drinkers" versus nondrinkers. An analysis of individual participant data using Cox regression was carried out using data from 11 case cohorts.

RESULTS

We identified 11 published studies suitable for inclusion in the meta-analysis. Moderate postdiagnosis alcohol consumption was not associated with overall survival [HR, 0.95; 95% confidence interval (CI), 0.85-1.05], but there was some evidence of better survival associated with prediagnosis consumption (HR, 0.80; 95% CI, 0.73-0.88). Individual data on alcohol consumption for 29,239 cases with 4,839 deaths were available from the 11 case cohorts, all of which had data on estrogen receptor (ER) status. For women with ER-positive disease, there was little evidence that pre- or postdiagnosis alcohol consumption is associated with breast cancer-specific mortality, with some evidence of a negative association with all-cause mortality. On the basis of a single study, moderate postdiagnosis alcohol intake was associated with a small reduction in breast cancer-specific mortality for women with ER-negative disease. There was no association with prediagnosis intake for women with ER-negative disease.

CONCLUSION

There was little evidence that pre- or post-diagnosis alcohol consumption is associated with breast cancer-specific mortality for women with ER-positive disease. There was weak evidence that moderate post-diagnosis alcohol intake is associated with a small reduction in breast cancer-specific mortality in ER-negative disease.

IMPACT

Considering the totality of the evidence, moderate postdiagnosis alcohol consumption is unlikely to have a major adverse effect on the survival of women with breast cancer.

Large-scale evaluation of common variation in regulatory T cell-related genes and ovarian cancer outcome

The presence of regulatory T cells (Treg) in solid tumors is known to play a role in patient survival in ovarian cancer and other malignancies. We assessed inherited genetic variations via 749 tag single-nucleotide polymorphisms (SNP) in 25 Treg-associated genes (CD28, CTLA4, FOXP3, IDO1, IL10, IL10RA, IL15, 1L17RA, IL23A, IL23R, IL2RA, IL6, IL6R, IL8, LGALS1, LGALS9, MAP3K8, STAT5A, STAT5B, TGFB1, TGFB2, TGFB3, TGFBR1, TGRBR2, and TGFBR3) in relation to ovarian cancer survival. We analyzed genotype and overall survival in 10,084 women with invasive epithelial ovarian cancer, including 5,248 high-grade serous, 1,452 endometrioid, 795 clear cell, and 661 mucinous carcinoma cases of European descent across 28 studies from the Ovarian Cancer Association Consortium (OCAC). The strongest associations were found for endometrioid carcinoma and IL2RA SNPs rs11256497 [HR, 1.42; 95% confidence interval (CI), 1.22-1.64; P = 5.7 × 10(-6)], rs791587 (HR, 1.36; 95% CI, 1.17-1.57; P = 6.2 × 10(-5)), rs2476491 (HR, = 1.40; 95% CI, 1.19-1.64; P = 5.6 × 10(-5)), and rs10795763 (HR, 1.35; 95% CI, 1.17-1.57; P = 7.9 × 10(-5)), and for clear cell carcinoma and CTLA4 SNP rs231775 (HR, 0.67; 95% CI, 0.54-0.82; P = 9.3 × 10(-5)) after adjustment for age, study site, population stratification, stage, grade, and oral contraceptive use. The rs231775 allele associated with improved survival in our study also results in an amino acid change in CTLA4 and previously has been reported to be associated with autoimmune conditions. Thus, we found evidence that SNPs in genes related to Tregs seem to play a role in ovarian cancer survival, particularly in patients with clear cell and endometrioid epithelial ovarian cancer.

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.

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.

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.

Enterolactone concentrations and prognosis after postmenopausal breast cancer: assessment of effect modification and meta-analysis

We previously reported that high concentrations of enterolactone, a lignan metabolite, are associated with lower mortality in 1,140 breast cancer patients from Germany. Using an extended set of 2,182 patients aged 50-74 years at diagnosis (2001-2005) and prospectively followed up until 2009, we investigated whether the association with mortality differs by lifestyle factors and tumor characteristics. Hazard ratios (HR) and 95% confidence intervals (CI) were estimated using multivariable Cox regression. Potential differential effects by tumor characteristics and lifestyle factors were assessed and a meta-analysis of five studies addressing lignan exposure and breast cancer prognosis was performed to summarize evidence. Median enterolactone concentrations were 17.4 (± 30.5 standard deviation) and 22.9 nmol L(-1) (± 44.8), respectively, for 269 deceased and 1,913 patients still alive. High enterolactone concentrations were significantly associated with lower all-cause mortality (per 10 nmol L(-1) : HR 0.94, 95% CI 0.90-0.98), breast cancer-specific mortality (HR 0.94, 0.89-0.99), and distant disease-free survival (HR 0.94, 0.90-0.98). Associations were found for stage 0-IIIA but not for stage IIIB-IV disease (p(het) = 0.01) and were stronger in patients with BMI <25 kg m(-2) than those with BMI ≥ 25 (p(het) = 0.04). In patients with healthy lifestyle (BMI <25, nonsmoker, physically active), the inverse association with all-cause mortality was still apparent (HR 0.92, 0.85-0.99). The meta-analysis yielded significant associations both for all-cause (HR 0.57, 0.42-0.78) and breast cancer-specific mortality (HR 0.54, 0.39-0.75). Our findings show that high lignan exposure is associated with reduced mortality in breast cancer patients. The inverse association observed in this study cannot be entirely explained by a healthy lifestyle.

Identification of new genetic susceptibility loci for breast cancer through consideration of gene-environment interactions

Genes that alter disease risk only in combination with certain environmental exposures may not be detected in genetic association analysis. By using methods accounting for gene-environment (G × E) interaction, we aimed to identify novel genetic loci associated with breast cancer risk. Up to 34,475 cases and 34,786 controls of European ancestry from up to 23 studies in the Breast Cancer Association Consortium were included. Overall, 71,527 single nucleotide polymorphisms (SNPs), enriched for association with breast cancer, were tested for interaction with 10 environmental risk factors using three recently proposed hybrid methods and a joint test of association and interaction. Analyses were adjusted for age, study, population stratification, and confounding factors as applicable. Three SNPs in two independent loci showed statistically significant association: SNPs rs10483028 and rs2242714 in perfect linkage disequilibrium on chromosome 21 and rs12197388 in ARID1B on chromosome 6. While rs12197388 was identified using the joint test with parity and with age at menarche (P-values = 3 × 10(-07)), the variants on chromosome 21 q22.12, which showed interaction with adult body mass index (BMI) in 8,891 postmenopausal women, were identified by all methods applied. SNP rs10483028 was associated with breast cancer in women with a BMI below 25 kg/m(2) (OR = 1.26, 95% CI 1.15-1.38) but not in women with a BMI of 30 kg/m(2) or higher (OR = 0.89, 95% CI 0.72-1.11, P for interaction = 3.2 × 10(-05)). Our findings confirm comparable power of the recent methods for detecting G × E interaction and the utility of using G × E interaction analyses to identify new susceptibility loci.

Circulating 25-hydroxyvitamin D and postmenopausal breast cancer survival: Influence of tumor characteristics and lifestyle factors?

We previously reported that lower post-diagnostic circulating 25-hydroxyvitamin D [25(OH)D] concentrations were associated with higher risk of overall mortality and distant disease in stage I-IV postmenopausal breast cancer survivors. This association was now re-examined in an extended dataset to investigate potential effect modification by tumor characteristics and lifestyle factors. A prospective cohort study was conducted in Germany including 2,177 incident stage I-IV postmenopausal breast cancer patients aged 50-74 years. Patients were diagnosed between 2001 and 2005 and median follow-up time was 5.3 years. Cox proportional hazards models were stratified by age at diagnosis, study center and season of blood collection and adjusted for other prognostic factors. A meta-analysis of studies on circulating 25(OH)D and mortality in breast cancer patients was performed to summarize evidence. Lower concentrations of 25(OH)D were significantly associated with higher risk of overall mortality [hazard ratio (HR) lowest vs. highest tertile = 1.86; 95% confidence interval (CI): 1.22, 2.82; p-trend = 0.002] and distant disease (HR = 1.76; 95% CI: 1.24, 2.49; p-trend = 0.003) in stage I-IIIa but not in stage IIIb-IV breast cancer patients. No significant interaction by lifestyle factors was observed (all p-interaction > 0.05). The meta-analysis yielded significant associations with overall and breast cancer-specific mortality (lowest vs. highest quantile: HR = 1.52; 95% CI: 1.22, 1.88 and HR = 1.74; 95% CI: 1.23, 2.40, respectively). In conclusion, post-diagnostic circulating 25(OH)D concentrations were associated with overall mortality and distant disease in stage I-IIIa postmenopausal breast cancer patients. This association was not strongly modified by lifestyle factors.

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.

Mortality and recurrence risk in relation to the use of lipid-lowering drugs in a prospective breast cancer patient cohort

Lipid-lowering drugs are used for the prevention of cardiovascular diseases. Statins are the most commonly used lipid-lowering drugs. Evidence from preclinical and observational studies suggests that statins might improve the prognosis of breast cancer patients. We analyzed data from the German MARIEplus study, a large prospective population-based cohort of patients aged 50 and older, who were diagnosed with breast cancer between 2001 and 2005. For overall mortality, breast-cancer specific mortality, and non-breast-cancer mortality, we included 3189 patients with invasive breast cancer stage I-IV, and for recurrence risk 3024 patients with breast cancer stage I-III. We used Cox proportional hazards models to assess the association with self-reported lipid-lowering drug use at recruitment. We stratified by study region, tumor grade, and estrogen/progesterone receptor status, and adjusted for age, tumor size, nodal status, metastases (stage I-IV only), menopausal hormone treatment, mode of detection, radiotherapy, and smoking. Mortality analyses were additionally adjusted for cardiovascular disease, diabetes mellitus and body-mass index. During a median follow-up of 5.3 years, 404 of 3189 stage I-IV patients died, and 286 deaths were attributed to breast cancer. Self-reported use of lipid-lowering drugs was non-significantly associated with increased non-breast cancer mortality (Hazard ratio (HR) 1.49, 95% confidence interval (CI) 0.88-2.52) and increased overall mortality (HR 1.21, 95% CI 0.87-1.69) whereas no association with breast cancer-specific mortality was found (HR 1.04, 0.67-1.60). Restricted to stage I-III breast cancer patients, 387 recurrences occurred during a median follow-up of 5.4 years. We found lipid-lowering drug use to be non-significantly associated with a reduced risk of recurrence (HR 0.83, 95% CI 0.54-1.24) and of breast cancer-specific mortality (HR 0.89, 95% CI 0.52-1.49). Although compatible with previous findings of an improved prognosis associated with statin use, our results do not provide clear supportive evidence for an association with lipid-lowering drug use due to imprecise estimates.

Intrinsic breast cancer subtypes defined by estrogen receptor signalling-prognostic relevance of progesterone receptor loss

The majority of luminal type breast carcinomas are slowly growing tumors with an overall favorable prognosis. However, a proportion of cases (luminal B tumors) are characterized by coactivation of growth factor receptors or non-canonical ER signaling and a poorer clinical outcome. The aim of our study was to evaluate whether the expression of proteins that are part of the ER signaling network may be used to distinguish low-risk from high-risk luminal tumors. Unsupervised hierarchical clustering of a set of proteins either involved in estrogen receptor signaling or associated with resistance to endocrine therapy was performed in a series of 443 postmenopausal breast carcinomas. Using this approach, we were able to reproduce the established classification with two distinct groups of luminal (estrogen receptor positive) tumors, one group of HER2-associated tumors and a group of triple-negative tumors. However, neither proliferation nor the expression of one or more of the ER-co-factors or resistance-associated factors, but PR-expression was identified as the most important stratifier distinguishing between the two luminal groups. In fact, not only the four identified clusters were shown to be significantly associated with patient outcome, PR-expression alone or in combination with Ki-67-stains stratified ER-positive tumors into a low-risk and a high-risk group. Our data indicate that defining luminal B tumors by the presence of high-risk criteria (loss of PR-expression or increased proliferation) provides a robust and highly significant stratification of ER-positive breast carcinomas into luminal A and B.

The UGT1A6_19_GG genotype is a breast cancer risk factor

Validation of an association between the UGT1A6_19_T>G (rs6759892) polymorphism and overall breast cancer risk. A pilot study included two population-based case-control studies from Germany (MARIE-GENICA). An independent validation study comprised four independent breast cancer case-control studies from Finland (KBCP, OBCS), Germany (BBCC), and Sweden (SASBAC). The pooled analysis included 7418 cases and 8720 controls from all six studies. Participants were of European descent. Genotyping was done by MALDI-TOF MS and statistical analysis was performed by logistic regression adjusted for age and study. The increased overall breast cancer risk for women with the UGT1A6_19_GG genotype which was observed in the pilot study was confirmed in the set of four independent study collections (OR 1.13, 95% CI 1.05–1.22; p = 0.001). The pooled study showed a similar effect (OR 1.09, 95% CI 1.04–1.14; p = 0.001). The risk effect on the basis of allele frequencies was highly significant, the pooled analysis showed an OR of 1.11 (95% CI 1.06–1.16; p = 5.8 × 10⁻⁶). We confirmed the association of UGT1A6_19_GG with increased overall breast cancer risk and conclude that our result from a well powered multi-stage study adds a novel candidate to the panel of validated breast cancer susceptibility loci.

Association of pre-diagnosis physical activity with recurrence and mortality among women with breast cancer

Evidence is emerging that physical activity (PA) may improve overall survival after breast cancer diagnosis. However, the effect of PA on breast cancer recurrence and on cause-specific mortality is less investigated. We assessed the association of pre-diagnosis PA with recurrence, overall and cause-specific survival in a prospective cohort study in Germany including 3,393 non-metastatic breast cancer patients aged 50-74 years. Cox proportional hazards models were calculated adjusted for relevant prognostic factors. During a median follow-up of 5.6 years, 367 patients deceased. Overall mortality was significantly inversely associated with pre-diagnosis recreational PA. However, this effect was mainly attributed to deaths due to causes other than breast cancer. Multiple fractional polynomial analyses yielded a nonlinear association with markedly increased non-breast cancer mortality for women who did not engage in any sports or cycling in the years before the breast cancer diagnosis with a hazard ratio (HR, none vs. any) of 1.71, 95% confidence interval (1.16, 2.52). There were no further risk reductions with increasing activity levels. The association with breast cancer-specific mortality showed a similar dose-response but was far less pronounced with HR (none vs. any) = 1.22 (0.91, 1.64). In contrast, regarding cancer recurrence the dose-response was linear. However, this association was restricted to estrogen/progesterone receptor-negative (ER-/PR-) cases (p interaction = 0.033) with HR (highest vs. no recreational PA) = 0.53 (0.24, 1.16), p trend = 0.0045. Thus, breast cancer patients with a physically inactive lifestyle pre-diagnosis may decease prematurely irrespective of their cancer prognosis. Higher levels of exercise may reduce the risk of recurrence of ER-/PR- breast tumors.