Skip to main content

ESR1, HK3 and BRSK1 gene variants are associated with both age at natural menopause and premature ovarian failure

Abstract

Background

Premature ovarian failure (POF) is a complex and heterogeneous disorder that is influenced by multiple genetic components. Numerous candidate gene studies designed to identify POF susceptibility loci have been published, but most positive findings have not been confirmed in follow up studies. We sought to determine if sequence variants previously associated with age at natural menopause (AANM) or early menopause (EM) contribute as well to genetic susceptibility to POF.

Methods

Our study was performed on 371 unrelated idiopathic women with POF and 800 women controls, all Chinese Han. Thirty six SNPs from previous genome-wide association studies (GWAS) responsible for AANM or EM and 3 additional SNPs in ESR1, and 2 additional SNPs in PTHB1 were tested using the Sequenom MassARRAY iPLEX platform for genotyping.

Results

Three SNPs - rs2278493 in HK3, rs2234693 in ESR1 and rs12611091 in BRSK1 - showed nominally significant association with POF. Thus, a plausible relationship could exist between ESR1, BRSK1, HK3 and POF.

Conclusions

This largest association study undertaken to determine correlation between POF and AANM/EM revealed three significant SNPs (rs2278493, rs2234693, and rs12611091). All are associated with not only AAWM and EM but also POF. Insights into shared genetic susceptibility between POF and AANM/EM will provide novel entry points for unraveling genetic mechanism involved in ovarian reserve and oocyte aging processes.

Background

Natural menopause is one of the most important physiological events, the cessation of ovarian function and the end of the reproductive lifespan. Twin and family studies indicate high heritability for age at natural menopause (AANM) (49-87%) [1]. The AANM in Chinese women is approximately 49 years [2]. Early onset of menopause increases risk for many postmenopausal health complications. Menopause occurring between 40 and 44 years has an incidence of 5% and is termed early menopause (EM) [3]. Premature ovarian failure (POF), also termed as primary ovarian insufficiency (POI), is defined as ovarian failure before the age 40 years [4, 5].

POF is characterized by secondary amenorrhea, hypoestrogenism and elevated gonadotropin serum levels (FSH > 40 IU/L). Approximately 1% of the population has POF before age 40 years; only 0.1% or less are affected prior to age 30. Although most cases of POF occur sporadically, approximately 10 - 15% of cases have an affected first-degree relative [6]. Based on either known roles in folliculogenesis or phenotypes of murine knockout models, many candidate genes have been interrogated to determine their role in POF. However, except for FSHR in Finnish women none explains more than 1-2% of POF in any given ethnic group [7]. Genes whose perturbations cause POF have nonetheless been identified in single patients or families, small patient groups or selected ethnic groups. This may reflect genetic heterogeneity, ethnic-specific genetic variation or failure to interrogate sufficient numbers of the candidate genes that could potentially account for POF.

Via genome-wide association studies (GWASs), it is possible to identify common genetic variants contributing to susceptibility to genetically complex or polygenic diseases. Despite the attractiveness of GWAS for POF, genome searches in a sufficiently large number of patients are practically absent [8]. Kang H.J. et al. showed association of two SNPs (rs3884597 and rs6944723) in the PTH-responsive B1 gene (PTHB1) in a small discovery set of 24 women and 24 controls, followed by confirmation in 101 cases and 87 controls [9]. The other GWAS showed association with ADAMTS19 in the discovery set of 99 Dutch women and 181 controls, but not in the replication set (60 POF cases and 90 controls) [10]. That 8% of the relatives of women with POF in Chinese suffer from either early menopause or POF (Table 1) suggests an interrelated genetic mechanism between POF and early menopause. Indeed, Tibiletti et al. proposed that POF and EM patients share similar genetic features and postulated that these conditions may be a variable expression of the same genetic disease [11]. Additionally, FMR1 is a single gene cause of POF, and it is clinically indicated to test for premutation in this gene. Genes like EIF2B5 and FMR1 that are known to be involved in POF have been reported to be associated also with AANM [12–14]. Few studies have explored the shared genetic associations among POF, AANM and EM. On the basis of four independent GWAS, 36 SNPs have been associated with AANM. These are located on chromosomes 2, 3, 5, 6, 7, 8, 9, 11, 12, 13, 15, 18, 19 and 20, respectively, with robust statistical evidence. Of these, 2 SNPs were related to EM [1, 15, 16]. Besides these 36 SNPs, we also assessed correlation between POF and three additional SNPs in the estrogen receptor 1 (ESR1) gene, given the key role of ER-α in folliculogenesis. ESR1 has also been considered a promising candidate gene for POF, with rs2234693, rs9340799 and rs1569788 of ESR1 reported to be associated with idiopathic POF [17–19]. ESR1 also has been one of the most popular candidates related with AANM in previous association studies [20, 21].

Table 1 Clinical features of cases and controls, which are not significantly different (p > 0.05).

Overall, we sought to determine the contribution to POF of the 36 SNPs from AANM GWAS, the 3 SNPs in the ESR1 gene and the 2 SNPs in the PTHB1 gene. Our association study involved 371 women with POF and 800 controls, all Chinese Han.

Methods

Participants and phenotypes

Our study was performed on 371 unrelated idiopathic women with POF and 800 women controls, all Chinese Han, collectively gathered from a number of high volume regional referral centers in China. Samples were collected from 18 provinces during 2002-2010. The age of secondary amenorrhea in our sample was young (mean age 23.85 ± 6.78 years) (Table 1). 8% (32/400) of cases had a positive family history, i.e., at least two first-or second-degree female family members with POF or early menopause including the index POF patient. Most of the affected relatives (27/32; thirteen < 40 yrs and fourteen < 45 yrs) were mothers of the index patient; the remaining five were sisters. Inclusion criteria for POF consisted of cessation of menstrual cycles before 40 years of age, with at least two serum follicle stimulating hormone (FSH) concentrations exceeding 40 IU/L. Serum FSH level was determined by the electrochemiluminescent immunoassay in the analyzer cobas e 601 (Roche, Switzerland). The analytical sensitivity of basal FSH was < 0.10 mIU/ml. All samples were run in the same assay. Early menopause is defined as menopause that occurs between ages 40 and 44. Women with known chromosomal abnormalities, previous ovarian surgery, chemotherapy, or radiotherapy were excluded. Controls attended the clinic because of tubal factor infertility. They were known to be menstruating regularly and confirmed by hormone assays (< 40 mIU/L) and ultrasound imaging to be free of POF-related and early menopause- related phenotypes. Informed consent was obtained from all subjects. The study was approved by the Institutional Review Board of Reproductive Medicine of Shandong University.

Genotyping and statistical analysis

SNPs were chosen for genotyping on the basis of published findings in association studies. A total of 41 SNPs are located in LHCGR, PPARG, SRD5A1, HK3, RAP80, SYCP2L, TNF, ESR1, IGF2R, PTHB1, NBN, TGFBR1, FSHB, PGR, ANKK1, IGF1, CYP19A1, POLG, SMAD7, BRSK1, TMEM150B, MCM8 or intergenic region (Table 2). We used the Sequenom MassARRAY iPLEX platform for genotyping. The Sequenom protocol involves a multiplex PCR reaction prior to a single-base primer extension reaction. The individual SNPs are identified by using matrix-assisted laser desorption/ionization time- of- fight mass spectrometry. SNPs were excluded if the minor allele frequency was < 1%. Ultimately, 36 SNPs were genotyped in the case control cohort. Allelic association analysis was conducted using SHEsis [22]. The genotype distributions for all SNPs were tested for Hardy-Weinberg equilibrium (HWE) using the Fisher's exact test. All SNPs were in Hardy- Weinberg equilibrium. The allelic odds ratio (OR) and 95% confidence interval (CI) were estimated assuming a multiplicative model. P-values were two-tailed. An alpha of 0.05 was used to claim statistical significance.

Table 2 Loci selected for association analysis in POF and controls.

Results

We analyzed 41 SNPs for association with POF in our collection of 400 cases and 800 controls. The analysis was nominally significant for three SNPs: rs2278493 in HK3, rs2234693 in ESR1 and rs12611091 in BRSK1, respectively. Results are summarized in Table 2. These findings demonstrate that certain SNPs associated with AANM also may contribute to POF, but need not necessarily play a major role in POF. No significant association with POF was observed for the remaining 38 SNPs.

Discussion

We conducted the largest association study reported to date in a POF cohort in order to test the contribution to POF of 36 SNPs, specifically known to be associated with AANM and EM. We also investigated the associations of 3 ESR1 SNPs and 2 PTHB1 SNPs in POF. We found contribution to POF for SNPs in ESR1, BRSK1 and HK3, 2 of these genes not previously considered as candidates for POF.

Estrogen regulates cyclic gonadotropin release at the hypothalamus- hypophysis- ovarian (HPO) axis by the estrogen receptor α (ER-α), which is encoded by ESR1 and enhances follicular development by the ER-β encoded by ESR2. Female ER-α knockout mice show anovulation and complete infertility. Rs2234693 lies in an intron of the ESR1gene and has been reported to be related to POF and the onset of natural menopause in Korean and Dutch [17–19], consistent with our results. Because this SNP does not change any amino acid in the ERS1 protein, it is presumed to be in linkage disequilibrium with a regulatory sequence(s) that could affect gene expression or function and consequently result in an alteration in estrogenic biological activity.

Rs12611091 is located in the intron region of BRSK1 (BR serine/threonine kinase 1), a gene that codes for an AMP-activated protein kinase (AMPK)-related kinase [23]. BRSK1 is highly expressed in human brain and moderately expressed in mammalian ovaries [24]. Mutation or overexpression of the BRSK1 gene deleteriously affects vesicle transport and release at the axon terminals [25]. BRSK1 might affect the secretion of gonadotropin- releasing hormone (GnRH) from the hypothalamus through this mechanism. Additionally, maternal embryonic leucine zipper kinase (MELK), which is highly expressed in spermatogonia and oocytes, is one of the downstream targets of BRSK1 [26, 27]. Plausible involvement in ovarian aging is therefore suggested for BRSK1.

Hexokinase 3 (white cell) (HK3), one of the four hexokinase family members, is involved in carbohydrate metabolism through phosphorylating glucose to produce glucose-6-phosphate [28]. HK3 is present in uterus, placenta, lung and adipose according to GeneAtlas [24]. Literature analysis for HK3 did not indicate an immediate functional explanation for the observed association. To elucidate the precise molecular mechanism, extensive sequence analysis and functional studies of the HK3gene variants are needed.

Idiopathic POF and EM differ in age of menopause onset and have long been considered to represent variable expression of the same genetic disease [11]. It is highly plausible that certain genes might contribute to both. That is, similar underlying genetic predisposition exists, possibly with different environmental factors or modifiers triggering disease progression and, hence, specific age of onset. Another possibility is that certain genes in pathways ordinarily predisposing only to EM interact or converge to cause even earlier age of onset, i.e, below age 40 and thus lead to clinical POF.

Conclusions

In summary, our data demonstrate that genetic susceptibility for both AANM/EM and POF is shared by the ESR1, HK3 and BRSK1 genes. This became evident by our showing that three SNPs previously related to AANM were associated with POF. At these newly identified loci fine mapping or sequencing might lead to identification of causal variants, and thus expand our knowledge of the underlying mechanism of POF. Insights into shared genetic susceptibility between POF and AANM/EM will further provide novel entry points for unraveling etiology involved in ovarian reserve and aging process. Identifying genetic causation may further be of diagnostic value to other family members, especially female offspring. If a younger relative has the same genetic variant, timeline for achieving pregnancy might be advanced and early child-bearing encouraged.

Abbreviations

AANM:

age at natural menopause

EM:

early menopause

POF:

premature ovarian failure

POI:

primary ovarian insufficiency

GWAS:

genome-wide association study

PTHB1:

PTH-responsive B1 gene

ESR1:

estrogen receptor 1

FSH:

follicle stimulating hormone

SNP:

single nucleotide polymorphism

HPO:

hypothalamus- hypophysis- ovarian

ER-α:

estrogen receptor α

AMPK:

AMP-activated protein kinase

MELK:

maternal embryonic leucine zipper kinase

BRSK1:

BR serine/threonine kinase 1

GnRH:

gonadotropin- releasing hormone

HK3:

Hexokinase 3 (white cell).

References

  1. Stolk L, Zhai G, van Meurs JB, Verbiest MM, Visser JA, Estrada K, Rivadeneira F, Williams FM, Cherkas L, Deloukas P, Soranzo N, de Keyzer JJ, Pop VJ, Lips P, Lebrun CE, van der Schouw YT, Grobbee DE, Witteman J, Hofman A, Pols HA, Laven JS, Spector TD, Uitterlinden AG: Loci at chromosomes 13, 19 and 20 influence age at natural menopause. Nat Genet. 2009, 41: 645-647. 10.1038/ng.387.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Nie GN, Wang XY, Yang HY, Aihua Ou: The investigation and analysis of the factors related with the menopausal age of urban women in China. China maternal and child healthcare (Chinsese). 2011, 8: 1191-1193.

    Google Scholar 

  3. Chang SH, Kim CS, Lee KS, Kim H, Yim SV, Lim YJ, Park SK: Premenopausal factors influencing premature ovarian failure and early menopause. Maturitas. 2007, 58: 19-30. 10.1016/j.maturitas.2007.04.001.

    Article  PubMed  Google Scholar 

  4. De Vos M, Devroey P, Fauser BC: Primary ovarian insufficiency. Lancet. 2010, 376: 911-921. 10.1016/S0140-6736(10)60355-8.

    Article  PubMed  Google Scholar 

  5. Nelson LM: Clinical practice. Primary ovarian insufficiency. N Engl J Med. 2009, 360: 606-614. 10.1056/NEJMcp0808697.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Simpson JL: Genetic and phenotypic heterogeneity in ovarian failure: overview of selected candidate genes. Ann N Y Acad Sci. 2008, 1135: 146-154. 10.1196/annals.1429.019.

    Article  PubMed  Google Scholar 

  7. Aittomaki K, Lucena JL, Pakarinen P, Sistonen P, Tapanainen J, Gromoll J, Kaskikari R, Sankila EM, Lehvaslaiho H, Engel AR, Nieschlag E, Huhtaniemi I, de la Chapelle A: Mutation in the follicle-stimulating hormone receptor gene causes hereditary hypergonadotropic ovarian failure. Cell. 1995, 82: 959-968. 10.1016/0092-8674(95)90275-9.

    Article  CAS  PubMed  Google Scholar 

  8. Christin-Maitre S, Tachdjian G: Genome-wide association study and premature ovarian failure. Ann Endocrinol (Paris). 2010, 71: 218-221. 10.1016/j.ando.2010.02.014.

    Article  CAS  Google Scholar 

  9. Kang H, Lee SK, Kim MH, Song J, Bae SJ, Kim NK, Lee SH, Kwack K: Parathyroid hormone-responsive B1 gene is associated with premature ovarian failure. Hum Reprod. 2008, 23: 1457-1465. 10.1093/humrep/den086.

    Article  CAS  PubMed  Google Scholar 

  10. Knauff EA, Franke L, van Es MA, van den Berg LH, van der Schouw YT, Laven JS, Lambalk CB, Hoek A, Goverde AJ, Christin-Maitre S, Hsueh AJ, Wijmenga C, Fauser BC: Genome-wide association study in premature ovarian failure patients suggests ADAMTS19 as a possible candidate gene. Hum Reprod. 2009, 24: 2372-2378. 10.1093/humrep/dep197.

    Article  CAS  PubMed  Google Scholar 

  11. Tibiletti MG, Testa G, Vegetti W, Alagna F, Taborelli M, Dalpra L, Bolis PF, Crosignani PG: The idiopathic forms of premature menopause and early menopause show the same genetic pattern. Hum Reprod. 1999, 14: 2731-2734. 10.1093/humrep/14.11.2731.

    Article  CAS  PubMed  Google Scholar 

  12. Fogli A, Gauthier-Barichard F, Schiffmann R, Vanderhoof VH, Bakalov VK, Nelson LM, Boespflug-Tanguy O: Screening for known mutations in EIF2B genes in a large panel of patients with premature ovarian failure. BMC Womens Health. 2004, 4: 8-10.1186/1472-6874-4-8.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Conway GS: Premature ovarian failure and FMR1 gene mutations: an update. Ann Endocrinol (Paris). 2010, 71: 215-217. 10.1016/j.ando.2010.02.009.

    Article  CAS  Google Scholar 

  14. He C, Kraft P, Chasman DI, Buring JE, Chen C, Hankinson SE, Pare G, Chanock S, Ridker PM, Hunter DJ: A large-scale candidate gene association study of age at menarche and age at natural menopause. Hum Genet. 2010, 128: 515-527. 10.1007/s00439-010-0878-4.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Murray A, Bennett CE, Perry JR, Weedon MN, Jacobs PA, Morris DH, Orr N, Schoemaker MJ, Jones M, Ashworth A, Swerdlow AJ: Common genetic variants are significant risk factors for early menopause: results from the Breakthrough Generations Study. Hum Mol Genet. 2011, 20: 186-192. 10.1093/hmg/ddq417.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. He C, Kraft P, Chen C, Buring JE, Pare G, Hankinson SE, Chanock SJ, Ridker PM, Hunter DJ, Chasman DI: Genome-wide association studies identify loci associated with age at menarche and age at natural menopause. Nat Genet. 2009, 41: 724-728. 10.1038/ng.385.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Kim S, Pyun JA, Cha DH, Ko JJ, Kwack K: Epistasis between FSHR and CYP19A1 polymorphisms is associated with premature ovarian failure. Fertil Steril. 2011

    Google Scholar 

  18. Yoon SH, Choi YM, Hong MA, Lee GH, Kim JJ, Im HJ, Min EG, Kang BM, Yoon BK, Moon SY: Estrogen receptor {alpha} gene polymorphisms in patients with idiopathic premature ovarian failure. Hum Reprod. 2010, 25: 283-287. 10.1093/humrep/dep375.

    Article  CAS  PubMed  Google Scholar 

  19. Bretherick KL, Hanna CW, Currie LM, Fluker MR, Hammond GL, Robinson WP: Estrogen receptor alpha gene polymorphisms are associated with idiopathic premature ovarian failure. Fertil Steril. 2008, 89: 318-324. 10.1016/j.fertnstert.2007.03.008.

    Article  CAS  PubMed  Google Scholar 

  20. Yang JJ, Cho LY, Lim YJ, Ko KP, Lee KS, Kim H, Yim SV, Chang SH, Park SK: Estrogen receptor-1 genetic polymorphisms for the risk of premature ovarian failure and early menopause. J Womens Health (Larchmt). 2010, 19: 297-304. 10.1089/jwh.2008.1317.

    Article  Google Scholar 

  21. 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: 3146-3150. 10.1210/jc.84.9.3146.

    CAS  PubMed  Google Scholar 

  22. Shi YY, He L: SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005, 15: 97-98. 10.1038/sj.cr.7290272.

    Article  CAS  PubMed  Google Scholar 

  23. Barnes AP, Lilley BN, Pan YA, Plummer LJ, Powell AW, Raines AN, Sanes JR, Polleux F: LKB1 and SAD kinases define a pathway required for the polarization of cortical neurons. Cell. 2007, 129: 549-563. 10.1016/j.cell.2007.03.025.

    Article  CAS  PubMed  Google Scholar 

  24. GeneAtlas. [http://genatlas.medecine.univ-paris5.fr/]

  25. Lu R, Niida H, Nakanishi M: Human SAD1 kinase is involved in UV-induced DNA damage checkPOFnt function. J Biol Chem. 2004, 279: 31164-31170. 10.1074/jbc.M404728200.

    Article  CAS  PubMed  Google Scholar 

  26. Kishi M, Pan YA, Crump JG, Sanes JR: Mammalian SAD kinases are required for neuronal polarization. Science. 2005, 307: 929-932. 10.1126/science.1107403.

    Article  CAS  PubMed  Google Scholar 

  27. Thelie A, Papillier P, Pennetier S, Perreau C, Traverso JM, Uzbekova S, Mermillod P, Joly C, Humblot P, Dalbies-Tran R: Differential regulation of abundance and deadenylation of maternal transcripts during bovine oocyte maturation in vitro and in vivo. BMC Dev Biol. 2007, 7: 125-10.1186/1471-213X-7-125.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Nakamura N, Shibata H, O'Brien DA, Mori C, Eddy EM: Spermatogenic cell-specific type 1 hexokinase is the predominant hexokinase in sperm. Mol Reprod Dev. 2008, 75: 632-640. 10.1002/mrd.20791.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank all the participants involved in this study. Supported by grants from the National Basic Research Program of China (973 Program-2012CB944700, 2010CB945002, 2007CB947403, 2007CB948103), A Foundation for the Author of National Excellent Doctoral Dissertation of PR China (201078), the National Natural Science Foundation of China (30973170, 81000236, 81000238), Science research foundation item of no-earnings health vocation from Ministry of Public Health (201002013).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Zi-Jiang Chen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

ZJC critically revised the manuscript, contributed to the study design, analysis, and supervised the research. YYQ drafted and critically revised the manuscript, performed the statistical analyses, and contributed to the interpretation of the data. MS, DYW, XYL, BZ and HJ coordinated clinic diagnosis, the blood collection and storage of cases and controls. YL and JLS performed genotyping and statistical analyses. JFX contributed to the bioinformatics and the interpretation of the data. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Qin, Y., Sun, M., You, L. et al. ESR1, HK3 and BRSK1 gene variants are associated with both age at natural menopause and premature ovarian failure. Orphanet J Rare Dis 7, 5 (2012). https://doi.org/10.1186/1750-1172-7-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1750-1172-7-5

Keywords