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

This corrects the article DOI: 10.1038/ncomms5999.

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.

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.

Validation of Intratumoral T-bet+ Lymphoid Cells as Predictors of Disease-Free Survival in Breast Cancer

We previously observed T-bet(+) lymphocytes to be associated with a good prognosis in a cohort of women with familial breast cancer. To validate this finding, we evaluated lymphocyte T-bet expression in an independent unselected prospectively accrued series of women with lymph node-negative breast carcinoma. T-bet and clinicopathologic data were available for 614 women. Hormone receptors, HER2, Ki-67, CK5, EGFR, p53, and T-bet status were determined using IHC and/or biochemical methods. Tumors were assigned to luminal A, luminal B, HER2, and basal subtypes based on the expression of IHC markers. Multiple cutpoints were examined in a univariate penalized Cox model to stratify tumors into T-bet(+/high) and T-bet(-/low). Fisher exact test was used to analyze T-bet associations with clinicopathologic variables, IHC markers, and molecular subtype. Survival analyses were by the Cox proportional hazards model. All tests were two sided. A test with a P value < 0.05 was considered statistically significant. T-bet(+/high) tumor status was significantly associated with large tumor size, high grade, hormone receptor negativity, CK5, EGFR and p53 positivity, high Ki-67, and basal subtype. With a median follow-up of 96.5 months, T-bet(-/low) tumor status was associated with a reduced disease-free survival compared with T-bet(+/high) tumor status in multivariate analysis (P = 0.0027; relative risk = 5.62; 95% confidence intervals, 1.48-50.19). Despite being associated with adverse clinicopathologic characteristics, T-bet(+) tumor-infiltrating lymphoid cells are associated with a favorable outcome. This supports their role in Th1-mediated antitumor activity and may provide insight for the development of new therapeutic strategies.

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.

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.

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.

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.

Identification of six new susceptibility loci for invasive epithelial ovarian cancer

Genome-wide association studies (GWAS) have identified 12 epithelial ovarian cancer (EOC) susceptibility alleles. The pattern of association at these loci is consistent in BRCA1 and BRCA2 mutation carriers who are at high risk of EOC. After imputation to 1000 Genomes Project data, we assessed associations of 11 million genetic variants with EOC risk from 15,437 cases unselected for family history and 30,845 controls and from 15,252 BRCA1 mutation carriers and 8,211 BRCA2 mutation carriers (3,096 with ovarian cancer), and we combined the results in a meta-analysis. This new study design yielded increased statistical power, leading to the discovery of six new EOC susceptibility loci. Variants at 1p36 (nearest gene, WNT4), 4q26 (SYNPO2), 9q34.2 (ABO) and 17q11.2 (ATAD5) were associated with EOC risk, and at 1p34.3 (RSPO1) and 6p22.1 (GPX6) variants were specifically associated with the serous EOC subtype, all with P < 5 × 10(-8). Incorporating these variants into risk assessment tools will improve clinical risk predictions for BRCA1 and BRCA2 mutation carriers.

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.

Candidate genetic modifiers for breast and ovarian cancer risk in BRCA1 and BRCA2 mutation carriers

BACKGROUND

BRCA1 and BRCA2 mutation carriers are at substantially increased risk for developing breast and ovarian cancer. The incomplete penetrance coupled with the variable age at diagnosis in carriers of the same mutation suggests the existence of genetic and nongenetic modifying factors. In this study, we evaluated the putative role of variants in many candidate modifier genes.

METHODS

Genotyping data from 15,252 BRCA1 and 8,211 BRCA2 mutation carriers, for known variants (n = 3,248) located within or around 445 candidate genes, were available through the iCOGS custom-designed array. Breast and ovarian cancer association analysis was performed within a retrospective cohort approach.

RESULTS

The observed P values of association ranged between 0.005 and 1.000. None of the variants was significantly associated with breast or ovarian cancer risk in either BRCA1 or BRCA2 mutation carriers, after multiple testing adjustments.

CONCLUSION

There is little evidence that any of the evaluated candidate variants act as modifiers of breast and/or ovarian cancer risk in BRCA1 or BRCA2 mutation carriers.

IMPACT

Genome-wide association studies have been more successful at identifying genetic modifiers of BRCA1/2 penetrance than candidate gene studies.

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.

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.

Abstract B110: Characterization of the role of TRAF2 and TNFRSF14 in breast cancer

Human breast cancer patients display significant diversity in terms of their survival, recurrence, metastasis as well as response to treatment. Our group recruited a multi-centre prospective cohort of over 1900 women with ANN breast cancer and followed them clinically. Tumors from 135 women in this group were selected for gene expression analysis on human cDNA microarrays to examine the molecular alterations in ANN patients that could be predictive of recurrence. From genes that demonstrated a statistically significant difference in expression in tumors from patients who have experienced an early recurrence versus who have remained disease free for over 10 years, TRAF2 (tumor necrosis factor (TNF) receptor-associated factor 2) was found to be highly expressed in early recurrence. TRAF2 is involved in immune system cytokine signaling pathways and interacts with many receptors in the TNF receptor superfamily and mediates the survival effects of these receptors including NF-κB, LTβR, JNK and MAPKs kinases. TRAF2 has been reported to be overexpressed or constitutively phosphorylated in some cancers, including lymphoma and pancreatic cancers.

The tumor necrosis factor receptor TNFRSF14, also called LIGHTR or HVEM is located on the short arm of chromosome 1p36. Deletion of 1p36 is especially frequent in breast tumors and is associated with progression and lymph node metastasis, poor prognosis, higher rate of recurrence, larger tumor size, and DNA aneuploidy7. Cheung and al. (2010) have described somatic TNFRSF14 mutations in Follicular Lymphoma (FL) which are associated with worse prognosis. They identified 46 cases (18.3%) with non-synonymous mutations affecting TNFRSF14 in a cohort of 251 FL patients. The significance of such acquired mutations on the structure-function relationship of TNFRSF14 is still unknown, but it appears that TNFRSF14 can be considered as a potential candidate gene that might contribute to FL development.

TNFRSF14 belongs to the TNF-receptor superfamily and contains multiple copies of a cysteine-rich motif in the extracellular region that is known to provide the ligand recognition motifs. Activation of TNF-receptor family members is induced by their cognate, trivalent, TNF-related ligand, clustering multiple receptors, which in turn recruits the cytosolic TNF receptor-associated factors (TRAFs) adaptors leading to serine kinase activation of the NFκB system. TNFRSF14 associated with TRAF2 and other members of the TRAF family and activates JNK as well as NF κB and AP-1, which control expression of multiple immune, inflammatory and acute phase action genes in response to infection or cellular stress. The short cytoplasmic regions are frequently capable of stimulating the transcription factors activator protein 1 (AP1) and NFKB1. This particular region of the TNFRSF14 gene interacts with members of the TNFR-associated factor (TRAF) family and mediates the signal transduction pathways that activate the immune response.

Since TNFRSF14 interacts with TRAF2 and/or other TRAFs to mediate the NF-κB activation and TRAF2 is highly expressed in breast cancer early recurrence, we hypothesized that altered TRAF2 and TNFRSF14 could contribute to recurrence in ANN breast cancer.

Single-strand conformation polymorphism (SSCP) and DNA Sequencing were performed to detect any small alteration in PCR-amplified product from 130 ANN breast cancer tumor samples and 48 germline breast cancer tumor samples. Known polymorphisms were excluded and direct sequencing revealed one new coding SNP (Pro262Leu; C&gt;T) located at exon 8. The genotypes CC were observed in 177 of the 178 samples in total. Only one ANN breast cancer tumor sample (936A) revealed a heterozygote genotype (CT). According to algorithms ‘SIFT’ and PolyPhen-2, the Pro262Leu amino acid exchange alters the secondary structure of TNFRSF14 and is predicted to affect the function of the protein. Due to the impact of TNFRSF14 mutations associated with poor prognosis during carcinogenesis, this variant may play a role in the therapy outcome in breast cancer.

Quantitative Real-time RT-PCR was performed to quantify TRAF2 expression in a subset of cell lines and tumor samples. Current studies include immunohistochemistry of tissue microarrays containing 888 ANN tumors to detect TRAF2 protein levels and localization and investigation of the association of TRAF2 expression with clinical outcome. In addition, the characterization of the importance of alterations in TRAF2 gene and its associated pathways is being performed using a combination of bioinformatics, to identify potential. Furthermore, the TNFRSF14 Pro262Leu_C&gt;T SNP is being analysed using TaqMan allelic discrimination to investigate a putative association of this variant with breast cancer risk or recurrence. These studies could lead to the discovery of new prognostic markers that may improve our understanding of the molecular pathways responsible for breast cancer development and progression and identify novel targets for breast cancer therapeutics.

Citation Format: Sandrine Tchatchou, Lucine Bosnoyan-Collins, Dushanthi Pinnaduwage, Shelley B. Bull, Irene L. Andrulis. Characterization of the role of TRAF2 and TNFRSF14 in breast cancer. [abstract]. In: Proceedings of the AACR Special Conference on Advances in Breast Cancer Research: Genetics, Biology, and Clinical Applications; Oct 3-6, 2013; San Diego, CA. Philadelphia (PA): AACR; Mol Cancer Res 2013;11(10 Suppl):Abstract nr B110.

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.

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.

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.

Molecular and clinical characterization of an in frame deletion of uncertain clinical significance in the BRCA2 gene

In this study, we analyzed a "variant of uncertain significance" (VUS) located in exon 23 of the BRCA2 gene exhibited by six members of five distinct families with hereditary breast cancer (BC). The variant was identified by DNA sequencing, and cDNA analysis revealed its co-expression with wild-type mRNA. We analyzed co-occurrence with other pathological mutations in BRCA1/2, performed a case-control study, looked for evolutionary data and used in-silico analyses to predict its potential clinical significance. Sequencing revealed an in frame deletion of 126 nucleotides in exon 23, leading to a deletion of 42 amino acids (c.9203_9328del126, p.Pro2992_Thr3033del). All of the VUS-carriers suffered from either BC or ovarian/pancreatic cancer. No other definite pathologic mutation of BRCA genes was found in the five families. The identified deletion could not be observed in a control cohort of 2,652 healthy individuals, but in 5 out of 916 (0.5%) tested BC families without a bona fide pathogenic BRCA1/2 mutation (P = 0.0011). According to these results, the in frame deletion c.9203_9328del126 is a rare mutation strongly associated with familial BC. In summary, our investigations indicate that this BRCA2 deletion is pathogenic.

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

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

Identification of a DMBT1 polymorphism associated with increased breast cancer risk and decreased promoter activity

According to present estimations, the unfavorable combination of alleles with low penetrance but high prevalence in the population might account for the major part of hereditary breast cancer risk. Deleted in Malignant Brain Tumors 1 (DMBT1) has been proposed as a tumor suppressor for breast cancer and other cancer types. Genomewide mapping in mice further identified Dmbt1 as a potential modulator of breast cancer risk. Here, we report the association of two frequent and linked single-nucleotide polymorphisms (SNPs) with increased breast cancer risk in women above the age of 60 years: DMBT1 c.-93C>T, rs2981745, located in the DMBT1 promoter; and DMBT1 c.124A>C, p.Thr42Pro, rs11523871(odds ratio [OR]=1.66, 95% confidence interval [CI]=1.21-2.29, P=0.0017; and OR=1.66; 95% CI=1.21-2.28, P=0.0016, respectively), based on 1,195 BRCA1/2 mutation-negative German breast cancer families and 1,466 unrelated German controls. Promoter studies in breast cancer cells demonstrate that the risk-increasing DMBT1 -93T allele displays significantly decreased promoter activity compared to the DMBT1 -93C allele, resulting in a loss of promoter activity. The data suggest that DMBT1 polymorphisms in the 5'-region are associated with increased breast cancer risk. In accordance with previous results, these data link decreased DMBT1 levels to breast cancer risk.

Nuclear receptor coregulator SNP discovery and impact on breast cancer risk

Background

Coregulator proteins are "master regulators", directing transcriptional and posttranscriptional regulation of many target genes, and are critical in many normal physiological processes, but also in hormone driven diseases, such as breast cancer. Little is known on how genetic changes in these genes impact disease development and progression. Thus, we set out to identify novel single nucleotide polymorphisms (SNPs) within SRC-1 (NCoA1), SRC-3 (NCoA3, AIB1), NCoR (NCoR1), and SMRT (NCoR2), and test the most promising SNPs for associations with breast cancer risk.

Methods

The identification of novel SNPs was accomplished by sequencing the coding regions of these genes in 96 apparently normal individuals (48 Caucasian Americans, 48 African Americans). To assess their association with breast cancer risk, five SNPs were genotyped in 1218 familial BRCA1/2-mutation negative breast cancer cases and 1509 controls (rs1804645, rs6094752, rs2230782, rs2076546, rs2229840).

Results

Through our resequencing effort, we identified 74 novel SNPs (30 in NCoR, 32 in SMRT, 10 in SRC-3, and 2 in SRC-1). Of these, 8 were found with minor allele frequency (MAF) >5% illustrating the large amount of genetic diversity yet to be discovered. The previously shown protective effect of rs2230782 in SRC-3 was strengthened (OR = 0.45 [0.21-0.98], p = 0.04). No significant associations were found with the other SNPs genotyped.

Conclusions

This data illustrates the importance of coregulators, especially SRC-3, in breast cancer development and suggests that more focused studies, including functional analyses, should be conducted.

Five polymorphisms and breast cancer risk: results from the Breast Cancer Association Consortium

Previous studies have suggested that minor alleles for ERCC4 rs744154, TNF rs361525, CASP10 rs13010627, PGR rs1042838, and BID rs8190315 may influence breast cancer risk, but the evidence is inconclusive due to their small sample size. These polymorphisms were genotyped in more than 30,000 breast cancer cases and 30,000 controls, primarily of European descent, from 30 studies in the Breast Cancer Association Consortium. We calculated odds ratios (OR) and 95% confidence intervals (95% CI) as a measure of association. We found that the minor alleles for these polymorphisms were not related to invasive breast cancer risk overall in women of European descent: ECCR4 per-allele OR (95% CI) = 0.99 (0.97-1.02), minor allele frequency = 27.5%; TNF 1.00 (0.95-1.06), 5.0%; CASP10 1.02 (0.98-1.07), 6.5%; PGR 1.02 (0.99-1.06), 15.3%; and BID 0.98 (0.86-1.12), 1.7%. However, we observed significant between-study heterogeneity for associations with risk for single-nucleotide polymorphisms (SNP) in CASP10, PGR, and BID. Estimates were imprecise for women of Asian and African descent due to small numbers and lower minor allele frequencies (with the exception of BID SNP). The ORs for each copy of the minor allele were not significantly different by estrogen or progesterone receptor status, nor were any significant interactions found between the polymorphisms and age or family history of breast cancer. In conclusion, our data provide persuasive evidence against an overall association between invasive breast cancer risk and ERCC4 rs744154, TNF rs361525, CASP10 rs13010627, PGR rs1042838, and BID rs8190315 genotypes among women of European descent.

Chromosome copy number variation and breast cancer risk

Breast cancer is the most frequent cancer among women. It is caused by genetic and environmental factors. Whereas mutations in high-penetrance susceptibility genes have been identified in familial breast cancer and several single nucleotide polymorphisms (SNPs) have been shown to be associated with both familial and sporadic breast cancer risk, the impact of genomic copy number variants (CNVs) on breast cancer risk has so far poorly been studied. An example of a CNV affecting the tumor suppressor gene MTUS1 that has been shown to be associated with familial breast cancer risk is given. Moreover, we discuss CNVs affecting detoxification genes like GSTM1 and GSTT1 whose association with breast cancer risk is controversial. Finally, the potential of array-based genome-wide CNV association studies is discussed.

A variant affecting a putative miRNA target site in estrogen receptor (ESR) 1 is associated with breast cancer risk in premenopausal women

MicroRNAs (miRNAs) negatively regulate expression of target transcripts by hybridization to complementary sites of their messenger RNA targets. Chen et al. have described several putative functional single nucleotide polymorphisms (SNPs) in miRNA target sites. Here, we selected 11 miRNA target site SNPs located in 3' untranslated regions of genes involved in cancer and breast cancer to analyze their impact on breast cancer risk using a large familial study population. Whereas no association was observed for 10 SNPs, a significant association was revealed for the variant affecting a miRNA target site in the estrogen receptor (ESR) 1. Age stratification showed that the association was stronger in premenopausal women [C versus T: odds ratio (OR) = 0.60, confidence interval (CI) = 0.41-0.89, P = 0.010]. Furthermore, the effect was stronger in high-risk familial cases (C versus T: OR = 0.42, CI = 0.25-0.71, P = 0.0009). Clinical studies have shown that elimination of ESR1 significantly reduces breast cancer risk. Thus, therapies that inhibit ESR1 are used for breast cancer treatment. According to in silico analysis, ESR1_rs2747648 affects the binding capacity of miR-453, which is stronger when the C allele is present. In contrast, the T allele attenuates the binding of miR-453, which might lead to a reduced miRNA-mediated ESR1 repression, in consequence higher ESR1 protein levels and an increased breast cancer risk. Thus, the breast cancer protective effect observed for the C allele in premenopausal women is biologically reasonable. The analysis of large study populations in multicentre collaboration will be needed to verify the association and answer questions regarding the possible impact of this variant on therapeutic and clinical outcome.

Aurora kinases A and B and familial breast cancer risk

Aurora genes play a crucial role in tumourigenesis and are overexpressed in many kinds of cancers. We investigated whether coding variants within the Aurora genes are associated with familial breast cancer risk. While AURKA Phe31Ile (1712T>A) and AURKB Thr298Met (893G>A) showed no association, the synonymous AURKB Ser295Ser (885A>G) polymorphism resulted in an increased breast cancer risk for carriers of the homozygous 885G genotype (OR=1.45, 95% CI=1.05-2.0, P=0.02). Due to the impact of aurora kinases in the loss of chromosomal integrity during carcinogenesis, this variant may also influence the therapy outcome in breast cancer.

Associations of genetic variants in the estrogen receptor coactivators PPARGC1A, PPARGC1B and EP300 with familial breast cancer

The mitogen effect of the ovarian steroid estrogen is a strong risk factor for breast cancer development. This effect is mainly mediated by the estrogen receptor alpha, a hormone inducible transcription factor, which activates gene expression through recruiting multiple coactivators, such as PPARGC1A, PPARGC1B and EP300. We tested the hypothesis that non-conservative, putative functional amino acid exchanges in PPARGC1A, PPARGC1B and EP300 act as low-penetrance familial breast cancer risk factors. The analysis of 816 BRCA1/2 mutation-negative familial breast cancer patients and 1012 controls revealed an association of the PPARGC1A Thr612Met polymorphism with familial breast cancer (OR = 1.35, 95% CI 1.00-1.81, P = 0.049), high-risk familial breast cancer (OR = 1.51, 95% CI 1.08-2.12, P = 0.017) and bilateral familial breast cancer (OR = 2.30, 95% CI 1.24-4.28, P = 0.009). Logistic regression analyses of the PPARGC1B Ala203Pro variant showed an increased familial breast cancer risk of heterozygous and homozygous variant allele carriers (OR = 1.48, 95% CI 1.15-1.91, P = 0.002). The genotype-combination analysis of the associated PPARGC1A Thr612Met variant and the associated PPARGC1B Ala203Pro variant suggests an allele dose-dependent breast cancer risk (P(trend) = 0.0004). Our results indicate for the first time the importance of inherited variants in the estrogen receptor coactivator genes PPARGC1A and PPARGC1B for familial breast cancer susceptibility. Owing to their impact on estrogen signaling, these polymorphisms might also influence adjuvant anti-estrogen therapy, using agents such as tamoxifen and raloxifen, and outcome of breast cancer patients.

Association of genetic variants in the Rho guanine nucleotide exchange factor AKAP13 with familial breast cancer

The A-kinase anchor protein 13 (AKAP13, alias BRX and lbc) tethers cAMP-dependent protein kinase to its subcellular environment and catalyses Rho GTPases activity as a guanine nucleotide exchange factor. The crucial role of members of the Rho family of GTPases in carcinogenesis is well established and targeting Rho proteins with antineoplastic compounds has become a major effort in the fight against cancer. Thus, genetic alterations within the candidate cancer susceptibility gene AKAP13 would be expected to provoke a constitutive Rho signalling, thereby facilitating the development of cancer. Here, we analysed the potential impact of four polymorphic non-conservative amino acid exchanges (Arg494Trp, Lys526Gln, Asn1086Asp and Gly2461Ser) in AKAP13 on familial breast cancer. We performed a case-control study using genomic DNA of BRCA1/2 mutation-negative German female index patients from 601 unrelated families, among a subset of 356 high-risk families, and 1053 German female unrelated controls. The newfound Lys526Gln polymorphism revealed a significant association with familial breast cancer (OR = 1.58, 95% CI = 1.07-2.35) and an even stronger association with high-risk familial breast cancer (OR = 1.85, 95% CI = 1.19-2.88). Haplotype analyses were in line with genotype results displaying a similar significance as analyses of individual polymorphisms. Due to the pivotal role of AKAP13 in the Rho GTPases signalling network, this variant might affect the susceptibility to other cancers as well.