- Research article
- Open Access
Genetic variation of ESR1 and its co-activator PPARGC1B is synergistic in augmenting the risk of estrogen receptor-positive breast cancer
Breast Cancer Research volume 13, Article number: R10 (2011)
Given the role of estrogen in breast carcinogenesis and the modification of estrogen receptor (ER) activity by its biochemical cofactors, we hypothesize that genetic variation within ER cofactor genes alters cellular response to estrogen exposure and consequently modifies the risk for ER-positive breast cancer.
We genotyped 790 tagging SNPs within 60 ER cofactor genes in 1,257 cases and 1,464 controls from Sweden and in 2,215 cases and 1,265 controls from Finland, and tested their associations with either ER-positive or ER-negative breast cancer.
Seven SNPs showed consistent association with ER-positive breast cancer in the two independent samples, and six of them were located within PPARGC1B, encoding an ER co-activator, with the strongest association at rs741581 (odds ratio = 1.41, P = 4.84 × 10-5) that survived Bonferroni correction for multiple testing in the combined ER-positive breast cancer sample (Pcorrected = 0.03). Moreover, we also observed significant synergistic interaction (P interaction = 0.008) between the genetic polymorphisms within PPARGC1B and ESR1 in ER-positive breast cancer. By contrast, no consistent association was observed in ER-negative breast cancer. Furthermore, we found that administration of estrogen in the MCF-7 cell line induced PPARGC1B expression and enhanced occupancies of ER and RNA polymerase II within the region of SNP association, suggesting the upregulation of PPARGC1B expression by ESR1 activation.
Our study revealed that DNA polymorphisms of PPARGC1B, coding a bona fide ER co-activator, are associated with ER-positive breast cancer risk. The feed-forward transcriptional regulatory loop between PPARGC1B and ESR1 further augments their protein interaction, which provides a plausible mechanistic explanation for the synergistic genetic interaction between PPARGC1B and ESR1 in ER-positive breast cancer. Our study also highlights that biochemically and genomically informed candidate gene studies can enhance the discovery of interactive disease susceptibility genes.
It is known that the risk of breast cancer is related to lifetime exposure to estrogen [1, 2]. Estrogen stimulates cell proliferation and increases the frequency of spontaneous mutations, leading to a malignant phenotype . Breast cells respond to estrogen via estrogen receptors (ERs) through a defined biochemical process: upon ligand binding, ERs undergo a conformational change that facilitates receptor dimerization, DNA binding, recruitment of ER cofactors, and modulation of target gene expression [4–6].
Endocrine therapy provides strong evidence that attenuation of ER (ESR1) activity can reduce breast cancer risk , and women with ER-positive tumor would be the most likely to benefit from these treatments [7, 8]. The genetic studies of ESR1, however, have had contradictory results. Only recently, through a very large genetic association study, has there been demonstrated a small but significant association of polymorphisms within ESR1 with the risk of breast cancer [9–11]. Two plausible explanations for the inconsistent results might be due to the small sample sizes and thus limited statistical power of these studies, or that the risk was not evaluated by stratifying breast cancer patients based on tumor ER status. However, there is at least one further possibility: ER cofactors can either enhance transcriptional activity of ER as co-activators or inhibit the activity as co-repressors. The genetic variants within ER cofactors have not been systematically investigated in term of association with breast cancer risk, although some coding variants within individual genes, such as NCOA3 and CCND1, have been investigated [12–15].
Given the modification of ER activity by its cofactors through their physical and functional interactions , the cofactor proteins that bind to ER may be as important as the receptor itself in mediating transcriptional response to estrogen exposure . We therefore hypothesized that genetic variation within ER cofactor genes may alter cellular response to estrogen exposure and consequently, alone or by interacting with genetic variations within ESR1, modify breast cancer risk in an ER status-dependent fashion. To assess this hypothesis, we investigated the association of common genetic variation, using a tagging SNP approach, within 60 cofactor genes in two large case-control samples of breast cancer from Sweden and Finland, and investigated their interaction with genetic variation within ESR1 in terms of influencing the risk of hormone-driven breast cancer.
Materials and methods
The Swedish sample was from a population-based case-control study that has been described in detail previously . Briefly, 1,322 cases were Swedish-born women diagnosed with incident primary invasive breast cancer between October 1993 and March 1995 who contributed blood samples. All cases were postmenopausal and between 50 and 74 years of age at diagnosis. All the cases were identified through the six regional cancer registries in Sweden. The controls (n = 1,524) were randomly selected from the Swedish Registry of Total Population with no previous breast cancer and were frequency-matched for age with the cases. Questionnaires were used to collect risk factor information.
The Finnish sample was from a hospital-based case-control study in which the cases consisted of two series of unselected breast cancer patients and additional familial patients diagnosed at the Helsinki University Central Hospital. The first set of cases were 884 patients collected in 1997/1998 and 2000, covering 79% of all newly diagnosed breast cancer cases during those periods [19, 20]. The second set of cases, consisting of 986 newly diagnosed breast cancer patients, were collected during 2001 to 2004 and covered 87% of all such patients during that period . An additional 538 familial breast cancer cases were also collected at the same hospital, as previously described [22, 23]. Women with a prior diagnosis of breast cancer in situ were excluded, leaving 2,215 invasive breast cancer cases for analysis. Healthy female population controls (n = 1,287) were collected from the same geographical regions of Finland as the cases.
Information on reproductive and hormonal risk factors was available for the Swedish sample and showed expected association patterns with breast cancer [24–26]. Such information was not available for the Finnish controls.
Hormone receptor status information was retrieved from medical records of all participating cases and was available for both the Swedish and Finnish cases.
Approval for the study was obtained from the Institutional Review Boards in Sweden, Finland and the National University of Singapore. All subjects provided written informed consent.
DNA was extracted from 4 ml whole blood using the QIAamp DNA Blood Maxi Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions.
Candidate gene and tagging SNP selection
In the present study, the keywords 'ER cofactor', 'ER coactivator' and 'ER corepressor' were used in a literature search to identify ER cofactor genes. Boolean searching ('AND' 'OR') was used to narrow or broaden the search in PubMed. Using this method, 60 ER cofactor genes were identified as candidate genes. Tagging SNPs within the 60 candidate genes were selected based on the HapMap CEU data (Rel #22/phase II Apr07, on NCBI B36 assembly, dbSNP b126) . In brief, for each gene, all common SNPs with a minor allele frequency >0.05 within the gene and 5 kb surrounding region were first identified from the HapMap database . Tagging SNPs were then selected in Haploview version 4.1  using a pair-wise SNP tagging approach with r2 > 0.8 used as the criterion for selection. A total of 806 tagging SNPs were selected within the 60 ER co-factor genes.
Illumina's GoldenGate assay was used for genotyping SNPs, following the manufacturers' instructions (Illumina, San Diego, CA, USA). In brief, all 806 tagging SNPs were subjected to genotyping assay design, out of which 790 SNPs were successfully designed and subjected to genotyping analysis. DNA samples were randomly assigned to the plates carrying positive and negative controls, and all genotyping results were generated and checked by laboratory staff unaware of the case-control status. SNPs with a call rate <96% (81 SNPs failed in the Swedish sample and 42 SNPs failed in the Finnish sample) and minor allele frequency <1% (18 SNPs in the Swedish sample and 40 SNPs in the Finnish sample) were excluded from further analysis. Deviation of genotype frequencies from those expected under Hardy-Weinberg Equilibrium were assessed in the control subjects. SNPs with Hardy-Weinberg Equilibrium P < 7.4 × 10-5 (0.05/675) were excluded (6 SNPs failed in the Swedish sample and 15 SNPs failed in the Finnish sample). In total, 685 SNPs from the Swedish sample and 693 SNPs from the Finnish sample were used for statistical analysis, and 675 shared SNPs between the Swedish and Finnish samples were used for analysis in the combined sample.
Genotyping was duplicated in 2% of samples (in both Swedish and Finnish samples) and there was concordance in >99% of the duplicated samples, suggesting high genotyping accuracy. With r2 > 0.8, the average coverage of common variation (minor allele frequency >5%) within the 60 candidate genes was 91%. Out of these, 51 genes had coverage over 80% (Additional file 1 Table S1).
Reverse transcriptase-quantitative PCR analysis
MCF-7 cells were cultured in DMEM (Invitrogen, Carlsbad, CA, USA) medium with 10% FBS (Invitrogen). Prior to hormone treatment, cells were maintained in phenol-red free DMEM F-12 containing 5% charcoal stripped serum for 72 hours for hormone depletion. Cells were treated with 10 nM 17β-estradiol (Sigma-Aldrich, St. Louis, MO, USA) for a period of 0 or 3 hours. Cells were harvested and total RNA and reverse transcriptase-quantitative PCR analysis was carried out as described previously . Dimethylsulfoxide (Sigma-Aldrich, St. Louis, MO, USA)/vehicle-treated cells were used as controls for the same time course. Real-time PCR analysis was performed in the ABI Prism 7700 sequence detection system using SYBR Green from ABI (Applied Biosystems, Foster City, CA,USA).
Primers were designed using the online Primer 3 program . All experiments were repeated at least twice. Two sets of primers were used for identifying different isoforms of PPARGC1B. The oligonucleotide sequences were as follows: PPARGC1B_1 isoform (NM_001172699.1) forward 5'-GAAGAGGAAGAAGGGGAGGA-3' and reverse 5'-CTCTGGTAGGGGCAGTGGT-3'; and PPARGC1B_2 isoform (NM_133263.3) forward 5'-CCTGAAGATGACGTGGGTCT-3' and reverse 5'-CCTTCCTTCTGGGTGTCAGA-3'. β-Actin specific primers (forward 5'-TCCCTGGAGAAGAGCTACGA-3' and reverse 5'-AGGAAGGAAGGCTGGAAGAG-3') were used as an internal control to normalize the amounts of reverse transcribed product used in the PCR reaction. Threshold cycle (Ct) values obtained for PPARGC1B isoforms were normalized to β-actin Ct values. The normalized Ct (ΔCt) values were then used to calculate the difference (ΔΔCt) between estradiol-treated and dimethylsulfoxide-treated samples. The fold change of PPARGC1B was calculated as 2-ΔΔCt.
To measure the magnitude of association between SNPs and breast cancer risk, per-allele odds ratios (ORs) (assuming a log-additive model) and 95% confidence intervals were estimated using logistic regression. As the controls were younger than cases in the Finnish samples, age at diagnosis/enrollment (as a continuous variable) was included in the regression models in the Finnish analysis for OR adjustment. The Cochran-Armitage trend test was used to calculate P values in the Swedish and Finnish sample sets, separately in subtypes, and in cases overall. Inverse variance weighting was used in a meta-analysis for two independent datasets. The individual OR was obtained from age-unadjusted analysis in the Swedish sample and age-adjusted analysis in the Finnish sample. To evaluate differences in ORs between studies, a test of homogeneity was carried out for each individual SNP analysis (data not shown).
To determine the model of inheritance, associations between SNPs within the PPARGC1B gene and ER-positive breast cancer risk were estimated by assuming dominant, recessive and additive models in the two sample sets. We then performed these analyses with meta-analysis using inverse variance weighting approach. Individual ORs from two independent studies followed-up age-unadjusted analysis in the Swedish sample and age-adjusted analysis in the Finnish sample.
Forward stepwise logistic regression was used to explore whether the associations at the six SNPs were independent of each other. The selection criterion was P < 0.2. The analysis was performed in ER-positive breast cancer risk in the two sample sets separately as well as in the combined ER-positive sample dataset. To account for different minor allele frequencies in the two populations, a binary indicator variable for study was included in the regression models as well as age in the combined data regression analysis.
Pair-wise interaction analysis was performed under a dominant mode of inheritance using logistic regression and likelihood ratio tests. To maximize the statistical power, we pooled sample sets from the Swedish and Finnish data. Age and study were included in the model as covariables. The full model included an interaction term between the two interacting variables for the risk of breast cancer. In this multivariate logistic regression analysis, each coefficient provided an estimate of the log OR whilst adjusting for all other variables included in the model. Likelihood ratio tests, comparing models with and without the interaction term, were used to generate P values.
All analyses were performed using STATA version 8.0 (StataCorp, College station, TX, USA). Linkage disequilibrium (LD) calculation was performed in Haploview version 4.1 . All statistical tests were two-sided.
Two independent case-control samples of breast cancer from Sweden and Finland were investigated in the present study, whose characteristics are summarized in Table 1. The cases and controls of the Swedish sample were frequency-matched on age, whereas the Finnish controls were younger than the Finnish cases (P < 0.0001). In the Swedish sample, there were significant differences between the cases and controls in terms of age at first birth (P = 0.0002), age at menopause (P = 0.0001), hormone replacement treatment use (P = 0.017), and parity (P = 0.0001), which is consistent with the well-established role of these reproductive factors in breast cancer development. The reproductive factor information was not available for the Finnish controls. In both the Swedish and Finnish cases, there were similar percentages of ER-positive (81.9% vs. 80.9%) and ER-negative (18.1% vs. 19.1%) cases.
SNP association analysis
First, single SNP association analyses were performed using trend tests in the Swedish and Finnish samples separately by stratifying the cases into ER-positive and ER-negative groups, with 685 SNPs being tested in the Swedish sample and 693 SNPs being tested in the Finnish sample. 48 SNPs (7.00%) in the Swedish sample and 50 SNPs (7.28%) in the Finnish sample showed association with ER-positive breast cancer risk with nominal P < 0.05. Seven SNPs showed consistent association between the two independent samples (Additional file 1 Table S2), and six of them were located within the PPARGC1B gene. In contrast, 21 and 50 SNPs showed association with ER-negative breast cancer with nominal P < 0.05 in the Swedish and Finnish samples, respectively, but no SNPs showed consistent associations between the two independent samples.
We then analyzed SNP associations in the combined Swedish and Finnish samples. In general, SNPs showed stronger evidence of association with ER-positive breast cancer than ER-negative breast cancer (Table 2; see also Additional file 1 Table S3). The most significant association was identified at rs741581 within the second intron of PPARGC1B (OR = 1.41, P = 4.84 × 10-5) in ER-positive breast cancer, which survived the Bonferroni correction for multiple testing (Pcorrected = 0.03). rs741581 was one of the seven SNPs that showed consistent associations between the Swedish and Finnish samples.
We also evaluated the SNP association with overall breast cancer risk and found 55 SNPs (8.03%) from the Swedish samples and 61 SNPs (8.80%) from the Finnish samples to show association with overall breast cancer risk with nominal P < 0.05. Only two SNPs, however, showed consistent association between the two independent samples (Additional file 1 Table S2), and none of the associations survived Bonferroni correction for multiple testing in the combined samples (smallest Pcorrected = 0.198).
Genotype association analysis of PPARGC1Bin ER-positive breast cancer
PPARGC1B is located on 5q33.1 and encodes for peroxisome proliferative activated receptor gamma coactivator beta protein (PGC-1β), a bona fide co-activator of ERα. In the present study, 40 tagging SNPs within PPARGC1B were successfully genotyped in both the Swedish and Finnish samples, which could capture 80% of common variants (131 out of 162 SNPs) within PPARGC1B with a minimal r2 value of 0.8 (mean r2 value = 0.95, according to HapMap CEU data).
To have a better understanding of the association within PPARGC1B, we performed genotype-based association analysis by assuming dominant, recessive and additive model of inheritance. We found that the top three SNPs yielding the most significant association evidence in the dominant model compared with other models (Additional file 1 Table S5). Under the dominant model, the same six SNPs (as for the trend tests) of the 40 SNPs within PPARGC1B showed consistent association with ER-positive breast cancer between the Swedish and Finnish samples (Table 3). The strength of the association (ORs) at the six SNPs was stronger in ER-positive breast cancer than in overall or ER-negative cancers, with the strongest association identified at rs741581 (P = 1.9 × 10-2 in the Swedish samples, P = 6.1 × 10-5 in the Finnish samples, and P = 1.8 × 10-5 in the combined samples).
The six SNPs showing consistent association with ER-positive breast cancer were located within two regions of high LD (Figure 1B), suggesting that the associations at those SNPs may not be completely independent. We therefore performed a forward stepwise logistic regression (cut-off P = 0.20) and revealed two independent associations with ER-positive breast cancer at rs741581 (P = 0.031) and rs6895698 (P = 0.014) in the combined sample. Under the dominant model, we found that rs741581, rs6895698, age and study sample were four independent variables associated with ER-positive breast cancer risk. Similarly, the same stepwise analysis of ER-positive breast cancer in the two individual samples also revealed two independent associations at rs741581 (P = 0.172) and rs2340621 (P = 0.036) in the Swedish sample and at rs741581 (P = 0.023) and rs6895698 (P = 0.053) in the Finnish sample. Notably, rs6895698 and rs2340621 lie within the same LD block and are highly correlated (r2 = 0.72, according to HapMap CEU data).
Association of ESR1variation with ER-positive breast cancer
Our previous study suggested an association between ESR1 polymorphisms and breast cancer risk [9, 32] in the same Swedish sample. The association was within a region flanked by rs988328 to rs3020318 and was manifested by three haplotypes. Using the haplotype information from the HapMap CEU data, we identified three common SNPs that were in high LD (r2 = 0.89) with one of the three haplotypes (TAG18~21) , while no SNPs were found with r2 > 0.5 for the other two haplotypes, based on the HapMap CEU data (Rel #22/phase II Apr07, on NCBI B36 assembly, dbSNP b126). Given that the three SNPs were in perfect LD (r2 = 1), we genotyped one of the three SNPs, rs7761846, in our Swedish and Finnish samples. Given that a large association study of ESR1 by the Breast Cancer Association Consortium also revealed a significant association within the same region under a dominant model , we searched for SNPs that were in high LD (r2 > 0.5) with the three haplotypes but were not genotyped in our previous study. Then we performed a genotype-based association analysis under a dominant model of inheritance. As expected, rs7761846 showed association with ER-positive breast cancer (OR = 1.28, P = 0.014) in the combined sample. The two independent Swedish and Finnish samples also revealed consistent association, although the association in the Finnish sample did not reach statistical significance (Table 4).
Genetic interaction between the polymorphisms of PPARGC1B (rs6895698, rs2340621 and rs741581) and ESR1(rs7761846)
Given the known modification of ER activity by PPARGC1B in cellular response to estrogen exposure, we investigated the genetic interaction between rs741581, rs2340621 and rs6895698 within PPARGC1B and rs7761846 within ESR1 in terms of modulating ER-positive breast cancer risk. The analysis in the combined sample identified a significant synergistic interaction between rs2340621 (representing PPARGC1B) and rs7761846 (representing ESR1) (Pinteraction = 0.008) (Table 5). Women carrying both PPARGC1B (rs2340621) and ESR1 (rs7761846) risk genotypes (GA/AA and CT/CC) had a much higher risk for breast cancer than noncarriers (GG and TT) (OR = 1.94, P = 2.03 × 10-6). Similar patterns of genetic interaction were also observed between the remaining SNPs rs741581 (PPARGC1B) and rs7761846 (ESR1) as well as rs6895698 (PPARGC1B) and rs7761846 (ESR1), although these interactions did not achieve statistical significance - probably due to the low minor allele frequencies of rs741581 and rs6895698. However, the significant genetic interactions could not be detected in overall or ER-negative breast cancer (Additional file 2, Tables S6 and S7).
Transcriptional regulation of PPARGC1Bby ERα
To understand the molecular mechanism underlying the observed genetic interaction, we investigated whether there was any transcriptional cross-talk between the two genes beyond the known ligand-dependent, co-activating interaction of the PGC-1β with ERα [33, 34], using the ER-responsive MCF7 breast cancer cell line.
First, we examined the expression of PPARGC1B in MCF7 and noted a twofold induction of PPARGC1B expression by ER activation after estradiol administration (Additional file 3 Figure S1). As a marker of transcriptional activity, chromatin immunoprecipitation-sequencing analysis in the same MCF7 cell line identified a significant peak of RNA polymerase II occupancy close to the transcriptional start site of PPARGC1B within the LD region of SNP association, and the RNA polymerase II occupancy was further enhanced by estradiol treatment. This observation confirms the transcriptional responsiveness of PPARGC1B to estradiol. Moreover, the chromatin immunoprecipitation-sequencing analysis also identified five ER binding sites in and around PPARGC1B (one site located approximately 50 kb 5' of the transcriptional start site, one in the second intron of the gene within 13 kb of the associated SNP rs741581, and the other three binding sites approximately 10, 31 and 57 kb 3' of the polyadenylation signal sequence) and within the LD region of significant association with ER-positive breast cancer (Figure 1F). Interestingly, the sites showing highest of ER occupancy were seen at two locations, one ~13 kb from the significant SNP rs741581 and the second within 31 kb 3' of the polyadenylation signal sequence.
We recently described the identification of all ER binding site interactions in the human genome [35, 36] and defined that genes engaged in chromatin loop formation by a transcription factor were definitively regulated by the factor. Our data indicated that all of the ER binding sites around PPARGC1B were engaged in chromatin loop formation centered on the PPARGC1B gene (Figure 1I), which indicates that ERα directly regulates PPARGC1B.
Taken together, these data strongly indicate that PPARGC1B expression could be directly regulated by ERα and - when coupled with the known enhancement of ERα transcriptional activity by the PGC-1β at the site of binding - suggest a feed-forward regulatory loop between the two genes that augments ER signaling when the two factors are present.
To our knowledge, this is the first comprehensive association analysis of common variation within ER cofactor genes in breast cancer where 36 ER co-activators and 24 ER co-repressors were investigated. The utilization of two independent case-control samples of northern European origin allowed us to identify an association based not only on the overall significance in the large combined sample, but also on the consistency of the SNP association between the two individual samples. We found significant associations between PPARGC1B polymorphisms and risk for ER-positive breast cancer, and, importantly, we revealed a synergistic effect between the genetic polymorphisms within PPARGC1B and ESR1.
Genetic association studies of ER cofactor genes have so far been limited. Burwinkel and colleagues reported a significant association of coding variants Q586 H and T960T of NCOA3 with familial breast cancer risk, and further suggested that familial breast cancer patients may condense the rare allele's contribution to the protective effect of breast cancer . Whilst two studies have reported an association of the variant Pro241Pro in CCND1 with breast cancer risk [37, 38], other studies have reported negative results for this variant [14, 39, 40]. In particular, Wirtenberger and colleagues investigated the coding variant Ala203Pro of PPARGC1B and found it to be associated with familial breast cancer susceptibility . In our study, we did not observe significant association between polymorphisms in NCOA3 and CCND1 with breast cancer risk. The Ala203Pro (rs7732671) variant of PPARGC1B, however, is 10 kb away and not correlated with PPARGC1B SNP rs741581 (r2 < 0.05 in HapMap CEU data), and thus would not have been detected by our tagging SNP approach. Nevertheless, both Wirtenberger and colleagues' study and our study support the association of genetic variation of PPARGC1B with particular subtypes of breast cancer.
Importantly, the association of PPARGC1B as well as its synergistic interaction with ESR1 was only observed in breast cancer patients with ER-positive tumors, as would be expected according to the biochemical mechanism of interaction. There is growing evidence that the impact of genetic risk factors on breast cancer varies by hormone receptor status. For example, recent studies by the Breast Cancer Association Consortium have led to the discovery of novel breast cancer susceptibility loci in FGFR2, TNRC9, 8q24, 2q35, and 5p12 that showed stronger association with ER-positive disease than with ER-negative disease [42–45], with fibroblast growth factor receptor also being a direct target of ER. These data suggest the risk of ER-positive tumors that has been shown to be driven by reproductive factors in epidemiologic studies also has a genomic basis based on the constituents of the ER gene regulatory network [46, 47]. In our study, although the sample sizes of two ER-positive datasets were smaller compared with the two overall datasets, the number of overlapping SNPs between the Swedish and Finnish studies was thus larger than that observed in the overall breast cancer analysis. Recently, we also demonstrated that genetic variation of the estrogen metabolism pathway - particularly the genes involved in the production of estrogen through androgen conversion - also influences the risk for the development of estrogen-sensitive breast cancer . As with this study, the effect size of the metabolism gene polymorphisms are relatively small but, taken together with PPARGC1B and fibroblast growth factor receptor, show that the estrogen receptor signaling axis that engages both upstream and downstream components may have, in the composite, a significant role in the genesis of the most common form of breast cancer.
The genetic interaction between PPARGC1B and ESR1 is biologically plausible. The PPARGC1B protein PGC-1β is a bona fide ER co-activator  that physically interacts with ERα and plays a role in amplifying ER signaling, which provides a convincing biological mechanism for the observed genetic interaction between the two genes. Furthermore, our series of transcriptional regulation analyses in the MCF7 ER-positive breast cancer cell line has demonstrated that PPARGC1B expression can be induced by estrogen treatment, and this transcriptional response of PPARGC1B is probably mediated by five functional ER binding sites around PPARGC1B that are all engaged in interlocking chromatin loops highly indicative of an ER regulated gene . PPARGC1B may thus be involved in a feed-forward control mechanism with ERα such that ER induction (for example, by estradiol treatment) heightens the expression of a co-activator PPARGC1B of ER, which in turn increases ER action at the DNA binding site. The feed-forward looping mechanism will therefore further augment the protein interaction between PPARGC1B and ESR1. This putative amplification effect, if confirmed, is another mechanistic model for epistatic interactions between genetic loci and may be one reason for the strength of its signal in the association study as compared with the other ER cofactors studied.
There are some limitations to our study. Coverage of common variation is not sufficient (< 80%) for some genes (Additional file 1, Table S1), so that some associations may have been missed. In addition, our tagging SNP selection provides a rather limited coverage of 5 kb surrounding sequences of the candidate genes, which may have contributed to some associations of regulatory SNPs being undetected, such as the one reported within ESR1 . The number of overlapping SNPs between the two datasets is small for both ER-positive and overall breast cancer analyses. The limited overlapping could be due to ethnic heterogeneity between the two population samples and their moderate sample sizes. On the one hand, the ethnic heterogeneity may partially explain the low overlapping SNPs between two datasets; on the other hand, the current sample size is not large enough to capture the moderate effect of associated SNPs. Some of the top SNPs for each individual sample set are therefore probably false positive, which causes the small overlap between the numbers of significant SNPs in both datasets. The sample size limitation in ER-negative patients also could lead to the nonsignificant results in ER-negative analysis, since we observed that some associations in ER-negative analysis are in the same direction with ER-positive analysis. ER cofactors are known to work as a multicomponent protein complex, but due to a sample size limitation we are unable to detect interaction among three or more genes simultaneously. It is also worth noting that the contribution of genetic variation to cancer risk is based on both their prevalence and penetrance, and thus the relative importance of individual SNPs may vary from population to population. Further confirmation of our findings in other populations is therefore warranted.
Our study has revealed an association of genetic variation within PPARGC1B with the risk of ER-positive breast cancer. Consistent with the known interaction of PPARGC1B and ER at the molecular level, where PPARGC1B modulates ER activity and thus ER signaling, our study revealed a synergistic effect between genetic variation within the PPARGC1B and ESR1 genes. PPARGC1B has been shown to alter responses to the selective ER modulator, tamoxifen . Kressler and colleagues also demonstrated that PPARGC1B indirectly co-activates tamoxifen-bound ERα, which cooperates with NCOA1 to enable tamoxifen agonism in kidney and osteosarcoma cell lines. Lastly, the synergism demonstrated in the present study also suggests that disrupting the interaction between an ER co-activator - such as PPARGC1B - and ERα, or blocking their mutual activation, may represent a sensitive and leveraged strategy for cancer prevention . Our study therefore provides new biological insight into the genetic basis of the more common ER-positive breast cancer and highlights that biochemically and genomically informed candidate gene study can enhance the discovery of interactive disease susceptibility genes.
- CCND1 :
Dulbecco's modified Eagle's medium
- ESR1 :
estrogen receptor 1
fetal bovine serum
- FGFR2 :
fibroblast growth factor receptor 2
- NCOA :
nuclear receptor coactivator
polymerase chain reaction
peroxisome proliferative activated receptor gamma coactivator beta protein
- PPARGC1B :
peroxisome proliferative activated receptor gamma coactivator beta
single nucleotide polymorphism
- TNRC9 :
TOX high mobility group box family member 3.
Key T, Appleby P, Barnes I, Reeves G: Endogenous sex hormones and breast cancer in postmenopausal women: reanalysis of nine prospective studies. J Natl Cancer Inst. 2002, 94: 606-616.
Missmer SA, Eliassen AH, Barbieri RL, Hankinson SE: Endogenous estrogen, androgen, and progesterone concentrations and breast cancer risk among postmenopausal women. J Natl Cancer Inst. 2004, 96: 1856-1865. 10.1093/jnci/djh336.
Dickson RB, Stancel GM: Estrogen receptor-mediated processes in normal and cancer cells. J Natl Cancer Inst Monogr. 2000, 27: 135-145.
Gao X, Loggie BW, Nawaz Z: The roles of sex steroid receptor coregulators in cancer. Mol Cancer. 2002, 1: 7-
Smith CL, O'Malley BW: Coregulator function: a key to understanding tissue specificity of selective receptor modulators. Endocr Rev. 2004, 25: 45-71. 10.1210/er.2003-0023.
Tremblay GB, Giguere V: Coregulators of estrogen receptor action. Crit Rev Eukaryot Gene Expr. 2002, 12: 1-22. 10.1615/CritRevEukaryotGeneExpr.v12.i1.10.
Osborne CK, Schiff R: Estrogen-receptor biology: continuing progress and therapeutic implications. J Clin Oncol. 2005, 23: 1616-1622. 10.1200/JCO.2005.10.036.
Johnston SR: New strategies in estrogen receptor-positive breast cancer. Clin Cancer Res. 2010, 16: 1979-1987. 10.1158/1078-0432.CCR-09-1823.
Einarsdottir K, Darabi H, Li Y, Low YL, Li YQ, Bonnard C, Sjolander A, Czene K, Wedren S, Liu ET, Hall P, Humphreys K, Liu J: ESR1 and EGF genetic variation in relation to breast cancer risk and survival. Breast Cancer Res. 2008, 10: R15-10.1186/bcr1861.
Dunning AM, Healey CS, Baynes C, Maia AT, Scollen S, Vega A, Rodriguez R, Barbosa-Morais NL, Ponder BA, Low YL, Bingham S, Haiman CA, Le Marchand L, Broeks A, Schmidt MK, Hopper J, Southey M, Beckmann MW, Fasching PA, Peto J, Johnson N, Bojesen SE, Nordestgaard B, Milne RL, Benitez J, Hamann U, Ko Y, Schmutzler RK, Burwinkel B, Schurmann P, et al: Association of ESR1 gene tagging SNPs with breast cancer risk. Hum Mol Genet. 2009, 18: 1131-1139. 10.1093/hmg/ddn429.
Zheng W, Long J, Gao YT, Li C, Zheng Y, Xiang YB, Wen W, Levy S, Deming SL, Haines JL, Gu K, Fair AM, Cai Q, Lu W, Shu XO: Genome-wide association study identifies a new breast cancer susceptibility locus at 6q25.1. Nat Genet. 2009, 41: 324-328. 10.1038/ng.318.
Burwinkel B, Wirtenberger M, Klaes R, Schmutzler RK, Grzybowska E, Forsti A, Frank B, Bermejo JL, Bugert P, Wappenschmidt B, Butkiewicz D, Pamula J, Pekala W, Zientek H, Mielzynska D, Siwinska E, Bartram CR, Hemminki K: Association of NCOA3 polymorphisms with breast cancer risk. Clin Cancer Res. 2005, 11: 2169-2174. 10.1158/1078-0432.CCR-04-1621.
Ceschi M, Sun CL, Van Den Berg D, Koh WP, Yu MC, Probst-Hensch N: The effect of cyclin D1 (CCND1) G870A-polymorphism on breast cancer risk is modified by oxidative stress among Chinese women in Singapore. Carcinogenesis. 2005, 26: 1457-1464. 10.1093/carcin/bgi093.
Haiman CA, Garcia RR, Hsu C, Xia L, Ha H, Sheng X, Le Marchand L, Kolonel LN, Henderson BE, Stallcup MR, Greene GL, Press MF: Screening and association testing of common coding variation in steroid hormone receptor co-activator and co-repressor genes in relation to breast cancer risk: the Multiethnic Cohort. BMC Cancer. 2009, 9: 43-10.1186/1471-2407-9-43.
Montgomery KG, Chang JH, Gertig DM, Dite GS, McCredie MR, Giles GG, Southey MC, Hopper JL, Campbell IG: The AIB1 glutamine repeat polymorphism is not associated with risk of breast cancer before age 40 years in Australian women. Breast Cancer Res. 2005, 7: R353-R356. 10.1186/bcr1009.
Hall JM, McDonnell DP: Coregulators in nuclear estrogen receptor action: from concept to therapeutic targeting. Mol Interv. 2005, 5: 343-357. 10.1124/mi.5.6.7.
McKenna NJ, Lanz RB, O'Malley BW: Nuclear receptor coregulators: cellular and molecular biology. Endocr Rev. 1999, 20: 321-344. 10.1210/er.20.3.321.
Wedren S, Lovmar L, Humphreys K, Magnusson C, Melhus H, Syvanen AC, Kindmark A, Landegren U, Fermer ML, Stiger F, Persson I, Baron J, Weiderpass E: Oestrogen receptor alpha gene haplotype and postmenopausal breast cancer risk: a case control study. Breast Cancer Res. 2004, 6: R437-R449. 10.1186/bcr811.
Kilpivaara O, Bartkova J, Eerola H, Syrjakoski K, Vahteristo P, Lukas J, Blomqvist C, Holli K, Heikkila P, Sauter G, Kallioniemi OP, Bartek J, Nevanlinna H: Correlation of CHEK2 protein expression and c.1100delC mutation status with tumor characteristics among unselected breast cancer patients. Int J Cancer. 2005, 113: 575-580. 10.1002/ijc.20638.
Syrjakoski K, Vahteristo P, Eerola H, Tamminen A, Kivinummi K, Sarantaus L, Holli K, Blomqvist C, Kallioniemi OP, Kainu T, Nevanlinna H: Population-based study of BRCA1 and BRCA2 mutations in 1035 unselected Finnish breast cancer patients. J Natl Cancer Inst. 2000, 92: 1529-1531. 10.1093/jnci/92.18.1529.
Fagerholm R, Hofstetter B, Tommiska J, Aaltonen K, Vrtel R, Syrjakoski K, Kallioniemi A, Kilpivaara O, Mannermaa A, Kosma VM, Uusitupa M, Eskelinen M, Kataja V, Aittomäki K, von Smitten K, Heikkila P, Lukas J, Holli K, Bartkova J, Blomqvist C, Bartek J, Nevanlinna H: NAD(P)H:quinone oxidoreductase 1 NQO1*2 genotype (P187S) is a strong prognostic and predictive factor in breast cancer. Nat Genet. 2008, 40: 844-853. 10.1038/ng.155.
Vahteristo P, Bartkova J, Eerola H, Syrjakoski K, Ojala S, Kilpivaara O, Tamminen A, Kononen J, Aittomäki K, Heikkila P, Holli K, Blomqvist C, Bartek J, Kallioniemi OP, Nevanlinna H: A CHEK2 genetic variant contributing to a substantial fraction of familial breast cancer. Am J Hum Genet. 2002, 71: 432-438. 10.1086/341943.
Eerola H, Blomqvist C, Pukkala E, Pyrhonen S, Nevanlinna H: Familial breast cancer in southern Finland: how prevalent are breast cancer families and can we trust the family history reported by patients?. Eur J Cancer. 2000, 36: 1143-1148. 10.1016/S0959-8049(00)00093-9.
Magnusson C, Baron JA, Correia N, Bergstrom R, Adami HO, Persson I: Breast-cancer risk following long-term oestrogen- and oestrogen-progestin-replacement therapy. Int J Cancer. 1999, 81: 339-344. 10.1002/(SICI)1097-0215(19990505)81:3<339::AID-IJC5>3.0.CO;2-6.
Magnusson C, Colditz G, Rosner B, Bergstrom R, Persson I: Association of family history and other risk factors with breast cancer risk (Sweden). Cancer Causes Control. 1998, 9: 259-267. 10.1023/A:1008817018942.
Magnusson CM, Persson IR, Baron JA, Ekbom A, Bergstrom R, Adami HO: The role of reproductive factors and use of oral contraceptives in the aetiology of breast cancer in women aged 50 to 74 years. Int J Cancer. 1999, 80: 231-236. 10.1002/(SICI)1097-0215(19990118)80:2<231::AID-IJC11>3.0.CO;2-R.
The International HapMap Consortium: A haplotype map of the human genome. Nature. 2005, 437: 1299-1320. 10.1038/nature04226.
Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA: Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet. 2004, 74: 106-120. 10.1086/381000.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.
Liu J, Vasant DK, Li Y, Banu S, Lee YK, Qu D, Heikkinen T, Aaltonen K, Muranen TA, Kajiji TS, Bonnard C, Aittomäki K, von Smitten K, Blomqvist C, Hopper JL, Southey MC, Brauch H, Consortium TG, Chenevix-Trench G, Beesley J, Spurdle AB, Chen X, Kathleen Cuningham Foundation Consortium for Research into Familial Breast Cancer AOCSG, Czene K, Hall P, Nevanlinna H, Liu ET: Germ-line variation at a functional p53 binding site increases susceptibility to breast cancer development. HUGO J. 2009, 3: 31-40. 10.1007/s11568-010-9138-x.
Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, NJ: Humana Press, 365-386.
Li Y, Sung WK, Liu JJ: Association mapping via regularized regression analysis of single-nucleotide-polymorphism haplotypes in variable-sized sliding windows. Am J Hum Genet. 2007, 80: 705-715. 10.1086/513205.
Kressler D, Hock MB, Kralli A: Coactivators PGC-1β and SRC-1 interact functionally to promote the agonist activity of the selective estrogen receptor modulator tamoxifen. J Biol Chem. 2007, 282: 26897-26907. 10.1074/jbc.M705596200.
Kressler D, Schreiber SN, Knutti D, Kralli A: The PGC-1-related protein PERC is a selective coactivator of estrogen receptor alpha. J Biol Chem. 2002, 277: 13918-13925. 10.1074/jbc.M201134200.
Fullwood MJ, Liu MH, Pan YF, Liu J, Xu H, Mohamed YB, Orlov YL, Velkov S, Ho A, Mei PH, Chew EG, Huang PY, Welboren WJ, Han Y, Ooi HS, Ariyaratne PN, Vega VB, Luo Y, Tan PY, Choy PY, Wansa KD, Zhao B, Lim KS, Leow SC, Yow JS, Joseph R, Li H, Desai KV, Thomsen JS, Lee YK, et al: An oestrogen-receptor-alpha-bound human chromatin interactome. Nature. 2009, 462: 58-64. 10.1038/nature08497.
Li G, Fullwood MJ, Xu H, Mulawadi FH, Velkov S, Vega V, Ariyaratne PN, Mohamed YB, Ooi HS, Tennakoon C, Wei CL, Ruan Y, Sung WK: ChIA-PET tool for comprehensive chromatin interaction analysis with paired-end tag sequencing. Genome Biol. 2010, 11: R22-10.1186/gb-2010-11-2-r22.
Onay UV, Aaltonen K, Briollais L, Knight JA, Pabalan N, Kilpivaara O, Andrulis IL, Blomqvist C, Nevanlinna H, Ozcelik H: Combined effect of CCND1 and COMT polymorphisms and increased breast cancer risk. BMC Cancer. 2008, 8: 6-10.1186/1471-2407-8-6.
Yu CP, Yu JC, Sun CA, Tzao C, Ho JY, Yen AM: Tumor susceptibility and prognosis of breast cancer associated with the G870A polymorphism of CCND1. Breast Cancer Res Treat. 2008, 107: 95-102. 10.1007/s10549-007-9522-y.
Driver KE, Song H, Lesueur F, Ahmed S, Barbosa-Morais NL, Tyrer JP, Ponder BA, Easton DF, Pharoah PD, Dunning AM: Association of single-nucleotide polymorphisms in the cell cycle genes with breast cancer in the British population. Carcinogenesis. 2008, 29: 333-341. 10.1093/carcin/bgm284.
Krippl P, Langsenlehner U, Renner W, Yazdani-Biuki B, Wolf G, Wascher TC, Paulweber B, Weitzer W, Leithner A, Samonigg H: The 870G > A polymorphism of the cyclin D1 gene is not associated with breast cancer. Breast Cancer Res Treat. 2003, 82: 165-168. 10.1023/B:BREA.0000004372.20461.33.
Wirtenberger M, Tchatchou S, Hemminki K, Schmutzhard J, Sutter C, Schmutzler RK, Meindl A, Wappenschmidt B, Kiechle M, Arnold N, Weber BH, Niederacher D, Bartram CR, Burwinkel B: Associations of genetic variants in the estrogen receptor coactivators PPARGC1A, PPARGC1B and EP300 with familial breast cancer. Carcinogenesis. 2006, 27: 2201-2208. 10.1093/carcin/bgl067.
Stacey SN, Manolescu A, Sulem P, Thorlacius S, Gudjonsson SA, Jonsson GF, Jakobsdottir M, Bergthorsson JT, Gudmundsson J, Aben KK, Strobbe LJ, Swinkels DW, van Engelenburg KC, Henderson BE, Kolonel LN, Le Marchand L, Millastre E, Andres R, Saez B, Lambea J, Godino J, Polo E, Tres A, Picelli S, Rantala J, Margolin S, Jonsson T, Sigurdsson H, Jonsdottir T, Hrafnkelsson J, et al: Common variants on chromosome 5p12 confer susceptibility to estrogen receptor-positive breast cancer. Nat Genet. 2008, 40: 703-706. 10.1038/ng.131.
Stacey SN, Manolescu A, Sulem P, Rafnar T, Gudmundsson J, Gudjonsson SA, Masson G, Jakobsdottir M, Thorlacius S, Helgason A, Aben KK, Strobbe LJ, Albers-Akkers MT, Swinkels DW, Henderson BE, Kolonel LN, Le Marchand L, Millastre E, Andres R, Godino J, Garcia-Prats MD, Polo E, Tres A, Mouy M, Saemundsdottir J, Backman VM, Gudmundsson L, Kristjansson K, Bergthorsson JT, Kostic J, et al: Common variants on chromosomes 2q35 and 16q12 confer susceptibility to estrogen receptor-positive breast cancer. Nat Genet. 2007, 39: 865-869. 10.1038/ng2064.
Milne RL, Benitez J, Nevanlinna H, Heikkinen T, Aittomäki K, Blomqvist C, Arias JI, Zamora MP, Burwinkel B, Bartram CR, Meindl A, Schmutzler RK, Cox A, Brock I, Elliott G, Reed MW, Southey MC, Smith L, Spurdle AB, Hopper JL, Couch FJ, Olson JE, Wang X, Fredericksen Z, Schurmann P, Bremer M, Hillemanns P, Dork T, Devilee P, van Asperen CJ, et al: Risk of estrogen receptor-positive and -negative breast cancer and single-nucleotide polymorphism 2q35-rs13387042. J Natl Cancer Inst. 2009, 101: 1012-1018. 10.1093/jnci/djp167.
Garcia-Closas M, Hall P, Nevanlinna H, Pooley K, Morrison J, Richesson DA, Bojesen SE, Nordestgaard BG, Axelsson CK, Arias JI, Milne RL, Ribas G, Gonzalez-Neira A, Benitez J, Zamora P, Brauch H, Justenhoven C, Hamann U, Ko YD, Bruening T, Haas S, Dork T, Schurmann P, Hillemanns P, Bogdanova N, Bremer M, Karstens JH, Fagerholm R, Aaltonen K, Aittomäki K, et al: Heterogeneity of breast cancer associations with five susceptibility loci by clinical and pathological characteristics. PLoS Genet. 2008, 4: e1000054-10.1371/journal.pgen.1000054.
Althuis MD, Fergenbaum JH, Garcia-Closas M, Brinton LA, Madigan MP, Sherman ME: Etiology of hormone receptor-defined breast cancer: a systematic review of the literature. Cancer Epidemiol Biomarkers Prev. 2004, 13: 1558-1568.
Garcia-Closas M, Chanock S: Genetic susceptibility loci for breast cancer by estrogen receptor status. Clin Cancer Res. 2008, 14: 8000-8009. 10.1158/1078-0432.CCR-08-0975.
Low YL, Li Y, Humphreys K, Thalamuthu A, Darabi H, Wedren S, Bonnard C, Czene K, Iles MM, Heikkinen T, Aittomäki K, Blomqvist C, Nevanlinna H, Hall P, Liu ET, Liu J: Multi-variant pathway association analysis reveals the importance of genetic determinants of estrogen metabolism in breast and endometrial cancer susceptibility. PLoS Genet. 2010, 6: e1001012-10.1371/journal.pgen.1001012.
The authors are grateful to all of the women who took the time and effort to participate in the present study. They would also like to thank Heng Khai Koon and Ong Eng Hu Jason for genotyping, Frans Verhoeff for processing the genotyping data and the Finnish Cancer Registry for providing cancer data. The present work was supported by funding from the Agency for Science, Technology and Research of Singapore (A_STAR), the Nordic Cancer Union (both for Sweden and Finland) and the Märit and Hans Rausing's Initiative Against Breast Cancer. In addition, the Finnish study has been financially supported by the Helsinki University Central Hospital Research Fund, Academy of Finland (132473), the Finnish Cancer Society, and the Sigrid Juselius Foundation.
The authors declare that they have no competing interests.
YQL, SW, KE, HN, KSC, PH, ETL and JJL initiated and designed the study. CB, KC and HN provided the study material and patient information. YQL, CB, TH, KA and SW collected and organized the data. YQL, YL, GLL, THC, DKV, JJL, KH and HD performed data analysis and interpreted results. YQL, SW, YL, DKV, KH, PH, ETL and JJL drafted the manuscript. CB, HD, KH, KE, HN, PH, ETL and JJL performed critical review and revised the manuscript. All authors read and approved the manuscript.
Electronic supplementary material
Additional file 1: Supplementary results of coverage evaluation of common variants and association analysis in ER cofactor genes. Table S1 presenting coverage evaluation of the common variant in 60 ER cofactor genes. Table S2 presenting ORs and P values of the consistent SNPs between the Swedish and Finnish samples from the analyses of ER-positive and overall breast cancer. Table S3 presenting the 25 most significant SNPs in ER-negative association analysis in Swedish and Finnish samples. Table S4 presenting the 25 most significant SNPs in overall association analysis in Swedish and Finnish samples. Table S5 presenting the comparison of P value among additive, dominant and recessive models in the analysis of ER-positive breast cancer in PPARGC1B in the combined Swedish and Finnish samples. (DOC 239 KB)
Additional file 2: Analysis of the pair-wise interaction effect between SNPs within PPARGC1B and ESR1 on the overall and ER-negative breast cancer in the combined Swedish and Finnish samples. Table S6 presenting analysis of the pair-wise interaction effect between SNPs within PPARGC1B and ESR1 on the overall breast cancer in the combined Swedish and Finnish samples. Table S7 presenting the analysis of pair-wise interaction effect between SNPs within PPARGC1B and ESR1 on the ER-negative breast cancer in the combined Swedish and Finnish samples. (DOC 55 KB)
Additional file 3: Relative expression of PPARGC1B gene in MCF7 cells 3 hrs post E2 treatment. Figure S1 presenting a relative expression of the PPARGC1B gene in MCF7 cells 3 hours post 17β-estradiol treatment. (DOC 34 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Li, Y., Li, Y., Wedrén, S. et al. Genetic variation of ESR1 and its co-activator PPARGC1B is synergistic in augmenting the risk of estrogen receptor-positive breast cancer. Breast Cancer Res 13, R10 (2011). https://doi.org/10.1186/bcr2817
- Breast Cancer
- Estrogen Receptor
- Breast Cancer Risk
- Familial Breast Cancer
- Swedish Sample