Skip to main content

Long-term outcomes of offspring from multiple gestations: a two-sample Mendelian randomization study on multi-system diseases using UK Biobank and FinnGen databases

Abstract

Background

Assisted reproductive technologies (ART) have increased the incidence of multiple births, which can have a negative impact on maternal and offspring health. The study aimed to investigate the association between genetically predicted multiple birth and the risk of 42 common diseases of the nervous, psychiatric, cardiovascular, respiratory, digestive, and endocrine systems.

Methods

The study utilized two-sample Mendelian randomization (MR) analysis to explore the potential causal relationship between genetically predicted multiple birth and the genetically predicted risk of diseases. The study used the FinnGen and UK Biobank datasets for analysis.

Results

The study found no significant causal relationship between multiple birth and psychiatric disorders. However, the lower limits of the 95% confidence intervals for bipolar affective disorder and anxiety disorders were not robust, indicating a need for further investigation. The study found that multiple birth may be a strong risk factor for infantile cerebral palsy, and caution is necessary in both natural and ART multiple births. The study revealed a potential causal relationship between multiple birth and coronary heart disease, ischemic heart disease, and deep vein thrombosis, which may be related to abnormal intrauterine environments in multiple pregnancies. Surprisingly, multiple birth appears to have a protective effect against some respiratory diseases, such as chronic obstructive pulmonary disease and asthma.

Conclusions

The study highlights the need for caution regarding the risk of infantile cerebral palsy, cardiovascular diseases, and psychiatric disorders in multiple birth. Our study can lead to the development of preventive strategies and improved clinical management for affected infants.

Background

The world is currently facing a critical issue of population decline, with several countries experiencing a decline in birth rates [1]. This trend has raised concerns about potential negative impacts on society, including labor shortages, decreased economic growth, and strain on healthcare and social systems [2]. Various causes contribute to negative population growth, with social pressures being a significant contributor to the incidence of infertility [3], identified as a leading cause of fertility decline in some developed countries [4]. Assisted reproductive technologies (ART), including in vitro fertilization (IVF), have revolutionized the field of reproductive medicine and have allowed many couples to conceive who would otherwise have been unable to do so [5]. While the development and maturity of ART is good news for infertile patients who are eager for children and has a positive effect on negative population growth, it has also led to a higher incidence of multiple birth [6]. Although multiple birth can further ameliorate negative population growth, they may increase the health risks of mothers and offspring [7].

Mothers with multiple pregnancies may be more vulnerable to a range of complications during pregnancy [8], including but not limited to gestational hypertension [9], gestational diabetes mellitus [10], antepartum hemorrhage [11], preterm birth [12], and anemia [13]. Pregnancy complications can cause adverse outcomes for mother and child, and even have long-term adverse effects on offspring [14, 15]. Several population-based studies have consistently shown that offspring of multiple pregnancies have worse long-term neurodevelopmental outcomes than singleton pregnancies [16, 17]. Therefore, as the use of ART continues to rise, it is becoming increasingly important to understand the potential long-term health consequences of different types of pregnancies. In particular, investigating the health outcomes of offspring of multiple pregnancies is crucial to ensure their long-term well-being and inform strategies to mitigate the negative impact of population decline.

However, it is not clear whether the susceptibility of twin offspring to certain diseases is due to innate genetic factors or acquired environmental factors, such as epigenetic changes. Furthermore, it is unclear whether twin pregnancy itself or complications caused by twin pregnancy cause the adverse long-term outcomes in offspring. This field is largely unexplored due to the length of follow-up required for such studies and the difficulty of intervention.

Mendelian randomization (MR), as an emerging epidemiological method in recent years, utilizes genetic variants as instrumental variables (IVs) to assess the causal effects of exposure factors on outcomes [18]. Here, "exposure" represents a risk factor, an intermediate phenotype, which may or may not serve as a potential causal biomarker which affects the "outcome" (for example, a certain disease) [19]. The possibility of reverse causation is minimized because genotype is established at conception (Mendel’s law of segregation). Thus, with the robust and reliable IVs provided by genome-wide association studies (GWAS), the MR analysis can be considered to be in accordance with the normal causal order [20].

This study aims to estimate the association between genetically predicted part of a multiple birth and the genetically predicted risk of several diseases. Using the FinnGen and UK Biobank databases, we assessed the susceptibility of offspring of multiple birth to 42 common diseases of the nervous, psychiatric, cardiovascular, respiratory, digestive, and endocrine systems.

Methods

Data sources

The study design is presented in Fig. 1 and adhered to the STROBE-MR guidelines [21]. We obtained genetic variants for multiple birth from the UK Biobank database [22] and limited out analysis to the largest European population, consisting of 414,217 individuals (including 9550 cases and 404,667 controls), in order to maintain consistency in ancestry between the exposure and outcome groups. The study collected responses to the question “Are you a twin, triplet or other multiple birth?” from all participants, except those who reported being adopted as a child. We obtained GWAS summary statistics for long-term outcomes from UK Biobank and FinnGen consortium R8 release data [23], and explored specific outcomes, including: attention deficit hyperactivity disorder, depression, autism, bipolar affective disorder, Alzheimer’s disease, cerebral palsy (or infantile cerebral palsy), epilepsy, cognitive impairment, schizophrenia, mood disorders, anxiety disorders, suicide or self-inflicted injury, migraine, mental retardation, stroke, subarachnoid hemorrhage, transient ischemic attack, pulmonary embolism, deep vein thrombosis, arterial hypertension, atrial fibrillation and flutter, ischemic heart disease, coronary heart disease (coronary atherosclerosis), cardiomyopathy, myocardial infarction, chronic obstructive pulmonary disease, asthma, bronchitis, tuberculosis, gastric ulcer, duodenal ulcer, Crohn's disease, fibrosis and cirrhosis of liver, gastritis (acute or chronic), thyrotoxicosis, hypothyroidism, thyroiditis, type 1 diabetes, type 2 diabetes, gout and obesity. Sex, age, principal components, and genotyping batch were corrected during the analysis.

Fig. 1
figure 1

Study design overview

Instrumental variable (IV)

In this study, we conducted a two-sample MR analysis [24] using single nucleotide polymorphisms (SNPs) as instrumental variables (IVs) to investigate the relationship between multiple birth and long-term outcomes [25]. As depicted in Fig. 1, our analysis satisfies three principal hypotheses of classical MR analysis: (1) IVs were strongly associated with the exposure; (2) IVs were independent of confounders; and (3) IVs influenced the risk of outcomes solely through the exposure.

In this study, we employed stringent criteria to select IVs for our two-sample MR analysis. All SNPs that were significantly associated with multiple birth (p < 5 × 10–6) were considered as IVs, and we applied LD pruning to ensure that the selected SNPs were independent within a 10 Mb window with an r2 < 0.001. We also searched for secondary phenotypes in PhenoScanner [26] and GWAS Catalog [27] to rule out any potential confounding effects, and SNPs corresponding to the outcomes of interest were excluded. We evaluated the strength of the IVs using variance (R2) and F-statistics to assess the extent of weak instrument bias [28]. R2 was calculated using the formula of 2 × MAF × (1-MAF) × β2, and F was calculated using the formula of R2 × (N − K − 1)/(1 − R2). If F > 10, we considered the correlation between IVs and exposure to be strong enough, and the MR analysis results would be less affected by weak instrument bias.

In total, 19 SNPs were selected as IVs to represent multiple birth (see Additional file 1: Table S1). Among these SNPs, rs605765 was associated with secondary phenotypes such as waist circumference, weight and number of operations, some of which are correlated with BMI. As BMI is a confounding factor in many diseases, we performed a one-leave-out test to assess the stability of our MR analysis results after removing this SNP. In further analyses, we deleted palindromic SNPs with a moderate allele frequency.

Statistical analysis

In this study, various MR methods, such as inverse variance weighted (IVW) [29], MR-Egger regression [30], weighted median [31], weighted mode [32], simple mode [33] and MR-PRESSO [34], were employed to investigate the potential causal relationship between multiple birth and long-term disease. The IVW was selected as the primary MR analysis, and the key features of each model were presented in Additional file 2: Table S2.

To ensure the robustness of our findings, we employed various sensitivity analyses. To assess the heterogeneity of the IVs, we conducted Cochran’s Q test, where a p value < 0.05 indicates heterogeneity [35]. We used MR-PRESSO with NbDistribution = 10,000 to identify any outliers and excluded them from the analysis. The MR-Egger intercept method was used to test for the presence of horizontal pleiotropic of IVs, and a p value < 0.05 suggests that the IVW estimate may be biased [30]. Moreover, we performed a leave-one-out sensitivity test to evaluate the impact of individual SNPs on the causal effect. Additionally, we generated funnel and forest plots to detect the presence of pleiotropy.

All statistical analyses were conducted using TwoSampleMR [36] packages in the R software (4.2.2), and all p values were two-tailed. A threshold of p < 0.05 was deemed statistically significant.

Results

Multiple birth and mental illness

In this part, we explored the association between multiple birth and nine common psychiatric disorders, including attention deficit hyperactivity disorder (ADHD), depression, autism, bipolar disorder, cognitive impairment, schizophrenia, mood disorders, anxiety disorders, and suicide or self-inflicted injury, using MR analysis. Figure 2 shows the results of the forest plot using the IVW method. The results suggest that there is no significant causal relationship between multiple birth and any of these psychiatric disorders, as observed in both FinnGen and UK Biobank datasets.

Fig. 2
figure 2

Association of genetically predicted multiple birth with risk of mental illness. Estimates were obtained from the inverse-variance weighted methods

Additional file 3: Table S3 shows the results of other MR methods, including MR Egger, weighted median, weighted mode, simple mode, and MR-PRESSO, all of which also did not show a significant causal relationship between multiple birth and these outcomes. Sensitivity analyses showed no heterogeneity except for the analysis of depression in the UK Biobank database. After removing outliers, still no significant causal relationship was observed between multiple birth and depression. The MR Egger intercept and MR-PRESSO outcomes indicate the absence of directional pleiotropy.

Additional file 9: Material S1 provides funnel plots, scatter plots, and one-leave-out test results for all outcomes to visualize heterogeneity more clearly. However, the results for bipolar disorder in the FinnGen database (OR 1.112, 95% CI 0.996–1.243) and anxiety disorders in the UK Biobank database (OR 1.092, 95% CI 0.996–1.198) were not robust, and the exclusion of certain SNPs led to the emergence of significant results that were previously non-significant, indicating a potential causal relationship between multiple birth and these two disorders.

Multiple birth and nervous system disease

In this section, we conducted an MR analysis to explore the association between multiple birth and eight neurological disorders, including stroke, subarachnoid hemorrhage, transient ischemic attack, epilepsy, migraine, Alzheimer's disease, adult and infantile cerebral palsy, and mental retardation. The IVW results, as shown in Fig. 3, suggest that there is no significant causal relationship between multiple birth and any of the neurological disorders investigated, based on outcome data from both FinnGen and UK Biobank.

Fig. 3
figure 3

Association of genetically predicted multiple birth with risk of nervous system disease. Estimates were obtained from the inverse-variance weighted methods

However, the MR Egger analysis results, presented in Additional file 4: Table S4, suggest that multiple birth is strongly associated with infantile cerebral palsy (OR 5.191, 95% CI 1.169–23.051), even though this association was not detected by the IVW method. It is worth noting that this result should be interpreted with caution due to the small number of cases in the UK Biobank data for infantile cerebral palsy (ncase = 166).

Despite the sensitivity analysis revealing heterogeneity in the results for epilepsy and cerebral palsy in the FinnGen database, MR-PRESSO did not detect any outliers. In addition, the association between multiple birth and transient ischemic attack from the FinnGen database was not reliable due to significant pleiotropy.

Additional file 10: Material S2 provides funnel plots, scatterplots, and one-leave-out test results for all analyses. It is noteworthy that the result for the association between multiple birth and epilepsy changed when rs72892862 was excluded. Multiple birth appeared to be a protective factor for epilepsy in this case (OR 0.859, 95% CI 0.743–0.993).

Overall, these findings suggest that multiple birth is not significantly associated with any of the neurological disorders investigated, except for a potential strong association with infantile cerebral palsy that requires further investigation.

Multiple birth and cardiovascular system disease

As shown in Fig. 4, MR analysis using IVW did not observe any significant causal relationship between multiple birth and arterial hypertension, atrial fibrillation and flutter, ischemic heart disease, coronary heart disease, cardiomyopathy, myocardial infarction, pulmonary embolism, or deep vein thrombosis.

Fig. 4
figure 4

Association of genetically predicted multiple birth with risk of cardiovascular system disease. Estimates were obtained from the inverse-variance weighted methods

Additional file 5: Table S5 provides the results of other analytical methods. The simple mode method suggested that multiple birth may be a risk factor for ischemic heart disease (OR 1.092, 95% CI 1.005–1.187, FinnGen). MR-PRESSO detected outliers in the results of multiple birth and coronary heart disease (FinnGen). When this outlier was removed, multiple birth appeared to be more strongly associated with coronary heart disease (OR 1.065, 95% CI 0.998–1.132) (Additional file 11: Material S3).

Although MR analysis found significant heterogeneity between multiple birth and arterial hypertension (FinnGen and UKBB), ischemic heart disease (FinnGen), and myocardial infarction (FinnGen), MR-PRESSO did not detect any single outlier. There seems to be a potential causal relationship between multiple birth and deep vein thrombosis, whether the data for deep vein thrombosis is from FinnGen (OR 1.101, 95% CI 0.992–1.223) or UKBB (OR 1.105, 95% CI 0.978–1.249). The results of the one-leave-out test also suggest that both of these results are susceptible to the influence of a single SNP, and the estimated causal relationship becomes more significant after removing a certain SNP (Additional file 11: Material S3).

Multiple birth and respiratory system disease

Figure 5 displays the correlation between multiple birth and four respiratory diseases: chronic obstructive pulmonary disease (COPD), asthma, bronchitis, and tuberculosis, as analyzed through MR analysis. The results suggest that multiple birth may be a protective factor against lifelong COPD (OR 0.933, 95% CI 0.873–0.998, FinnGen). Regardless of the MR analysis method used (as shown in Additional file 6: Table S6), no significant causal relationship was found between multiple birth and other respiratory diseases, and all results showed no significant heterogeneity or pleiotropy. As shown in Additional file 12: Material S4, when rs143496908 was removed, the MR analysis result showed a reduced probability of multiple birth being associated with asthma (OR 0.953, 95% CI 0.909–0.999, FinnGen).

Fig. 5
figure 5

Association of genetically predicted multiple birth with risk of respiratory system disease. Estimates were obtained from the inverse-variance weighted methods

Multiple birth and digestive system disease

Next, we investigated the association between multiple birth and several common digestive system disorders, including gastric ulcer, duodenal ulcer, Crohn's disease, fibrosis and cirrhosis of liver, and gastritis (chronic and acute). As shown in Fig. 6, the probability of multiple birth being diagnosed with duodenal ulcer was decreased (OR 0.853, 95% CI 0.738–0.986, FinnGen). However, this result should be interpreted with caution because this outcome was only observed in the FinnGen database, and the trend was actually reversed in the UK Biobank (OR 1.126, 95% CI 0.966–1.314). In addition, multiple birth seemed to be more susceptible to chronic gastritis (OR 1.101, 95% CI 1.008–1.204). The MR analysis results from other methods and sensitivity analyses are presented in Additional file 7: Table S7. Although MR Egger suggested a potential association between multiple birth and gastric ulcer (OR 0.693, 95% CI 0.530–0.906), the result was not reliable due to significant pleiotropy (p < 0.05). Furthermore, the analysis of fibrosis and cirrhosis of liver (UK Biobank) showed significant heterogeneity, and even after removing outliers using MR-PRESSO, no significant association was found between multiple birth and these diseases.

Fig. 6
figure 6

Association of genetically predicted multiple birth with risk of digestive system disease. Estimates were obtained from the inverse-variance weighted methods

Additional file 13: Material S5 provides the funnel plot, scatter plot, and one-leave-out test results. Previously discovered associations between multiple birth and duodenal ulcer (FinnGen) and chronic gastritis (UKBB) were found to change upon removal of a single SNP, indicating that the results are not robust.

Multiple birth and endocrine system disease

As shown in Fig. 7 and Additional file 8: Table S8, the IVW analysis indicated that the likelihood of developing seven common endocrine disorders, including thyrotoxicosis, hypothyroidism, thyroiditis, type 1 diabetes, type 2 diabetes, obesity, and gout, did not significantly differ between multiple birth and singleton birth. Only the MR Egger analysis suggested that multiple birth may be protective against thyroiditis (OR 0.631, 95% CI 0.409–0.971, FinnGen). None of the analyses showed significant pleiotropy, and some results with heterogeneity did not change after removing outliers. Moreover, all results were robust to the exclusion of individual SNPs (Additional file 14: Material S6).

Fig. 7
figure 7

Association of genetically predicted multiple birth with risk of endocrine system disease. Estimates were obtained from the inverse-variance weighted methods

Discussion

Assisted reproductive technologies have raised multiple birth rates, potentially heightening health risks for both mothers and offspring. This study aimed to investigate the association between genetically predicted part of a multiple birth and the genetically predicted risk of 42 common diseases of the nervous, psychiatric, cardiovascular, respiratory, digestive, and endocrine systems.

Regarding psychiatric disorders, we found no significant causal relationship between multiple birth and any investigated psychiatric disorders, whether the data came from FinnGen or UK Biobank. This result may differ from some clinical studies [37, 38]. Possible explanations could be that childhood psychiatric disorders might be influenced by environmental or epigenetic factors [39]. The upbringing and surroundings of multiple offspring differ from those of singletons, and certain family and social factors, such as strict parenting, could contribute to these differences in environmental and epigenetic factors [40]. Additionally, while our study suggests that multiple pregnancies themselves may not increase the risk of psychiatric disorders in offspring, pregnant women with multiple pregnancies may be more susceptible to various pregnancy complications [41, 42], some of which, such as preterm birth [43], have been associated with psychiatric disorders in offspring [44]. Besides, one-leave-out tests revealed that the results for bipolar affective disorder (OR 1.112, 95% CI 0.996–1.243, FinnGen) and anxiety disorders (OR 1.092, 95% CI 0.996–1.198, UKBB) lacked robustness, with the lower limits of their 95% confidence intervals approaching non-significant values. Removing specific individual SNPs could render these correlations significant. In fact, in 2019, Amrhein and other statisticians cautioned researchers about avoiding the "significance trap": science researcher should never conclude that there is "no difference" or "no association" just because a p value is larger than a threshold such as 0.05 or, equivalently, because a confidence interval includes zero [45]. Therefore, we still need to be cautious about the risk of multiple birth offspring developing bipolar affective disorder and anxiety disorders in the future. It is important to acknowledge that the upbringing environment within families with multiple birth can be more intricate than anticipated. A child in a twin family may experience enhanced happiness due to having a companion during their growth, but they might also encounter dissatisfaction from the need to share their belongings. Furthermore, parental attitudes towards their children could also play a role, particularly evident in dizygotic twins (especially when of different genders). These factors could potentially increase or decrease a child's susceptibility to certain mental disorders. However, no current observational studies compare the lifetime prevalence of these two disorders between multiple birth and singleton offspring, which could be a starting point for future research.

IVW analysis found no significant causal relationship between multiple birth and several common neurological diseases. However, according to MR Egger, multiple birth may be a strong risk factor for infantile cerebral palsy (OR 5.191, 95% CI 1.169–23.051). Unlike the IVW method, which assumes that all genetic variants affecting the outcome variable come from the same causal pathway, the MR Egger method is more flexible and can draw causal conclusions from non-causal genetic correlations as long as a weak bias assumption is met [46]. Therefore, although MR Egger's analysis not being the primary approach, the substantial OR value demands careful consideration. Previous clinical studies have also shown an increased risk of infantile cerebral palsy in multiple birth [47,48,49,50]. Spangmose et al. found no significant rise in ORs of cerebral palsy in ART twins compared to naturally conceived twins, providing strong evidence against the continued use of multiple embryo transfer in most ART settings [51]. The primary risk factors for infantile cerebral palsy are prematurity and low birth weight [52, 53], occurrences more prevalent in multiple pregnancies than in singleton ones [12, 54], which may be the reason for the increased incidence of cerebral palsy in multiple births. Our study results indicate that caution should be exercised regarding the risk of infantile cerebral palsy in both natural and ART multiple births. Additionally, since the phenotype of infantile cerebral palsy is only provided in the UK Biobank and the number of cases is small (ncase = 166), the non-significant IVW results may be due to the small sample size. Besides, it is worth noting that the lifetime risk of epilepsy in multiple birth appears to be lower. Further studies are needed to explore the potential mechanisms underlying the association between multiple birth and epilepsy, which may ultimately lead to the development of preventive strategies and improved clinical management for affected infants.

After excluding outliers, IVW analysis revealed a potential causal relationship between multiple birth and coronary heart disease (OR 1.065, 95% CI 0.998–1.132, FinnGen). Simple mode also suggested an increased risk of ischemic heart disease in multiple births (OR 1.092, 95% CI 1.005–1.187, FinnGen). In addition, GWAS data from FinnGen (OR 1.101, 95% CI 0.992–1.223) and UK Biobank (OR 1.105, 95% CI 0.978–1.249) suggested a potential causal relationship between multiple birth and deep vein thrombosis, and both results could become significant after the exclusion of a single SNP. Previous studies have suggested that abnormalities in the early intrauterine environment may lead to adverse long-term effects in offspring, recognized as metabolic memory (MM). Multiple clinical and experimental studies have shown that MM may cause persistent endothelial dysfunction in offspring, amplifying the vulnerability to long-term cardiovascular complications [14, 15]. Our study suggests an augmented susceptibility to certain cardiovascular diseases in multiple birth, potentially tied to aberrations in the intrauterine environment in multiple pregnancies. Previous studies have shown that pregnant women with multiple pregnancies face a heightened risk of metabolic disorders such as gestational diabetes [55], which may lead to a high glucose environment in the uterus and result in metabolic memory in offspring during early life. Therefore, pregnant women with multiple pregnancies should diligently monitor their blood glucose levels throughout gestation, and individuals with multiple births should be vigilant about the occurrence of cardiovascular diseases in adulthood. This conclusion needs to be confirmed by large-scale clinical studies and further investigation into changes in the intrauterine environment in multiple births and their impact on the development of lifelong cardiovascular diseases in offspring.

Surprisingly, multiple birth appears to shield against some respiratory diseases, such as COPD (OR 0.933, 95% CI 0.873–0.998, FinnGen) and asthma (OR 0.953, 95% CI 0.909–0.999, FinnGen), excluding single SNP. One possible explanation is that twins could acquire early-life immune safeguards via shared placenta and fetal membranes, potentially curbing future asthma and COPD risk. Furthermore, in the Finnish population, multiple births also exhibit reduced susceptibility thyroiditis (OR 0.631, 95% CI 0.409–0.971, FinnGen). Non-identical twins sharing the uterine environment are exposed to non-identical antigens from each other, which may induce a higher level of tolerance. This partially elucidates multiple births’ conceivable role as a protective factor against asthma and thyroiditis, necessitating validation through subsequent clinical and basic research.

MR analysis regarding the link between multiple birth and duodenal ulcer, based on data from FinnGen and UK Biobank, exhibited divergent trends. FinnGen’s data suggests that multiple birth may be a protective factor for duodenal ulcer (OR 0.853, 95% CI 0.738–0.986), while UK Biobank’s data suggests a potential causal relationship between multiple birth and increased risk of duodenal ulcer (OR 1.126, 95% CI 0.966–1.314). Notably, most results obtained from the two databases do not match completely, possibly stemming from differences in population and age. UK Biobank focused on individuals aged 40–69 within the UK, while FinnGen's study population may include any Finnish residents. In addition, FinnGen's data also indicates an increased risk of chronic gastritis in multiple birth (OR 1.101, 95% CI 1.008–1.204). These results need to be interpreted with caution, and studies on the impact of twins on disease incidence rates among different populations are necessary.

Although multiple births can ameliorate negative population growth, they may increase the health risks of mothers and offspring. Pregnancy complications can cause adverse outcomes for mother and child, and even have long-term adverse effects on offspring. Thus, as the use of ART gain prominence, it is crucial to investigate the health prospects of offspring from multiple pregnancies, safeguarding their long-term well-being and guiding strategies to counterbalance the adverse consequences of population decline. Moreover, ART-conceived and naturally conceived multiple births may exhibit some differences in the genetic component. ART, encompassing in vitro fertilization and cultivation, involves additional laboratory procedures and interventions that could potentially influence genetic component. Processes like embryo screening within ART might lead to differences in the frequency of certain gene types or variants in ART-conceived individuals. Similarly, the process of embryo cultivation and implantation could also impact their genetic composition. However, it should be noted that in most genetic characteristics, ART-conceived and naturally conceived children should be similar since the genetic information primarily originates from the genetic contributions of their parents. The genetic information of ART-conceived infants still derives from their biological parents' genetic material, and the in vitro fertilization process does not involve modifications to the genome.

The study’s results have implications for clinicians and policymakers. They highlight the potential for varied disease risk among offspring of multiple pregnancies. Vigilant monitoring of the long-term health trajectories of these offspring, particularly those born from ART, becomes imperative. Further research is needed to investigate the potential associations between multiple birth and other diseases and conditions.

The study's strengths include its large sample size, utilization of MR analysis, and the availability of data from two large databases. By using "a more realistic p value", Smith et al. suggest that compared to conventional observational epidemiology, Mendelian randomization studies can capitalize on the limited associations between a specific genetic variant and other genetic variants, as well as between genetic and nongenetic variables. This enables them to provide an unbiased estimation of the relationship between the factors directly influenced by the genetic variant and disease outcomes, making them less susceptible to confounding [56].

However, this study has several limitations. Firstly, the absence of pertinent data hindered the ability to perform separate analyses for monozygotic and dizygotic multiple births. Specifically, in dizygotic multiple births, there may be significant differences in genotypes among individual fetuses. Similarly, comparing distinctions between ART and non-ART twins proved challenging, given potential genetic variations introduced by the ART process. Secondly, the study population included only individuals of European ancestry, and the results need to be validated in more diverse populations. Thirdly, while two-sample Mendelian randomization is a powerful discovery tool, one-sample MR would confirm these findings and enable a direct evaluation of potential confounding factors and covariates that may contribute to pleiotropic effects. Lastly, the results obtained through Mendelian randomization in this study lack direct support from relevant clinical randomized controlled trials, warranting cautious interpretation. Despite these limitations, the study provides valuable insights into the potential health risks associated with multiple births.

Conclusion

In summary, this research adds to the growing body of evidence on the health risks associated with multiple birth. Although some of the findings lack robustness, the study suggests that individuals who are born from multiple pregnancies should be aware of the increased lifetime risk for bipolar affective disorder, anxiety disorders, infantile cerebral palsy, coronary heart disease, ischemic heart disease, deep vein thrombosis, and chronic gastritis. This study emphasizes the importance of investigating the health outcomes of individuals who are born from multiple pregnancies and highlights the need for further research on the long-term health consequences of various types of pregnancies. The study's findings have important clinical implications for managing multiple pregnancies and the long-term health outcomes of offspring.

Availability of data and materials

The summary data of Pan-UK Biobank can be downloaded from the website https://pan.ukbb.broadinstitute.org/downloads/index.html. The summary data of FinnGen can be downloaded from the website https://www.finngen.fi/en/access_results. The other datasets generated and/or analysed during the current study are publicly available and included in this published article and its supplementary information files.

Abbreviations

ADHD:

Attention deficit hyperactivity disorder

ART:

Assisted reproductive technologies

COPD:

Chronic obstructive pulmonary disease

GWAS:

Genome-wide association studies

IVF:

In vitro fertilization

IVs:

Instrumental variables

MM:

Metabolic memory

MR:

Mendelian randomization

SNP:

Single nucleotide polymorphism

UKBB:

UK Biobank

References

  1. Lesthaeghe R. The second demographic transition, 1986–2020: sub-replacement fertility and rising cohabitation—a global update. Genus. 2020;76(1):10.

    Article  Google Scholar 

  2. Ogura S, Jakovljevic MM. Editorial: Global population aging—health care, social and economic consequences. Front Public Health. 2018;6:335.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Ferraretti AP, Goossens V, de Mouzon J, et al. Assisted reproductive technology in Europe, 2008: results generated from European registers by ESHRE. Hum Reprod. 2012;27(9):2571–84.

    Article  CAS  PubMed  Google Scholar 

  4. Boivin J, Bunting L, Collins JA, et al. International estimates of infertility prevalence and treatment-seeking: potential need and demand for infertility medical care. Hum Reprod. 2007;22(6):1506–12.

    Article  PubMed  Google Scholar 

  5. Scott RT Jr, de Ziegler D. Introduction: Key performance indicators in assisted reproductive technologies. Fertil Steril. 2020;114(1):4–5.

    Article  PubMed  Google Scholar 

  6. Sunderam S, Kissin DM, Zhang Y, et al. Assisted reproductive technology surveillance—United States, 2018. MMWR Surveill Summ. 2022;71(4):1–19.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Wang Y, Shi H, Chen L, et al. Absolute risk of adverse obstetric outcomes among twin pregnancies after in vitro fertilization by maternal age. JAMA Netw Open. 2021;4(9):e2123634.

    Article  PubMed  PubMed Central  Google Scholar 

  8. ACOG Practice Bulletin No. 144. Multifetal gestations: twin, triplet, and higher-order multifetal pregnancies. Obstet Gynecol. 2014;123(5):1118–32.

    Article  Google Scholar 

  9. Narang K, Szymanski LM. Multiple gestations and hypertensive disorders of pregnancy: what do we know? Curr Hypertens Rep. 2020;23(1):1.

    Article  PubMed  Google Scholar 

  10. Hiersch L, Berger H, Okby R, et al. Gestational diabetes mellitus is associated with adverse outcomes in twin pregnancies. Am J Obstet Gynecol. 2019;220(1):1021e01-1102.

    Article  Google Scholar 

  11. Loussert L, Schmitz T, Korb D, et al. Risk of severe postpartum hemorrhage in twin pregnancies according to the sum of birth weights. Obstet Gynecol. 2022;140(6):958–64.

    Article  PubMed  Google Scholar 

  12. Roman A, Ramirez A, Fox NS. Prevention of preterm birth in twin pregnancies. Am J Obstet Gynecol MFM. 2022;4(2S):100551.

    Article  PubMed  Google Scholar 

  13. Giorgione V, D’Antonio F, Manji A, et al. Perinatal outcome of pregnancy complicated by twin anemia-polycythemia sequence: systematic review and meta-analysis. Ultrasound Obstet Gynecol. 2021;58(6):813–23.

    Article  CAS  PubMed  Google Scholar 

  14. Zhang E, Guo Q, Gao H, et al. Metformin and resveratrol inhibited high glucose-induced metabolic memory of endothelial senescence through SIRT1/p300/p53/p21 pathway. PLoS ONE. 2015;10(12):e0143814.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Singh KK, Mantella LE, Pan Y, et al. A global profile of glucose-sensitive endothelial-expressed long non-coding RNAs. Can J Physiol Pharmacol. 2016;94(9):1007–14.

    Article  CAS  PubMed  Google Scholar 

  16. Lorenz JM. Neurodevelopmental outcomes of twins. Semin Perinatol. 2012;36(3):201–12.

    Article  PubMed  Google Scholar 

  17. Wadhawan R, Oh W, Perritt RL, et al. Twin gestation and neurodevelopmental outcome in extremely low birth weight infants. Pediatrics. 2009;123(2):e220-227.

    Article  PubMed  Google Scholar 

  18. Smith GD, Ebrahim S. “Mendelian randomization”: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.

    Article  PubMed  Google Scholar 

  19. Lee K, Lim CY. Mendelian randomization analysis in observational epidemiology. J Lipid Atheroscler. 2019;8(2):67–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  21. Skrivankova VW, Richmond RC, Woolf BAR, et al. Strengthening the Reporting of Observational Studies in Epidemiology Using Mendelian Randomization: The STROBE-MR Statement. JAMA. 2021;326(16):1614–21.

    Article  PubMed  Google Scholar 

  22. UK Biobank GWAS Results.

  23. Kurki MI, Karjalainen J, Palta P, et al. FinnGen: Unique genetic insights from combining isolated population and national health register data. medRxiv, 2022:2022.2003.2003.22271360.

  24. Hemani G, Tilling K, Davey SG. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13(11):e1007081.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Lawlor DA, Harbord RM, Sterne JA, et al. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27(8):1133–63.

    Article  PubMed  Google Scholar 

  26. Kamat MA, Blackshaw JA, Young R, et al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. 2019;35(22):4851–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Buniello A, MacArthur JAL, Cerezo M, et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019;47(D1):D1005–12.

    Article  CAS  PubMed  Google Scholar 

  28. Burgess S, Thompson SG, Collaboration CCG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64.

    Article  PubMed  Google Scholar 

  29. Burgess S, Bowden J, Fall T, et al. Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants. Epidemiology. 2017;28(1):30–42.

    Article  PubMed  Google Scholar 

  30. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Bowden J, Davey Smith G, Haycock PC, et al. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 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(6):1985–98.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Milne RL, Kuchenbaecker KB, Michailidou K, et al. Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer. Nat Genet. 2017;49(12):1767–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Verbanck M, Chen CY, Neale B, et al. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Greco MF, Minelli C, Sheehan NA, et al. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med. 2015;34(21):2926–40.

    Article  Google Scholar 

  36. Hemani G, Zheng J, Elsworth B, et al. The MR-base platform supports systematic causal inference across the human phenome. Elife. 2018. https://doi.org/10.7554/eLife.34408.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Taylor CL, Rice ML, Christensen D, et al. Prenatal and perinatal risks for late language emergence in a population-level sample of twins at age 2. BMC Pediatr. 2018;18(1):41.

    Article  PubMed  PubMed Central  Google Scholar 

  38. You J, Shamsi BH, Hao MC, et al. A study on the neurodevelopment outcomes of late preterm infants. BMC Neurol. 2019;19(1):108.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Panariello F, Fanelli G, Fabbri C, et al. Epigenetic basis of psychiatric disorders: a narrative review. CNS Neurol Disord Drug Targets. 2022;21(4):302–15.

    Article  CAS  PubMed  Google Scholar 

  40. Chubar V, Van Leeuwen K, Bijttebier P, et al. Gene-environment interaction: New insights into perceived parenting and social anxiety among adolescents. Eur Psychiatry. 2020;63(1):e64.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Cheong-See F, Schuit E, Arroyo-Manzano D, et al. Prospective risk of stillbirth and neonatal complications in twin pregnancies: systematic review and meta-analysis. BMJ. 2016;354:i4353.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Francisco C, Wright D, Benko Z, et al. Competing-risks model in screening for pre-eclampsia in twin pregnancy according to maternal factors and biomarkers at 11–13 weeks’ gestation. Ultrasound Obstet Gynecol. 2017;50(5):589–95.

    Article  CAS  PubMed  Google Scholar 

  43. Herlihy N, Odom E, Cohen N, et al. Long-term outcomes of small for gestational age twins born at 34 weeks or later. Am J Perinatol. 2018;35(3):254–61.

    Article  PubMed  Google Scholar 

  44. Crump C, Sundquist J, Sundquist K. Preterm or early term birth and risk of autism. Pediatrics. 2021. https://doi.org/10.1542/peds.2020-032300.

    Article  PubMed  Google Scholar 

  45. Amrhein V, Greenland S, McShane B. Scientists rise up against statistical significance. Nature. 2019;567(7748):305–7.

    Article  CAS  PubMed  Google Scholar 

  46. Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Williams K, Hennessy E, Alberman E. Cerebral palsy: effects of twinning, birthweight, and gestational age. Arch Dis Child Fetal Neonatal Ed. 1996;75(3):F178-182.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Scher AI, Petterson B, Blair E, et al. The risk of mortality or cerebral palsy in twins: a collaborative population-based study. Pediatr Res. 2002;52(5):671–81.

    Article  PubMed  Google Scholar 

  49. Sellier E, Goldsmith S, McIntyre S, et al. Cerebral palsy in twins and higher multiple births: a Europe-Australia population-based study. Dev Med Child Neurol. 2021;63(6):712–20.

    Article  PubMed  Google Scholar 

  50. Tronnes H, Wilcox AJ, Lie RT, et al. Risk of cerebral palsy in relation to pregnancy disorders and preterm birth: a national cohort study. Dev Med Child Neurol. 2014;56(8):779–85.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Spangmose AL, Christensen LH, Henningsen AA, et al. Cerebral palsy in ART children has declined substantially over time: a Nordic study from the CoNARTaS group. Hum Reprod. 2021;36(8):2358–70.

    Article  PubMed  Google Scholar 

  52. Graham HK, Rosenbaum P, Paneth N, et al. Cerebral palsy. Nat Rev Dis Primers. 2016;2:15082.

    Article  PubMed  Google Scholar 

  53. Paul S, Nahar A, Bhagawati M, et al. A review on recent advances of cerebral palsy. Oxid Med Cell Longev. 2022;2022:2622310.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Townsend R, Khalil A. Fetal growth restriction in twins. Best Pract Res Clin Obstet Gynaecol. 2018;49:79–88.

    Article  CAS  PubMed  Google Scholar 

  55. Monteiro SS, Fonseca L, Santos TS, et al. Gestational diabetes in twin pregnancy: a predictor of adverse fetomaternal outcomes? Acta Diabetol. 2022;59(6):811–8.

    Article  CAS  PubMed  Google Scholar 

  56. Smith GD, Lawlor DA, Harbord R, et al. Clustered environments and randomized genes: a fundamental distinction between conventional and genetic epidemiology. PLoS Med. 2007;4(12):e352.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank the participants of all GWAS cohorts included in the present work and the investigators of the UK Biobank, PhenoScanner, GWAS Catalog and FinnGen for sharing the GWAS summary statistics.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualisation, YJ; methodology, YJ, YD, RS, and XZ; data Analysis (Psychiatric disorders), XZ and LW; data Analysis (Neurological Diseases), JZ and SZ; data Analysis (Cardiovascular Diseases), HZ and CF; data Analysis (Respiratory Diseases), YC and PG; data Analysis (Digestive Diseases), LZ and SW; data Analysis (Endocrine Diseases), JY and MH; writing—original draft, YJ; writing—review and editing, WD and LF. All the authors participated in planning, execution, and analysis and have read and approved the final submitted version.

Corresponding authors

Correspondence to Wencheng Ding or Ling Feng.

Ethics declarations

Ethics approval and consent to participate

Our analysis used publicly available genome-wide association study (GWAS) summary statistics. No new data were collected, and no new ethical approval was required.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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: Table S1

. SNPs that were used as instrumental variables met the conditions: p value < 5 × 10–6 and F > 10. * The secondary phenotypes of rs605765 obtained by PhenoScanner include waist circumference, weight and number of operations, some of which are associated with BMI, and BMI was identified as a confounder for lots of diseases. In the subsequent one-leave-out test, we focused on whether the MR Results were stable after removing this SNP.

Additional file 2: Table S2

. Characteristics of each method used in this study.

Additional file 3: Table S3

. Two-sample Mendelian randomization estimations showing the effects, heterogeneity and horizontal pleiotropy of multiple birth on the risk of mental illness.

Additional file 4: Table S4

. Two-sample Mendelian randomization estimations showing the effects, heterogeneity and horizontal pleiotropy of multiple birth on the risk of nervous system disease.

Additional file 5: Table S5

. Two-sample Mendelian randomization estimations showing the effects, heterogeneity and horizontal pleiotropy of multiple birth on the risk of cardiovascular system disease.

Additional file 6: Table S6

. Two-sample Mendelian randomization estimations showing the effects, heterogeneity and horizontal pleiotropy of multiple birth on the risk of respiratory system disease.

Additional file 7: Table S7

. Two-sample Mendelian randomization estimations showing the effects, heterogeneity and horizontal pleiotropy of multiple birth on the risk of digestive system disease.

Additional file 8: Table S8

. Two-sample Mendelian randomization estimations showing the effects, heterogeneity and horizontal pleiotropy of multiple birth on the risk of endocrine system disease.

Additional file 9: Material S1

. The scatter plot, funnel plot and leave-one-out plot for the MR analysis of multiple birth and mental illness.

Additional file 10: Material S2

. The scatter plot, funnel plot and leave-one-out plot for the MR analysis of multiple birth and nervous system disease.

Additional file 11: Material S3

. The scatter plot, funnel plot and leave-one-out plot for the MR analysis of multiple birth and cardiovascular system disease.

Additional file 12: Material S4

. The scatter plot, funnel plot and leave-one-out plot for the MR analysis of multiple birth and respiratory system disease.

Additional file 13: Material S5

. The scatter plot, funnel plot and leave-one-out plot for the MR analysis of multiple birth and digestive system disease.

Additional file 14: Material S6

. The scatter plot, funnel plot and leave-one-out plot for the MR analysis of multiple birth and endocrine system disease.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jiang, Y., Du, Y., Su, R. et al. Long-term outcomes of offspring from multiple gestations: a two-sample Mendelian randomization study on multi-system diseases using UK Biobank and FinnGen databases. J Transl Med 21, 608 (2023). https://doi.org/10.1186/s12967-023-04423-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-023-04423-w

Keywords