- Research article
- Open Access
Estrogen receptor alpha gene polymorphism and endometrial cancer risk – a case-control study
BMC Cancer volume 8, Article number: 322 (2008)
Estrogen is an established endometrial carcinogen. One of the most important mediators of estrogenic action is the estrogen receptor alpha. We have investigated whether polymorphic variation in the estrogen receptor alpha gene (ESR1) is associated with endometrial cancer risk.
In 702 cases with invasive endometrial cancer and 1563 controls, we genotyped five markers in ESR1 and used logistic regression models to estimate odds ratios (OR) and 95 percent confidence intervals (CI).
We found an association between rs2234670, rs2234693, as well as rs9340799, markers in strong linkage disequilibrium (LD), and endometrial cancer risk. The association with rs9340799 was the strongest, OR 0.75 (CI 0.60–0.93) for heterozygous and OR 0.53 (CI 0.37–0.77) for homozygous rare compared to those homozygous for the most common allele. Haplotype models did not fit better to the data than single marker models.
We found that intronic variation in ESR1 was associated with endometrial cancer risk.
Estrogens, whether from endogenous or exogenous sources, are potent mitogens in the endometrium and thus constitute a major carcinogen in this tissue . Endometrial cancer is therefore a good model for investigating clinical effects of estrogen signaling. Estrogen receptor alpha is the principal estrogen receptor expressed in the endometrium [2, 3] and it is considered to be crucial in the development of endometroid endometrial carcinoma, the most common histological subtype among endometrial neoplasms . Since altered ER function could have impact on the risk of endometrial cancer, the estrogen receptor alpha gene (ESR1) is a plausible endometrial cancer candidate gene which has been investigated in a few studies [4–6]. We have investigated whether polymorphic variation in the ESR1 is associated with invasive endometrial cancer risk, overall and in subgroups defined by other hormonally related factors.
This nation-wide population-based case-control study encompassed all cases of incident, primary histopathologically confirmed invasive endometrial cancer among women 50 to 74 years of age resident in Sweden between January 1, 1994 and December 31, 1995 as previously described in detail . Endometrial cancer patients were identified at diagnosis through a notification system organized within the six Swedish regional cancer registries, to which reporting of all malignant tumors is mandatory. Controls were women randomly selected in 5-year age strata (to match the expected age frequency distribution among the cases) from the Swedish Total Population Register. Only women with an intact uterus and no previous history of endometrial cancer were eligible for this study.
Endometrial cancer patients were asked to participate by their respective physicians. After their informed consent, we mailed self-administered questionnaires asking for detailed information about use of menopausal hormones and oral contraceptives, weight, height, reproductive history, medical history, and other lifestyle factors. Controls were contacted directly with the questionnaire. Seventy-six percent of eligible cases (n = 802) and 84 percent of the controls (n = 3 550) ultimately participated in the parent study by answering the questionnaire. Among these controls, 491 who failed to return the mailed questionnaire were interviewed by phone. Slides with tumor tissue were collected from all cases and reviewed by one study pathologist to confirm the diagnosis, and to document tumor grade and myometrial invasion. The specimens were reclassified as: endometroid adenocarcinoma (n = 648), seropapillary carcinoma (n = 36), clear cell carcinoma (n = 10), adenoacanthoma (n = 3), adenosquamous carcinoma (n = 12), endometrial atypical hyperplasia (n = 80) (adenomatous hyperplasia with slight, moderate or severely pronounced atypia), and anaplastic carcinoma (n = 13). Results from the parent study have been published [7–10].
Selection of present study population
In the analysis, we included all women with endometrial cancer (n = 802) and 802 randomly selected controls (frequency-matched by age) from the parent study. In order to increase statistical power in subgroup analyses, we additionally selected all remaining eligible controls who had taken menopausal hormone treatment, either preparations containing medium potency estrogen only (mainly estradiol or conjugated estrogens) or medium potency estrogen in combination with progestin for at least 2 years (277 controls) and all remaining control women with self-reported diabetes mellitus (124 controls). For a detailed description of the hormone treatment variables, please see . In total, 802 cases and 1 203 controls were selected. In addition, we included 871 controls from the parent study who had also participated and been genotyped in a parallel breast cancer study, who fulfilled the inclusion criteria (including having an intact uterus). The Institutional Review Board at Karolinska Institutet, Stockholm, Sweden and the six other Swedish regional review boards approved this study.
Collection of biological material
All selected living women were contacted by mail and those who gave informed consent received a blood sampling kit by mail. Whole blood samples were drawn at a primary health care facility close to the woman's home and were sent to us by standard mail. A majority of the samples arrived at Karolinska Institutet within 24 hours after blood draw. All blood samples were immediately stored at -20°C. Endometrial cancer cases who declined to donate a blood sample were asked to permit to our use of normal tissue from archived paraffin embedded tissue taken at cancer surgery (consisting of, for example, cancer-free lymph nodes, uterine tube or myometrium). We also attempted to retrieve archived tissue samples from all deceased endometrial cancer cases. Samples were coded and transferred to the laboratories. Blood samples or archived tissue samples were obtained for 603 and 104 endometrial cancer patients, respectively, and blood samples from in total 1 574 control women (922 selected for this study + 652 genotyped in the breast cancer study), yielding participation rates of 88 percent for cases (707 of 802 eligible) and 76 percent for controls (1 574 of 2 074 eligible). Considering also the participation rate in the parent study, overall participation rates were 67 percent and 64 percent among cases and controls, respectively.
Selection of SNPs and microsatellite markers
We selected ESR1 polymorphisms to be analyzed from the literature [11–14]. The polymorphisms were not intended to capture variation across the entire gene. The 5' promoter microsatellite (TAn) rs2234670 is located -1174-base pairs upstream of exon 1, while the SNPs rs2234693 (also known as PvuII) and rs9340799 (also known as XbaI) are found in intron 1, and the synonymous coding rs4986934 (codon 243 CGC>CGT) and rs1801132 SNPs (codon 325 CCC>CCG) are located in exon 3 and 4, respectively.
DNA was isolated from 3 ml whole blood using the Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA) according to the manufacturer's instructions. From non-malignant cells in paraffin embedded tissue, we extracted DNA using a standard phenol/chloroform/isoamylalcohol protocol .
SNPs and microsatellite markers
All primers and probes for the 5' promoter microsatellite rs2234670, the intron 1 SNPs rs2234693 and rs9340799, and the synonymous exonic rs4986934 and rs1801132 SNPs, were designed based on the reference sequences AF082876, AF326912 and NM_000125, respectively. For the fluorescence polarization and the minisequencing assays , we designed minisequencing primers complementary to the sequence immediately adjacent to the SNPs. For the Molecular Beacon assay , we designed two fluorescently labeled allele specific probes for each SNP, carrying the variable position in the middle of the loop region . We used a web-based DNA folding program (mfold)  to estimate the stability of the stem and loop structure of the MB probes. Supplementary tables 1 and 2 (Additional file 1) give the primer and probe sequences and their modifications.
Minisequencing assay using fluorescence polarization detection
The region containing the intron 1 SNPs, rs2234693 and rs9340799, was amplified in one PCR fragment. PCR was performed using 20 ng of genomic DNA, 0.02 U/μl Taq Gold DNA polymerase (Applied Biosystems, Foster City, CA, USA), 100 μM dNTPs (Amersham Pharmacia Biotech, Uppsala, Sweden) and 0.1 μM of both PCR primers in 30 μl of GeneAmp PCR buffer supplied with the enzyme at 95°C for 7 minutes, followed by 40 cycles at 92°C for 15 seconds, 56°C for 1 minute and 72°C for 45 seconds. Five microliters of PCR product was aliquoted to 384-well microtiter plates and 5 μl of a PCR cleanup mixture containing 0.1 U/μl shrimp alkaline phosphatase (USB Corporation, Cleveland, OH) and 0.1 U/μl Exonuclease I (USB Corporation) in 20 mM Tris-HCl pH 8.0 and 10 mM MgCl2 was added.
The plate was incubated at 37°C for 1 hour, followed by 15 minutes at 95°C to inactivate the enzymes. Five microliters of extension mixture containing 5 μM extension primers, and 0.08 U/μl Thermosequenase (Amersham Pharmacia Biotech) in 50 mM Tris-HCl pH 9, 50 μM KCl, 5 mM NaCl, 5 mM MgCl2, 8% glycerol, and fluorescently labeled and unlabeled ddNTPs was added to each reaction. The two ddNTPs relevant for the particular SNP were labeled and included at a 1:5 ratio relative to unlabeled ddNTPs. The total concentration of ddNTP was 0.125 μM. The fluorescent dye rhodamine 110 (R110) was used to label ddGTP and ddCTP, and 6-carboxytetramethylrhodamine (TAMRA) were used to label ddATP and ddUTP (PerkinElmer Life Sciences, Boston, MA, USA). The cyclic extension reaction was performed at 92°C for 1 minute followed by 40 cycles of 92°C for 10 seconds and 55° for 30 seconds. The fluorescence signals in the 384 well plates were read on an Analyst AD™ (Molecular Devices Corporation, Sunnyvale, CA, USA). Genotypes were assigned using the software AlleleCaller ™ supplied with the instrument, and a custom-made Excel macro. In each run, positive controls for the three genotypes and negative controls were included. Both DNA polarities were analyzed and the results were concordant in all samples. Approximately 30 percent of the assays were repeated and the results were identical. In addition, 3% of the genotypes were validated by solid-phase minisequencing.
Molecular Beacon assay
To analyze the exonic SNPs rs4986934 and rs1801132, 10 ng of genomic DNA was amplified in a 25 μl reaction containing 0.02 U/μl of AmpliTaq Gold DNA polymerase (Applied Biosystems), 3.5 mM MgCl2, 250 μM of dNTPs (Amersham Pharmacia Biotech) and 1 μM of PCR primers in TaqMan Buffer A containing the ROX dye (Applied Biosystems). The two MB probes were included at 0.34 μM concentration. The cycling conditions were 95°C for 10 minutes, followed by 40 cycles at 95°C for 30 seconds, 55°C (rs1801132) or 57°C (rs4986934) for 1 minute and 72°C for 45 seconds and was performed using an ABI Prism 7700 Sequence Detection System (Applied Biosystems). The increase in fluorescent signal was registered during the annealing step of the reaction and the end-point signals were used to assign the genotypes as previously described . In each run, positive controls for the three genotypes and negative controls were included. Both possible nucleotides of the SNP were interrogated in the same reaction. Approximately 3% of the assays were repeated and the results were identical. Two percent of the assays were validated by solid-phase minisequencing. The Molecular Beacon assay has previously been quantitatively validated in our laboratory .
Solid-phase minisequencing assay
Solid-phase minisequencing in microtiter plates  was used in part to genotype rs4986934. The assay also served as a reference method for the other SNP assays.
The dinucleotide repeat rs2234670 region was amplified using an ABI-877 Integrated Thermal cycler PCR robot with standard reagents (Applied Biosystems) . The cycling conditions were 96°C for 10 minutes followed by 36 cycles of 30 seconds at 96°C, 30 seconds at 57°C, and 1 minute at 72°C with a final extension step for 7 minutes. The HEX-labeled fluorescent PCR products were mixed with an internal-lane standard (GS-500 TAMRA) and separated on a 6% polyacrylamide gel using a 96-well ABI 377 automatic sequenator. Gel data were analyzed using Genescan Analysis 3.1 and allele sizes were determined using Genotyper 2.0 (all Applied Biosystems). In each run, 2–3 negative controls were included. Approximately 1% of the assays were repeated and gave concordant results. After calculation of actual number of TA repeats there was no difference between repeated runs of the same sample. The assay was validated by control sequencing of three different repeat lengths.
We verified that genotype frequencies among controls were in accordance with expectations under Hardy-Weinberg equilibrium (HWE) using χ2 statistics. We considered the rs2234670 in its original form, as 8 categories, or dichotomized, into long and short (≤ 14 or > 14 repeats), as the TAn lengths were bimodally distributed with peaks at 11 and 18 repeats, with a dip at 14 repeats. In addition, we estimated pairwise linkage disequilibrium values |D'| and r2 [22, 23].
We used conditional logistic regression to calculate odds ratios (OR) and 95 percent confidence intervals (CI) for association between single loci genotypes and endometrial cancer risk. We conditioned on the variables used for selection, i.e. age in 5-year categories, use of menopausal hormones for two years or more, and self-reported diabetes mellitus.
We tested for associations between ESR1 polymorphisms and other covariates by scrutinizing 3 × k tables and calculating χ2 statistics. Covariates that were approximately normally distributed (body mass index (BMI, weight in kg/height in m2) and age at menopause) were tested by means of two-sample t test over genotype categories. The following covariates were considered: age (in 5-year age groups), use of any preparation containing medium potency estrogen with or without added progestin (never, < 2 years, ≥ 2 years), use of preparation containing only medium potency estrogen without progestin (never, < 2 years, ≥ 2 years), BMI (< 25, 25–27.99, ≥ 28 or as a continuous variable), age at menopause (in tertiles among controls or as continuous variable), parity (nulliparous, 1 childbirth, ≥ 2 childbirths), smoking at least 100 cigarettes in their lifetime or regularly for at least 1 year (ever, never), use of combined oral contraceptives (ever, never), age at last birth (tertiles or as a continuous variable, only including women who ever had given birth), and age at menarche (tertiles or as a continuous variable). For detailed descriptions of how the above variables were defined, please see . All covariates were introduced into the logistic regression model to detect any indications of confounding or other associations potentially affecting the primary association between genotypes and endometrial cancer. We used likelihood ratio (LR) tests test to compare successive models, all of which included the same set of endometrial cancer cases and controls. We performed these analyses using the SAS system (Release 8.02, SAS Institute Inc. Cary, NC, USA).
We tested for disease-haplotype association using likelihood-based approaches contained in the software EH plus  as well as in routines written in the S-PLUS (Insightful) programming language. Most available software (including EH plus) for testing association between haplotypes and case/control status assume multiplicative penetrance, that is, the OR comparing two haplotype copies against none is assumed to be the square of the OR comparing one against none.
We successfully genotyped 695 cases and 1564 controls for the rs2234670, 698 cases and 1561 controls for the rs2234693, 698 cases and 1562 controls for the rs9340799, 701 cases and 1563 controls for the rs4986934, and 702 cases and 1563 controls for the rs1801132 marker. In controls, SNP genotype frequencies were in HWE but the rs2234670 genotype frequencies were not, whether in dichotomised form (p = 0.003) or when all alleles were considered (p < 0.0001).
The two intron 1 markers were in strong LD while there was weaker pairwise LD between the remaining markers (Table 1).
Risk factor exposure differences between cases and controls reflected established associations (Table 2). Ninety-three percent of the tumours in the following analyses were of endometroid histological type. Myometrial invasion was available for 79 percent of the cases, out of which 31 percent had invasion through more than half of the myometrial thickness. Thirty-eight, 45 and 17 percent were of histopathological grade 1, grade 2, and grade 3, respectively. We found no convincing associations between genotypes and any of the covariates (data not shown). There was no difference in genotype frequencies between cases who donated blood samples and those who participated via tissue samples (P = 0.96, 0.66, 0.48, 0.56, and 0.51 for rs2234670, rs2234693, rs9340799, rs4986934, and rs1801132, respectively) or between primary selected controls and controls included after they had been genotyped for the parallel breast cancer study (P = 0.52, 0.98, 0.77, 0.31, 0.29 for rs2234670, rs2234693, rs9340799, rs4986934, and rs1801132, respectively). Further, genotype frequencies did not differ between those included in multivariate models and those excluded due to missing information about one or more of the included covariates (data not shown). Cancers in women participating via tissue samples were more often of more advanced stage as measured by degree of myometrial invasion (no invasion 2.8 vs 10.5 percent, less than half of the myometrial thickness 65.3 vs 61.7 percent, more than half of the myometrial thickness 26.4 vs 25.9 percent, through the serosa 5.6 vs 1.9 percent, p = 0.05).
Non-participants, i.e. those that were eligible for this study but that either were not selected or did not consent to participate, were on average older than participants. Non-participating cases more often had cancer with more advanced myometrial invasion; 50 percent had cancer invasion to more than half of the myometrial thickness as compared to 31 percent among participants (p < 0.0001). Participants were also more likely to have been oral contraceptives users. Due to the over-sampling scheme, participants were also more likely to have used menopausal hormones. Despite our over-sampling of women with diabetes mellitus, participants were less likely to have this disease than non-participants. Participants and non-participants were similar with regard to smoking habits, family history of endometrial cancer, BMI, age at last birth and age at menopause (data not shown).
Association with endometrial cancer
In univariate analyses conditioned on the selection variables the rs2234670, the rs2234693 and the rs9340799 were associated with endometrial cancer risk so that subjects homozygous for the rarer alleles (denoted 11 in Table 3) were at decreased risk as compared to homozygotes for the most common allele (denoted 00). None of the covariates appeared to be confounders. However, including information about parity, BMI, smoking status and use of oral contraceptives significantly improved the model fit. Age at last birth and age at menopause further bettered model fit but since adjustment left point estimates unchanged these variables were excluded due to missing information in many individuals. Multivariate analyses yielded slightly stronger results compared to those from univariate analyses (Table 3). The difference could not be attributed to any one variable. Heterozygotes were at intermediate risk. Excluding cases with DNA derived from tissue left estimates largely unaltered and associations using cases with tissue derived DNA did not speak against pooling all successfully genotyped cases in the further analysis (Additional file 1, Supplementary table 3).
We examined whether including more than one marker would improve model fit, indicating influence on risk by combinations of alleles at different loci (i.e. haplotype effects). The association with rs9340799 was the strongest. The addition of either rs2234670 or rs2234693 to a model with rs9340799 did not significantly improve explanation of disease status (LR test p = 0.70 and p = 0.50, respectively). Conversely, starting with rs2234670 or rs2234693 and adding rs9340799 significantly improved model fit (LR test p = 0.02 and p = 0.03, respectively). Adding any of the last two markers, rs4986934 or rs1801132, neither improved fit nor changed estimates nor confidence intervals for the first three markers considered one at the time. When adding to the model with rs9340799 also rs2234693 and the interactions between these loci, the fit was not improved (LR test, p = 0.48). Neither did adding the rs1801132 (LR test p = 0.31) or the rs2234670 (LR test p = 0.83) and the respective interactions improve fit. Only three out of the four possible haplotypes of rs9340799 and rs2234693 existed in our study population (Table 4), so that we were unable to compare the model containing rs9340799 and rs2234693 main effects and interactions with the haplotype model. Other haplotype models (under the assumption of multiplicative penetrance), using rs9340799 and one of rs2234670 or rs1801132, did not provide a significant improvement in fit over the single locus (rs9340799) model.
We investigated the association between rs9340799 and endometrial cancer stratified by the covariates (Table 5). We were primarily interested in potential interactions between the marker and use of menopausal hormones and BMI, since these factors are the most important determinants of serum estrogen levels in postmenopausal women, though we also examined other hormonally related variables. We could not estimate main effects of variables that we used in our sampling scheme (see Methods section) and thus the results presented (first in Table 5) from the analyses stratified by these factors convey only genetic effects. In subdividing the cases and controls by use of menopausal hormones, the protective effect of rs9340799 emerged among never users while the pattern was unclear among users, possibly due too small numbers of observations (Table 5). P-values for interaction were 0.007, 0.35, and 0.08 for compounds containing estrogen only, estrogen plus progestins cyclically and estrogen plus progestins continuously, respectively. Despite the small p-values, there were no allele-dose response relations. Association analyses stratified by various other risk factors for endometrial cancer are given in Table 6. BMI did modify the association between rs9340799 and endometrial cancer (p = 0.91). This result was not sensitive to BMI category cutpoints (p = 0.75 and p = 0.45 using BMI>30 or BMI>31 as the highest category, respectively). The protective effect of rs9340799 was clearer among non-smokers and among non-users of combined oral contraceptives although tests for interaction were non-significant (p = 0.17 and p = 0.19, respectively). There appeared to be no modification of the association between rs9340799 and endometrial cancer risk by use of oral low potency estrogen (p = 0.50), diabetes mellitus (p = 0.91), parity (p = 0.88), age at menopause (in tertiles, p = 0.72), age at last birth (in tertiles, p = 0.46), or family history of endometrial cancer (p = 0.51).
Only 45 cases were of non-endometroid histotype. Excluding these from the analyses did not change the results of any of the above analyses (data not shown).
Analyses where cases were stratified by degree of myometrial invasion (Additional file 1, Supplementary table 3) indicated a more pronounced protective effect of the rs9340799 variant allele among those with invasion through more than 50 percent of the myometrium (adjusted OR 0.67 CI 0.45–0.98 for heterozygotes and OR 0.35 CI 0.16–0.75 for variant allele homozygotes). This latter analysis was based on only 139 cases. Stratifying on histopathological grade did not reveal any clear pattern to indicate a differential effect of genotype on risk for endometrial cancers of various degree of differentiation (Additional file 1, Supplementary table 4).
We show that non-coding variation in ESR1 is associated with endometrial cancer risk. With regard to the rs9340799 locus, homozygotes for the rare G allele had an almost halved risk for endometrial cancer compared to common allele homozygotes. This protection seemed more pronounced when only cases with invasion through more than 50 percent of the myometrium were considered. The latter analysis was based on few cases and needs further confirmation. We found no convincing modification of this association by other endometrial cancer risk factors or histopathological grade.
We found no evidence of additional effects of other markers than rs9340799. However, we cannot determine whether this SNP is functional or in LD with some other functional locus. Furthermore, we cannot be completely certain that rs9340799 is more influential than rs2234693 because of the strong LD between these markers, separated by only 46 bp.
The association between rs9340799 and endometrial cancer was not convincingly modified by exposures that alter the availability of estradiol or estrone such as use of menopausal hormone preparations or high BMI. In our parallel breast cancer study  the ESR1-disease risk association was stronger in women with high BMI but this appeared not to be the case in endometrial cancer. The effect of rs9340799 was largely the same in any stratum with sufficient sample size.
Our study is population-based and the largest investigating the relation between ESR1 and endometrial cancer to date. Furthermore, the study population is genetically homogenous, which limits concern about population stratification and which also preserves power as no additional stratifications based on ethnicity need to be made. The vast majority of cases in this study were of the strongly estrogen related endometroid type, and thus constitute a sample in which functional ESR1 alterations are expected to become evident.
The markers that we have studied here were selected because they were known to be polymorphic according to the literature available when the study was planned (around 1996). They were not selected, as is usually the case with current studies, in an attempt to represent the entire gene. This is a weakness but it should not discredit the reported association.
It appears as if the ESR1 gene is composed of three different LD blocks . All markers that are typed in the present study are located within the first region, towards the 5' end of the gene. We did not find evidence of any interactions between markers, which would have indicated that functional combinations of alleles exist in the gene. We initially intended to investigate if we could establish whether any SNP was functional or instead in LD with a functional locus, by comparing a model with genotype main effects and interactions, to a haplotype model, using SNPs rs9340799 and rs2234693. If the haplotype model had exhibited a significantly better fit then there would have been indication that other, linked, un-typed, SNPs were responsible for the association . However, only three out of the four possible haplotypes of rs9340799 and rs2234693 exist in our data, so that we could not make this comparison. In contrast to our findings in the parallel breast cancer study , a haplotype model with rs9340799 and rs1801132 did not provide any additional explanatory information.
Our main results replicate our own findings in a separate but equivalent Swedish study population where we observed associations of similar strength and magnitude . They also corroborate results of a small Japanese study where rs9340799 GG and rs2234693 CC genotypes were associated to a lower risk for endometrial cancer compared to those homozygous for the common allele ; OR 0.26 (CI 0.09–0.79) and OR 0.23 (CI 0.07–0.82), respectively. Another Japanese study found no such association  but a decreased risk with allele-dose effect with the rare allele of a ESR1 codon 10 SNP.
Weel and colleagues found the rs2234693 CC genotype to be associated with an earlier onset of natural menopause , which would entail a lower risk of endometrial cancer, but we did not find any similar association, nor any interaction between ESR1 genotype and age at menopause.
There is no evidence in the literature that rs9340799 affects the amount or function of the estrogen receptor protein. Herrington and colleagues found that the adjacent rs2234693 C allele produced an additional binding site for the myb family of transcription factors . Binding of B-myb to this site appeared to have the capability of enhancing ESR1 transcription. Since B-myb is in itself estrogen responsive there is a possibility of a positive feedback loop that could amplify estrogenic response. It is possible that in fact rs2234693 is the functional variant and that our finding was due to the strong LD between the markers in addition to random variation in the data. Ongoing studies that include dense mapping of ESR1 variants will give guidance as to what variants need to be further functionally validated.
Although our results, which are in line with a priori theories, might indicate a true influence of ESR1 variation on estrogen dependent phenotypes they could conceivably be due to selection bias, a potential problem in case-control studies. Our participation rates, 88 and 76 percent among cases and controls, respectively (67 and 64 percent using those eligible for the parent study in the denominator) could lead to bias if participation were related to genotype. Non-participating cases had a more advanced disease as measured by myometrial invasion. If, for example, a genetic variant is associated to severe cancer but not to less severe cancer and the more severe cases are less likely to participate because they have died, any association with cancer overall would be biased towards the null. Furthermore, genotype frequencies were not different between participants who donated blood and those who took part via a tissue sample, the latter of which were more often deceased. Our results cannot be generalized to include all endometrial cancer because there is data that indicates differing etiologies between endometroid and non-endometroid endometrial cancer .
Another possible mechanism for selection bias is if an ESR1 variant influences a phenotype that differentially affects ascertainment or recruitment probability among cases and controls. In this case one would expect to find case-control associations with multiple otherwise unrelated phenotypes. ESR1 variation has indeed been associated with a host of other outcomes such as coronary heart disease , height , bone mineral density [32, 33], multiple sclerosis , cognitive impairment , and age at menarche . However, despite the diversity among these outcomes there is reasonable evidence of etiological roles for estrogen. Our comparisons of characteristics among participants and non-participants do not point to any plausible mechanism for selection bias as an explanation to the observed association.
In conclusion, we found an association between genetic variation in ESR1 and endometrial cancer risk; the rs9340799 GG genotype was associated to an almost 50 percent decreased risk for endometrial cancer compared to the AA genotype. Our results strengthen belief in this a priori strong candidate gene but we cannot, based on our data, establish the importance of any single genetic locus within the gene.
Lax SF: Molecular genetic pathways in various types of endometrial carcinoma: from a phenotypical to a molecular-based classification. Virchows Arch. 2004, 444 (3): 213-223. 10.1007/s00428-003-0947-3. Epub 2004 Jan 2028.
Lecce G, Meduri G, Ancelin M, Bergeron C, Perrot-Applanat M: Presence of estrogen receptor beta in the human endometrium through the cycle: expression in glandular, stromal, and vascular cells. J Clin Endocrinol Metab. 2001, 86 (3): 1379-1386. 10.1210/jc.86.3.1379.
Matsuzaki S, Fukaya T, Suzuki T, Murakami T, Sasano H, Yajima A: Oestrogen receptor alpha and beta mRNA expression in human endometrium throughout the menstrual cycle. Mol Hum Reprod. 1999, 5 (6): 559-564. 10.1093/molehr/5.6.559.
Iwamoto I, Fujino T, Douchi T, Nagata Y: Association of estrogen receptor alpha and beta3-adrenergic receptor polymorphisms with endometrial cancer. Obstet Gynecol. 2003, 102 (3): 506-511. 10.1016/S0029-7844(03)00578-7.
Sasaki M, Tanaka Y, Kaneuchi M, Sakuragi N, Dahiya R: Polymorphisms of estrogen receptor alpha gene in endometrial cancer. Biochem Biophys Res Commun. 2002, 297 (3): 558-564. 10.1016/S0006-291X(02)02248-9.
Weiderpass E, Persson I, Melhus H, Wedren S, Kindmark A, Baron JA: Estrogen receptor alpha gene polymorphisms and endometrial cancer risk. Carcinogenesis. 2000, 21 (4): 623-627. 10.1093/carcin/21.4.623.
Weiderpass E, Baron JA, Adami HO, Magnusson C, Lindgren A, Bergstrom R, Correia N, Persson I: Low-potency oestrogen and risk of endometrial cancer: a case-control study. Lancet. 1999, 353 (9167): 1824-1828. 10.1016/S0140-6736(98)10233-7.
Weiderpass E, Adami HO, Baron JA, Magnusson C, Bergstrom R, Lindgren A, Correia N, Persson I: Risk of endometrial cancer following estrogen replacement with and without progestins. J Natl Cancer Inst. 1999, 91 (13): 1131-1137. 10.1093/jnci/91.13.1131.
Weiderpass E, Adami HO, Baron JA, Magnusson C, Lindgren A, Persson I: Use of oral contraceptives and endometrial cancer risk (Sweden). Cancer Causes Control. 1999, 10 (4): 277-284. 10.1023/A:1008945721786.
Weiderpass E, Persson I, Adami HO, Magnusson C, Lindgren A, Baron JA: Body size in different periods of life, diabetes mellitus, hypertension, and risk of postmenopausal endometrial cancer (Sweden). Cancer Causes Control. 2000, 11 (2): 185-192. 10.1023/A:1008946825313.
McGuire WL, Chamness GC, Fuqua SA: Abnormal estrogen receptor in clinical breast cancer. J Steroid Biochem Mol Biol. 1992, 43 (1–3): 243-247. 10.1016/0960-0760(92)90214-4.
Roodi N, Bailey LR, Kao WY, Verrier CS, Yee CJ, Dupont WD, Parl FF: Estrogen receptor gene analysis in estrogen receptor-positive and receptor-negative primary breast cancer. J Natl Cancer Inst. 1995, 87 (6): 446-451. 10.1093/jnci/87.6.446.
Sano M, Inoue S, Hosoi T, Ouchi Y, Emi M, Shiraki M, Orimo H: Association of estrogen receptor dinucleotide repeat polymorphism with osteoporosis. Biochem Biophys Res Commun. 1995, 217 (1): 378-383. 10.1006/bbrc.1995.2787.
Yaich L, Dupont WD, Cavener DR, Parl FF: Analysis of the PvuII restriction fragment-length polymorphism and exon structure of the estrogen receptor gene in breast cancer and peripheral blood. Cancer Res. 1992, 52 (1): 77-83.
Isola J, DeVries S, Chu L, Ghazvini S, Waldman F: Analysis of changes in DNA sequence copy number by comparative genomic hybridization in archival paraffin-embedded tumor samples. Am J Pathol. 1994, 145 (6): 1301-1308.
Chen X, Levine L, Kwok PY: Fluorescence polarization in homogeneous nucleic acid analysis. Genome Res. 1999, 9 (5): 492-498.
Tyagi S, Bratu DP, Kramer FR: Multicolor molecular beacons for allele discrimination. Nat Biotechnol. 1998, 16 (1): 49-53. 10.1038/nbt0198-49.
Tapp I, Malmberg L, Rennel E, Wik M, Syvanen AC: Homogeneous scoring of single-nucleotide polymorphisms: comparison of the 5'-nuclease TaqMan assay and Molecular Beacon probes. Biotechniques. 2000, 28 (4): 732-738.
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.
Syvanen AC, Sajantila A, Lukka M: Identification of individuals by analysis of biallelic DNA markers, using PCR and solid-phase minisequencing. Am J Hum Genet. 1993, 52 (1): 46-59.
Holgersson S, Karlsson JA, Kihlgren A, Rosen B, Savolainen P, Gyllensten U: Fluorescent-based typing of the two short tandem repeat loci HUMTH01 and HUMACTBP2: reproducibility of size measurements and genetic variation in the Swedish population. Electrophoresis. 1994, 15 (7): 890-895. 10.1002/elps.11501501127.
Lewontin RC, Feldman MW: A general asymptotic property of two-locus selection models. Theor Popul Biol. 1988, 34 (2): 177-193. 10.1016/0040-5809(88)90041-X.
Wall JD, Pritchard JK: Haplotype blocks and linkage disequilibrium in the human genome. Nat Rev Genet. 2003, 4 (8): 587-597. 10.1038/nrg1123.
Zhao JH, Curtis D, Sham PC: Model-free analysis and permutation tests for allelic associations. Hum Hered. 2000, 50 (2): 133-139. 10.1159/000022901.
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 (4): R437-449. 10.1186/bcr811.
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 (1): R15-10.1186/bcr1861.
Cordell HJ, Clayton DG: A unified stepwise regression procedure for evaluating the relative effects of polymorphisms within a gene using case/control or family data: application to HLA in type 1 diabetes. Am J Hum Genet. 2002, 70 (1): 124-141. 10.1086/338007.
Weel AE, Uitterlinden AG, Westendorp IC, Burger H, Schuit SC, Hofman A, Helmerhorst TJ, van Leeuwen JP, Pols HA: Estrogen receptor polymorphism predicts the onset of natural and surgical menopause. J Clin Endocrinol Metab. 1999, 84 (9): 3146-3150. 10.1210/jc.84.9.3146.
Herrington DM, Howard TD, Brosnihan KB, McDonnell DP, Li X, Hawkins GA, Reboussin DM, Xu J, Zheng SL, Meyers DA, Bleecker ER: Common estrogen receptor polymorphism augments effects of hormone replacement therapy on E-selectin but not C-reactive protein. Circulation. 2002, 105 (16): 1879-1882. 10.1161/01.CIR.0000016173.98826.88.
Schuit SC, Oei HH, Witteman JC, Geurts van Kessel CH, van Meurs JB, Nijhuis RL, van Leeuwen JP, de Jong FH, Zillikens MC, Hofman A, Pols HA, Uitterlinden AG: Estrogen receptor alpha gene polymorphisms and risk of myocardial infarction. Jama. 2004, 291 (24): 2969-2977. 10.1001/jama.291.24.2969.
Schuit SC, van Meurs JB, Bergink AP, Klift van der M, Fang Y, Leusink G, Hofman A, van Leeuwen JP, Uitterlinden AG, Pols HA: Height in pre- and postmenopausal women is influenced by estrogen receptor alpha gene polymorphisms. J Clin Endocrinol Metab. 2004, 89 (1): 303-309. 10.1210/jc.2003-031095.
van Meurs JB, Schuit SC, Weel AE, Klift van der M, Bergink AP, Arp PP, Colin EM, Fang Y, Hofman A, van Duijn CM, van Leeuwen JP, Pols HA, Uitterlinden AG: Association of 5' estrogen receptor alpha gene polymorphisms with bone mineral density, vertebral bone area and fracture risk. Hum Mol Genet. 2003, 12 (14): 1745-1754. 10.1093/hmg/ddg176.
Yamada Y, Ando F, Niino N, Ohta S, Shimokata H: Association of polymorphisms of the estrogen receptor alpha gene with bone mineral density of the femoral neck in elderly Japanese women. J Mol Med. 2002, 80 (7): 452-460. 10.1007/s00109-002-0348-0.
Kikuchi S, Fukazawa T, Niino M, Yabe I, Miyagishi R, Hamada T, Tashiro K: Estrogen receptor gene polymorphism and multiple sclerosis in Japanese patients: interaction with HLA-DRB1*1501 and disease modulation. J Neuroimmunol. 2002, 128 (1–2): 77-81. 10.1016/S0165-5728(02)00140-6.
Yaffe K, Lui LY, Grady D, Stone K, Morin P: Estrogen receptor 1 polymorphisms and risk of cognitive impairment in older women. Biol Psychiatry. 2002, 51 (8): 677-682. 10.1016/S0006-3223(01)01289-6.
Stavrou I, Zois C, Ioannidis JP, Tsatsoulis A: Association of polymorphisms of the oestrogen receptor alpha gene with the age of menarche. Hum Reprod. 2002, 17 (4): 1101-1105. 10.1093/humrep/17.4.1101.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/8/322/prepub
The authors are indebted to all the women who have contributed to this study by answering the questionnaire and donating biological samples. Dr Anders Lindgren performed the pathological reviews. Anna Christensson, Boel Bissmarck and Anders Holmberg provided invaluable assistance in collecting samples from the participants. Maria Branting, Anders Westermark, Erika Svensson, Kristina Larsson, and Birgitta Sundelin provided excellent technical support. We also wish to express our gratitude towards primary health care centers and pathology departments all over Sweden for their unselfish cooperation. Hans-Olov Adami contributed with skills, experience and enthusiasm to all stages of the project. This study was supported by the National Institutes of Health, grant number 5 RO1 CA 77973-03, by the Swedish Cancer Society, and by the K&A Wallenberg foundation (Wallenberg Consortium North). K. Humphreys was supported by a grant from the Wallenberg Consortium North.
The authors declare that they have no competing interests.
SW planned and coordinated the collecting of biological samples, participated in the genotyping, performed most of the biostatistical analyses, interpreted results, and wrote the manuscript. LL performed the bulk of the genotyping and took part in interpreting results and writing the manuscript. KH oversaw all and performed some of the biostatistical analyses, interpreted the results. CM, IP, JB and EW planned the study and interpreted the results. HM, ACS, AK, UL, MLF planned, coordinated and oversaw all DNA extraction and genotyping. FS performed genotyping. All authors critically read and took part in finalizing the manuscript.
Electronic supplementary material
Additional file 1: Supplementary material Wedren ESR1 endometrial cancer. Tables of primers and probes used for genotyping and tables of main results stratified by DNA source, myometrial invasion, and histological grade. (PDF 20 KB)
Authors’ original submitted files for images
About this article
Cite this article
Wedrén, S., Lovmar, L., Humphreys, K. et al. Estrogen receptor alpha gene polymorphism and endometrial cancer risk – a case-control study. BMC Cancer 8, 322 (2008). https://doi.org/10.1186/1471-2407-8-322
- Endometrial Cancer
- Myometrial Invasion
- Endometrial Cancer Patient
- Endometrial Cancer Risk
- Endometrial Cancer Case