- 1Department of Reproductive Medicine, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2National Research Institute for Family Planning, Beijing, China
- 3Graduate School, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
Background: Endometriosis is a chronic inflammatory disease of women during their reproductive years. The relationship between the severity and location of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery remains unclear.
Methods: We explored the association between the various phenotypes of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery, using two-sample Mendelian randomization (MR) and summary data on endometriosis stages and locations from the FinnGen consortium and women’s menstruation, ovulation, reproductive function, and mode of delivery from OpenGWAS and ReproGen. Inverse-variance weighting was used for the primary MR analysis. In addition, a series of sensitivity analyses, confounding analyses, co-localization analyses, and multivariate MR analyses were performed.
Results: MR analysis showed a negative effect of moderate to severe endometriosis on age at last live birth (OR = 0.973, 95% CI: 0.960–0.986) and normal delivery (OR = 0.999, 95% CI: 0.998–1.000; values for endpoint were excluded), ovarian endometriosis on age at last live birth (OR = 0.976, 95% CI: 0.965–0.988) and normal delivery (OR = 0.999, 95% CI: 0.998–1.000; values for endpoint were excluded), and fallopian tubal endometriosis on excessive irregular menstruation (OR = 0.966, 95% CI: 0.942–0.990). Bidirectional MR analysis showed that age at menarche had a negative causal effect on intestinal endometriosis (OR = 0.417, 95% CI: 0.216–0.804). All MR analyses were confirmed by sensitivity analyses, and only the genetic effects of moderate to severe endometriosis on normal delivery and age at last live birth were supported by co-localization evidence.
Conclusion: Our findings deepen the understanding of the relationship between various types of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery and clarify the important role of moderate to severe endometriosis.
Introduction
Endometriosis is a chronic inflammatory disease that primarily affects women of reproductive age, with a prevalence of up to 10% (1, 2). It is characterized by the presence of endometrial tissue outside the body of the uterus and is associated with pelvic pain, dysmenorrhea, and infertility and its treatment is limited to hormonal therapy or surgical removal of disease (1, 3). Endometriosis is typically staged according to the revised American Society of Reproductive Medicine (ASRM) criteria, with milder lesions in stages 1–2 and more severe lesions in stages 3–4. In addition, ectopic endometrium can invade almost any part of the body, including the lungs and pleura, but most commonly the pelvic organs and parietal peritoneum (4).
The etiology of endometriosis remains unknown, but several risk factors have been reported. To date, it is known that earlier age at menarche and shorter menstrual cycles are associated with an increased risk of endometriosis, whereas higher parity is associated with a lower risk (5). However, these studies have lacked associations between various phenotypic features of endometriosis, including ASRM stages and locations, and menstruation, ovulation, reproductive function, and mode of delivery, consistently focusing on the broad category of endometriosis, with far less information on its subcategories. Moreover, most of these studies only included risk factors for endometriosis and lacked findings on the possible consequences of endometriosis.
In addition, studies on factors such as menstruation, ovulation, reproductive function, and mode of delivery in relation to endometriosis have been mostly limited to epidemiologic observations and cohort studies. Therefore, it is difficult to draw causal conclusions from these studies due to possible biases and potential confounders. In this manuscript, we investigated the risk factors and consequences of endometriosis based on GWAS using Mendelian randomization (MR). MR is an analytical method for assessing causal inference in epidemiologic studies in which endometriosis can be either an exposure or an outcome, and uses genetic variants that are strongly associated with exposure as instrumental variables (IVs) to assess the causal relationship between exposure and outcome.
Here, we reported a comprehensive GWAS-based MR analysis of the association between different stages and locations of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery, elucidating the causal effects of endometriosis sub-phenotypes leading to various components of the complex pathophysiology and improving our understanding of the risk factors and consequences of endometriosis.
Methods
We selected genetic variants that were strongly associated with different ASRM stages and locations of endometriosis from the FinnGen consortium with a sample size of 210,870 women, including endometriosis ASRM stages 1–2 and ASRM stages 3–4, and endometriosis of ovarian, fallopian tube and pelvic peritoneum, rectovaginal septum and vagina, and intestine. At the same time, we collected the summary level data from OpenGWAS and ReproGen study for GWAS on menstruation (excessive irregular menstruation, and menstrual cycle length), ovulation (age at menarche, age at natural menopause, years ovulating), reproductive function (age at first live birth, age at last live birth, number of live births, spontaneous abortion), and mode of delivery (normal delivery, caesarean section) (6, 7). In this study, we used two-sample MR analysis to investigate the genetic basis between different stages and locations of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery and assess the causal effects among them.
Data information
Phenotypes were described in detail as follows: (1) different ASRM stages and locations of endometriosis; (2) menstruation (excessive irregular menstruation and menstrual cycle length); (3) ovulation (age at menarche, age at natural menopause, and years ovulating); (4) reproductive function (age at first live birth, age at last live birth, number of live births, and spontaneous abortion); (5) mode of delivery (normal delivery and caesarean section). All genetic association estimates were derived from GWAS studies conducted in European populations, and the design and analysis protocols (such as adjustment for covariates) were available for query in publications (7–9). An overview of the data is shown in Figure 1, and detailed information can be found in Additional File 2: Table 1.
Figure 1 Data overview. It shows different ASRM stages and locations of endometriosis, menstruation, ovulation, reproductive function, and mode of delivery. Each phenotype labels the data source and the number of cases and controls.
Endometriosis
In FinnGen, N14 codes for endometriosis. The subitems of N14 were obtained, which were ASRM stages 1–2 and ASRM stages 3–4, as well as endometriosis of ovary, fallopian tube, pelvic peritoneum, rectovaginal septum and vagina, and intestine. Controls for a particular phenotype of endometriosis were those individuals who were not cases.
Menstruation
In OpenGWAS, ICD10: N92.1 was included in the analysis, and women with this code were diagnosed with excessive and frequent menstruation with irregular cycles (N = 2,048). It is referred to in this manuscript as excessive irregular menstruation. The length of menstrual cycle refers to the time since the last menstrual period (N = 50,309).
Ovulation
A total of 243,944 individuals of European ancestry who self-reported their age at menarche, ranging from 5 to 25 years, were included in OpenGWAS. Summary data on the age at natural menopause were reported by Ruth et al., who reported an age of 40 to 60 years, and were stored in the ReproGen database (N = 201,323) (8). Years ovulating was defined as the time between age at menarche and menopause, calculated by D’Urso et al. using UK Biobank data, as detailed in the publication (10).
Reproductive function
Age at first and last live birth was collected from women who reported having more than one child. The number of live births ranged from 0 to 22, with a median of 2. ICD10: O03 coded for spontaneous abortion.
Mode of delivery
GWAS summary data on normal deliveries and caesarean sections were available in OpenGWAS, which included 7,946 normal deliveries and 23,610 caesarean sections.
Selection of genetic instruments and data harmonization
We selected all single-nucleotide polymorphisms (SNPs) strongly associated with exposure as instrumental variables (IVs), with a genome-wide significance threshold of p < 5 × 10−8, which was p < 5 × 10−6 for endometriosis of the fallopian tube and intestine. Moreover, the strength of each IV was assessed, and when F >10, the SNP was included (10). All SNPs in linkage disequilibrium (LD) were trimmed with a window of r2 < 0·001 and size <10,000 kb, and at the same time, SNPs with minor allele frequency (MAF) >0.5% were retained. IVs that met the above quality control are shown in Additional File 1: Table 1. Then, SNP exposure and SNP outcome were harmonized with the use of the “harmonise_data” function (11). SNPs for incompatible alleles and being palindromic with intermediate allele frequencies were removed during harmonization, if present.
Mendelian randomization analysis
Univariate MR analyses were performed to estimate the causal effect of each type of endometriosis on (1) menstruation, (2) ovulation, (3) reproductive function, and (4) mode of delivery. Inverse variance weighted (IVW) MR was used as the primary analysis method unless only one SNP was available (in which case the Wald ratio was used). IVW can provide unbiased estimates of outcomes in the absence of horizontal pleiotropy in IVs (12). In addition, the IVW method can account for the heterogeneity of causal estimates obtained from individual variants (13, 14).
Bidirectional Mendelian randomization analysis
To investigate the possibility of reverse causality between exposure and outcome, bidirectional MR analyses were performed (exposures and outcomes were reversed in the analyses). We used a p < 5 × 10−8 to select the IVs associated with exposure, and the p value threshold was appropriately relaxed to 5 × 10−6 for excessive irregular menstruation, menstrual cycle length, spontaneous abortion, normal delivery, and cesarean section. Other analysis details were the same as for the MR analysis. The available IVs are shown in Additional File 1: Table 2.
Sensitivity analyses
For each causal effect detected, we used a series of sensitivity analyses to assess the robustness of the MR analysis. First, Cochran’s Q statistic, including both IVW and MR–Egger methods, was used to evaluate the heterogeneity of IVs in the causal effects. Second, MR–Egger regression was used to evaluate the presence of horizontal pleiotropy; if the intercept term equals zero, it indicates that horizontal pleiotropy does not exist. Third, the MR‐pleiotropy residual sum and outlier (MR-PRESSO) analysis obtained horizontal pleiotropy by detecting abnormal IVs in causal effects. If abnormal IVs were present, outliers were removed, and MR analyses and sensitivity analyses were repeated. When <50% of the instruments showed horizontal pleiotropy, the MR-PRESSO test was most appropriate (15). Finally, we performed a leave-one-out sensitivity analysis by stepwise elimination of each IV to see if the results changed after each SNP was removed.
Confounding analysis
We used the “phenoscanner” package (version 1.0) for confounding analysis (16). Diseases/physical conditions significantly associated with IV were explored using p < 5 × 10−8 as the threshold.
Colocalization analysis
The coloc R package (version 5.2.2) was used to perform a Bayesian test for co-localization of the two traits to estimate the posterior probability of shared variation (17). We retrieved all SNPs within the upper and lower 50 kb of each IV for co-localization analysis to analyze the posterior probability of H4 (PP.H4), with PP.H4 >0.95 indicating co-localization of the two traits.
Multivariable Mendelian randomization
The exposures included in the multivariable MR (MVMR) analysis were ASRM stages 3–4 and ovarian endometriosis. Univariate MR analysis suggested a causal relationship between these exposures and age at last live birth and normal delivery. We used MVMR to account for potential horizontal pleiotropy (18). The inclusion of IVs for MVMR was the same as for MR analysis.
Statistical analysis
The results of MR analysis were presented as odds ratio (OR) and 95% confidence interval (CI). Due to the presence of multiple comparisons with causal effects, we corrected the p-values with the false discovery rate (FDR), with a threshold of an adjusted p-value < 0.1. This work was performed under R packages “TwoSampleMR” (version 0.5.7), “MR-PRESSO” (version 1.0), and “MVMR” (version 0.4) (11, 15).
Results
Causal effects of endometriosis stage and location on menstruation, ovulation, reproductive function, and mode of delivery
As shown in Figure 2, MR analysis showed a negative effect of ASRM stages 3–4 endometriosis on age at last live birth (OR = 0.973, 95% CI: 0.960–0.986) and normal delivery (OR = 0.999, 95% CI: 0.998–1.000, excluding the value of the endpoint 1.000). We also assessed the effect of different locations of endometriosis on various outcomes. IVW estimates showed a negative effect of ovarian endometriosis on age at last live birth (OR = 0.976, 95% CI: 0.965–0.988) and normal delivery (OR = 0.999, 95% CI: 0.998–1.000, values for endpoint 1.000 were excluded). Endometriosis in the fallopian tube was also observed to have a negative effect on excessive irregular menstruation (OR = 0.966, 95% CI: 0.942–0.990) (Figure 2). All MR estimates can be viewed in Additional File 1: Table 3.
Figure 2 MR estimates for the association between endometriosis stages and locations and menstruation, ovulation, reproductive function, and mode of delivery.
Horizontal pleiotropy and heterogeneity analyses showed that the causal effects of ASRM stages 3–4 endometriosis and ovarian endometriosis on normal labor and tubal on excessive menstruation were robust (Additional File 2: Tables 2, 4, 6). MR-PRESSO analyses showed significant horizontal pleiotropy between IVs in the causal effect of ASRM stages 3–4 endometriosis and ovarian endometriosis on age at last live birth, identifying three abnormal IVs, respectively (Additional File 2: Tables 3, 5). After removing the outliers, the MR analysis was re-performed, and horizontal pleiotropy and heterogeneity between IVs were corrected (Additional File 2: Tables 3, 5). The leave-one-out test did not identify SNPs with abnormal effects (Additional File 2: Figure 1). This manuscript only listed the corrected MR Results (the same below).
Bidirectional Mendelian randomization analysis
Bidirectional MR analysis showed that age at menarche had a negative causal effect on intestinal endometriosis (OR = 0.417, 95% CI: 0.216–0.804), with an increase in age at menarche and a decrease in risk of the disease (Figure 3). No other causal effects were detected. Additional File 1: Table 4 demonstrated all MR estimates.
There was horizontal pleiotropy in the IVs for age at menarche; MR-PRESSO analysis excluded one SNP (Additional File 2: Table 7). After removal of outliers, horizontal pleiotropy of IVs was corrected, whereas heterogeneity among IVs remained (Additional File 2: Table 7). However, our results were generated based on the random-effects model (IVW), which remained robust in the presence of heterogeneity (13, 14). The leave-one-out test showed no outliers in the IVs (Additional File 2: Figure 1).
Confounding analysis
IVs in causal effects were summarized and analyzed with various confounders by Phenoscanner. We found that potential confounders of ASRM stages 3–4 and ovarian endometriosis mainly included trunk predicted mass, birth weight, coagulation dysfunction, height, heel bone mineral density, atopic dermatitis, potential confounders of fallopian tubal endometriosis mainly including mean plasma volume, mean hemoglobin, and erythrocyte distribution width, and potential confounders of age at menarche mainly including body mass index, birth weight, and triglycerides (Additional File 2: Figures 2, 3, 4, 5).
Colocalization analysis
By co-localization analysis (PP.H4 >0.95), we found evidence to support the existence of shared causal variants between ASRM stages 3 and 4 endometriosis with normal delivery and age at last live birth (Additional File 2: Table 8). There was no evidence of co-localization for the remaining four causal effects.
Multivariable Mendelian randomization
We further assessed the association between ASRM stages 3–4 and ovarian endometriosis with normal delivery and age at last live birth. Both ASRM stages 3–4 and ovarian endometriosis were not associated with normal delivery and age at last live birth in the MVMR framework (Additional File 2: Tables 9, 10).
Discussion
In this study, we performed a two-sample MR analysis using large genome-wide association datasets to evaluate putative causal relationships between seven types of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery. We found negative causal effects of ASRM stage 3–4 endometriosis on age at last live birth and on normal delivery, ovarian endometriosis on age at last live birth, and normal delivery, and fallopian tubal endometriosis on excessive irregular menstruation. In addition, bidirectional MR showed that age at menarche had a negative causal effect on intestinal endometriosis. All MR analyses were confirmed by sensitivity analyses, and further, the genetic effects of ASRM stage 3–4 endometriosis on normal delivery and age at last live birth were supported by co-localization evidence. Our findings highlight the genetic association between stages and locations of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery.
Previous efforts, including cohort studies and observational studies, have focused on the effect of the severity of endometriosis on the ability of assisted reproductive technology (ART) to achieve fertility (19–22). We reported an association between endometriosis severity and normal delivery and age at last live birth, with moderate-to-severe endometriosis impairing normal delivery and advancing the age at last live birth. Moreover, our findings supported that menstruation, ovulation, reproductive function, and mode of delivery were not associated with mild endometriosis.
We also investigated the causal effects of endometriosis locations on menstruation, ovulation, reproductive function, and mode of delivery. In this study, our analyses found that ovarian endometriosis showed a greater impact on women, followed by fallopian tube, and endometriosis of pelvic peritoneum, rectovaginal septum and vagina, and intestine were unrelated exposures. In brief, ovarian endometriosis brought forward the age at last live birth and impaired normal delivery, whereas fallopian tubal endometriosis reduced the risk of excessive irregular menstruation.
The association between ovarian endometriosis and age at first live birth was noted by Tuominen et al. They found that the age at first live birth was greatest in women with ovarian endometriosis compared with those with peritoneal and deep endometriosis (23). This is inconsistent with our finding that we did not find such a correlation. First, we compared ovarian endometriosis with non-endometriosis and endometriosis elsewhere, whereas Tuominen et al. compared ovarian endometriosis with peritoneal and deep endometriosis. Second, our study included a more general population, whereas the study by Tuominen et al. was limited to the first live birth before surgical confirmation of endometriosis and excluded women with mild symptoms who were conservatively diagnosed and treated for endometriosis. In addition, Parazzini et al. also suggested that age at first live birth was not associated with endometriosis (24).
Menstrual cycle length has been widely reported to be associated with endometriosis (5). Short menstrual cycles increase the risk of endometriosis (25, 26). However, we did not find an association between menstrual cycle length and endometriosis. In fact, the effect of the menstrual cycle length on endometriosis is inconsistent across studies, and several studies suggested that the menstrual cycle length was not related to endometriosis (27, 28). Similarly, we did not find an association between endometriosis and the number of live births, although it was previously thought that the risk of endometriosis decreased with increasing number of births (24, 28, 29).
In addition to menstrual cycle length, age at menarche has also been reported. The earlier the age at menarche, the higher the risk of endometriosis (25, 30, 31). However, it has also been shown that there was no association between endometriosis and age at menarche (24, 28, 32). We found that the effect of age at menarche on endometriosis may be limited. The earlier the age at menarche, the higher the risk of intestinal endometriosis. This finding may explain the discrepancy between studies.
In addition, negative effects on age at last live birth and normal delivery were observed for both moderate-to-severe endometriosis and ovarian endometriosis. These two types of endometriosis can lead to abnormal delivery and an earlier age at last live birth, shortening reproductive life. This may suggest that clinical attention needs to be paid to the early detection and timely intervention of these two types of endometriosis to ameliorate their adverse effects on women’s labor and reproductive life span.
When both exposures (moderate-to-severe endometriosis and ovarian endometriosis) were considered in the same multivariate model, neither was associated with outcomes, indicating collinearity of moderate-to-severe endometriosis and ovarian endometriosis. In fact, there was substantial overlap between IVs in moderate-to-severe endometriosis and ovarian endometriosis. We speculated that there might be a causal effect between moderate-to-severe endometriosis and ovarian endometriosis, i.e., ovarian endometriosis tends to be more severe. This is beyond the scope of this manuscript, and therefore we did not explore it further, but it could be explored in the future.
The main strength of this study is the utilization of summary-level data on non-overlapping exposures and outcomes in the two-sample MR framework for causal inference. This is a comprehensive analysis of the genetic association between the stages and locations of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery. In addition, we conduct a series of sensitivity analyses, confounding analyses, and co-localization analyses to explore the bias introduced by violating the MR assumption. We also used MVMR to explore specific biases caused by horizontal pleiotropy via moderate-to-severe endometriosis and ovarian endometriosis.
Several limitations of this study should be considered when interpreting the results. First, our findings may have been affected by cohort selection bias, with the population in the sample representing a subset of the total population. In addition, because all genetic variants were derived from genomic studies of European populations, it is not known whether the findings of this study apply to other ethnic populations. Second, endometriosis may occur not only at the five locations described in this manuscript but also at other locations, such as lungs and pleura. We only performed the analyses of the five common locations of endometriosis, and the reason for not analyzing other locations was the lack of corresponding data. Next, due to the lack of IVs, SNPs used for several exposures in the analysis did not meet the traditional GWAS significance threshold (p < 5 × 10−8), which was relaxed to <5 × 10−6. This did not affect the reliability of the study’s conclusions, and our main conclusions are based on strict thresholds and a series of subsequent analyses.
In conclusion, this two-sample MR study provided a comprehensive analysis of the genetic associations between stages and locations of endometriosis and menstruation, ovulation, reproductive function, and mode of delivery. Our study elucidates that moderate-to-severe endometriosis and ovarian endometriosis can lead to abnormal deliveries and a shortened reproductive life span. Both types of endometriosis should be detected promptly and intervened early in clinical practice.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
Only summary level data were used in this study. Ethical approval is available in the publications (6–8).
Author contributions
LS: Data curation, Investigation, Validation, Writing – original draft. JL: Methodology, Visualization, Writing – review & editing. HZ: Conceptualization, Writing – review & editing. YZ: Conceptualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the Basic Research Project for Young and Middle-Aged Physicians of the Beijing Health Promotion Association (BJHPA-2022-SHZHYX ZHQNYJ-JICH-002) and Scientific Research Starting Foundation for Returned Overseas Scholars (2020HGRY004) of Tongji Hospital.
Acknowledgments
We want to acknowledge the participants and investigators of the FinnGen, OpenGWAS, ReproGen, and UK Biobank study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2024.1328403/full#supplementary-material
References
1. Zondervan KT, Becker CM, Missmer SA. Endometriosis. N Engl J Med. (2020) 382:1244–56. doi: 10.1056/NEJMra1810764
2. Rowlands IJ, Abbott JA, Montgomery GW, Hockey R, Rogers P, Mishra GD. Prevalence and incidence of endometriosis in Australian women: A data linkage cohort study. BJOG. (2021) 128:657–65. doi: 10.1111/1471-0528.16447
3. Marquardt RM, Tran DN, Lessey BA, Rahman MS, Jeong J-W. Epigenetic dysregulation in endometriosis: implications for pathophysiology and therapeutics. Endocr Rev. (2023) 44:1074–95. doi: 10.1210/endrev/bnad020
4. Arafah M, Rashid S, Akhtar M. Endometriosis: A comprehensive review. Adv Anat Pathol. (2021) 28:30–43. doi: 10.1097/PAP.0000000000000288
5. Shafrir AL, Farland LV, Shah DK, Harris HR, Kvaskoff M, Zondervan K, et al. Risk for and consequences of endometriosis: A critical epidemiologic review. Best Pract Res Clin Obstet Gynaecol. (2018) 51:1–15. doi: 10.1016/j.bpobgyn.2018.06.001
6. Lyon MS, Andrews SJ, Elsworth B, Gaunt TR, Hemani G, Marcora E. The variant call format provides efficient and robust storage of gwas summary statistics. Genome Biol. (2021) 22:32. doi: 10.1186/s13059-020-02248-0
7. Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner KM, et al. Finngen provides genetic insights from a well-phenotyped isolated population. Nature. (2023) 613:508–18. doi: 10.1038/s41586-022-05473-8
8. Ruth KS, Day FR, Hussain J, Martínez-Marchal A, Aiken CE, Azad A, et al. Genetic insights into biological mechanisms governing human ovarian ageing. Nature. (2021) 596:393–7. doi: 10.1038/s41586-021-03779-7
9. Available online at: http://www.nealelab.is/uk-biobank/.
10. D'Urso S, Arumugam P, Weider T, Hwang L-D, Bond TA, Kemp JP, et al. Mendelian randomization analysis of factors related to ovulation and reproductive function and endometrial cancer risk. BMC Med. (2022) 20:419. doi: 10.1186/s12916-022-02585-w
11. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The mr-base platform supports systematic causal inference across the human phenome. Elife. (2018) 7:e34408. doi: 10.7554/eLife.34408
12. Burgess S, Dudbridge F, Thompson SG. Combining information on multiple instrumental variables in mendelian randomization: comparison of allele score and summarized data methods. Stat Med. (2016) 35:1880–906. doi: 10.1002/sim.6835
13. Burgess S, Davey Smith G, Davies NM, Dudbridge F, Gill D, Glymour MM, et al. Guidelines for performing mendelian randomization investigations. Wellcome Open Res. (2019) 4:186. doi: 10.12688/wellcomeopenres.15555.2
14. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.21758
15. Verbanck M, Chen C-Y, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
16. Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. Phenoscanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. (2019) 35:4851–3. doi: 10.1093/bioinformatics/btz469
17. Giambartolomei C, Vukcevic D, SChadt EE, Franke L, Hingorani AD, Wallace C, et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PloS Genet. (2014) 10:e1004383. doi: 10.1371/journal.pgen.1004383
18. Sanderson E, Spiller W, Bowden J. Testing and correcting for weak and pleiotropic instruments in two-sample multivariable mendelian randomization. Stat Med. (2021) 40:5434–52. doi: 10.1002/sim.9133
19. Horton J, Sterrenburg M, Lane S, Maheshwari A, Li TC, Cheong Y. Reproductive, obstetric, and perinatal outcomes of women with adenomyosis and endometriosis: A systematic review and meta-analysis. Hum Reprod Update. (2019) 25:592–632. doi: 10.1093/humupd/dmz012
20. Soriano D, Adler I, Bouaziz J, Zolti M, Eisenberg VH, Goldenberg M, et al. Fertility outcome of laparoscopic treatment in patients with severe endometriosis and repeated in vitro fertilization failures. Fertil Steril. (2016) 106:1264–9. doi: 10.1016/j.fertnstert.2016.06.003
21. Zimmermann A, Faust C, Miquel L, Berbis J, Perrin J, Courbiere B. Impact of moderate-to-severe endometriosis on ivf cumulative live birth rate: A retrospective matched cohort study. Reprod BioMed Online. (2023) 47:103186. doi: 10.1016/j.rbmo.2023.02.012
22. Huang X-d, Sun Z-Y, Yu Q. Impact of endometriosis on embryo quality and endometrial receptivity in women undergoing assisted reproductive technology. Reprod Biol. (2023) 23:100733. doi: 10.1016/j.repbio.2023.100733
23. Tuominen A, Saavalainen L, Niinimäki M, Gissler M, But A, Härkki P, et al. First live birth before surgical verification of endometriosis-a nationwide register study of 18 324 women. Hum Reprod. (2023) 38:1520–8. doi: 10.1093/humrep/dead120
24. Parazzini F, Ferraroni M, Fedele L, Bocciolone L, Rubessa S, Riccardi A. Pelvic endometriosis: reproductive and menstrual risk factors at different stages in lombardy, northern Italy. J Epidemiol Community Health. (1995) 49:61–4. doi: 10.1136/jech.49.1.61
25. Missmer SA, Hankinson SE, Spiegelman D, Barbieri RL, Malspeis S, Willett WC, et al. Reproductive history and endometriosis among premenopausal women. Obstet Gynecol. (2004) 104:965–74. doi: 10.1097/01.AOG.0000142714.54857.f8
26. 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:389–93. doi: 10.1007/s00404-007-0479-1
27. Blum S, Fasching PA, Hildebrandt T, Lermann J, Heindl F, Born T, et al. Comprehensive characterization of endometriosis patients and disease patterns in a large clinical cohort. Arch Gynecol Obstet. (2022) 305:977–84. doi: 10.1007/s00404-021-06200-w
28. Peterson CM, Johnstone EB, Hammoud AO, Stanford JB, Varner MW, Kennedy A, et al. Risk factors associated with endometriosis: importance of study population for characterizing disease in the endo study. Am J Obstet Gynecol. (2013) 208:451.e1–.11. doi: 10.1016/j.ajog.2013.02.040
29. Sangi-Haghpeykar H, Poindexter AN. Epidemiology of endometriosis among parous women. Obstet Gynecol. (1995) 85:983–92. doi: 10.1016/0029-7844(95)00074-2
30. Nnoaham KE, Webster P, Kumbang J, Kennedy SH, Zondervan KT. Is Early Age at Menarche a Risk Factor for Endometriosis? A Systematic Review and Meta-Analysis Of case-Control Studies. Fertil Steril. (2012) 98:702–12. doi: 10.1016/j.fertnstert.2012.05.035
31. Kim HJ, Lee HS, Kazmi SZ, Hann HJ, Kang T, Cha J, et al. Familial Risk for Endometriosis and Its Interaction with Smoking, Age at Menarche and Body Mass Index: A Population-Based Cohort Study among Siblings. BJOG. (2021) 128:1938–48. doi: 10.1111/1471-0528.16769
Keywords: endometriosis, menstruation, ovulation, reproductive function, mode of delivery, Mendelian randomization
Citation: Shen L, Li J, Zhang H and Zhao Y (2024) Causal effects of endometriosis stages and locations on menstruation, ovulation, reproductive function, and delivery modes: a two-sample Mendelian randomization study. Front. Endocrinol. 15:1328403. doi: 10.3389/fendo.2024.1328403
Received: 26 October 2023; Accepted: 09 July 2024;
Published: 02 August 2024.
Edited by:
Achmad Kemal Harzif, University of Indonesia, IndonesiaReviewed by:
Maria Stamou, Massachusetts General Hospital and Harvard Medical School, United StatesZhibo Tan, Heyou International Hospital, China
Copyright © 2024 Shen, Li, Zhang and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yiqing Zhao, Yiqing_Zhao99@hotmail.com
†These authors have contributed equally to this work