- Research article
- Open Access
Gene polymorphisms in RANKL/RANK/OPG pathway are associated with ages at menarche and natural menopause in Chinese women
BMC Women's Health volume 15, Article number: 32 (2015)
Age at menarche (AAM) and age at natural menopause (AANM) have been shown intimately associated with woman’s health later in life. Previous studies have indicated that AAM and AANM are highly heritable. RANKL/RANK/OPG signaling pathway is essential for mammary gland development, which is also found associated with post-menopausal and hormone-related diseases. The aim of this study was to evaluate associations between the polymorphisms in the TNFSF11, TNFRSF11A and TNFRSF11B genes in the RANKL/RANK/OPG pathway with AAM and AANM in Chinese women.
Post-menopausal Chinese women (n = 845) aged from 42 to 89 years were recruited in the study. Information about AAM and AANM were obtained through questionnaires and the genomic DNA was isolated from peripheral blood from the participants. Total 21 tagging single nucleotide polymorphisms (SNPs) of TNFSF11, TNFRSF11A and TNFRSF11B were genotyped.
Three SNPs of TNFRSF11A (rs4500848, rs6567270 and rs1805034) showed significant association with AAM (P < 0.01, P = 0.02 and P = 0.01, respectively), and one SNP (rs9962159) was significantly associated with AANM (P = 0.03). Haplotypes TC and AT (rs6567270-rs1805034) of TNFRSF11A were found to be significantly associated with AAM (P = 0.01 and P = 0.02, respectively), and haplotypes GC and AC (rs9962159-rs4603673) of TNFRSF11A showed significant association with AANM (P = 0.03 and P < 0.01, respectively). No significant association between TNFSF11 or TNFRSF11B gene with AAM or AANM was found.
The present study suggests that TNFRSF11A but not TNFSF11 and TNFRSF11B genetic polymorphisms are associated with AAM and AANM in Chinese women. The findings provide evidence that genetic variations in RANKL/RANK/OPG pathway may be associated with the onset and cessation of the menstruation cycle.
Age at menarche (AAM) and age at natural menopause (AANM) have been shown intimately associated with woman’s health later in life. Women with early menarche have high risks of breast cancer , ovarian cancer , type 2 diabetes  or metabolic syndrome , whereas late menarche can increase the risk of osteoporosis . On the other hand, early AANM is associated with increased risk of cardiovascular diseases  and osteoporosis . Recent data have shown that AAM and AANM were associated with all-cause mortality . However, the factors that affect AAM and AANM are not entirely clear.
AAM and AANM are complex traits which are influenced by both genetic and environmental factors and their interactions . Twin and familial studies have indicated that AAM and AANM are highly heritable, ranging from 45% to 74% for AAM [10,11] and from 49% to 87% for AANM [12,13]. Genes involved in hormone biosynthesis and metabolic pathways were found to be associated with AAM and AANM [14,15], however, no specific genes have been identified yet.
The receptor activator of nuclear factor-kappa B ligand (RANKL), its receptor RANK and the decoy receptor osteoprotegerin (OPG) belong to the tumor necrosis factor superfamily and they are encoded by genes TNFSF11, TNFRSF11A and TNFRSF11B, respectively. RANKL/RANK/OPG signaling pathway plays important roles in bone modeling and remodeling , cell death and proliferation, inflammation, and immunity [17,18]. RANKL/RANK/OPG pathway is also found associated with post-menopausal and hormone-related diseases, such as osteoporosis  and reproductive cancer . Furthermore, RANKL is found to be essential for mammary gland development in mice by promoting proliferation and maintaining survival of mammary epithelial cells . Mammary gland changes are one of the hallmarks during menarche and menopause [22,23]. Therefore, RANKL/RANK/OPG pathway may involve in modulating the onset and cessation of the menstrual cycle. The present study investigated the associations of single nucleotide polymorphisms and haplotypes in TNFSF11, TNFRSF11A and TNFRSF11B genes in RANKL/RANK/OPG pathway with AAM and AANM in Chinese females.
A total 1026 post-menopausal women from ten community centers in Nanchang from December 2011 to December 2012 were enrolled in the study. All the participants were from Han Chinese ethnic group. Age at interview, AAM, AANM, detailed medical history, birth history (number of live delivery), and abortion information (number of abortions) were obtained through a self-designed questionnaire, all the information collected in the study was self-reported. AAM was defined as the age at the first menstrual period. AANM was defined as one year without menstruation after the age at the last menstrual period. For each participant, height (cm) and weight (kg) were measured. The body mass index (BMI) was calculated as weight/height2.
All of the participants were subjected to blood counts, liver and kidney function tests, fasting plasma glucose tests. Subjects included in the study had normal blood counts, normal liver and kidney functions and blood glucose levels. Subjects were excluded from the study if they suffered from diseases and surgeries that could affect menstruation, such as severe chronic diseases, rheumatic diseases (e.g. systemic lupus erythematosus, rheumatoid arthritis), severe endocrine and metabolic diseases (e.g. diabetes, hyperparathyroidism, pituitary or adrenal diseases), malabsorption diseases (e.g. chronic diarrhea, anorexia nervosa), cancer, and uterine or ovarian resection. Participants who had taken glucocorticosteroid or sex hormone within the past 3 months were also excluded. Finally, 845 subjects were included in the study. The study was approved by the Ethics Committee of The Third Hospital of Nanchang. Written informed consent was obtained from every participant.
Tagging SNPs of the three genes were selected from the software program Haploview version 4.2  (http://www.broad.mit.edu/mpg/haploview/) with minor allele frequencies (MAF) > 10% in the Chinese Han population in HapMap (http://www.hapmap.org/), and the pairwise linkage disequilibrium (LD) was greater than a threshold of r2 (r2 = 0.8). In addition to, SNPs reported in previous studies or potentially functional SNPs in three candidate genes were forced into the SNP selection process. Finally, a total of 21 SNPs were selected in three genes (9 in TNFRSF11A gene, 6 in TNFSF11 gene, and 6 in TNFRSF11B gene). Of these, 18 SNPs are located in the introns of the three genes, one in 5′-UTR, two in the exonic region. All of these SNPs were authenticated using the NCBI (http://www.ncbi.nlm.nih.gov/SNP/) and HapMap databases.
Approximately 5 mL of venous blood was collected from all of the participants after a minimum of 10 h fasting and stored in tubes containing 100 μL of 10% ethylene diaminetetraacetic acid (EDTA). Genomic DNA was extracted from whole blood samples using the QIAamp DNA Mini Kit (Qiagen Inc., Hilden, Germany). DNA samples concentration and quality were detected spectrophotometrically at 260/280 nm and stored at −80°C until analysed. Genotyping was performed using the high-throughput Sequenom genotyping platform (MassARRAY MALDI-TOF MS system, Sequenom Inc., San Diego, CA). For quality control, 5% of the samples were repeatedly genotyped, and the results were found to be 100% concordant.
Genotype frequencies and concordance of the SNPs were analyzed for the Hardy-Weinberg equilibrium (HWE) using the χ2 test. Data were expressed as mean ± standard deviation. The stepwise multiple regression analysis was used to analyze the relationships between the SNPs and AAM and AANM, subsequently, each SNP with different genotype was analyzed independently using one-way univariate analysis of variance (ANOVA), BMI, age at interview, number of deliveries and abortions were considered as covariates and were adjusted during analysis. Bonferroni correction was used to adjust the P values for multiple comparisons. The statistical analyses were performed using SPSS version 13.0 for Windows (SPSS Inc., Chicago, IL, USA). The linkage disequilibrium structure and allele frequencies were examined using Haploview 4.2 software . The significance of each haplotype within the defined blocks was analyzed by PLINK software  (http://pngu.mgh.harvard.edu/~purcell/plink/). All analyses were two-tailed, and P -value < 0.05 was considered statistically significant.
Characteristics of the study participants
The basic characteristics of the 845 participants aged from 42 to 89 years were shown in Table 1. The mean age at interview was 60.88 ± 8.72 years, the mean AAM was 14.97 ± 2.00 years and AANM was 48.77 ± 4.16 years. No statistically significant association was observed between AAM and AANM (P = 0.15).
SNP genotyping and linkage disequilibrium
The basic characteristics of the SNPs are listed in Table 2. All study SNPs had a minor allele frequency of at least 0.1 and were in agreement with Hardy-Weinberg equilibrium (P > 0.05). Linkage disequilibrium between alleles at polymorphic loci was shown in Figure 1. Four haplotype blocks and seventeen of the most common haplotypes (frequency > 5%) were further analyzed for the association of haplotype with AAM and AANM.
Association analyses of the SNP and haplotypes with AAM and AANM
Three SNPs in TNFRSF11A, i.e. rs4500848, rs6567270 and rs1805034, showed significant association with AAM (P < 0.01, P = 0.02 and P = 0.01, respectively), whereas only rs9962159 in TNFRSF11A was significantly associated with AANM (P = 0.03) (Table 2). After correction of age at interview, BMI, number of deliveries and abortions, the associations between those SNPs and AAM or AANM were found significant. After the Bonferroni correction, the rs4500848 was still significantly associated with AAM (P = 0.04), however, the associations between the others SNPs with AAM or AANM were no longer statistically significant (Table 2). Individuals with the T/T genotype of SNP rs4500848 had an earlier onset of menarche by 0.59 years than did those with the C/C genotype. Likewise, women with the G/G genotype of SNP rs9962159 had an earlier menopause by 0.79 years than those with the A/A genotype (Table 3).
Two haplotypes (TC and AT) of block rs6567270-rs1805034 of TNFRSF11A were found significantly associated with AAM (P = 0.01 and P = 0.02, respectively) (Table 4). Haplotypes GC and AC of block rs9962159-rs4603673 of TNFRSF11A were significantly associated with AANM (P = 0.03 and P < 0.01, respectively). Haplotype TAGCGT of block rs9525641-rs2277439-rs2324851-rs2875459-rs2200287-rs9533166 of TNFSF11 showed marginally significant association with AAM (P = 0.06). Notably, all the significantly associated SNPs and haplotypes were observed in TNFRSF11A. SNPs and haplotypes in TNFSF11 and TNFRSF11B genes did not show significant association with either AAM or AANM.
According to the previous studies, there was a direct relationship between AAM and AANM, women with earlier menarche had earlier menopause in Poland . However, other studies had reported no association between AAM and AANM . In this tudy, no statistically significant association was observed between AAM and AANM. The present study revealed that three SNPs (rs4500848, rs6567270 and rs1805034) and two haplotypes of TNFRSF11A showed significant association with AAM in Chinese women. These findings are in line to a previous report by Pan et al. , the authors found five SNPs (rs7239261, rs8094884, rs3826620, rs8089829, and rs9956850) and seven haplotypes of TNFRSF11A significantly associated with AAM in Chinese women. Thus, polymorphisms in TNFRSF11A are highly associated with AAM in Chinese women. In contrast to TNFRSF11A, SNPs of TNFSF11 did not show association with AAM in our study. Noticeably, two SNPs (rs9525641 and rs2200287) of TNFSF11 displayed a strong association with AAM in white women , but no significant association was observed between the two SNPs and AAM in our study. The inconsistency between the results of the present study and the other  may due to the different ethnic populations used, different sample sizes and statistical approaches. The inconsistent results were also observed in the association of TNFSF11 gene polymorphisms and AANM. In the present study, no significant association between polymorphisms of TNFSF11 and AANM was found in Chinese women, however, such relationship was reported in white women, two SNPs (rs346578 and rs9525641) of TNFSF11 showed association with AANM . Regardless of the discrepancy in TNFSF11, we and the others  both found a strong association between polymorphisms of TNFRSF11A and AANM.
The associations between polymorphisms of TNFRSF11A with AAM and AANM found in the present study can be explained by its possible roles in mammary gland development and menstruation. First, TNFRSF11A belongs to RANKL/RANK/OPG signaling pathway. RANKL plays an important role in mammary gland development, indicating its potential role in regulating or responding to sex hormone fluctuation and subsequently influencing menstrual cycles. Studies have shown that gonadotropin-releasing hormone (GnRH) can modulate RANKL expression in breast cancer cells , and expressions of RANK and RANKL in different cell lines are controlled by estrogen , follicle-stimulating hormone , and dehydroepiandrosterone . Estrogen is also found to regulate gene expression and ratio of the RANKL/OPG . Second, RANKL signaling pathway can stimulate ductal side-branching and alveologenesis in the mammary gland in mouse . RANKL can be induced in mammary epithelium and can regulate the proliferation of cells . Therefore, RANK signaling pathway may influence onset of puberty and menstrual cycle by regulating mammary gland development. Third, genome-wide association (GWA) studies have identified some novel genetic loci associated with AAM and AANM [37,38]. Gene set enrichment pathway analyses using the GWA dataset found that nuclear factor-kappa B (NF-κB) signaling pathway may be associated with timing of menopause . Recent studies have revealed that NF-κB pathway plays an important role in mammary ductal morphogenesis , and ovarian cell function in animals . It was well established that the RANKL/RANK/OPG pathway can activate NF-κB and its downstream players . Thus, genes (e.g. TNFRSF11A) in the RANKL/RANK/OPG pathway may have role in the onset and cessation of the menstruation cycle.
The present study has some limitations. First, beside genetic other factors can influence timing of menarche and menopause, e.g. environment and socioeconomic status. We studied only the relationship between genetic variations and AAM and AANM. Furthermore, gene-environment interactions may also play a role in causing variation in the AAM and AANM. Second, The data of AAM and AANM were collected through retrospective self-report, which may cause recall bias. The participants in the study were aged from 42 to 89 years, with long interval periods, which might potentially incur recall error. It is reported that the accuracy of long-term recall of AAM and AANM varied from 70% to 84% [43,44]. In fact, it was found that some participants could not remember the exact age at the first menstrual period, and those subjects were excluded from the study. Large-scale studies are needed to confirm current findings, and the precise mechanisms underlying the observed associations in our study remain to be determined.
The present study, for the first time, demonstrated that TNFRSF11A but not TNFSF11 and TNFRSF11B genetic polymorphisms are associated with AAM and AANM in Chinese women. The findings provide evidence that genetic variations in RANKL/RANK/OPG pathway may be associated with the onset and cessation of the menstruation cycle.
Age at menarche
Age at natural menopause
Receptor activator of nuclear factor-kappa B
Receptor activator of nuclear factor-kappa B ligand
Single nucleotide polymorphisms
Body mass index
Nuclear factor-kappa B
Ritte R, Lukanova A, Tjønneland A, Olsen A, Overvad K, Mesrine S, et al. Height, age at menarche and risk of hormone receptor -positive and -negative breast cancer: a cohort study. Int J Cancer. 2013;132:2619–29.
Gong TT, Wu QJ, Vogtmann E, Lin B, Wang YL. Age at menarche and risk of ovarian cancer: a meta-analysis of epidemiological studies. Int J Cancer. 2013;132:2894–900.
Stöckl D, Döring A, Peters A, Thorand B, Heier M, Huth C, et al. Age at menarche is associated with prediabetes and diabetes in women (aged 32–81 years) from the general population: the KORA F4 study. Diabetologia. 2012;55:681–8.
Akter S, Jesmin S, Islam M, Sultana SN, Okazaki O, Hiroe M, et al. Association of age at menarche with metabolic syndrome and its components in rural Bangladeshi women. Nutr Metab (Lond). 2012;9:99.
Eastell R. Role of oestrogen in the regulation of bone turnover at the menarche. J Endocrinol. 2005;185:223–34.
Cui R, Iso H, Toyoshima H, Date C, Yamamoto A, Kikuchi S, et al. Relationships of age at menarche and menopause, and reproductive year with mortality from cardiovascular disease in Japanese postmenopausal women: the JACC study. J Epidemiol. 2006;16:177–84.
Sioka C, Fotopoulos A, Georgiou A, Xourgia X, Papadopoulos A, Kalef-Ezra JA. Age at menarche, age at menopause and duration of fertility as risk factors for osteoporosis. Climacteric. 2010;13:63–71.
Li S, Rosenberg L, Wise LA, Boggs DA, LaValley M, Palmer JR. Age at natural menopause in relation to all-cause and cause-specific mortality in a follow-up study of US black women. Maturitas. 2013;75:246–52.
Dvornyk V. Waqar-ul-Haq: genetics of age at menarche: a systematic review. Hum Reprod Update. 2012;18:198–210.
Snieder H, MacGregor AJ, Spector TD. Genes control the cessation of a woman’s reproductive life: a twin study of hysterectomy and age at menopause. J Clin Endocrinol Metab. 1998;83:1875–80.
Sharma K. Genetic basis of human female pelvic morphology: a twin study. Am J Phys Anthropol. 2002;117:327–33.
Murabito JM, Yang Q, Fox C, Wilson PW, Cupples LA. Heritability of age at natural menopause in the Framingham heart study. J Clin Endocrinol Metab. 2005;90:3427–30.
de Bruin JP, Bovenhuis H, van Noord PA, Pearson PL, van Arendonk JA, te Velde ER, et al. The role of genetic factors in age at natural menopause. Hum Reprod. 2001;16:2014–8.
He C, Kraft P, Chen C, Buring JE, Paré G, Hankinson SE, et al. Genome-wide association studies identify loci associated with age at menarche and age at natural menopause. Nat Genet. 2009;41:724–8.
Perry JR, Hsu YH, Chasman DI, Johnson AD, Elks C, Albrecht E, et al. DNA mismatch repair gene MSH6 implicated in determining age at natural menopause. Hum Mol Genet. 2014;23:2490–7.
Wada T, Nakashima T, Hiroshi N, Penninger JM. RANKL-RANK signaling in osteoclastogenesis and bone disease. Trends Mol Med. 2006;12:17–25.
Boyce BF, Xing L. Biology of RANK, RANKL, and osteoprotegerin. Arthritis Res Ther. 2007;9 Suppl 1:S1.
Hanada R, Hanada T, Sigl V, Schramek D, Penninger JM. RANKL/RANK- beyond bones. J Mol Med (Berl). 2011;89:647–56.
Rinotas V, Niti A, Dacquin R, Bonnet N, Stolina M, Han CY, et al. Novel genetic models of osteoporosis by overexpression of human RANKL in transgenic mice. J Bone Miner Res. 2014;29:1158–69.
Lane D, Matte I, Laplante C, Garde-Granger P, Rancourt C, Piché A. Osteoprotegerin (OPG) activates integrin, focal adhesion kinase (FAK), and Akt signaling in ovarian cancer cells to attenuate TRAIL-induced apoptosis. J Ovarian Res. 2013;6:82.
Cao Y, Bonizzi G, Seagroves TN, Greten FR, Johnson R, Schmidt EV, et al. IKKalpha provides an essential link between RANK signaling and cyclin D1 expression during mammary gland development. Cell. 2001;107:763–75.
Karapanou O, Papadimitriou A. Determinants of menarche. Reprod Biol Endocrinol. 2010;8:115.
Russo J, Russo IH. Development of the human breast. Maturitas. 2004;49:2–15.
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Kaczmarek M. The timing of natural menopause in Poland and associated factors. Maturitas. 2007;57:139–53.
Otero UB, Chor D, Carvalho MS, Faerstein E, Lopes Cde S, Werneck GL. Lack of association between age at menarche and age at menopause: Pró-Saúde Study, Rio de Janeiro, Brazil. Maturitas. 2010;67:245–50.
Pan R, Liu YZ, Deng HW, Dvornyk V. Association analyses suggest the effects of RANK and RANKL on age at menarche in Chinese women. Climacteric. 2012;15:75–81.
Lu Y, Liu P, Recker RR, Deng HW, Dvornyk V. TNFRSF11A and TNFSF11 are associated with age at menarche and natural menopause in white women. Menopause. 2010;17:1048–54.
Schubert A, Schulz H, Emons G, Grundker C. Expression of osteoprotegerin and receptor activator of nuclear factor-κB ligand (RANKL) in HCC70 breast cancer cells and effects of treatment with gonadotropin-releasing hormone on RANKL expression. Gynecol Endocrinol. 2008;24:331–8.
Bord S, Ireland DC, Beavan SR, Compston JE. The effects of estrogen on osteoprotegerin, RANKL, and estrogen receptor expression in human osteoblasts. Bone. 2003;32:136–41.
Cannon JG, Kraj B, Sloan G. Follicle-stimulating hormone promotes RANK expression on human monocytes. Cytokine. 2011;53:141–4.
Wang YD, Tao MF, Wang L, Cheng WW, Wan XP. Selective regulation of osteoblastic OPG and RANKL by dehydroepiandrosterone through activation of the estrogen receptor β-mediated MAPK signaling pathway. Horm Metab Res. 2012;44:494–500.
Bashir A, Mak YT, Sankaralingam S, Cheung J, McGowan NW, Grigoriadis AE, et al. Changes in RANKL/OPG/RANK gene expression in peripheral mononuclear cells following treatment with estrogen or raloxifene. Steroids. 2005;70:847–55.
Fernandez-Valdivia R, Mukherjee A, Ying Y, Li J, Paquet M, DeMayo FJ, et al. The RANKL signaling axis is sufficient to elicit ductal side-branching and alveologenesis in the mammary gland of the virgin mouse. Dev Biol. 2009;328:127–39.
Kim NS, Kim HJ, Koo BK, Kwon MC, Kim YW, Cho Y, et al. Receptor activator of NF-kappaB ligand regulates the proliferation of mammary epithelial cells via Id2. Mol Cell Biol. 2006;26:1002–13.
Delahanty RJ, Beeghly-Fadiel A, Long JR, Gao YT, Lu W, Xiang YB, et al. Evaluation of GWAS-identified genetic variants for age at menarche among Chinese women. Hum Reprod. 2013;28:1135–43.
Pyun JA, Kim S, Cho NH, Koh I, Lee JY, Shin C, et al. Genome-wide association studies and epistasis analyses of candidate genes related to age at menarche and age at natural menopause in a Korean population. Menopause. 2014;21:522–9.
Stolk L, Perry JR, Chasman DI, He C, Mangino M, Sulem P, et al. Meta-analyses identify 13 loci associated with age at menopause and highlight DNA repair and immune pathways. Nat Genet. 2012;44:260–8.
Voutilainen M, Lindfors PH, Lefebvre S, Ahtiainen L, Fliniaux I, Rysti E, et al. Ectodysplasin regulates hormone-independent mammary ductal morphogenesis via NF-κB. Proc Natl Acad Sci U S A. 2012;109:5744–9.
Pavlová S, Klucska K, Vašíček D, Ryban L, Harrath AH, Alwasel SH, et al. The involvement of SIRT1 and transcription factor NF-κB (p50/p65) in regulation of porcine ovarian cell function. Anim Reprod Sci. 2013;140:180–8.
Hanada R, Hanada T, Penninger JM. Physiology and pathophysiology of the RANKL/RANK system. Biol Chem. 2010;391:1365–70.
Must A, Phillips SM, Naumova EN, Blum M, Harris S, Dawson-Hughes B, et al. Recall of early menstrual history and menarcheal body size: after 30 years, how well do women remember? Am J Epidemiol. 2002;155:672–9.
den Tonkelaar I. Validity and reproducibility of self-reported age at menopause in women participating in the DOM-project. Maturitas. 1997;27:117–23.
The authors would like to thank Prof. Gao Xin Yuan and Dr. Zhang Zeng for technical assistance. We thank all the participants in this study. This study was supported by grants from the National Natural Science Foundation of China (no. 81260133) and Key Projects of Health Department of Jiangxi province, China (no. 20114030).
The authors declare that they have no competing interests.
PD participated in the molecular genetic studies, conducted statistical analyses and drafted the manuscript. ZMW assisted in genetic studies and drafted the manuscript. JL and LNW collected all samples and questionnaire information. ZY helped in recruiting subjects and statistical analyses. PT conceived of the study, participated in its design, recruited participants, revised the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Duan, P., Wang, ZM., Liu, J. et al. Gene polymorphisms in RANKL/RANK/OPG pathway are associated with ages at menarche and natural menopause in Chinese women. BMC Women's Health 15, 32 (2015). https://doi.org/10.1186/s12905-015-0192-3
- Age at menarche
- Age at natural menopause