Xinyu Han
Tianqiang Wu1- 1Department of First Clinical Medical College, Heilongjiang University of Chinese Medicine, Harbin, China
- 2Department of Endocrinology, The First Affiliated Hospital of Heilongjiang University of Chinese Medicine, Harbin, China
Background: Numerous observational studies have revealed a correlation between telomere length (TL) and adverse pregnancy outcomes (APOs). However, the impacts of TL on APOs are still unclear.
Methods: Mendelian randomization (MR) was carried out using summary data from genome-wide association studies (GWAS). Inverse variance weighted (IVW) was employed as the primary analysis to explore the causal relationship between TL and APOs. The exposure data came from a GWAS dataset of IEU analysis of the United Kingdom Biobank phenotypes consisting of 472,174 European participants. Summary-level data for five APOs were obtained from the GWAS datasets of the FinnGen consortium. We also performed multivariate MR (MVMR), adjusting for smoking, alcohol intake, body mass index (BMI), and number of live births. In addition, we conducted a series of rigorous analyses to further examine the validity of our MR findings.
Results: After Bonferroni correction and rigorous quality control, univariable MR (UVMR) demonstrated that a shorter TL was significantly associated with an increased risk of spontaneous abortion (SA) (odds ratio [OR]: 0.815; 95% confidence interval [CI]: 0.714–0.930; P = 0.002) and preterm birth (PTB) (OR: 0.758; 95% CI: 0.632-0.908; P = 0.003) in the IVW model. There was a nominally significant relationship between TL and preeclampsia (PE) in the IVW model (OR: 0.799; 95% CI: 0.651-0.979; P = 0.031). However, no significant association was found between TL and gestational diabetes mellitus (GDM) (OR: 0.950; 95% CI: 0.804-1.122; P = 0.543) or fetal growth restriction (FGR) (OR: 1.187; 95% CI: 0.901-1.565; P = 0.223) among the five statistical models. Furthermore, we did not find a significant causal effect of APOs on TL in the reverse MR analysis. MVMR analysis showed that the causal effects of TL on SA remained significant after accounting for smoking, alcohol intake, BMI, and number of live births.
Conclusion: Our MR study provides robust evidence that shorter telomeres were associated with an increased risk of SA. Further work is necessary to investigate the potential mechanisms. UVMR and MVMR findings showed limited evidence that TL affects the risk of PTB, PE, GDM, and FGR, illustrating that the outcomes of previous observational studies may have been confounded.
1 Introduction
Located at the extremities of chromosomes, telomeres are DNA-protein complexes that safeguard genomic stability and integrity by shielding chromosome ends from erosion and fusion (1). With each cell division, telomeres become shorter, which is regarded as an indicator of cellular aging (2). Multiple prior investigations have shown that genetics is the primary determinant of TL (3, 4). Telomere attrition can cause DNA damage, induce cellular senescence and apoptosis, and is a well-known risk factor for a variety of age-related disorders, such as premature ovarian failure, infertility (5), diabetes (6), cardiovascular diseases (7), and neurodegenerative diseases (8). APOs include SA, PTB, PE, GDM, and FGR, which pose a grave threat to maternal, fetal, and neonatal health and significantly cause maternal and infant mortality (9). Given their negative immediate and long-term effects, it is imperative to address APOs. The shortening of telomeres, which can influence the aging of placental and fetal membrane cells, has been associated with aberrant placental aging and APOs (10). Increasing evidence indicates that TL in placentas or fetal membrane may function as a potential biomarker for predicting APOs (11). The length of telomeres can be measured in numerous types of tissue, with minor variations depending on the cell type. Currently, peripheral blood leukocyte TL is the most common type in TL research, possibly due to the convenience of obtaining blood samples. In addition, a recent study evaluated the feasibility of using peripheral blood leukocyte TL to represent other types of tissue TL by detecting the relative telomere length (RTL) of various types of tissues from 952 Genotype-Tissue Expression (GTEx) donors. In 23 categories of tissues, whole blood RTL correlated positively with 15 tissue-specific RTL measurements, demonstrating that peripheral blood leukocyte RTL can represent RTL in the majority of tissues (12). Despite the fact that no study has directly confirmed the correlation between peripheral blood leukocyte TL and placental TL, some studies have found a positive correlation between the TERC mRNA level in term placenta of PE patients at 36-38 weeks and the TERC mRNA level in their peripheral blood Monocyte PBMC (r=0.21, p=0.024) (13). TERC is a gene for telomere maintenance that restores telomerase activity and prolongs cell life. Insufficient expression of this gene can contribute to telomeric lesions, accelerate the accumulation of aging cells, and reduce TL (14). The article examines the possible correlation between TL in peripheral blood and placenta through an indirect approach. Two additional studies also imply the synchronicity of TL in peripheral blood and fetal membranes in an indirect way. Alrefaei et al. discovered that the telomere length of human fetal mesenchymal stem cells (hFM MSCs) in mothers aged 30-39 and 40 was substantially shorter compared to that of young mothers aged 20-29 (15). A study conducted by Nsereko et al. on 297 expectant women from Rwanda revealed a significant negative correlation between maternal age and the TL of peripheral blood leukocytes, consistent with the findings of studies on fetal membrane TL (16). Moreover, similar to TL of placental or fetal membrane, multiple studies have also found an association between maternal peripheral blood leukocyte TL and APOs. A case-control study revealed that couples with idiopathic recurrent pregnancy loss (iRPL) had substantially shorter peripheral blood leukocyte TL than healthy control couples (17). Another prospective study involving 100 expectant Mexican women found that peripheral blood TL during late gestation was substantially lower in women with PTB compared to those with full-term delivery (P=0.02), supporting the association between maternal TL and the incidence of PTB (18). In a case-control study comparing relative TL of genomic DNA extracted from peripheral blood leukocytes in 113 GDM patients and 396 normal pregnant women, GDM patients demonstrated significantly shorter RTL than normal controls (P = 0.046), suggesting that TL is negatively correlated with the incidence of GDM (19). However, a study involving 1,228 women of reproductive age attempting natural conception with a history of 1-2 pregnancy losses revealed that preconception leukocyte TL is unrelated to fecundability, pregnancy loss, or live birth (20). This inconsistency in observational study findings can be attributed to sample size bias and residual confounders. Therefore, the causal relationships between TL and APOs remain ambiguous and require more convincing evidence for validation. Using genetic variations that are closely linked to exposure as instrumental variables (IVs), MR is a reliable method for determining the causality between exposure and outcome (21). Similar to a randomized controlled trial, the MR design is less susceptible to bias from confounding factors and reverse causality due to the random assignment of genetic variants at conception (22, 23). MVMR is an emerging technique that allows simultaneous assessment of relevant exposures by incorporating genetic variants of multiple risk factors into the same model to minimize the impact of confounding variables (24). Therefore, UVMR and MVMR analyses were conducted to ascertain the causal relationships between TL and the five prevalent APOs.
2 Materials and methods
2.1 Study design
We performed a two-sample MR analysis to assess the causal relationship of TL with APOs using publicly available data. Figure 1 illustrates the study’s methodology. The genetic variants selected to estimate the causal effect must satisfy three key assumptions (25): Assumption 1, genetic variants should be strongly associated with TL (P < 5 × 10-8); Assumption 2, The association between genetic variants of TL and APOs is independent of confounding factors; Assumption 3, Genetic variants affect the risk of APOs directly through TL, not through other pathways. Evidence from previous observational clinical trials and MR studies indicated that smoking, alcohol intake, BMI, and number of live births are risk factors for the development of APOs (26–31). Therefore, we further conducted MVMR adjusting for genetic liability to smoking, alcohol intake, BMI, and number of live births. Since our data were derived from publicly accessible GWAS summary statistics, no ethical approval was necessary.
Figure 1 Assumptions of the Mendelian randomization (MR) analysis for TL and the risk of APOs.
2.2 Data sources
The summary statistics for TL were derived from the largest GWAS on TL to date (dataset ID: ieu-b-4879), which included 472,174 European participants from the UK Biobank (UKB), with a comparable distribution of males (45.8%) and females (54.2%) (32). Quantitative polymerase chain reaction (PCR) analyses were performed to obtain Leukocyte TL measurements for these participants from the UKB (33). GWAS summary data for SA (9,113 cases and 89,340 controls), PTB (5,480 cases and 98,626 controls), PE (3,903 cases and 114,735 controls), GDM (5,687 cases and 123,579 controls), and FGR (2,579 cases and 171,167 controls) were obtained from the FinnGen consortium (34). All participants in this study are of European descent. In addition, we gained aggregate data regarding smoking, alcohol intake, BMI, and number of live births from the Neale Lab or MRC-IEU consortium. All GWAS data samples originate from the European populations. Table 1 provides a summary of all datasets included in this investigation.
Table 1 Details of studies included in Mendelian randomization (MR) analyses.
2.3 Selection and evaluation of instrumental variable
To satisfy the three assumptions of the MR analysis, we selected the IVs through the following procedure. Initially, all SNPs substantially associated with TL (P < 5 × 10 -8) were chosen as IVs. When SNPs that strongly predicted APOs were extracted at the genome-wide significance level (P < 5 × 10 -8), the number of available SNPs was low or even absent, consequently, a cut off (P < 5 × 10 -6) was adopted to obtain SNPs that predicting APOs in the reverse MR analysis. The corresponding linkage disequilibrium was then evaluated to confirm that there were SNPs in a linkage disequilibrium state and that the SNPs were independent by removing SNPs from a 10,000-kb window with an r2 < 0.001 thresholds. Thirdly, putative pleiotropic effects were eliminated by retrieving the secondary phenotype of each SNP from PhenoScan V2 (35). SNPs corresponding to the phenotypes associated with the results were excluded from further analysis, while the remaining SNPs were utilized.
Variance (R2) and F-statistic were employed to assess the robustness of IVs to avoid weak tool bias. We calculated the F-statistic for each SNP by adopting the formula: F = R2/(1-R2) [(N-K-1)/K], where N is the sample size, k is the total number of SNPs selected for MR analysis, and R2 is the total proportion of phenotypic variations explained by all the SNPs in our MR model (36). The following formula was used to calculate the R2 for each SNP: R2 = Σ [2 × (1 – MAF) × MAF × β2/(SE2 × N)], where SE and β are the standard error and β coefficient for effect size, MAF is the minor allele frequency for each SNP (36). An F-statistic larger than 10 was deemed significant enough for the association between IVs and exposure to prevent the MR analyses from being influenced by weak tool bias (25). Statistical power for each outcome was calculated by utilizing the online tool (https://shiny.cnsgenomics.com/mRnd/) (37). A sufficient power of over 80% was recommended.
2.4 Statistical analysis
The IVW method was employed as the primary analytical method for estimating causal effects between TL and APOs. Additional analysis methodologies included MR-Egger, weighted median, weighted mode, and simple mode. The IVW method is an extension of the Wald ratio estimator based on meta-analytic principles that can provide an accurate estimate in an ideal state where all included SNPs are presumed to be valid IVs without pleiotropy (38). MR-Egger permits certain SNPs to influence the outcome through mechanisms other than exposure modification (39). In addition, the MR-Egger intercept can detect and adjust pleiotropy (39). To employ the weighted median method, at least 50 percent of intravenous solutions must satisfy the assumption that they are valid intravenous solutions (40). The weighted model method clusters SNPs and calculates estimates based on the cluster containing the greatest number of SNPs (39). Lastly, even though the simple model method is less effective than IVW, it provides robustness for pleiotropy (41). Based on prior research (26–31), we adjusted for smoking, alcohol intake, BMI, and number of live births in multivariate MR to demonstrate a causal effect of TL on APOs. The methods we applied to execute MVMR included IVW and MR-Lasso (42).
Various sensitivity analyses were implemented in this study to ensure the stability and reliability of the MR results. First, Cochran’s Q test was applied to evaluate the heterogeneity between SNPs, where a p-value greater than 0.05 indicated no heterogeneity. Secondly, the MR-Egger intercept was utilized to quantify the horizontal pleiotropy of IVs. Thirdly, we performed the leave‐one‐out analysis to check whether the MR results were driven by any single SNP. Finally, we conducted the MR-PRESSO to detect potential outlier SNPs (42). A two-sided p-value < 0.05 was considered statistically significant. We modified the p-value by Bonferroni correction for the number of outcomes. For the primary analysis (association of TL with five APOs), the association with a two-sided p-value < 0.01 (where α = 0.05/5 outcomes) was deemed statistically significant and a two-sided p-value < 0.05 was thought suggestive. MR analyses were conducted using the TwoSampleMR (version 0.5.6) and MVMR (version 0.3) in R. All data analyses were performed with R version 4.2.2.
3 Result
3.1 Genetic instruments
In the present study, 117, 107, 109, 108, and 93 SNPs were eventually obtained as the IVs for TL to assess the associations between TL and SA, PTB, PE, GDM, and FGR, respectively (Supplementary Tables S1–S5). The F statistic for all these genetic variants was above the threshold of 100, indicating a low likelihood of weak instrumental bias.
The statistical power of SA, PTB, and PE was above 80 percent, which is sufficient to prove the reliability of the results. However, the statistical power of GDM and FGR is less than 80%, which may lead to false negatives. The detailed results of the selected IV’s strength and statistical power are shown in Table 2. For the MR analysis of APOs on TL, 15, 9, 7, 10, and 10 SNPs associated with SA, PTB, PE, GDM, and FGR were selected respectively (Supplementary Tables S6–S10). However, the F statistic of most IVs was less than 10, suggesting the possibility of weak instrumental bias.
Table 2 The results of the selected IVs strength and statistical power.
3.2 Estimated causal effect of TL on APOs
After Bonferroni correction, we observed that a shorter TL was significantly associated with a higher risk of SA (OR: 0.815; 95%CI: 0.714-0.930; P = 0.002) and PTB (OR: 0.758; 95% CI: 0.632-0.908; P = 0.003) in IVW model. There was a nominally significant relationship between TL and PE in the IVW model (OR: 0.799; 95% CI: 0.651-0.979; P = 0.031). However, no significant association was found between TL and GDM (OR: 0.950; 95% CI: 0.804-1.122; P = 0.543) or FGR (OR: 1.187; 95% CI: 0.901-1.565; P = 0.223) among the five statistical models. The causal associations of genetically predicted TL and the risk of APOs are presented in Figure 2. The scatter plots of the association between TL and APOs are shown in Figures 3, 4. The Cochran’s Q test showed no heterogeneity, and the MR-Egger intercept test found no evidence of horizontal pleiotropy in the MR analysis results. Additionally, the MR-PRESSO results showed no outlier SNPs. Table 3 provides details on the sensitivity analysis results. The leave-one-out plots further support the robustness of our results and suggest that the effects of any single SNP were unlikely to influence causal estimates (Supplementary Figures S1–S5).
Figure 2 Association of genetically predicted telomere length (TL) and the risk of adverse pregnancy outcomes (APOs). *Indicates that the relationship has nominal statistical significance; **Indicates that the P-value meets the Bonferroni correction threshold.
Figure 3 Scatter plots for Mendelian randomization (MR) analyses of the correlation between telomere length (TL) and adverse pregnancy outcomes (APOs). (A) SA; (B) PE; (C) PTB.
Figure 4 Scatter plots for Mendelian randomization (MR) analyses of the correlation between telomere length (TL) and adverse pregnancy outcomes (APOs). (A) GDM; (B) FGR.
Table 3 Heterogeneity, horizontal pleiotropy, and MR-PRESSO tests of the associations between TL and APOs.
In the multivariable MR analysis, adjusting for smoking, alcohol intake, BMI, and number of live births, the magnitude of associations between genetic liability to TL and SA diminished, whereas the causal relationship between TL and SA (OR: 0.867; 95%CI: 0.763-0.985; P = 0.028) remained significant. After removing heterogeneous SNPs, the results of the MR-Lasso test remained constant. The F-statistics were 42.50, indicating that there is no potential mild instrument bias. The association between TL and PTB or PE did not persist, however, after adjusting for smoking, alcohol intake, BMI, and number of live births. Table 4 presents the MVMR results in detail.
Table 4 Causal estimates of TL on SA, PTB and PE in MVMR.
3.3 Estimated causal effect of APOs on TL
Genetic susceptibility to SA (OR: 1.008; 95% CI: 0.987-1.028, P = 0.418), PTB(OR: 0.996; 95% CI: 0.983-1.009, P = 0.521), PE(OR: 1.004; 95% CI: 0.992-1.017, P = 0.478), GDM(OR: 1.000; 95% CI: 0.987-1.013; P = 0.980), or FGR(OR: 1.002; 95% CI: 0.993-1.011; P = 0.710) was all not causally related to TL in all five statistical models (Figure 5). In the sensitivity analyses, no evidence for directional pleiotropy was found when we reanalyzed the five results using MR-Egger regression. Heterogeneity was assessed by Cochrane’s Q test, and the results of IVW and MR Egger analyses revealed heterogeneity only for SA exposure and not for any other APOs. As a result, a random effects model was employed to evaluate the causal relationship between SA and TL. For five APOs, however, no outliers were identified in the MR-PRESSO model. Table 5 provides details on the sensitivity analysis results for the reverse MR study. In addition, leave-one-out plots demonstrated that individual SNPs were not anticipated to influence causal estimates (Supplementary Figures S6–S10). In addition, scatter plots (Supplementary Figures S11, S12) displaying the effect size of each SNP on TL were provided.
Figure 5 The causal relationship of genetically predicted APOs and TL.
Table 5 Heterogeneity, horizontal pleiotropy, and MR-PRESSO tests of the associations between APOs and TL.
4 Discussion
In this investigation, we examined the causal effect of TL on five APOs utilizing MR analyses. We discovered that genetically predicted a shorter TL was associated with an increased risk of certain APOs. The significant correlations were found between TL and SA or PTB After the Bonferroni correction. There was a suggestive association between TL and PE. The casual associations between TL and SA remained significant after controlling for smoking, alcohol intake, BMI, and number of live births. In the MVMR models, the correlation between TL and PTB or PE did not persist.
Our findings that a shortened TL was correlated with an increased risk of SA were consistent with those of prior research. For instance, a recent case-control study found that women with iRPL (n=40) had significantly shorter peripheral blood leukocyte TL in early follicular stages than healthy women (n=41) with a comparable age and body mass index (P=0.006) (43). In studies of infertile women undergoing in vitro fertilization (IVF), shorter peripheral blood leukocyte TL was linked with a higher incidence of embryonic aneuploidy, a major contributor to pregnancy loss (P=0.01) (44). All of these findings support the protective function of a longer TL in SA and suggest that special attention should be paid to patients with a shorter TL along with certain measures should be taken to prevent the development of SA.
A recently published study comparing maternal blood and placental tissue telomere length in PTB and uncomplicated full-term pregnancy (NTP) revealed that maternal blood TL was significantly shorter in PTB (n=11) than in NTP (n=19) and that the proportion of trophoblast cells with shortened telomeres was significantly higher in PTB placental samples (P<0.001) (45). In addition, another cohort study found that shorter peripheral blood TL in early pregnancy was significantly associated with an earlier gestational age at delivery (r = 0.35, P = 0.02). This association persisted after accounting for maternal age in a linear regression model (OR=0.30, 95% CI: 0.04-0.56, P=0.03) (10). Our univariate MR results indicated that shortened TL was significantly associated with the incidence of PTB; however, after adjusting for smoking, alcohol intake, BMI, and number of live births, the correlation became nonsignificant. The discrepancy between the results of MR analysis and observational studies may be due to interference by confounding factors.
Prior investigations have reached contradictory conclusions regarding the relationship between TL and PE, with one study reporting a shorter TL in placenta samples from expectant women with both TL and PE (46). The proportion of trophoblast cells with short TL was significantly higher in early (n=7) and late (n=6) pre-eclamptic placenta samples than in healthy controls (n=13; P=0.03), according to a separate study (47). In contrast, a case-control study from China demonstrated that a lengthier leukocyte TL in peripheral blood was associated with PE (48). A study found no statistically significant difference in telomere length between maternal blood and placental samples from patients with PE at 31 to 40 weeks of gestation (n=31) and normal pregnancy controls (n=30) (49). This is consistent with our findings, in which a nominally significant association was found between TL shortening and PE incidence in our univariate MR results, but this relationship became nonsignificant when smoking, alcohol intake, BMI, and number of live births were considered in a multivariable model. These findings suggest that the relationship between TL length and PE is divergent or even contradictory, requiring further investigation.
Our MR investigation identified no correlation between peripheral blood TL and either GDM or FGR. The relative TL of genomic DNA extracted from peripheral blood leukocytes was substantially shorter in patients with GDM than in controls (p = 0.046) in a previous case-control study that did not adjust for confounding variables (19). After adjusting for age, another prospective cohort study found no significant correlation between TL and GDM risk (50). This corresponds with our findings. Numerous studies have found a link between shortened placental trophoblast telomere length and FGR (51–53), but none have explored the relation between peripheral blood TL and FGR. Our study found no correlation between peripheral blood TL and FGR. The discrepancy between this result and those of observational studies may be due to the influence of confounding variables on the one hand, and the inability of peripheral blood leukocyte TL to fully reflect placental TL on the other hand, despite the fact that TL is typically consistent across tissues (12). The relationship between leukocyte TL and placental TL and its association with FGR requires additional study.
The results of our reverse MR analysis did not provide genetic evidence of a causal relationship between APOs and TL. There is currently no observational study examining the relationship between APOs and TL; however, studies have found that women with APOs, such as SA, PE, and PTB, have a higher risk of cardiovascular diseases closely related to short telomeres (7, 54, 55), suggesting that APOs may be correlated with TL shortening. This inference is not confirmed by the results of our reverse MR analysis. However, there is a weak instrumental bias in this result, which may contribute to an underestimation of the correlation between exposure and outcome in the two-sample MR study. Insufficient sample size is the primary contributor to this bias. Therefore, we should explain this result with caution and employ larger sample sizes of SA-related GWAS data for verification in the future.
The precise mechanism underlying the connection between TL and APOs is still unknown. Despite this, several studies have pointed out that accelerated TL shortening leads to increased oxidative stress and immune-inflammatory responses, both of which are known to play crucial roles in the development of APOs (56, 57). Telomere attrition may promote placental and fetal membrane cell senescence, which may contribute to the development of APOs (58). Moreover, an earlier study reveals that TL can indicate ovarian reserve function in women. In particular, the age of natural menopause is delayed by approximately 10.2 months for every 1 kilobase (kb) increase in leukocyte TL (59). Reduced ovarian reserve function is a risk factor for miscarriage (60).
The following is a discussion on the potential mechanism of a causal relationship between TL and SA at the genetic level. TL is primarily determined by genes, and previous research has demonstrated that genetic factors can account for 36% to 84% of TL variation (61). Certain gene mutations or functional abnormalities may result in telomere shortening. For example, Poly ADP Ribose Polymerase1 (PARP1), Poly ADP Ribose Polymerase1 (PARP2), and Telomeric Repeat Binding Factor 2 (TERF2/TRF2) (62). PARP1 and PARP2 are ADP ribose transferase (ART) family members. ART is a multifunctional protein Post-translational modification enzyme that participates in numerous cellular processes, such as DNA damage repair, lipid metabolism, immune response, transcriptional regulation, and cell death and plays a crucial role in maintaining gene stability and TL (62). In addition, PARP can prevent telomere degradation and lengthen telomeres by regulating telomerase reverse transcriptase gene transcription (63). Studies have demonstrated that catalytic inhibition or gene deletion of PARP-1 and PARP-2 in the mouse uterus can increase p53 signaling and the population of senescent decidual cells, causing decidualization failure and pregnancy loss in rodents (64). Additionally, PARP-2 was found to be upregulated in the endometrial stromal area (implantation area) during the receptive stage but downregulated during pregnancy failure and pseudopregnancy. Experiments in vitro have demonstrated that PARP-2 regulates endometrial receptivity by taking part in cysteine protease 8-dependent inflammation and apoptosis (65). TERF2/TRF2 is a multifunctional telomeric protein that protects the extremities of chromosomes by binding to repetitive sequences on telomere DNA and forming structurally stable telomere caps. The expression of TERF2/TRF2 in the abortion tissue of women with iRPL is substantially reduced compared to that of induced abortion patients, and the TL of iRPL is significantly shorter. This suggests that the downregulation of TERF2/TRF2 expression results in telomere uncapping, which plays a crucial role in pregnancy loss (66). Other studies have demonstrated that TERF2 can regulate autophagy by binding to High Mobility Group Box 1 (HMGB1), a non-histone Chromatin-related protein (67). In recent years, a large number of studies have suggested that abnormalities in autophagy may be associated with SA by modulating decidualization, trophoblast cells, and immune cells at the maternal fetal interface (68). In our study, SNPs used as IVs include rs139795227, rs3093888, and rs3785074, which correspond to the PARP1, PARP2, and TERF2/TRF2 genes, respectively. Therefore, we conclude that genes associated with telomere shortening may be involved in the pathogenesis of SA via inflammatory response, apoptosis, autophagy, and other mechanisms.
According to our knowledge, this is the first time the MR framework was employed for estimating the genetic causality between TL and APOs. This MR study has several strengths that should be noted. First, we eliminated genetic variants linked with potential confounders frequently observed in epidemiological studies and selected only SNPs strongly associated with TL. Second, the large sample size of our MR analysis boosted our statistical power and provided solid evidence of the existence of relationships. Thirdly, we conducted multiple sensitivity analyses to validate the reliability of these results. Lastly, we utilized MVMR to investigate the direct impact of TL on APOs after adjusting for smoking, alcohol intake, BMI, and number of live birth. Despite the benefit, there are some restrictions. First, as a result of using summary-level data from the GWAS database, we were unable to evaluate the nonlinear correlations between TLs and APOs. Second, the GWAS dataset for TL contained 45.8% male and 54.2% female participants, which may introduce bias. Due to a lack of information for the female subgroup, we were unable to examine our results stratified according to the sex classification of TL. We anticipate that sex-disaggregated TL data will be utilized in the future to investigate the relationship between TL and APOs in greater depth. Aside from that, the majority of participants in this study were of European descent, reducing population stratification bias but limiting the applicability of our findings to other populations. Finally, Even though the F-statistic confirmed that no weak IVs existed, we observed low statistical power (below 80 percent) for certain phenotypes, which could result in false negatives.
5 Conclusion
Our MR study provides strong evidence that genetically predicted shortened telomeres are associated with a higher risk of SA. The potential mechanisms require further examination. UVMR and MVMR findings showed limited evidence that TL affects the risk of PTB, PE, GDM, and FGR, illustrating that the outcomes of previous observational studies may have been confounded.
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
The study used summary-level data from publicly available datasets, which were not collected at the individual level. All of the participants provided written informed consent in each of the contributing studies. Therefore, ethical approval was not obtained.
Author contributions
XH: Conceptualization, Investigation, and Writing – original draft. TW: Investigation and Resources – original draft. CL: Revising manuscript, responding to comments, and polishing language. All authors contributed to the article and approved the submitted version.
Acknowledgments
We would like to thank the researchers and study participants for their contributions.
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.2023.1225600/full#supplementary-material
References
1. Chakravarti D, LaBella KA, DePinho RA. Telomeres: history, health, and hallmarks of aging. Cell (2021) 184:306–22. doi: 10.1016/j.cell.2020.12.028
2. Srinivas N, Rachakonda S, Kumar R. Telomeres and telomere length: A general overview. Cancers (Basel) (2020) 12:558. doi: 10.3390/cancers12030558
3. Andrew T, Aviv A, Falchi M, Surdulescu GL, Gardner JP, Lu X, et al. Mapping genetic loci that determine leukocyte telomere length in a large sample of unselected female sibling pairs. Am J Hum Genet (2006) 78:480–6. doi: 10.1086/500052
4. Jeanclos E, Schork NJ, Kyvik KO, Kimura M, Skurnick JH, Aviv A. Telomere length inversely correlates with pulse pressure and is highly familial. Hypertension (2000) 36:195–200. doi: 10.1161/01.hyp.36.2.195
5. Fattet A-J, Toupance S, Thornton SN, Monnin N, Guéant J-L, Benetos A, et al. Telomere length in granulosa cells and leukocytes: a potential marker of female fertility? A systematic review of the literature. J Ovarian Res (2020) 13:96. doi: 10.1186/s13048-020-00702-y
6. Cheng F, Luk AO, Shi M, Huang C, Jiang G, Yang A, et al. Shortened leukocyte telomere length is associated with glycemic progression in type 2 diabetes: A prospective and mendelian randomization analysis. Diabetes Care (2022) 45:701–9. doi: 10.2337/dc21-1609
7. Deng Y, Li Q, Zhou F, Li G, Liu J, Lv J, et al. Telomere length and the risk of cardiovascular diseases: A mendelian randomization study. Front Cardiovasc Med (2022) 9:1012615. doi: 10.3389/fcvm.2022.1012615
8. Rodríguez-Fernández B, Gispert JD, Guigo R, Navarro A, Vilor-Tejedor N, Crous-Bou M. Genetically predicted telomere length and its relationship with neurodegenerative diseases and life expectancy. Comput Struct Biotechnol J (2022) 20:4251–6. doi: 10.1016/j.csbj.2022.08.006
9. Rana S, Lemoine E, Granger JP, Karumanchi SA. Preeclampsia: pathophysiology, challenges, and perspectives. Circ Res (2019) 124:1094–112. doi: 10.1161/CIRCRESAHA.118.313276
10. Panelli DM, Diwan M, Cruz GI, Leonard SA, Chueh J, Gotlib IH, et al. An exploratory analysis of leukocyte telomere length among pregnant and non-pregnant people. Brain Behav Immun Health (2022) 25:100506. doi: 10.1016/j.bbih.2022.100506
11. Kohlrausch FB, Keefe DL. Telomere erosion as a placental clock: From placental pathologies to adverse pregnancy outcomes. Placenta (2020) 97:101–7. doi: 10.1016/j.placenta.2020.06.022
12. Demanelis K, Jasmine F, Chen LS, Chernoff M, Tong L, Delgado D, et al. Determinants of telomere length across human tissues. Science (2020) 369:eaaz6876. doi: 10.1126/science.aaz6876
13. Lekva T, Roland MCP, Estensen ME, Norwitz ER, Tilburgs T, Henriksen T, et al. Dysregulated non-coding telomerase RNA component and associated exonuclease XRN1 in leucocytes from women developing preeclampsia-possible link to enhanced senescence. Sci Rep (2021) 11:19735. doi: 10.1038/s41598-021-99140-z
14. Lu W, Zhang Y, Liu D, Songyang Z, Wan M. Telomeres-structure, function, and regulation. Exp Cell Res (2013) 319:133–41. doi: 10.1016/j.yexcr.2012.09.005
15. Alrefaei GI, Alkarim SA, Abduljabbar HS. Impact of mothers' age on telomere length and human telomerase reverse transcriptase expression in human fetal membrane-derived mesenchymal stem cells. Stem Cells Dev (2019) 28:1632–45. doi: 10.1089/scd.2019.0144
16. Nsereko E, Uwase A, Muvunyi CM, Rulisa S, Ntirushwa D, Moreland P, et al. Association between micronutrients and maternal leukocyte telomere length in early pregnancy in rwanda. BMC Pregnancy Childbirth (2020) 20:692. doi: 10.1186/s12884-020-03330-y
17. Thilagavathi J, Mishra SS, Kumar M, Vemprala K, Deka D, Dhadwal V, et al. Analysis of telomere length in couples experiencing idiopathic recurrent pregnancy loss. J Assist Reprod Genet (2013) 30:793–8. doi: 10.1007/s10815-013-9993-1
18. Page RL, Han G, Akinlotan M, Patron MP, Gandhi H, Kochan KJ. Telomere length and preterm birth in pregnant mexican-origin women. Matern Child Health J (2021) 25:1798–805. doi: 10.1007/s10995-021-03209-0
19. Weng Q, Deng K, Wu F, Gan M, Li J, Dai Y, et al. Leukocyte telomere length, lipid parameters and gestational diabetes risk: a case-control study in a Chinese population. Sci Rep (2019) 9:8483. doi: 10.1038/s41598-019-44968-9
20. Purdue-Smithe AC, Kim K, Andriessen VC, Pollack AZ, Sjaarda LA, Silver RM, et al. Preconception leukocyte telomere length and pregnancy outcomes among women with demonstrated fecundity. Hum Reprod (2021) 36:3122–30. doi: 10.1093/humrep/deab201
21. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet (2014) 23:R89–98. doi: 10.1093/hmg/ddu328
22. Zheng J, Baird D, Borges M-C, Bowden J, Hemani G, Haycock P, et al. Recent developments in mendelian randomization studies. Curr Epidemiol Rep (2017) 4:330–45. doi: 10.1007/s40471-017-0128-6
23. Sheng J, Liu J, Chan KHK. Evaluating the causal effects of gestational diabetes mellitus, heart disease, and high body mass index on maternal alzheimer’s disease and dementia: multivariable mendelian randomization. Front Genet (2022) 13:833734. doi: 10.3389/fgene.2022.833734
24. Sanderson E. Multivariable mendelian randomization and mediation. Cold Spring Harb Perspect Med (2021) 11:a038984. doi: 10.1101/cshperspect.a038984
25. Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ (2018) 362:k601. doi: 10.1136/bmj.k601
26. Yuan S, Liu J, Larsson SC. Smoking, alcohol and coffee consumption and pregnancy loss: a Mendelian randomization investigation. Fertil Steril (2021) 116:1061–7. doi: 10.1016/j.fertnstert.2021.05.103
27. Nanninga EK, Weiland S, Berger MY, Feijen-de Jong EI, Erwich JJHM, Peters LL. Adverse maternal and infant outcomes of women who differ in smoking status: E-cigarette and tobacco cigarette users. Int J Environ Res Public Health (2023) 20:2632. doi: 10.3390/ijerph20032632
28. O’Leary CM, Nassar N, Kurinczuk JJ, Bower C. The effect of maternal alcohol consumption on fetal growth and preterm birth. BJOG (2009) 116:390–400. doi: 10.1111/j.1471-0528.2008.02058.x
29. Bar-Zeev Y, Haile ZT, Chertok IA. Association between prenatal smoking and gestational diabetes mellitus. Obstet Gynecol (2020) 135:91–9. doi: 10.1097/AOG.0000000000003602
30. Sun Y, Shen Z, Zhan Y, Wang Y, Ma S, Zhang S, et al. Effects of pre-pregnancy body mass index and gestational weight gain on maternal and infant complications. BMC Pregnancy Childbirth (2020) 20:390. doi: 10.1186/s12884-020-03071-y
31. Qiu J, Du T, Li W, Zhao M, Zhao D, Wang Y, et al. Impact of recurrent pregnancy loss history on reproductive outcomes in women undergoing fertility treatment. Am J Obstet Gynecol (2023) 228:66.e1–9. doi: 10.1016/j.ajog.2022.08.014
32. Codd V, Denniff M, Swinfield C, Warner SC, Papakonstantinou M, Sheth S, et al. Measurement and initial characterization of leukocyte telomere length in 474,074 participants in UK biobank. Nat Aging (2022) 2:170–9. doi: 10.1038/s43587-021-00166-9
33. Welsh S, Peakman T, Sheard S, Almond R. Comparison of DNA quantification methodology used in the DNA extraction protocol for the UK biobank cohort. BMC Genomics (2017) 18:26. doi: 10.1186/s12864-016-3391-x
34. 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
35. 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
36. Papadimitriou N, Dimou N, Tsilidis KK, Banbury B, Martin RM, Lewis SJ, et al. Physical activity and risks of breast and colorectal cancer: a mendelian randomisation analysis. Nat Commun (2020) 11:597. doi: 10.1038/s41467-020-14389-8
37. Brion M-JA, Shakhbazov K, Visscher PM. Calculating statistical power in Mendelian randomization studies. Int J Epidemiol (2013) 42:1497–501. doi: 10.1093/ije/dyt179
38. Pagoni P, Dimou NL, Murphy N, Stergiakouli E. Using Mendelian randomisation to assess causality in observational studies. Evid Based Ment Health (2019) 22:67–71. doi: 10.1136/ebmental-2019-300085
39. Hemani G, Bowden J, Davey Smith G. Evaluating the potential role of pleiotropy in mendelian randomization studies. Hum Mol Genet (2018) 27:R195–R208. doi: 10.1093/hmg/ddy163
40. 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:304–14. doi: 10.1002/gepi.21965
41. Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol (2017) 46:1985–98. doi: 10.1093/ije/dyx102
42. 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
43. Onat T, Çaltekin MD, Inandiklioglu N, Baser E, Kirmizi DA, Kara M, et al. Telomere length in idiopathic recurrent pregnancy loss. Z Geburtshilfe Neonatol (2021) 225:119–24. doi: 10.1055/a-1345-9821
44. Hanson BM, Tao X, Zhan Y, Kim JG, Klimczak AM, Herlihy NS, et al. Shorter telomere length of white blood cells is associated with higher rates of aneuploidy among infertile women undergoing in vitro fertilization. Fertil Steril (2021) 115:957–65. doi: 10.1016/j.fertnstert.2020.09.164
45. Farladansky-Gershnabel S, Dekel N, Biron-Shental T, Shechter-Maor G, Amiel A, Weisz A, et al. Spontaneous preterm birth: Elevated galectin-3 and telomere shortening may reflect a common pathway of enhanced inflammation and senescence. Reprod Sci (2023) 30:487–93. doi: 10.1007/s43032-022-00972-5
46. Biron-Shental T, Sukenik-Halevy R, Sharon Y, Goldberg-Bittman L, Kidron D, Fejgin MD, et al. Short telomeres may play a role in placental dysfunction in preeclampsia and intrauterine growth restriction. Am J Obstet Gynecol (2010) 202:381.e1–7. doi: 10.1016/j.ajog.2010.01.036
47. Farladansky-Gershnabel S, Gal H, Kidron D, Krizhanovsky V, Amiel A, Sukenik-Halevy R, et al. Telomere homeostasis and senescence markers are differently expressed in placentas from pregnancies with early- versus late-onset preeclampsia. Reprod Sci (2019) 26:1203–9. doi: 10.1177/1933719118811644
48. Zhang R, Du J, Xiao Z, Jiang Y, Jin L, Weng Q. Association between the peripartum maternal and fetal telomere lengths and mitochondrial DNA copy numbers and preeclampsia: a prospective case-control study. BMC Pregnancy Childbirth (2022) 22:483. doi: 10.1186/s12884-022-04801-0
49. Broady AJ, Loichinger MH, Ahn HJ, Davy PMC, Allsopp RC, Bryant-Greenwood GD. Protective proteins and telomere length in placentas from patients with pre-eclampsia in the last trimester of gestation. Placenta (2017) 50:44–52. doi: 10.1016/j.placenta.2016.12.018
50. Lin Y, Zhu Y, Wu J, Hinkle SN, Rawal S, Han J, et al. A prospective study of leukocyte telomere length and risk of gestational diabetes in a multiracial cohort. Epidemiology (2019) 30 Suppl 2:S10–6. doi: 10.1097/EDE.0000000000001081
51. Niu Z, Li K, Xie C, Wen X. Adverse birth outcomes and birth telomere length: A systematic review and meta-analysis. J Pediatr (2019) 215:64–74.e6. doi: 10.1016/j.jpeds.2019.08.040
52. Toutain J, Prochazkova-Carlotti M, Cappellen D, Jarne A, Chevret E, Ferrer J, et al. Reduced placental telomere length during pregnancies complicated by intrauterine growth restriction. PloS One (2013) 8:e54013. doi: 10.1371/journal.pone.0054013
53. Biron-Shental T, Sukenik-Halevy R, Sharon Y, Laish I, Fejgin MD, Amiel A. Telomere shortening in intra uterine growth restriction placentas. Early Hum Dev (2014) 90:465–9. doi: 10.1016/j.earlhumdev.2014.06.003
54. Kyriacou H, Al-Mohammad A, Muehlschlegel C, Foster-Davies L, Bruco MEF, Legard C, et al. The risk of cardiovascular diseases after miscarriage, stillbirth, and induced abortion: a systematic review and meta-analysis. Eur Heart J Open (2022) 2:oeac065. doi: 10.1093/ehjopen/oeac065
55. Cho GJ, Um JS, Kim SJ, Han SW, Lee SB, Oh M-J, et al. Prior pregnancy complications and maternal cardiovascular disease in young Korean women within 10 years after pregnancy. BMC Pregnancy Childbirth (2022) 22:229. doi: 10.1186/s12884-022-04578-2
56. Barnes RP, Fouquerel E, Opresko PL. The impact of oxidative DNA damage and stress on telomere homeostasis. Mech Ageing Dev (2019) 177:37–45. doi: 10.1016/j.mad.2018.03.013
57. Córdova-Oriz I, Chico-Sordo L, Varela E. Telomeres, aging and reproduction. Curr Opin Obstet Gynecol (2022) 34:151–8. doi: 10.1097/GCO.0000000000000779
58. Polettini J, da Silva MG. Telomere-related disorders in fetal membranes associated with birth and adverse pregnancy outcomes. Front Physiol (2020) 11:561771. doi: 10.3389/fphys.2020.561771
59. Gray KE, Schiff MA, Fitzpatrick AL, Kimura M, Aviv A, Starr JR. Leukocyte telomere length and age at menopause. Epidemiology (2014) 25:139–46. doi: 10.1097/EDE.0000000000000017
60. Bunnewell SJ, Honess ER, Karia AM, Keay SD, Al Wattar BH, Quenby S. Diminished ovarian reserve in recurrent pregnancy loss: a systematic review and meta-analysis. Fertil Steril (2020) 113:818–827.e3. doi: 10.1016/j.fertnstert.2019.11.014
61. Aviv A. Genetics of leukocyte telomere length and its role in atherosclerosis. Mutat Res (2012) 730:68–74. doi: 10.1016/j.mrfmmm.2011.05.001
62. Muoio D, Laspata N, Fouquerel E. Functions of ADP-ribose transferases in the maintenance of telomere integrity. Cell Mol Life Sci (2022) 79:215. doi: 10.1007/s00018-022-04235-z
63. Ghosh U, Das N, Bhattacharyya NP. Inhibition of telomerase activity by reduction of poly(ADP-ribosyl)ation of TERT and TEP1/TP1 expression in HeLa cells with knocked down poly(ADP-ribose) polymerase-1 (PARP-1) gene. Mutat Res (2007) 615:66–74. doi: 10.1016/j.mrfmmm.2006.10.002
64. Kelleher AM, Setlem R, Dantzer F, DeMayo FJ, Lydon JP, Kraus WL. Deficiency of PARP-1 and PARP-2 in the mouse uterus results in decidualization failure and pregnancy loss. Proc Natl Acad Sci U.S.A. (2021) 118:e2109252118. doi: 10.1073/pnas.2109252118
65. Soni UK, Chadchan SB, Joshi A, Kumar V, Maurya VK, Verma RK, et al. Poly(ADP-ribose) polymerase-2 is essential for endometrial receptivity and blastocyst implantation, and regulated by caspase-8. Mol Cell Endocrinol (2020) 518:110946. doi: 10.1016/j.mce.2020.110946
66. Pirzada RH, Orun O, Erzik C, Cagsin H, Serakinci N. Role of TRF2 and TPP1 regulation in idiopathic recurrent pregnancy loss. Int J Biol Macromol (2019) 127:306–10. doi: 10.1016/j.ijbiomac.2019.01.011
67. Iachettini S, Ciccarone F, Maresca C, D’ Angelo C, Petti E, Di Vito S, et al. The telomeric protein TERF2/TRF2 impairs HMGB1-driven autophagy. Autophagy (2023) 19:1479–90. doi: 10.1080/15548627.2022.2138687
Keywords: Mendelian randomization, telomere length, spontaneous abortion, preterm birth, adverse pregnancy outcomes, causality
Citation: Han X, Wu T and Liu Cy (2023) Univariable and multivariable Mendelian randomization investigating the effects of telomere length on the risk of adverse pregnancy outcomes. Front. Endocrinol. 14:1225600. doi: 10.3389/fendo.2023.1225600
Received: 19 May 2023; Accepted: 17 July 2023;
Published: 03 August 2023.
Edited by:
Eytan R. Barnea, BioIncept, LLC, United StatesCopyright © 2023 Han, Wu and Liu. 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: Chun yan Liu, 624031871@qq.com