Skip to main content
Download PDF

Age at menarche, age at natural menopause, and risk of rheumatoid arthritis — a Mendelian randomization study

Arthritis Research & Therapy volume 23, Article number: 108 (2021) Cite this article

Abstract

Background

Hormonal reproductive factors have been suggested to play an important role in the etiology of rheumatoid arthritis (RA), an autoimmune inflammatory disorder affecting primarily women. We conducted a two-sample Mendelian randomization (MR) study examining three relevant exposures, age at menarche (AAM), age at natural menopause (ANM), and age at first birth (AFB) with the risk of RA.

Methods

We collected summary statistics from the hitherto largest GWAS conducted in AAM (N = 329,345), ANM (N = 69,360), AFB (N = 251,151), and RA (Ncase = 14,361, Ncontrol = 43,923), all of European ancestry. We constructed strong instruments using hundreds of exposure-associated genetic variants and estimated causal relationship through different MR approaches including an inverse-variance weighted method, an MR-Egger regression and a weighted median method. We conducted a multivariable MR to control for pleiotropic effect acting in particular through obesity and socioeconomic status. We also performed important sensitivity analyses to verify model assumptions.

Results

We did not find any evidence in support for a causal association between genetically predicted reproductive factors and risk of RA (ORper-SD increment in AAM = 1.06 [0.98–1.15]; ORper-SD increment in ANM = 1.05 [0.98–1.11], OR per-SD increment in AFB = 0.85 [0.65–1.10]). Results remained consistent after removing palindromic SNPs (ORper-SD increment in AAM = 1.06 [0.97–1.15], ORper-SD increment in ANM = 1.05 [0.98–1.13], ORper-SD increment in AFB = 0.81 [0.61–1.07]) or excluding SNPs associated with potential confounding traits (ORper-SD increment in AAM = 1.03 [0.94–1.12], ORper-SD increment in ANM = 1.04 [0.95–1.14]). No outlying instrument was identified through the leave-one-out analysis.

Conclusions

Our MR study does not convincingly support a casual effect of reproductive factors, as reflected by age at menarche, age at menopause, and age at first birth, on the development of RA. Despite the largely augmented set of instruments we used, these instruments only explained a modest proportion of phenotypic variance of exposures. Our knowledge regarding this topic is still insufficient and future studies with larger sample size should be designed to replicate or dispute our findings.

Background

Rheumatoid arthritis (RA) is a chronic autoimmune inflammatory disorder affecting mainly women of reproductive age and often leads to disabling outcomes as well as shortened life expectancy if left untreated or not properly controlled. The sex difference in the prevalence of RA has been well documented where the disease strikes women more frequently and severe. For example, the incidence of RA has been estimated to be 4–5 times higher in women than in men below age 50 and twice higher during age 60–70 [1]. In addition, observational studies have suggested that in general, female RA patients do worse than male patients [2].

The reasons for this overrepresentation of women remain unclear but X-linked genetic factors and hormonal components are likely to be involved. Some women develop RA at transitional periods when sex hormones are shifting, for example after pregnancy and/or before menopause [3]. Medications that modulate hormone levels including long-term oral contraceptive use [4] and/or postmenopausal hormone therapy [5] have been found to be associated with a reduced risk of RA. These observations highlight a role of hormonal and reproductive factors in the disease etiology.

Several large-scale epidemiological studies have investigated the relationship between female reproductive factors and RA using three most readily available measures — age at menarche (AAM), age at natural menopause (ANM) and age at first birth (AFB) — yet results remain controversial. For example, the longitudinal Nurses’ Health Study enrolling 121,700 women during 1976–2002 identified an association of age at menarche ≤ 10 years with an increased risk of seropositive RA (RR 1.6, 95% CI 1.1–2.4) but not a significant association with risk of overall RA [6]. The study (NHS, 1976–2010; NHSII 1989–2011) also revealed that early age at natural menopause (≤ 44 years) was associated with an increased risk of seronegative RA (pooled HR 2.4, 95% CI 1.5–4.0) [7]. Data from the Swedish EIRA study, a population-based case-control study of female incident RA cases (2035 cases and 2911 controls, aged 18–70 years) showed an increased risk of ACPA (antibodies to citrullinated peptide antigens)-negative RA in those who were at a young age at first birth (< 23 years) (OR 2.5, 95% CI 1.5–4.1) compared to nulliparous women [8]. An analysis using cross-sectional data from 1892 participants in the Third National Health and Nutrition Examination Survey did not find any association between age at menarche or pregnancy history with RA after menopause [9]. These discrepancies are perhaps not surprising since conventional epidemiological studies generally rely on environmental information and results are likely to be influenced by measurement error, confounding, and reverse causality.

Hormonal reproductive factors including puberty and fertility are influenced by genetic, nutritional, socioeconomic, and environmental factors and can be highly heterogeneous among women. Nevertheless, the genetic regulations in AAM, ANM, and AFB have been recently highlighted by discoveries from large-scale genome-wide association studies (GWAS) leveraging millions of women of European ancestry. These results provide a valuable opportunity to utilize a novel statistical approach Mendelian randomization (MR) to make causal inference — an approach that uses genetic variants (single nucleotide polymorphisms, SNPs) as instrumental variables (IVs) to assess a causal effect of a risk factor on an outcome from observational data. Since SNPs are randomly assigned at conception and always precede disease onset, results from MR are less susceptible to confounding and reverse causation, which are the major limitations of conventional observational studies [10]. To the best of our knowledge, no MR analysis has been performed to examine a potential causal association between hormonal reproductive factors and development of RA, of which findings may help address patient concerns in topics of puberty, fertility, motherhood, and RA as well as improve our knowledge on the biological mechanisms underlying RA.

Therefore, we aim to conduct the first and also the largest two-sample MR on three female reproductive factors (AAM, ANM, and AFB) with the risk of RA. Genetic variants associated with each reproductive event were used as instrumental variable (IVs). IV-exposure associations were extracted from the recently published and also the largest GWAS(s) conducted in AAM (N = 329,345), ANM (N = 69,360), and AFB (N = 251,151) [11,12,13]. IV-outcome associations were extracted from the largest GWAS conducted in RA (NRA = 14,361, Ncontrol = 43,923) [14].

Methods

We conducted current study applying a standard two-sample framework where IV-exposure associations and IV-outcome associations come from two sets of independent non-overlapping individuals. To reduce population stratification, we included only individuals of European ancestry.

Exposure

Three reproductive exposures demonstrated by previous GWAS(s) as having polygenic components were involved. Age at menarche, a milestone in female pubertal development, varies markedly among females. The genetic regulation in AAM has been highlighted by a recent meta-GWAS incorporating 329,345 women collected by the ReproGen consortium, 23andMe, and UK Biobank and identified 389 independent AAM-associated signals spreading over 10 biological pathways [11]. Age at natural menopause poses a substantial impact on infertility and disease risk including breast cancer and cardiovascular events. The genetic architecture of ANM has been examined by a recent GWAS of 69,360 women identifying 54 independent signals located in 44 genomic regions, most of which were found to contain one or more DNA damage response pathway genes [12]. Reproductive behavior such as age at first birth is known to be partly driven by biological processes. A recent GWAS has examined the genetic architecture of reproductive tempo defined by AFB in 251,151 women and identified 10 AFB-associated loci [13].

In all three GWAS(s), independent signals (our IVs) were defined as the following: A list of index variants was first defined using a distance-based metric, by which any SNPs passing the two-tailed threshold of significance (P < 5 × 10− 8) within 1 Mb of another significant SNP were considered to be located in the same locus. This list of signals was then augmented using approximate conditional analysis in GCTA, using an LD reference panel from the UK Biobank study. Only secondary signals that were uncorrelated (r2 < 0.05) were included in the final list. All IVs passed quality control procedures under minor allele frequency > 0.001 and Hardy-Weinberg Equilibrium > 1 × 10− 6. IV-exposure associations were extracted from each GWAS [11,12,13].

Outcome

IV-outcome associations were obtained from a meta-GWAS involving 18 participating cohorts totaling 14,361 RA cases and 43,923 controls of European ancestry. To the best of our knowledge, none of the participants in these 18 studies overlapped with participants in the exposure GWAS(s) [14].

Statistical analysis

Mendelian randomization analysis

We applied several MR methods including an inverse-variance weighted approach (IVW) [15], a maximum likelihood-based method [16], an MR-Egger regression [17] and a weighted median approach [18].

Briefly, IVW represents the main conventional approach which only gives consistent estimates if all genetic variants in the analysis are valid instrumental variables. When the IVs are weak, IVW tends to underestimate the true variation of the estimate, while the likelihood method gives appropriately estimated confidence intervals. MR-Egger evaluates the directional pleiotropic effect of instrumental variables, of which the intercept term can be interpreted as an estimate of the average pleiotropy of genetic variations. The weighted median method is robust to outliers and provides consistent estimates even when 50% of the genetic variants are invalid IVs; and is considered as relatively more robust to horizontal pleiotropy.

Sensitivity analysis

A valid MR analysis is defined by three key model assumptions — the IVs are strongly associated with the risk factor of interest (relevance), share no common cause with the outcome (independence), and affect outcome solely through the exposure (exclusion restriction) [10]. Upon the satisfaction of all three assumptions, causal inferences between exposure(s) and outcome(s) can be made based on observational data.

We performed several important sensitivity analyses to verify MR model assumptions. For each index SNP, we searched for its potential association with confounding traits in GWAS catalog and conducted analysis excluding pleiotropic SNPs. Moreover, we used a robust adjusted profile score (MR-RAPS) approach which is robust to both systemic and idiosyncratic pleiotropy and performed excellently in all the numerical examples [19]. Educational attainment and obesity are two important confounders affecting both reproductive traits and risk of RA [20]. We further integrated GWAS summary statistics and additional IVs on education and BMI, and conducted an IVW-based multivariable MR (MVMR) to estimate the direct effect of reproductive factors controlling for the effect of BMI and education [21, 22]. Finally, we excluded one SNP at-a-time and performed IVW on the remaining SNPs to identify outlying IVs.

We calculated statistical power using the non-centrality parameter of the test statistic as suggested by Brion et al. (http://cnsgenomics.com/shiny/mRnd/). All analyses were conducted with packages “TwoSampleMR”, “MRInstruments”, and “Mendelian Randomization” in R v3.6.3.

Ethics/consent statement

Our study is a secondary analysis of existing, de-identified, summary-level GWAS data. Specific ethics and consent statement for each GWAS examined in this study can be found in the original GWAS publications.

Results

Overall, we did not find convincing evidence in support for a causal relationship between the three hormonal related exposures and risk of RA. Specifically, genetically predicted AAM did not significantly influence the risk of RA using IVW approach (ORper-SD increment in AAM [95% CI], 1.06 [0.98–1.15]). Estimates remained consistent across different methods (ORper-SD increment in AAM [95% CI] for maximum likelihood 1.07 [1.00–1.14], for MR-Egger regression 1.11 [0.90–1.36], for weighted median 1.08 [0.97–1.21]). We did not observe apparent sign of pleiotropy (P for MR-Egger intercept = 0.69). Similarly, we did not find any compelling evidence supporting a casual association of genetically instrumented ANM with RA either using IVW (ORper-SD increment in ANM [95% CI], 1.05 [0.98–1.11]) or MR-Egger regression (ORper-SD increment in ANM [95% CI], 1.04 [0.90–1.20]). Null finding was identified using the weighted median approach (ORper-SD increment in ANM [95% CI], 1.03 [0.98–1.08]). Consistently, for genetically predicted AFB, we did not observe any significant association with RA using IVW (ORper-SD increment in AFB [95% CI], 0.85 [0.65–1.10]), MR-Egger regression (ORper-SD increment in AFB [95% CI], 3.32 [0.36–30.81]), or weighted median approach (ORper-SD increment in AFB [95% CI], 0.90 [0.73–1.10]). Indeed, for both ANM and AFB, significant results appeared using the maximum likelihood approach, yet this method is known to provide better power by neglecting horizontal pleiotropy and results were not supported by other methods. We did not observe apparent signs of horizontal pleiotropy (P for MR-Egger intercept = 0.92 for ANM and = 0.26 for AFB) (Table 1).

Table 1 Genetically predicted age at menarche, age at menopause and, age at first birth with the risk of rheumatoid arthritis. Results from primary Mendelian randomization analysis as well as sensitivity analyses based on a subset of instruments
Full size table

Results remained consistent after removing palindromic SNPs (OR [95% CI], 1.06 [0.97–1.14] for AAM; 1.05 [0.98–1.13] for ANM; 0.81 [0.61–1.07] for AFB). As shown in Supplementary Tables 1 and 2, the AAM and ANM IVs were also found to be associated with potential confounders while none of the 10 AFB-associated IVs was cited by the NHGRI-EBI Catalog (Supplementary Table 3). Using 283 AAM-associated IVs and 42 ANM-associated IVs excluding pleiotropic SNPs, we did not detect a causal effect of AAM or ANM on RA risk (OR [95% CI], 1.03 [0.94–1.12] for AAM, 1.04 [0.95–1.14] for ANM), corroborating our primary findings (Table 1).

To effectively control for pleiotropy, we next looked into the Robust Adjusted Profile Score (RAPS) approach which is robust to both systemic and idiosyncratic pleiotropy [19]. We performed MR-RAPS estimator and found that results remained largely consistent with our primary findings (Table 2).

Table 2 Genetically predicted age at menarche, age at menopause, and age at first birth with the risk of rheumatoid arthritis. A sensitivity analysis using MR-RAPS assuming over-dispersion
Full size table

Education and BMI are two modifiable risk factors, both of which play an important role in the etiology of RA and shape the reproductive exposures. We next conducted an IVW-based MVMR to estimate a direct effect of reproductive factors on RA accounting for the confounding effect from obesity and socioeconomic status. The results of MVMR remained consistent with our primary findings. The effect of AAM with RA did not alter substantially after adjusting for BMI (OR [95% CI], 0.97 [0.83–1.13]) or education (OR [95% CI], 1.07 [0.98–1.16]). Similarly, for ANM, we did not observe any significant effect with RA after adjusting for BMI (OR [95% CI], 1.06 [0.99–1.27]) or education (OR [95% CI], 1.04 [0.98–1.11]). For AFB, similar null effect was found after adjusting for BMI (OR [95% CI], 0.85 [0.57–1.24]) or education (OR [95% CI], 1.10 [0.58–2.11]) (Table 3).

Table 3 Genetically predicted age at menarche, age at menopause, and age at first birth with the risk of rheumatoid arthritis. Multivariable analysis adjusting for the effect of body mass index and year of education
Full size table

In the leave-one-out analysis where we iteratively removed one SNP at a time and performed IVW using the remaining SNPs, we did not observe apparent outlying SNPs and the odds ratios were in accordance with our primary findings, aggregating closely around the expected value of estimation (Fig. 1).

Fig. 1
figure 1

Sensitivity analysis leaving one SNP out at a time for the association between reproductive factors and RA risk. a The distribution of odds ratios from 389 leave-one-out analysis conducted for age at menarche and RA risk. b The distribution of odds ratios from 54 leave-one-out analysis conducted for age at menopause and RA risk. c The distribution of odds ratios from 10 leave-one-out analysis conducted for age at first birth and RA risk

Full size image

Finally, we calculated the power of our analysis. As shown in Table 4, the sample size of the RA GWAS was 58,284 with 24.64% cases. According to the three exposure GWAS(s), 7.4% of phenotypic variance of AAM could be explained by the 389 index SNPs, 5.7% of ANM phenotypic variance could be explained by the 54 index SNPs, and 0.2% of AFB phenotypic variance could be explained by the 10 SNPs. Under current situation, for AAM, our study had 80% power to detect a causal effect of a 10.4% (i.e., ORs of 1.104) increase in RA risk. For ANM, the minimal detectable effect was 12% increase (i.e., ORs of 1.12). For AFB, the minimal detectable effect was 70% increase (i.e., ORs of 1.70). We presented a range of power estimations in Table 4.

Table 4 Power calculation of the current analysis
Full size table

Discussion

In this study, we examined a putative causal relationship between three hormonal reproductive traits (AAM, ANM, and AFB) and an autoimmune inflammatory disease RA which affects mainly women. We capitalized on the summary statistics of the largest GWAS(s) conducted for these traits in European ancestry populations and constructed strong instruments using hundreds of SNPs associated with the exposures (F-statistic for AAM 67.6, for ANM 77.6, for AFB 50.3). We did not find convincing evidence in support for a causal effect of reproductive factors on RA using univariable MR analyses. Consistent null associations were identified by sensitivity analysis and multivariable MR analysis, demonstrating the robustness of our findings.

Current results from conventional epidemiological studies on this topic remain controversial, yet many studies point towards a positive association. For example, a study enrolling 121,700 female nurses found that age at menarche ≤ 10 years was associated with an increased risk of seropositive RA (RR 1.60, 95% CI 1.10–2.40) [6]. A community-based health survey including 30,447 subjects (18,326 women) between 1991 and 1996 found an association between early age at menopause (≤ 45 years) and subsequent development of RA (OR 2.42, 95% CI 1.32–4.45), which remained significant after adjusting for smoking, level of education and length of breastfeeding (OR 1.92, 95% CI 1.02–3.64) [23]. A prospective cohort study of 31,336 North America women reported similar findings (RRmenopause >51 vs. menopause <45 0.64, 95%CI 0.41–1.00) [24].

Our large-scale MR, however, did not identify a putative causal link between the three well-defined hormonal exposures and risk of RA. Several reasons underlie such a discrepancy. First of all, reproductive factors are highly complicated and heterogenous traits shaped by both genetic and environmental factors and genetics alone does not fully capture the phenotypic variance of these traits. For example, age at first birth is a human behavioral trait influenced largely by psychosocial, cultural, and financial factors rather than the genetics. Secondly, results from previous epidemiological studies are likely to be impaired by confounding factors. For example, obesity is an important confounder affecting both the exposure and the outcome. An MR study demonstrated that a 1-year delay in age at menarche reduced adult BMI by 0.38 kg/m2 (95% CI 0.25–0.51 kg/m2) [25]. Global adiposity is a robust causal risk factor for RA as demonstrated by our recently published MR [26]. It is likely that traditional epidemiological investigations did not adequately control for the confounding effects from obesity. The protective effect of education on RA has been reported by observational studies [27, 28]. An MR study identified that a 1-year later in age at menarche increased 0.14 years (53 days) of time spent in education [29]. We performed a MVMR to control for the effect of adiposity and education, and the negative results corroborating our main findings on a null association. Finally, it is also likely that the true causal effect of reproductive factors on RA is modest, which our study is underpowered to identify.

Biological mechanisms underlying hormonal factors and the development of RA remain unclear. The effect of sex hormones on the immune system and their interaction with environmental and genetic factors may partly explain the higher prevalence of RA observed among women. Estrogen is a complex modulator to the immune system exerting both a stimulatory and an inhibitory effect [30]. For example, estrogens at periovulatory to pregnancy levels stimulate B cells and the Th2 response and support the survival of auto-reactive T and B cell clones. On other hand, estrogens could inhibit cell-mediated responses such as the differentiation to Th17 cells [30,31,32]. A reduced risk of RA onset during pregnancy compared to an increased risk postpartum suggests a role the hormonal changes or the exposure to fetus paternal HLA in RA onset [33].

Our study has several strengths. To the best of our knowledge, no MR has been performed to assess the relationship between reproductive factors and RA. We incorporated three different reproductive traits (age at menarche, age at natural menopause, and age at first birth) reflecting the length of reproductive period and complementing each other well. Moreover, we conducted important sensitivity analyses to verify MR model assumptions. We selected the most significant independent SNPs identified by the largest GWAS, so all were robustly associated with exposure of interest, guaranteeing “relevance” assumption. We excluded SNPs associated with potential confounders on the exposure-outcome relationship to satisfy “exclusion restriction” assumption. The consistent results observed across different approaches, further lend support to our findings.

We have to acknowledge several limitations. Firstly, our analysis was performed using the European populations which restricted its generalizability. Secondly, the genetic instruments of three exposures (AAM, ANM, and AFB) we used as proxies for hormonal reproductive characteristics captured only a modest proportion of phenotypic variance. Reproductive factors are complex traits influenced by different components such as genetic, environmental, and socioeconomic factors as well as their complex interactions. The design of our study disables us to take into account environmental impacts. Thirdly, the association between genetically predicted age at each of the reproductive events and risk of RA was evaluated fitting the exposure as a continuous variable — we can still not exclude a non-linear effect which was not captured by our study with the current availability of data. Future work on such topics may be focused on categorized age of reproductive events. Fourthly, our study was conducted using overall RA (a majority of which are seropositive RA, > 85%) without specifying disease subsets characterized by the presence/absence of antibodies to citrullinated peptides or rheumatic factors. It is possible that hormonal factors influence different RA subsets via a distinct way. It is also likely that other factors such as hormone use and health conditions confound our results, in addition to the only two confounders (obesity and education) considered in the current study. However, it is difficult to control for the effect of hormone therapy due to limited availability of genetic data underlying this trait. Finally, power calculations showed that potential weak effects were difficult to be detected in our analysis.

Conclusions

In summary, using both univariable and multivariable Mendelian randomization approaches, we could not provide evidence supporting a casual effect of reproductive factors as reflected by age at menarche, age at menopause, or age at first birth in the development of RA. Our result is not so surprising considering the relatively weak genetic instruments and power. The findings represent a preliminary but important step towards the identification of causal associations between female hormonal reproductive factors and a female disease RA. As some hormonal factors are potentially modifiable, understanding their precise role is essential for future preventive interventions focusing on women at high risk. Our knowledge regarding this topic is still insufficient and future studies with larger sample size and better power should be designed to increase our knowledge in this field.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

RA:

Rheumatoid arthritis

AAM:

Age at menarche

ANM:

Age at natural menopause

AFB:

Age at first birth

ACPA:

Antibodies to citrullinated peptide antigens

MR:

Mendelian randomization

GWAS:

Genome-wide association studies

SNPs:

Single nucleotide polymorphisms

IV:

Instrumental variable

IVW:

Inverse-variance weighted method

BMI:

Body mass index

MVMR:

Multivariable Mendelian randomization

References

  1. 1.

    Kvien TK, Uhlig T, Ødegård S, Heiberg MS. Epidemiological aspects of rheumatoid arthritis: the sex ratio. Ann N Y Acad Sci. 2006;1069(1):212–22. https://doi.org/10.1196/annals.1351.019.

    Article  PubMed  Google Scholar 

  2. 2.

    van Vollenhoven RF. Sex differences in rheumatoid arthritis: more than meets the eye. BMC Med. 2009;7(1):12. https://doi.org/10.1186/1741-7015-7-12.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Desai MK, Brinton RD. Autoimmune disease in women: endocrine transition and risk across the lifespan. Front Endocrinol. 2019;10:265. https://doi.org/10.3389/fendo.2019.00265.

    Article  Google Scholar 

  4. 4.

    Orellana C, Saevarsdottir S, Klareskog L, Karlson EW, Alfredsson L, Bengtsson C. Oral contraceptives, breastfeeding and the risk of developing rheumatoid arthritis: results from the Swedish EIRA study. Ann Rheum Dis. 2017;76(11):1845–52. https://doi.org/10.1136/annrheumdis-2017-211620.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Orellana C, Saevarsdottir S, Klareskog L, Karlson EW, Alfredsson L, Bengtsson C. Postmenopausal hormone therapy and the risk of rheumatoid arthritis: results from the Swedish EIRA population-based case-control study. Eur J Epidemiol. 2015;30(5):449–57. https://doi.org/10.1007/s10654-015-0004-y.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Karlson EW, Mandl LA, Hankinson SE, Grodstein F. Do breast-feeding and other reproductive factors influence future risk of rheumatoid arthritis? Results from the Nurses' health study. Arthritis Rheum. 2004;50(11):3458–67. https://doi.org/10.1002/art.20621.

    Article  PubMed  Google Scholar 

  7. 7.

    Bengtsson C, Malspeis S, Orellana C, Sparks JA, Costenbader KH, Karlson EW. Association between menopausal factors and the risk of seronegative and seropositive rheumatoid arthritis: results from the nurses’ health studies. Arthritis Care Res (Hoboken). 2017;69(11):1676–84. https://doi.org/10.1002/acr.23194.

    CAS  Article  Google Scholar 

  8. 8.

    Orellana C, Wedren S, Kallberg H, Holmqvist M, Karlson EW, Alfredsson L, Bengtsson C, Group ES. Parity and the risk of developing rheumatoid arthritis: results from the Swedish Epidemiological Investigation of Rheumatoid Arthritis study. Ann Rheum Dis. 2014;73(4):752–5. https://doi.org/10.1136/annrheumdis-2013-203567.

    Article  PubMed  Google Scholar 

  9. 9.

    Beydoun HA, el-Amin R, McNeal M, Perry C, Archer DF. Reproductive history and postmenopausal rheumatoid arthritis among women 60 years or older: Third National Health and Nutrition Examination Survey. Menopause. 2013;20(9):930–5. https://doi.org/10.1097/GME.0b013e3182a14372.

    Article  PubMed  Google Scholar 

  10. 10.

    Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362:k601.

    Article  Google Scholar 

  11. 11.

    Day FR, Thompson DJ, Helgason H, Chasman DI, Finucane H, Sulem P, Ruth KS, Whalen S, Sarkar AK, Albrecht E, et al. Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk. Nat Genet. 2017;49(6):834–41. https://doi.org/10.1038/ng.3841.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Day FR, Ruth KS, Thompson DJ, Lunetta KL, Pervjakova N, Chasman DI, Stolk L, Finucane HK, Sulem P, Bulik-Sullivan B, et al. Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair. Nat Genet. 2015;47(11):1294–303. https://doi.org/10.1038/ng.3412.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Barban N, Jansen R, de Vlaming R, Vaez A, Mandemakers JJ, Tropf FC, Shen X, Wilson JF, Chasman DI, Nolte IM, et al. Genome-wide analysis identifies 12 loci influencing human reproductive behavior. Nat Genet. 2016;48(12):1462–72. https://doi.org/10.1038/ng.3698.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikari K, Kochi Y, Ohmura K, Suzuki A, Yoshida S, et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature. 2014;506(7488):376–81. https://doi.org/10.1038/nature12873.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG, Consortium E-I. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52. https://doi.org/10.1007/s10654-015-0011-z.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65. https://doi.org/10.1002/gepi.21758.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25. https://doi.org/10.1093/ije/dyv080.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14. https://doi.org/10.1002/gepi.21965.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Zhao Q, Wang J, Hemani G, Bowden J, Small D. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Ann Stat. 2018;48(3):1–59.

  20. 20.

    Sanderson E. Multivariable Mendelian randomization and mediation. Cold Spring Harb Perspect Med. 2021;11(2):a038984.

  21. 21.

    Lee JJ, Wedow R, Okbay A, Kong E, Maghzian O, Zacher M, Nguyen-Viet TA, Bowers P, Sidorenko J, Karlsson Linnér R, et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat Genet. 2018;50(8):1112–21. https://doi.org/10.1038/s41588-018-0147-3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Yengo L, Sidorenko J, Kemper KE, Zheng Z, Wood AR, Weedon MN, Frayling TM, Hirschhorn J, Yang J, Visscher PM. Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry. Hum Mol Genet. 2018;27(20):3641–9. https://doi.org/10.1093/hmg/ddy271.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Pikwer M, Bergstrom U, Nilsson JA, Jacobsson L, Turesson C. Early menopause is an independent predictor of rheumatoid arthritis. Ann Rheum Dis. 2012;71(3):378–81. https://doi.org/10.1136/ard.2011.200059.

    Article  PubMed  Google Scholar 

  24. 24.

    Merlino LA, Cerhan JR, Criswell LA, Mikuls TR, Saag KG. Estrogen and other female reproductive risk factors are not strongly associated with the development of rheumatoid arthritis in elderly women. Semin Arthritis Rheum. 2003;33(2):72–82. https://doi.org/10.1016/S0049-0172(03)00084-2.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Gill D, Brewer CF, Del Greco MF, Sivakumaran P, Bowden J, Sheehan NA, Minelli C. Age at menarche and adult body mass index: a Mendelian randomization study. Int J Obes. 2018;42(9):1574–81. https://doi.org/10.1038/s41366-018-0048-7.

    Article  Google Scholar 

  26. 26.

    Tang B, Shi H, Alfredsson L, Klareskog L, Padyukov L, Jiang X. Obesity-related traits and the development of rheumatoid arthritis - evidence from genetic data. Arthritis Rheumatol 2021;73(2):203–11.

  27. 27.

    Kwon JM, Rhee J, Ku H, Lee EK. Socioeconomic and employment status of patients with rheumatoid arthritis in Korea. Epidemiol Health. 2012;34:e2012003. https://doi.org/10.4178/epih/e2012003.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Bergstrom U, Jacobsson LT, Nilsson JA, Wirfalt E, Turesson C. Smoking, low formal level of education, alcohol consumption, and the risk of rheumatoid arthritis. Scand J Rheumatol. 2013;42(2):123–30. https://doi.org/10.3109/03009742.2012.723744.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Gill D, Del Greco MF, Rawson TM, Sivakumaran P, Brown A, Sheehan NA, Minelli C. Age at menarche and time spent in education: a Mendelian randomization study. Behav Genet. 2017;47(5):480–5. https://doi.org/10.1007/s10519-017-9862-2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Straub RH. The complex role of estrogens in inflammation. Endocr Rev. 2007;28(5):521–74. https://doi.org/10.1210/er.2007-0001.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Sthoeger ZM, Chiorazzi N, Lahita RG. Regulation of the immune response by sex hormones. I. In vitro effects of estradiol and testosterone on pokeweed mitogen-induced human B cell differentiation. J Immunol (Baltimore, Md : 1950). 1988;141(1):91–8.

    CAS  Google Scholar 

  32. 32.

    Hughes GC, Choubey D. Modulation of autoimmune rheumatic diseases by oestrogen and progesterone. Nat Rev Rheumatol. 2014;10(12):740–51. https://doi.org/10.1038/nrrheum.2014.144.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Silman A, Kay A, Brennan P. Timing of pregnancy in relation to the onset of rheumatoid arthritis. Arthritis Rheum. 1992;35(2):152–5. https://doi.org/10.1002/art.1780350205.

    CAS  Article  PubMed  Google Scholar 

Download references

Acknowledgements

Summary statistics for the genetic associations with AAM, ANM, AFB, BMI, education, and RA were obtained from GWAS respectively by Day FR et al., Day FR et al., Barban N et al., Yengo L et al., Lee JJ et al., and the RA meta-GWAS conducted by Okada Y et al. We thank all investigators for sharing the genome-wide summary statistics.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or nonprofit sectors. Open Access funding provided by Karolinska Institute.

Author information

Affiliations

  1. The First Affiliated Hospital of Zhejiang Chinese Medical University, Hangzhou, China

    Jingjing Zhu

  2. Department of Gynecology, Affiliated Hangzhou First People’s Hospital, Zhejiang University School of Medicine, Hangzhou, China

    Zheng Niu

  3. Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institute, Tomtebodavägen 5, 17 177, Stockholm, Sweden

    Lars Alfredsson & Xia Jiang

  4. Center for Molecular Medicine, Karolinska Institute, Stockholm, Sweden

    Lars Klareskog & Leonid Padyukov

  5. Program in Genetic Epidemiology and Statistical Genetics, Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, USA

    Xia Jiang

  6. West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu, China

    Xia Jiang

Authors
  1. Jingjing Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Zheng Niu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Lars Alfredsson
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Lars Klareskog
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Leonid Padyukov
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xia Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

JJZ and XJ designed the study, analyzed the data, interpreted the results, and wrote the manuscript. ZN, LA, LK, and LP interpreted the results and revised the manuscript. All authors provided the corresponding author with permission to be named in the manuscript. XJ is the guarantor of this study. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Xia Jiang.

Ethics declarations

Ethics approval and consent to participate

The summary statistics used in the present study are aggregated level of data which do not contain any personal information. The original GWAS have obtained ethical approval from relevant ethics review committees.

Consent for publication

Not applicable

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table 1.

The characteristic of age at menarche associated index SNPs, their effect sizes with exposure and outcome, as well as their associations with potential confounders. Supplementary Table 2. The characteristic of age at natural menopause associated index SNPs, their effect sizes with exposure and outcome, as well as their associations with potential confounders. Supplementary Table 3. The characteristic of age at first birth associated index SNPs, their effect sizes with exposure and outcome.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhu, J., Niu, Z., Alfredsson, L. et al. Age at menarche, age at natural menopause, and risk of rheumatoid arthritis — a Mendelian randomization study. Arthritis Res Ther 23, 108 (2021). https://doi.org/10.1186/s13075-021-02495-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13075-021-02495-x

Keywords

  • Mendelian randomization
  • Rheumatoid arthritis
  • Age at menarche
  • Age at natural menopause
  • Age at first birth

Arthritis Research & Therapy

ISSN: 1478-6362

Contact us