Type II transmembrane serine proteases (TTSPs) are related to tumor growth, invasion, and metastasis in cancer. Genetic variants in these genes may alter their function, leading to cancer onset and progression, and affect patient outcome. Here, 464 breast cancer cases and 370 controls were genotyped for 82 single-nucleotide polymorphisms covering eight genes. Association of the genotypes was estimated against breast cancer risk, breast cancer–specific survival, and survival in different treatment groups, and clinicopathological variables. SNPs in TMPRSS3 (rs3814903 and rs11203200), TMPRSS7 (rs1844925), and HGF (rs5745752) associated significantly with breast cancer risk (Ptrend = 0.008–0.042). SNPs in TMPRSS1 (rs12151195 and rs12461158), TMPRSS2 (rs2276205), TMPRSS3 (rs3814903), and TMPRSS7 (rs2399403) associated with prognosis (P = 0.004–0.046). When estimating the combined effect of the variants, the risk of breast cancer was higher with 4–5 alleles present compared to 0–2 alleles (P = 0.0001; OR, 2.34; 95% CI, 1.39–3.94). Women with 6–8 survival-associating alleles had a 3.3 times higher risk of dying of breast cancer compared to women with 1–3 alleles (P = 0.001; HR, 3.30; 95% CI, 1.58–6.88). The results demonstrate the combined effect of variants in TTSPs and their related genes in breast cancer risk and patient outcome. Functional analysis of these variants will lead to further understanding of this gene family, which may improve individualized risk estimation and development of new strategies for treatment of breast cancer.
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/m2 (OR = 1.26, 95% CI 1.15–1.38) but not in women
with a BMI of 30 kg/m2 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.
breast cancer risk; gene-environment interaction; polymorphisms; body mass index; case-control study
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.
Invasive lobular breast cancer (ILC) accounts for 10–15% of invasive breast cancer and is generally ER positive (ER+). To date, none of the genome-wide association studies that have identified loci that predispose to breast cancer in general or to ER+ or ER-negative breast cancer have focused on lobular breast cancer. In this lobular breast cancer study we identified a new variant that appears to be specific to this morphological subtype. We also ascertained which of the known variants predisposes specifically to lobular breast cancer and show for the first time that some of these loci are also associated with lobular carcinoma in situ, a non-obligate precursor of breast cancer and also a risk factor for contralateral breast cancer. Our study shows that the genetic pathways of invasive lobular cancer and ER+ ductal carcinoma mostly overlap, but there are important differences that are likely to provide insights into the biology of lobular breast tumors.
Breast cancer is the most common cancer among women. Common variants at 27 loci have been identified as associated with susceptibility to breast cancer, and these account for ~9% of the familial risk of the disease. We report here a meta-analysis of 9 genome-wide association studies, including 10,052 breast cancer cases and 12,575 controls of European ancestry, from which we selected 29,807 SNPs for further genotyping. These SNPs were genotyped in 45,290 cases and 41,880 controls of European ancestry from 41 studies in the Breast Cancer Association Consortium (BCAC). The SNPs were genotyped as part of a collaborative genotyping experiment involving four consortia (Collaborative Oncological Gene-environment Study, COGS) and used a custom Illumina iSelect genotyping array, iCOGS, comprising more than 200,000 SNPs. We identified SNPs at 41 new breast cancer susceptibility loci at genome-wide significance (P < 5 × 10−8). Further analyses suggest that more than 1,000 additional loci are involved in breast cancer susceptibility.
Estrogen receptor (ER)-negative tumors represent 20–30% of all breast cancers, with a higher proportion occurring in younger women and women of African ancestry1. The etiology2 and clinical behavior3 of ER-negative tumors are different from those of tumors expressing ER (ER positive), including differences in genetic predisposition4. To identify susceptibility loci specific to ER-negative disease, we combined in a meta-analysis 3 genome-wide association studies of 4,193 ER-negative breast cancer cases and 35,194 controls with a series of 40 follow-up studies (6,514 cases and 41,455 controls), genotyped using a custom Illumina array, iCOGS, developed by the Collaborative Oncological Gene-environment Study (COGS). SNPs at four loci, 1q32.1 (MDM4, P = 2.1 × 10−12 and LGR6, P = 1.4 × 10−8), 2p24.1 (P = 4.6 × 10−8) and 16q12.2 (FTO, P = 4.0 × 10−8), were associated with ER-negative but not ER-positive breast cancer (P > 0.05). These findings provide further evidence for distinct etiological pathways associated with invasive ER-positive and ER-negative breast cancers.
Menopausal hormone therapy (MHT) is associated with an elevated risk of breast cancer in postmenopausal women. To identify genetic loci that modify breast cancer risk related to MHT use in postmenopausal women, we conducted a two-stage genome-wide association study (GWAS) with replication. In stage I, we performed a case-only GWAS in 731 invasive breast cancer cases from the German case-control study Mammary Carcinoma Risk Factor Investigation (MARIE). The 1,200 single nucleotide polymorphisms (SNPs) showing the lowest P values for interaction with current MHT use (within 6 months prior to breast cancer diagnosis), were carried forward to stage II, involving pooled case-control analyses including additional MARIE subjects (1,375 cases, 1,974 controls) as well as 795 cases and 764 controls of a Swedish case-control study. A joint P value was calculated for a combined analysis of stages I and II. Replication of the most significant interaction of the combined stage I and II was performed using 5,795 cases and 5,390 controls from nine studies of the Breast Cancer Association Consortium (BCAC). The combined stage I and II yielded five SNPs on chromosomes 2, 7, and 18 with joint P values <6 × 10−6 for effect modification of current MHT use. The most significant interaction was observed for rs6707272 (P = 3 × 10−7) on chromosome 2 but was not replicated in the BCAC studies (P = 0.21). The potentially modifying SNPs are in strong linkage disequilibrium with SNPs in TRIP12 and DNER on chromosome 2 and SETBP1 on chromosome 18, previously linked to carcinogenesis. However, none of the interaction effects reached genome-wide significance. The inability to replicate the top SNP × MHT interaction may be due to limited power of the replication phase. Our study, however, suggests that there are unlikely to be SNPs that interact strongly enough with MHT use to be clinically significant in European women.
Postmenopausal breast cancer risk; Menopausal hormone therapy; Polymorphisms; Gene-environment interaction; Genome-wide association study; Case-only study
Breast cancer is the most common cancer in women in developed countries and has a well-established genetic component. Germline mutations in a network of genes encoding BRCA1, BRCA2, and their interacting partners confer hereditary susceptibility to breast cancer. Abraxas directly interacts with the BRCA1 BRCT (BRCA1 carboxyl-terminal) repeats and contributes to BRCA1-dependent DNA damage responses, making Abraxas a candidate for yet unexplained disease susceptibility. Here, we have screened 125 Northern Finnish breast cancer families for coding region and splice-site Abraxas mutations and genotyped three tagging single-nucleotide polymorphisms within the gene from 991 unselected breast cancer cases and 868 female controls for common cancer-associated variants. A novel heterozygous alteration, c.1082G>A (Arg361Gln), that results in abrogated nuclear localization and DNA response activities was identified in three breast cancer families and in one additional familial case from an unselected breast cancer cohort, but not in healthy controls (P = 0.002). On the basis of its exclusive occurrence in familial cancers, disease cosegregation, evolutionary conservation, and disruption of critical BRCA1 functions, the recurrent Abraxas c.1082G>A mutation connects to cancer predisposition. These findings contribute to the concept of a BRCA-centered tumor suppressor network and provide the identity of Abraxas as a new breast cancer susceptibility gene.
We tested the hypotheses that CHEK2*1100delC heterozygosity is associated with increased risk of early death, breast cancer–specific death, and risk of a second breast cancer in women with a first breast cancer.
Patients and Methods
From 22 studies participating in the Breast Cancer Association Consortium, 25,571 white women with invasive breast cancer were genotyped for CHEK2*1100delC and observed for up to 20 years (median, 6.6 years). We examined risk of early death and breast cancer–specific death by estrogen receptor status and risk of a second breast cancer after a first breast cancer in prospective studies.
CHEK2*1100delC heterozygosity was found in 459 patients (1.8%). In women with estrogen receptor–positive breast cancer, multifactorially adjusted hazard ratios for heterozygotes versus noncarriers were 1.43 (95% CI, 1.12 to 1.82; log-rank P = .004) for early death and 1.63 (95% CI, 1.24 to 2.15; log-rank P < .001) for breast cancer–specific death. In all women, hazard ratio for a second breast cancer was 2.77 (95% CI, 2.00 to 3.83; log-rank P < .001) increasing to 3.52 (95% CI, 2.35 to 5.27; log-rank P < .001) in women with estrogen receptor–positive first breast cancer only.
Among women with estrogen receptor–positive breast cancer, CHEK2*1100delC heterozygosity was associated with a 1.4-fold risk of early death, a 1.6-fold risk of breast cancer–specific death, and a 3.5-fold risk of a second breast cancer. This is one of the few examples of a genetic factor that influences long-term prognosis being documented in an extensive series of women with breast cancer.
γ-secretase is a large ubiquitously expressed protease complex composed of four core subunits: presenilin, Aph1, PEN-2, and nicastrin. The function of γ-secretase in the cells is to proteolytically cleave various proteins within their transmembrane domains. Presenilin and Aph1 occur as alternative variants belonging to mutually exclusive γ-secretase complexes and providing the complexes with heterogeneous biochemical and physiological properties. γ-secretase is proposed to have a role in the development and progression of cancer and γ-secretase inhibitors are intensively studied for their probable anti-tumor effects in various types of cancer models. Here, we for the first time determined mRNA expression levels of presenilin-1, presenilin-2, Aph1a, Aph1b, PEN-2, and nicastrin in a set of breast cancer tissue samples (N = 55) by quantitative real-time PCR in order to clarify the clinical significance of the expression of different γ-secretase complex components in breast cancer. We found a high positive correlation between the subunit expression levels implying a common regulation of transcription. Our univariate Kaplan-Meier survival analyses established low expression level of γ-secretase complex as a risk factor for breast cancer specific mortality. The tumors expressing low levels of γ-secretase complex were characterized by high histopathological tumor grade, low or no expression of estrogen and progesterone receptors and consequently high probability to fall into the class of triple negative breast cancer tumors. These results may provide novel tools to further categorize breast cancer tumors, especially the highly aggressive and poorly treatable breast cancer type of triple negative cases, and suggest a significant role for γ-secretase in breast cancer.
TERT-locus single nucleotide polymorphisms (SNPs) and leucocyte telomere measures are reportedly associated with risks of multiple cancers. Using the iCOGs chip, we analysed ~480 TERT-locus SNPs in breast (n=103,991), ovarian (n=39,774) and BRCA1 mutation carrier (11,705) cancer cases and controls. 53,724 participants have leucocyte telomere measures. Most associations cluster into three independent peaks. Peak 1 SNP rs2736108 minor allele associates with longer telomeres (P=5.8×10−7), reduced estrogen receptor negative (ER-negative) (P=1.0×10−8) and BRCA1 mutation carrier (P=1.1×10−5) breast cancer risks, and altered promoter-assay signal. Peak 2 SNP rs7705526 minor allele associates with longer telomeres (P=2.3×10−14), increased low malignant potential ovarian cancer risk (P=1.3×10−15) and increased promoter activity. Peak 3 SNPs rs10069690 and rs2242652 minor alleles increase ER-negative (P=1.2×10−12) and BRCA1 mutation carrier (P=1.6×10−14) breast and invasive ovarian (P=1.3×10−11) cancer risks, but not via altered telomere length. The cancer-risk alleles of rs2242652 and rs10069690 respectively increase silencing and generate a truncated TERT splice-variant.
Our recent genome-wide association study identified a novel breast cancer susceptibility locus at 9q31.2 (rs865686).
To further investigate the rs865686–breast cancer association, we conducted a replication study within the Breast Cancer Association Consortium, which comprises 37 case–control studies (48,394 cases, 50,836 controls).
This replication study provides additional strong evidence of an inverse association between rs865686 and breast cancer risk [study-adjusted per G-allele OR, 0.90; 95% confidence interval (CI), 0.88; 0.91, P = 2.01 × 10–29] among women of European ancestry. There were ethnic differences in the estimated minor (G)-allele frequency among controls [0.09, 0.30, and 0.38 among, respectively, Asians, Eastern Europeans, and other Europeans; P for heterogeneity (Phet) = 1.3 × 10–143], but no evidence of ethnic differences in per allele OR (Phet = 0.43). rs865686 was associated with estrogen receptor–positive (ER+) disease (per G-allele OR, 0.89; 95% CI, 0.86–0.91; P = 3.13 × 10–22) but less strongly, if at all, with ER-negative (ER–) disease (OR, 0.98; 95% CI, 0.94–1.02; P = 0.26; Phet = 1.16 × 10–6), with no evidence of independent heterogeneity by progesterone receptor or HER2 status. The strength of the breast cancer association decreased with increasing age at diagnosis, with case-only analysis showing a trend in the number of copies of the G allele with increasing age at diagnosis (P for linear trend = 0.0095), but only among women with ER+ tumors.
This study is the first to show that rs865686 is a susceptibility marker for ER+ breast cancer.
The findings further support the view that genetic susceptibility varies according to tumor subtype.
Recent genome-wide association studies identified 11 single nucleotide polymorphisms (SNPs) associated with breast cancer (BC) risk. We investigated these and 62 other SNPs for their prognostic relevance. Confirmed BC risk SNPs rs17468277 (CASP8), rs1982073 (TGFB1), rs2981582 (FGFR2), rs13281615 (8q24), rs3817198 (LSP1), rs889312 (MAP3K1), rs3803662 (TOX3), rs13387042 (2q35), rs4973768 (SLC4A7), rs6504950 (COX11) and rs10941679 (5p12) were genotyped for 25 853 BC patients with the available follow-up; 62 other SNPs, which have been suggested as BC risk SNPs by a GWAS or as candidate SNPs from individual studies, were genotyped for replication purposes in subsets of these patients. Cox proportional hazard models were used to test the association of these SNPs with overall survival (OS) and BC-specific survival (BCS). For the confirmed loci, we performed an accessory analysis of publicly available gene expression data and the prognosis in a different patient group. One of the 11 SNPs, rs3803662 (TOX3) and none of the 62 candidate/GWAS SNPs were associated with OS and/or BCS at P<0.01. The genotypic-specific survival for rs3803662 suggested a recessive mode of action [hazard ratio (HR) of rare homozygous carriers=1.21; 95% CI: 1.09–1.35, P=0.0002 and HR=1.29; 95% CI: 1.12–1.47, P=0.0003 for OS and BCS, respectively]. This association was seen similarly in all analyzed tumor subgroups defined by nodal status, tumor size, grade and estrogen receptor. Breast tumor expression of these genes was not associated with prognosis. With the exception of rs3803662 (TOX3), there was no evidence that any of the SNPs associated with BC susceptibility were associated with the BC survival. Survival may be influenced by a distinct set of germline variants from those influencing susceptibility.
A recent two-stage genome-wide association study (GWAS) identified five novel breast cancer susceptibility loci on chromosomes 9, 10 and 11. To provide more reliable estimates of the relative risk associated with these loci and investigate possible heterogeneity by subtype of breast cancer, we genotyped the variants rs2380205, rs1011970, rs704010, rs614367, rs10995190 in 39 studies from the Breast Cancer Association Consortium (BCAC), involving 49,608 cases and 48,772 controls of predominantly European ancestry. Four of the variants showed clear evidence of association (P ≤ 3 × 10−9) and weak evidence was observed for rs2380205 (P = 0.06). The strongest evidence was obtained for rs614367, located on 11q13 (per-allele odds ratio 1.21, P = 4 × 10−39). The association for rs614367 was specific to estrogen receptor (ER)-positive disease and strongest for ER plus progesterone receptor (PR)-positive breast cancer, whereas the associations for the other three loci did not differ by tumor subtype.
breast cancer susceptibility; polymorphisms; genome wide association; risk factors; hormone receptor status; 11q13
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−6). 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.
UGT1A6; polymorphism; breast cancer risk; validation; metabolism
Breast cancer is the most common cancer among women. To date, 22 common breast cancer susceptibility loci have been identified accounting for ~ 8% of the heritability of the disease. We followed up 72 promising associations from two independent Genome Wide Association Studies (GWAS) in ~70,000 cases and ~68,000 controls from 41 case-control studies and nine breast cancer GWAS. We identified three new breast cancer risk loci on 12p11 (rs10771399; P=2.7 × 10−35), 12q24 (rs1292011; P=4.3×10−19) and 21q21 (rs2823093; P=1.1×10−12). SNP rs10771399 was associated with similar relative risks for both estrogen receptor (ER)-negative and ER-positive breast cancer, whereas the other two loci were associated only with ER-positive disease. Two of the loci lie in regions that contain strong plausible candidate genes: PTHLH (12p11) plays a crucial role in mammary gland development and the establishment of bone metastasis in breast cancer, while NRIP1 (21q21) encodes an ER co-factor and has a role in the regulation of breast cancer cell growth.
Common genetic variants contribute to the observed variation in breast cancer risk for BRCA2 mutation carriers; those known to date have all been found through population-based genome-wide association studies (GWAS). To comprehensively identify breast cancer risk modifying loci for BRCA2 mutation carriers, we conducted a deep replication of an ongoing GWAS discovery study. Using the ranked P-values of the breast cancer associations with the imputed genotype of 1.4 M SNPs, 19,029 SNPs were selected and designed for inclusion on a custom Illumina array that included a total of 211,155 SNPs as part of a multi-consortial project. DNA samples from 3,881 breast cancer affected and 4,330 unaffected BRCA2 mutation carriers from 47 studies belonging to the Consortium of Investigators of Modifiers of BRCA1/2 were genotyped and available for analysis. We replicated previously reported breast cancer susceptibility alleles in these BRCA2 mutation carriers and for several regions (including FGFR2, MAP3K1, CDKN2A/B, and PTHLH) identified SNPs that have stronger evidence of association than those previously published. We also identified a novel susceptibility allele at 6p24 that was inversely associated with risk in BRCA2 mutation carriers (rs9348512; per allele HR = 0.85, 95% CI 0.80–0.90, P = 3.9×10−8). This SNP was not associated with breast cancer risk either in the general population or in BRCA1 mutation carriers. The locus lies within a region containing TFAP2A, which encodes a transcriptional activation protein that interacts with several tumor suppressor genes. This report identifies the first breast cancer risk locus specific to a BRCA2 mutation background. This comprehensive update of novel and previously reported breast cancer susceptibility loci contributes to the establishment of a panel of SNPs that modify breast cancer risk in BRCA2 mutation carriers. This panel may have clinical utility for women with BRCA2 mutations weighing options for medical prevention of breast cancer.
Women who carry BRCA2 mutations have an increased risk of breast cancer that varies widely. To identify common genetic variants that modify the breast cancer risk associated with BRCA2 mutations, we have built upon our previous work in which we examined genetic variants across the genome in relation to breast cancer risk among BRCA2 mutation carriers. Using a custom genotyping platform with 211,155 genetic variants known as single nucleotide polymorphisms (SNPs), we genotyped 3,881 women who had breast cancer and 4,330 women without breast cancer, which represents the largest possible, international collection of BRCA2 mutation carriers. We identified that a SNP located at 6p24 in the genome was associated with lower risk of breast cancer. Importantly, this SNP was not associated with breast cancer in BRCA1 mutation carriers or in a general population of women, indicating that the breast cancer association with this SNP might be specific to BRCA2 mutation carriers. Combining this BRCA2-specific SNP with 13 other breast cancer risk SNPs also known to modify risk in BRCA2 mutation carriers, we were able to derive a risk prediction model that could be useful in helping women with BRCA2 mutations weigh their risk-reduction strategy options.
Various common genetic susceptibility loci have been identified for breast cancer; however, it is unclear how they combine with lifestyle/environmental risk factors to influence risk. We undertook an international collaborative study to assess gene-environment interaction for risk of breast cancer. Data from 24 studies of the Breast Cancer Association Consortium were pooled. Using up to 34,793 invasive breast cancers and 41,099 controls, we examined whether the relative risks associated with 23 single nucleotide polymorphisms were modified by 10 established environmental risk factors (age at menarche, parity, breastfeeding, body mass index, height, oral contraceptive use, menopausal hormone therapy use, alcohol consumption, cigarette smoking, physical activity) in women of European ancestry. We used logistic regression models stratified by study and adjusted for age and performed likelihood ratio tests to assess gene–environment interactions. All statistical tests were two-sided. We replicated previously reported potential interactions between LSP1-rs3817198 and parity (Pinteraction = 2.4×10−6) and between CASP8-rs17468277 and alcohol consumption (Pinteraction = 3.1×10−4). Overall, the per-allele odds ratio (95% confidence interval) for LSP1-rs3817198 was 1.08 (1.01–1.16) in nulliparous women and ranged from 1.03 (0.96–1.10) in parous women with one birth to 1.26 (1.16–1.37) in women with at least four births. For CASP8-rs17468277, the per-allele OR was 0.91 (0.85–0.98) in those with an alcohol intake of <20 g/day and 1.45 (1.14–1.85) in those who drank ≥20 g/day. Additionally, interaction was found between 1p11.2-rs11249433 and ever being parous (Pinteraction = 5.3×10−5), with a per-allele OR of 1.14 (1.11–1.17) in parous women and 0.98 (0.92–1.05) in nulliparous women. These data provide first strong evidence that the risk of breast cancer associated with some common genetic variants may vary with environmental risk factors.
Breast cancer involves combined effects of numerous genetic, environmental, and behavioral risk factors that are unique to each individual. High risk genes, such as BRCA1 and BRCA2, account for only a small proportion of disease occurrence. Recent genome-wide research has identified more than 20 common genetic variants, which individually alter breast cancer risk very moderately. We undertook an international collaborative study to determine whether the effect of these genetic variants vary with environmental factors, such as parity, body mass index (BMI), height, oral contraceptive use, menopausal hormone therapy use, alcohol consumption, cigarette smoking, and physical activity, which are known to affect risk of developing breast cancer. Using pooled data from 24 studies of the Breast Cancer Association Consortium (BCAC), we provide first convincing evidence that the breast cancer risk associated with a genetic variant in LSP1 differs with the number of births and that the risk associated with a CASP8 variant is altered by high alcohol consumption. The effect of an additional genetic variant might also be modified by reproductive factors. This knowledge will stimulate new research towards a better understanding of breast cancer development.
A genome-wide association study (GWAS) identified single-nucleotide polymorphisms (SNPs) at 1p11.2 and 14q24.1 (RAD51L1) as breast cancer susceptibility loci. The initial GWAS suggested stronger effects for both loci for estrogen receptor (ER)-positive tumors. Using data from the Breast Cancer Association Consortium (BCAC), we sought to determine whether risks differ by ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), grade, node status, tumor size, and ductal or lobular morphology. We genotyped rs11249433 at 1p.11.2, and two highly correlated SNPs rs999737 and rs10483813 (r2= 0.98) at 14q24.1 (RAD51L1), for up to 46 036 invasive breast cancer cases and 46 930 controls from 39 studies. Analyses by tumor characteristics focused on subjects reporting to be white women of European ancestry and were based on 25 458 cases, of which 87% had ER data. The SNP at 1p11.2 showed significantly stronger associations with ER-positive tumors [per-allele odds ratio (OR) for ER-positive tumors was 1.13, 95% CI = 1.10–1.16 and, for ER-negative tumors, OR was 1.03, 95% CI = 0.98–1.07, case-only P-heterogeneity = 7.6 × 10−5]. The association with ER-positive tumors was stronger for tumors of lower grade (case-only P= 6.7 × 10−3) and lobular histology (case-only P= 0.01). SNPs at 14q24.1 were associated with risk for most tumor subtypes evaluated, including triple-negative breast cancers, which has not been described previously. Our results underscore the need for large pooling efforts with tumor pathology data to help refine risk estimates for SNP associations with susceptibility to different subtypes of breast cancer.
Estrogen receptor (ER)-negative breast cancer shows a higher incidence in women of African ancestry compared to women of European ancestry. In search of common risk alleles for ER-negative breast cancer, we combined genome-wide association study (GWAS) data from women of African ancestry (1,004 ER-negative cases and 2,745 controls) and European ancestry (1,718 ER-negative cases and 3,670 controls), with replication testing conducted in an additional 2,292 ER-negative cases and 16,901 controls of European ancestry. We identified a common risk variant for ER-negative breast cancer at the TERT-CLPTM1L locus on chromosome 5p15 (rs10069690: per-allele odds ratio (OR) = 1.18 per allele, P = 1.0 × 10−10). The variant was also significantly associated with triple-negative (ER-negative, progesterone receptor (PR)-negative and human epidermal growth factor-2 (HER2)-negative) breast cancer (OR = 1.25, P = 1.1 × 10−9), particularly in younger women (<50 years of age) (OR = 1.48, P = 1.9 × 10−9). Our results identify a genetic locus associated with estrogen receptor negative breast cancer subtypes in multiple populations.
The histone demethylase GASC1 (JMJD2C) is an epigenetic factor suspected of involvement in development of different cancers, including breast cancer. It is thought to be overexpressed in the more aggressive breast cancer types based on mRNA expression studies on cell lines and meta analysis of human breast cancer sets. This study aimed to evaluate the prognostic and predictive value of GASC1 for women with invasive breast cancer.
All the 355 cases were selected from a cohort enrolled in the Kuopio Breast Cancer Project between April 1990 and December 1995. The expression of GASC1 was studied by immunohistochemistry (IHC) on tissue microarrays. Additionally relative GASC1 mRNA expression was measured from available 57 cases.
In our material, 56% of the cases were GASC1 negative and 44% positive in IHC staining. Women with GASC1 negative tumors had two years shorter breast cancer specific survival and time to relapse than the women with GASC1 positive tumors (p=0.017 and p=0.034 respectively). The majority of GASC1 negative tumors were ductal cases (72%) of higher histological grade (84% of grade II and III altogether). When we evaluated estrogen receptor negative and progesterone receptor negative cases separately, there was 2 times more GASC1 negative than GASC1 positive tumors in each group (chi2, p= 0.033 and 0.001 respectively). In the HER2 positive cases, there was 3 times more GASC1 negative cases than GASC1 positives (chi2, p= 0.029). Patients treated with radiotherapy (n=206) and hormonal treatment (n=62) had better breast cancer specific survival, when they were GASC1 positive (Cox regression: HR=0.49, p=0.007 and HR=0.33, p=0.015, respectively). The expression of GASC1 mRNA was in agreement with the protein analysis.
This study indicates that the GASC1 is both a prognostic and a predictive factor for women with invasive breast cancer. GASC1 negativity is associated with tumors of more aggressive histopathological types (ductal type, grade II and III, ER negative, PR negative). Patients with GASC1 positive tumors have better breast cancer specific survival and respond better to radiotherapy and hormonal treatment.
Epigenetics; GASC1; Breast cancer; Survival; Tissue microarrays
The 19p13.1 breast cancer susceptibility locus is a modifier of breast cancer risk in BRCA1 mutation carriers and is also associated with risk of ovarian cancer. Here we investigated 19p13.1 variation and risk of breast cancer subtypes, defined by estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor-2 (HER2) status, using 48,869 breast cancer cases and 49,787 controls from the Breast Cancer Association Consortium (BCAC). Variants from 19p13.1 were not associated with breast cancer overall or with ER-positive breast cancer but were significantly associated with ER-negative breast cancer risk [rs8170 Odds Ratio (OR)=1.10, 95% Confidence Interval (CI) 1.05 – 1.15, p=3.49 × 10-5] and triple negative (TN) (ER, PR and HER2 negative) breast cancer [rs8170 OR=1.22, 95% CI 1.13 – 1.31, p=2.22 × 10-7]. However, rs8170 was no longer associated with ER-negative breast cancer risk when TN cases were excluded [OR=0.98, 95% CI 0.89 – 1.07, p=0.62]. In addition, a combined analysis of TN cases from BCAC and the Triple Negative Breast Cancer Consortium (TNBCC) (n=3,566) identified a genome-wide significant association between rs8170 and TN breast cancer risk [OR=1.25, 95% CI 1.18 – 1.33, p=3.31 × 10-13]. Thus, 19p13.1 is the first triple negative-specific breast cancer risk locus and the first locus specific to a histological subtype defined by ER, PR, and HER2 to be identified. These findings provide convincing evidence that genetic susceptibility to breast cancer varies by tumor subtype and that triple negative tumors and other subtypes likely arise through distinct etiologic pathways.
genetic susceptibility; association study; subtype; neoplasms; common variant
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.
The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER−) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26–1.48), p = 7.6×10−14 in Asians and 1.09 (95% CI 1.07–1.11), p = 6.8×10−18 in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19–1.41), p = 1.2×10−9 in Asians and 1.12 (95% CI 1.08–1.17), p = 3.8×10−9 in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER− than ER+ tumours in Europeans [OR (ER−) = 1.20 (95% CI 1.15–1.25), p = 1.8×10−17 versus OR (ER+) = 1.07 (95% CI 1.04–1.1), p = 1.3×10−7, pheterogeneity = 5.1×10−6]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER− tumours.
Recently, a locus on chromosome 6q22.33 (rs2180341) was reported to be associated with increased breast cancer risk in the Ashkenazi Jewish (AJ) population, and this association was also observed in populations of non-AJ European ancestry. In the present study, we performed a large replication analysis of rs2180341 using data from 31,428 invasive breast cancer cases and 34,700 controls collected from 25 studies in the Breast Cancer Association Consortium (BCAC). In addition, we evaluated whether rs2180341 modifies breast cancer risk in 3,361 BRCA1 and 2,020 BRCA2 carriers from 11 centers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Based on the BCAC data from women of European ancestry, we found evidence for a weak association with breast cancer risk for rs2180341 (per-allele odds ratio (OR) = 1.03, 95% CI 1.00–1.06, p = 0.023). There was evidence for heterogeneity in the ORs among studies (I2 = 49.3%; p = <0.004). In CIMBA, we observed an inverse association with the minor allele of rs2180341 and breast cancer risk in BRCA1 mutation carriers (per-allele OR = 0.89, 95%CI 0.80–1.00, p = 0.048), indicating a potential protective effect of this allele. These data suggest that that 6q22.33 confers a weak effect on breast cancer risk.