Skip to main content
  • Short communication
  • Published:

Failure to detect significant association between estrogen receptor-alpha gene polymorphisms and endometriosis in Japanese women

Abstract

Objectives

The aim of the study was to test whether estrogen receptor 1 (ESR1) gene polymorphisms are correlated with the risk of the development of endometriosis in Japanese women, as a preliminary study.

Methods

To compare allelic frequencies and genotype distributions, a case-control study of 100 affected women and 143 women with no evidence of disease was performed using 10 microsatellite repeat markers and 66 single-nucleotide polymorphisms (SNPs) in the ESR1 gene region.

Results

Although our results might be insufficient to detect genetic susceptibility, owing to the small sample size and low genetic power, statistical analysis of the differences in allelic frequency between the cases and controls at each microsatellite locus demonstrated that no microsatellite locus in the ESR1 gene displayed a significant association with the disease when multiple testing was taken into account. Also, there were no statistically significant differences in the SNP allele frequencies and genotypes between the cases and controls when multiple testing was taken into account.

Conclusion

The findings in our pilot study suggest that ESR1 polymorphisms do not contribute to endometriosis susceptibility.

Introduction

Endometriosis (MIM131200) is a common gynecological disorder and reportedly affects 6–10% of women of reproductive age [1], with substantial annual health costs [2] and health burdens for individuals [3]. The condition, in which endometrial cells are found in sites outside the uterine cavity and which is defined as the occurrence of ectopic steroid hormone-dependent endometrium-like tissue consisting of glands and stroma [4, 5]. While the exact etiology and pathogenesis of endometriosis is unclear, the susceptibility to this disease depends on a complex interaction of environmental and genetic factors. As for genetic factor(s) of susceptibility to endometriosis, the incidence of the disease is approximately seven times higher in relatives of women with endometriosis than in those without a family history [6].

A relationship between exposure to dioxin and the risk of developing endometriosis was previously reported in rhesus monkeys. That study indicated that dioxin could induce the development of endometriosis, with the incidence and severity of the condition being proportional to the dioxin dose [7]. Moreover, Ohtake et al. showed that dioxin-activated arylhydrocarbon receptors (AhRs) and estrogen receptors (ESRs) interacted functionally and the binding of dioxins to AhRs controlled the transcriptional activation of estrogen-responsive elements and estrogenic effects [8]. There is also growing evidence that exposure to environmental contaminants could contribute to the pathogenesis of endometriosis [9, 10].

Estradiol-17β (E2) stimulates uterine and vaginal epithelial proliferation in vivo [11]. The predominant biological effects of estrogen are mediated through two distinct intracellular receptors, ESR α and -β, that are encoded by the ESR1 and ESR2 genes, respectively, which are both expressed in normal endometrium and in endometriotic lesions [12]. E2 elicits its effects via the ESR, which functions as a ligand-activated transcription factor to turn on target genes in E2-responsive tissues. Both the ESR1 and ESR2 genes have been previously disrupted by targeted mutagenesis [13], which resulted in pleiotropic effects with distinct phenotype.

The ESR1 gene, which consists of eight exons that are transcribed into a 6.8-kb mRNA, spans a genomic sequence of over 27 kb located on chromosome 6 at position 6q25.1. ESR1s are members of the nuclear receptor superfamily and act as ligand-regulated transcription factors [14]. Some previous studies have reported on the existence of a positive association between endometriosis and ESR1 polymorphisms and splicing variants [1517]. However, the results of other studies are inconsistent, cannot be replicated, and may vary because of differences between cohorts and populations, diagnostic errors, differences in sample sizes and in choice and type of genetic markers, and errors in laboratory methodology [5, 18, 19].

In order to gain a better understanding of the association between ESR1 polymorphisms and endometriosis, we investigated and genotyped Japanese females with endometriosis (cases) and controls, using microsatellite and single-nucleotide polymorphism (SNP) markers within the ESR1 DNA sequence.

Subjects, materials, and methods

Subjects

All patients and controls were females and we obtained their written informed consent to participate in this study. The patients with endometriosis were diagnosed by laparoscopy and histological examination and the severity of the disease was staged according to the revised classification of the American Society of Reproductive Medicine (r-ASRM) [20]. Of the 100 endometriosis cases (age 38 ± 4.8 years), 27, 10, 31, and 32 cases were classified as stage I, II, III, and IV disease, respectively [21]. The 143 healthy controls (age 43 ± 12.3 years) were patients with unexplained infertility who had received diagnostic laparoscopy and had no evidence of endometriosis (Table 1). Patients with apparent disease(s) were excluded from the control group. Ethics approval for this study was obtained from the Institutional Review Board at Tokai University School of Medicine.

Table 1 Characteristics of endometriosis cases and controls

DNA extraction

Genomic DNA samples were isolated from peripheral blood leukocytes in blood drawn from each subject into tubes containing heparin. The genomic DNA samples were isolated by using a Qiagen DNA extraction Kit (Qiagen, Tokyo, Japan) after lysis with proteinase K and 0.5% sodium dodecyl sulfate (SDS) at 37°C for 1 h.

Analysis of microsatellite polymorphism in the ESR1 gene region

To determine the number of repeat units of the microsatellite loci exhibiting polymorphisms in the ESR gene, unilateral primers were synthesized by labeling at the 5′ end with the fluorescent reagent, 6-FAM (PE Biosystems, Foster City, CA, USA). Thirteen polymerase chain reaction (PCR) primer sets were used for the amplification of 13 microsatellite loci in the ESRα gene region, respectively (Supplementary Table 1). PCR amplification and GenScan (Applied Biosystems, Tokyo, Japan) detection of these microsatellites were carried out as previously described [22].

Genotyping

To examine the distribution of alleles and the genotype frequency of each SNP, we performed direct sequencing analysis of the ESR1 gene using 43 pairs of oligonucleotide PCR primers (Supplementary Table 2). The reaction mixture (20 μl) contained 2 μl of dNTP mixture (each 2.5 mM of dATP, dCTP, dGTP, and dTTP), genomic DNA (5 μl, 2 ng/μl), 2 μl of 10× buffer (100 mM Tris–HCl, pH 8.3, 500 mM KCl, 15 mM MgCl2), 20 pmol of forward and reverse primers, and 0.5 U of rTaq (TaKARa Bio. Inc., Shiga, Japan). PCR amplification was performed in a GeneAmp PCR system 9700 automated thermal cycler (Applied Biosystems, Tokyo, Japan). PCR reaction conditions consisted of initial denaturation for 5 min at 96°C, followed by 30 amplification cycles of 45 s at 96°C, 45 s at various temperatures depending on the primers used (Supplementary Table 2) for annealing, and 2 min at 72°C with a final extension of 7 min at 72°C. Each PCR product was purified by exonuclease I and then sequenced using an ABI 3100 automated sequencer (Applied Biosystems, Tokyo, Japan).

Statistical analyses

Allele and genotype frequencies were determined by direct counting. The significance of differences in the distribution of alleles and genotypes between the patients and controls was tested, using a case-control design, by Fisher’s exact probability test (P value test). A finding of P c < 0.05; that is, a significant correction for multiple testing, using the Bonferroni method, was considered as statistically significant. The odds ratio (OR) and 95% confidence interval (CI) were also calculated for all SNPs. The Haploview program (Dr. Mark Daly's lab, MIT/Harvard Broad Institute) was used to estimate pairwise linkage disequilibrium (LD) and haplotype frequency [23]. The CaTS algorithm (MIT/Harvard Broad Institute, Massachusetts, USA) was used to calculate the genetic power for an association study (http://www.sph.umich.edu/csg/abecasis/CaTS/).

Results

To investigate genetic variations in the ESR1 gene associated with endometriosis, a total of 100 Japanese patients with endometriosis and 143 healthy controls were enrolled for the association analysis, which was carried out using 13 microsatellites found in the region of the ESR1 gene. Three microsatellites with the dinucleotide repeats (CA) n , (AAAC) n , and (GTTT) n , located in introns 4and 6 and the 3′ flanking region, respectively, were not polymorphic (data not shown). As shown in Table 2, none of the microsatellites in the ESR1 gene were significantly associated with endometriosis after the P values were corrected (P c) for multiple testing (10 tests).

Table 2 Results of association tests for ten microsatellite markers in the ESR1 gene

A total of 66 SNPs in the ESR1 gene region were genotyped in the 100 patients with endometriosis and 143 controls. None of the SNPs in the ESR1 gene were significantly (P c > 0.05) associated with endometriosis (Supplementary Table 3) after the P values were corrected for 66 tests. The SNP site (C/G) in exon 4 was located at the nucleotide position 186,086 of the gene and was not involved in an amino acid change, i.e., Pro (rs1801132). The prevalences of the T/T, C/C, and C/C genotypes of the rs1884049, rs1884053, and rs1884054 polymorphisms were significantly higher in the patients than in the control subjects (P = 0.022, 0.014, and 0.028, Table 3) without correction for multiple testing. This significance was lost (P > 0.05) after correcting for 66 statistical tests.

Table 3 Genotype frequencies of the human ESR1 gene and the significance of their association with endometriosis

We evaluated the LD extension of approximately 280 kb of the ESR1 genomic region with 64 SNPs. For the LD block evaluation, we included SNPs with a minor allele frequency of >0.2, genotype success rate of >0.8, and P > 0.001 in the Hardy–Weinberg equilibrium test. Of these 64 SNPs, 58 met the criteria, and the pairwise LD index (∆) was calculated and plotted (Supplementary Fig. 1). We identified five LD blocks with a threshold of ∆ = 0.6. We numbered these blocks 1–5 from ATG to the stop codon in the ESR1 gene. Blocks 1–5 spanned roughly 49.2, 3.0, 87.7, 33.8, and 79,7 kb, respectively. We further analyzed the haplotype constitution with the three SNPs, rs1884049, rs1884053, and rs1884054. Four haplotypes were estimated to have a frequency of >0.01, using Haploview. Haplotype 1 (TCC) was more frequently observed and over-represented in the patients, with a significant difference in frequency between the case and the control group (P c = 0.084, Table 4).

Table 4 Association between ESR1 gene haplotypes and endometriosis

Discussion

Endometriosis is considered to be an estrogen-dependent disorder in women, with estrogen and ESR playing major roles in its pathogenesis, although genetic defects, immunity, and environmental and other hormonal factors might also contribute to the disease [4, 7]. Some studies observed a positive association between endometriosis and ESR1 polymorphisms and splicing variants [1517], whereas others reported equivocal or negative results [10, 19]. Polymorphisms and ER-related genotypes may also determine the function of the sex-steroid system both at the receptor level and at the level of hormone synthesis [24, 25]. Therefore, to elucidate these questions about the risk of endometriosis, we performed an association analysis to determine the relationship between polymorphisms within the ESR1 gene region and endometriosis in Japanese women, using 10 microsatellite repeats and 66 SNPs distributed from the promoter region to the 3′ flanking region. Statistical analysis of the differences in allelic frequencies between the cases and controls at each microsatellite locus demonstrated that two microsatellite loci in introns 4 and 5 of the ESR1 gene displayed a significant association with the disease (P = 0.025 and 0.011, respectively, Table 1). There were also statistically significant differences (uncorrected for multiple testing) between the cases and controls for nine allele frequencies and three genotypes for the SNP markers in intron 4 of the ESR1 gene (Supplementary Table 3). However, when we corrected our results by Bonferroni’s method for multiple testing, we did not obtain any positive markers at a significance level of P c < 0.05. Therefore, we conclude that the microsatellite or SNP allelic genotypic markers that we investigated in the ESR1 gene region are not significantly associated with endometriosis.

Despite our conclusion on the absence of an association between the ESR1 gene polymorphisms and endometriosis, it remains problematic that the statistical corrections using the Bonferroni test may be too stringent for association studies and that we may have missed a biologically important or genetically interesting positive association by creating false-negatives or a type II error with the corrections [26, 27].

For multifactorial or complex diseases, including endometriosis, association studies with a high density of genetic markers have identified many statistically valid associations for the susceptibility genes, as well as several quantitative traits based on the ‘common disease, common variant hypothesis’, i.e., that the etiology of a common disease is mediated by commonly occurring genomic variants in a population. These studies can be divided into two types, that is direct and indirect approach. The former is an association study to analyze a particular functional variant that alters amino acids in coding regions or expression levels in regulatory regions of candidate genes for a disease. On the other hand, in an indirect study, a particular set of genetic markers is investigated to expect the disease locus by LD analysis between disease allele and genetic marker. Of interest, in a multifactorial disease such as endometriosis, it is important to determine the level of the population-attributable risk (PAR) in order to detect the susceptibility genes. As an example, if a causal variant is rare and the disease is common, the PAR and consequently the OR will be low, and the genetic power would be decreased. In a model proposed by Madsen and Browning [28], using 1,000 cases and 1,000 controls, and PAR values of 0.02, 0.1, and 0.25, the power of the test at these PARs was estimated to be 0, 0.68, and 1, respectively. Thus, the detection of susceptibility genes by the indirect approach is dependent not only on the degree of linkage disequilibrium (LD) between the disease variant and the genetic marker but also on the genotype relative risk, whose factors influence the genetic power for an association study. However, it is widely known that allele frequencies and LD strengths vary considerably in different populations and chromosome regions [29]. As for the number of genetic markers used in an association study, a large number of genetic markers with various allele frequencies is needed, because the frequency of the disease variant is unknown. However, if a large number of markers were to be used for an association study, there would be a possibility that the P value may not show a significant association due to the Bonferroni adjustment.

Also, the calculation of the power and sample size required for association studies is important to ensure that the study is sufficiently powered to detect the subtle genetic effects that contribute to most complex diseases. The factors that affect statistical power in the detection of genetic associations are the disease prevalence, effect size (genotype relative risk), allele frequency, and genotype errors [30, 31]. Therefore, we calculated the power in the present association study using the following factors: disease prevalence 0.06–0.10, disease allele frequency 0.2–0.5, and effect size 1.10–1.50. The calculations for the power showed values between 0.07 and 0.68 (Supplementary Table 4). However, in a recent report of a genome-wide association study of endometriosis, the OR for the identified SNP was approximately 1.20 [32]. Given the results of that report, the power in the present study would be 0.14–0.20 (Supplementary Table 4). To increase the power to 0.90, this study would need 1,000 cases and 1,000 controls. In this regard, future follow-up studies using greater sample numbers might help to determine whether the significant differences (P < 0.05) in the frequencies of the microsatellite or SNP markers between the control and disease groups that we identified without the statistical corrections could still be replicated and recognized as genetically significant. Also, it was reported that body habitus, personal habits, and menstrual characteristics were significantly associated with the development of endometriosis [33]. Therefore, for genetic analysis, it might be important to identify the causative factors for endometriosis.

In conclusion, our results suggest that ESR1 gene polymorphisms are not significantly associated with the development of endometriosis in Japanese women. However, this lack of association might be due to our study’s small sample size and low genetic power to detect disease susceptibility. The mechanism by which these polymorphisms might affect the function of the estrogen receptor still needs to be clarified in order to better understand the pathogenesis of endometriosis and other estrogen-dependent diseases.

References

  1. Giudice LC, Kao LC. Endometriosis. Lancet. 2004;364(9447):1789–99.

    Article  PubMed  Google Scholar 

  2. Simoens S, Hummelshoj L, D’Hooghe T. Endometriosis: cost estimates and methodological perspective. Hum Reprod Update. 2007;13(4):395–404.

    Article  PubMed  CAS  Google Scholar 

  3. Jones GL, Kennedy SH, Jenkinson C. Health-related quality of life measurement in women with common benign gynecologic conditions: a systematic review. Am J Obstet Gynecol. 2002;187(2):501–11.

    Article  PubMed  Google Scholar 

  4. Seli E, Arici A. Endometriosis: interaction of immune and endocrine systems. Semin Reprod Med. 2003;21:135–44.

    Article  PubMed  CAS  Google Scholar 

  5. Bulun SE. Endometriosis. N Engl J Med. 2009;360(3):268–79.

    Article  PubMed  CAS  Google Scholar 

  6. Matalliotakis IM, Arici A, Cakmak H, Goumenou AG, Koumantakis G, Mahutte NG. Familial aggregation of endometriosis in the Yale Series. Arch Gynecol Obstet. 2008;278(6):507–11.

    Article  PubMed  Google Scholar 

  7. Rier SE, Martin DC, Bowman RE, Dmowski WP, Becker JL. Endometriosis in rhesus monkey (Macaca mulatta) following chronic exposure to 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin. Fundam Appl Toxicol. 1993;21(4):433–41.

    Article  PubMed  CAS  Google Scholar 

  8. Ohtake F, Takeyama K, Matsumoto T, Kitagawa H, Yamamoto Y, Nohara K, et al. Modulation of oestrogen receptor signalling by association with the activated dioxin receptor. Nature. 2003;423:545–50.

    Article  PubMed  CAS  Google Scholar 

  9. Bruner-Tran KL, Yeaman GR, Crispens MA, Igarashi TM, Osteen KG. Dioxin may promote inflammation-related development of endometriosis. Fertil Steril. 2008;89(5 Suppl):1287–98.

    Article  PubMed  CAS  Google Scholar 

  10. Guo SW, Simsa P, Kyama CM, Mihályi A, Fülöp V, Othman EE, et al. Reassessing the evidence for the link between dioxin and endometriosis: from molecular biology to clinical epidemiology. Mol Hum Reprod. 2009;15(10):609–24.

    Article  PubMed  CAS  Google Scholar 

  11. Galand P, Leroy F, Chrétien J. Effect of oestradiol on cell proliferation and histological changes in the uterus and vagina of mice. J Endocrinol. 1971;49:243–5.

    Article  PubMed  CAS  Google Scholar 

  12. Matsuzaki S, Murakami T, Uehara S, Canis M, Sasano H, Okamura K. Expression of estrogen receptor alpha and beta in peritoneal and ovarian endometriosis. Fertil Steril. 2001;75:1198–205.

    Article  PubMed  CAS  Google Scholar 

  13. Couse JF, Korach KS. Reproductive phenotypes in the estrogen receptor-alpha knockout mouse. Ann Endocrinol. 1999;60:143–8.

    CAS  Google Scholar 

  14. Beato M, Klug J. Steroid hormone receptors: an update. Hum Reprod Update. 2000;6:225–36.

    Article  PubMed  CAS  Google Scholar 

  15. Fujimoto J, Ichigo S, Hirose R, Sakaguchi H, Tamaya T. Expression of estrogen receptor wild type and exon 5 splicing variant mRNAs in normal and endometriotic endometria during the menstrual cycle. Gynecol Endocrinol. 1997;11:11–6.

    Article  PubMed  CAS  Google Scholar 

  16. Georgiou I, Syrrou M, Bouba I, Dalkalitsis N, Paschopoulos M, Navrozoglou I, et al. Association of estrogen receptor gene polymorphisms with endometriosis. Fertil Steril. 1999;72:164–6.

    Article  PubMed  CAS  Google Scholar 

  17. Kitawaki J, Obayashi H, Ishihara H, Koshiba H, Kusuki I, Kado N, et al. Oestrogen receptor-alpha gene polymorphism is associated with endometriosis, adenomyosis and leiomyomata. Hum Reprod. 2001;16:51–5.

    Article  PubMed  CAS  Google Scholar 

  18. Guo SW. Association of endometriosis risk and genetic polymorphisms involving sex steroid biosynthesis and their receptors: a meta-analysis. Gynecol Obstet Invest. 2006;61:90–105.

    Article  PubMed  CAS  Google Scholar 

  19. Tempfer CB, Simoni M, Destenaves B, Fauser BC. Functional genetic polymorphisms and female reproductive disorders:part II—endometriosis. Hum Reprod Update. 2009;15:97–118.

    Article  PubMed  CAS  Google Scholar 

  20. American Society for Reproductive Medicine. Revised American Society for Reproductive Medicine classification of endometriosis: 1996. Fertil Steril. 1997;67:817–21.

    Article  Google Scholar 

  21. American Fertility Society. Revised American Fertility Society classification of endometriosis: 1985. Fertil Steril. 1985;43:351–2.

    Google Scholar 

  22. Matsuzaka Y, Makino S, Okamoto K, Oka A, Tsujimura A, Matsumiya K, et al. Susceptibility locus for non-obstructive azoospermia is localized within the HLA-DR/DQ subregion: primary role of DQB1*0604. Tissue Antigens. 2002;60:53–63.

    Article  PubMed  CAS  Google Scholar 

  23. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5.

    Article  PubMed  CAS  Google Scholar 

  24. Sano M, Inoue S, Hosoi T, Ouchi Y, Emi M, Shiraki M, et al. Association of estrogen receptor dinucleotide repeat polymorphism with osteoporosis. Biochem Biophys Res Commun. 1995;217:378–83.

    Article  PubMed  CAS  Google Scholar 

  25. Zofková I, Zajícková K, Hill M. The estrogen receptor alpha gene determines serum androstenedione levels in postmenopausal women. Steroids. 2002;67:815–9.

    Article  PubMed  Google Scholar 

  26. Sankoh AJ, Huque MF, Dubey SD. Some comments on frequently used multiple endpoint adjustment methods in clinical trials. Stat Med. 1997;6(22):2529–42.

    Article  Google Scholar 

  27. Perneger TV. What’s wrong with Bonferroni adjustments. Br Med J. 1998;316(7139):1236–8.

    Article  CAS  Google Scholar 

  28. Madsen BE, Browning SR. A groupwise association test for rare mutations using a weighted sum statistic. PLoS Genet. 2009;5(2):e1000384.

    Article  PubMed  Google Scholar 

  29. Goddard KA, Hopkins PJ, Hall JM, Witte JS. Linkage disequilibrium and allele-frequency distributions for 114 single-nucleotide polymorphisms in five populations. Am J Hum Genet. 2000;66(1):216–34.

    Article  PubMed  CAS  Google Scholar 

  30. Gordon D, Finch SJ. Factors affecting statistical power in the detection of genetic association. J Clin Invest. 2005;115(6):1408–18.

    Article  PubMed  CAS  Google Scholar 

  31. Fisher SA, Lewis CM. Power of genetic association studies in the presence of linkage disequilibrium and allelic heterogeneity. Hum Hered. 2008;66(4):210–22.

    Article  PubMed  Google Scholar 

  32. Painter JN, Anderson CA, Nyholt DR, Macgregor S, Lin J, Lee SH. Genome-wide association study identifies a locus at 7p15.2 associated with endometriosis. Nat Genet. 2011;43(1):51–4.

    Article  PubMed  CAS  Google Scholar 

  33. Matalliotakis IM, Cakmak H, Fragouli YG, Goumenou AG, Mahutte NG, Arici A. Epidemiological characteristics in women with and without endometriosis in the Yale series. Arch Gynecol Obstet. 2008;277(5):389–93.

    Article  PubMed  Google Scholar 

Download references

Acknowledgments

This study was partly supported by Health Sciences Research Grants from the Ministry of Health, Labor and Welfare of Japan. This work is supported by a Grant from the Ministry of Health, Labor and Welfare of Japan and we thank the members of the Teaching and Research Support Center, Tokai University, for DNA sequencing.

Conflict of interest

The authors declare that they have no conflict of interest.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Yasunari Matsuzaka or Minoru Kimura.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary material 1 (XLS 170 kb)

Rights and permissions

Reprints and permissions

About this article

Cite this article

Matsuzaka, Y., Kikuti, Y.Y., Izumi, Si. et al. Failure to detect significant association between estrogen receptor-alpha gene polymorphisms and endometriosis in Japanese women. Environ Health Prev Med 17, 423–428 (2012). https://doi.org/10.1007/s12199-011-0259-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s12199-011-0259-1

Keywords