Skip to main content

Chorionicity-associated variation in metabolic phenotype of cord blood in twin

Abstract

Background

Monochorionic (MC) twins present a higher incidence of unfavorable clinical perinatal outcomes than dichorionic (DC) twins, often in association with placental vascular anastomosis. In this study, we profiled the umbilical cord plasma metabolomes of uncomplicated MC and DC twin pregnancies and related these to several offspring outcomes, previously associated with birthweight.

Methods

Umbilical vein blood samples were collected at birth from 25 pairs of uncomplicated MC twins and 24 pairs of uncomplicated DC twins. The samples were subjected to gas chromatography-mass spectrometry-based metabolomics. 152 metabolites were identified from the cord plasma samples of MC and DC twins. Partial least squares discriminant analysis and pathway analysis were performed to compare within DC/MC twin pairs and between DC and MC twins. A generalized estimating equation (GEE) model was utilized to explore the correlation between metabolic differences and birthweight discordance within and between twin pairs.

Results

Our study revealed clear differences between the metabolite profiles of umbilical cord plasma of MC and DC twins. Metabolite profiles in MC within twin pairs and DC within twin pairs were characterized by the differences in 2 − hydroxyglutaramic acid levels and nicotinamide levels, respectively. The metabolic pathways of GSH, tryptophan, and fatty acid metabolism, were significantly downregulated in MC twins compared to DC twins. In addition, the concentration of caffeine and decamethyl-cyclopentasiloxane (D5) was positively correlated with birthweight in MC and DC twins.

Conclusion

This study demonstrated that the altered metabolites in umbilical plasma made contributions to the different chorionicities between uncomplicated MC twins and DC twins. The chorionicity of twins seems to affect the metabolic cross-talk between co-twin pairs and be related to birthweight discordance of twins.

Background

In recent decades, advanced maternal age and the increasing use of assisted reproductive technologies have resulted in a significant increase in multiple births [1]. Twin pregnancies accounted for 0.78% of all births in China in 1989 [2] but more than doubled to 1.88% from 2007 to 2014 [3]. Twin pregnancies are associated with a higher incidence of maternal and fetal complications [4, 5] compared to singletons. Rates of stillbirth and neonatal mortality for twin gestations are 14.1 and 26.1 per 1000 total births, respectively [3]. Compared with singleton pregnancies, monochorionic (MC) twin pregnancies showed a thirteenfold increase, while dichorionic (DC) twins showed a fivefold increase, in rates of stillbirth [6,7,8]. MC twins have a significantly greater incidence of perinatal death (11.6% in MC twins versus 5.0% in DC twins), necrotising enterocolitis (OR 4.05, 95% CI 1.97–8.35), and neurological injury compared to DC twins [9]. Unequal sharing of the placental territory and vascular communications between twins made a remarkable contribution to differentiating fetal development [9, 10]. Although multiple studies have aimed to address the issue of clinical treatment and delivery outcomes of twins through retrospective analysis [11,12,13], the underlying metabolic differences in twin pregnancies, associated with adverse perinatal outcomes, have not been explored in detail. This is particularly true for MC relative to DC twin pairs.

Metabolomics enables the investigation of both normal physiology and the pathophysiology of many diseases by using advanced analytical chemistry techniques [14]. Identification of specific metabolites associated with twin pregnancies generally, or distinct to MC relative to DC pregnancies, using global untargeted approaches, has the potential to provide insights into clinical management. A growing number of studies have attempted to identify metabolic variations associated with the pathophysiology of MC-specific complications [15,16,17] and most have found disrupted amino acid and fatty acid metabolism in umbilical cord blood and/or placental tissue [15, 16]. However, these studies failed to investigate the metabolism of uncomplicated DC twin pregnancies and to differentiate between DC and MC pairs. Comparing uncomplicated MC twins with DC twins provide a favorable comparative study to investigate the effect of chorionicity between siblings on placental metabolite allocation.

Here, we hypothesized that the umbilical cord blood of MC and DC twins from uncomplicated pregnancies would show differing metabolite profiles, potentially in association with placenta-specific factors. Thus, the present study aimed to enrich the understanding of chorionicity and enable comparison with pathological twin gestations in future studies.

Methods

Study design

Women with twin pregnancies were recruited from the Peking University Third Hospital, and included in the study as part of the University Hospital Advanced Age Pregnant Cohort (clinicaltrials.gov Identifier: NCT03220750). The research protocols involving human participants were approved by The Human Research Ethics Committee of Peking University Third Hospital (IRB00006761-2016145). All participants have written informed consent prior to study recruitment. A total number of 258 twin pregnancies were enrolled in the cohort between September 2017 and December 2018. Briefly, exclusion criteria included maternal chronic diseases, malnutrition, smoking, drug abuse, fetal congenital and genetic anomalies, intrauterine fetal death (IUFD), twin-to-twin transfusion syndrome (TTTS), twin anemia polycythemia sequence (TAPS), selective intrauterine growth restriction (sIUGR), other adverse twin pregnancy outcomes, and those lost to follow-up at birth (Supplementary Fig. 1). Gestational age was determined by last menstrual period and confirmed with the ultrasound measurement of the crown-rump length (CRL) of the larger twin during 10–14 weeks’ gestation by first-trimester ultrasound examination for spontaneous conception according to the guideline of the International Society of Ultrasound in Obstetrics and Gynecology (ISUOG) [18]. After these exclusions, uncomplicated twin pregnancies with intertwin estimated fetal weight (EFW) discordances less than 25% remained in our study were 25 MC and 24 DC twin pairs. The fetal birth weight was measured by a clean electronic balance at birth. Thus, we take the following three groups for comparison: Comparison 1 compared the larger and smaller twin of DC twins (DC-L/DC-S). Comparison 2 compared the larger and smaller twin of MC twins (MC-L/MC-S). Comparison 3 compared the twins of DC and MC. All three comparisons are shown in Fig. 1.

Fig. 1
figure 1

Overview of study design. 49 twin pregnancies were included in this study, including 24 DC and 25 MC pregnancies. Comparison 1 identifies within-pair metabolite variation between larger and smaller twins in DC twin pairs. Comparison 2 identified similar differences in MC twin Pairs. Comparison 3 is a comparison between DC and MC pregnancies

Ultrasound assessment

Three different obstetric sonographers were assigned to perform all ultrasound examinations using a Voluson E10 (GE Healthcare, Zipf, Austria) ultrasound machine equipped with a C5-2 convex array probe and a Rm6c volume probe on participants whose gestation age was prior to 16 weeks. After recruitment at 16 weeks of gestation, ultrasound examinations were performed by one certified twin specialist obstetric sonographer. Moreover, fetal biometry/Doppler indices of all recruited participants were reported by the same sonographer at least once every two weeks based on the recommendations of the ISUOG [19].

Diagnosis of MC/DC pregnancy

The diagnosis of chorionicity and amnionicity was determined by a fixed obstetric sonographer. Twin pregnancy was identified based on ultrasonography conducted in the first trimester of gestation. Chorionicity was identified by (a) assessment of the membrane thickness at the insertion site of the amniotic membrane into the placenta, (b) determination of the T-sign or lambda sign between 11 and 14 gestational weeks, and (c) the number of placentas. Besides, the chorionicity of twin pregnancies was further evaluated after delivery via the pathological characteristics of the placenta.

Sample collection

Umbilical cord blood samples were collected from each of the umbilical veins into EDTA-coated blood collection tubes immediately after delivery and processed within twelve hours. Plasma was obtained by standard density gradient centrifugation, which was centrifuged twice at 3,000 rpm for 10 min at 4 °C. And then, plasma was transferred into cryopreservation tubes (Micronic, Holland) and stored at − 80 °C until analysis.

Metabolite extraction from cord plasma

Plasma supernatants were isolated after centrifugation at 12,000 rpm for 15 min at 4 °C. The supernatants were then dehydrated in a Speed Vac (Labconco, USA) at room temperature for 7 h. Dried supernatants were kept at − 80 °C prior to chemical derivatization.

Methyl Chloroform Derivatization and Gas Chromatography-Mass Spectrometry (GS-MS) analysis

The extracted samples were chemically derivatized via the methyl chloroformate (MCF) method based on the recommendations published by Smart et al. [20]. All cord plasma samples were analyzed in a single batch and derivatized compounds were examined by an Agilent GC7890B system using a ZB-1701 GC capillary column. An MSD5977A mass selective detector with the electron impact voltage set to 70 eV was applied to analyze the compounds. The GC column used for metabolite separation was the ZB-1701 GC capillary column (30 m × 250 μm id × 0.15 μm with a 5 m guard column, Phenomenex). The GC temperature was set up according to the protocol of Han et al. [21].

Metabolite identification, GS-MS data mining, and data normalization

The automated Mass Spectral Deconvolution and Identification System (AMDIS) was implemented for deconvolution. Metabolite identifications were determined based on the MS fragmentation pattern and respective chromatographic retention time via the in-house MCF mass spectral library established by Silas Villas Boas’s metabolomics laboratory in New Zealand. The peak height of the most abundant fragmented ion mass was selected using the MassOmics R scripts to extract the relative concentrations of the identified metabolites. The identified compound’s abundances were then normalized by the relative level of the internal standard (D4-alanine, D5-phenylalanine, or D2-tyrosine) in the corresponding sample. Median centering was performed to remove batch variation via nine QC samples (three QC samples per batch), and the dilution effects of plasma were corrected by a total ion chromatogram. Then, blank samples were used to subtract background contamination and carryover effects from identified metabolites.

Statistical analysis

Prior to statistical analysis, the cord plasma metabolite levels were adjusted by log transformation and Pareto scaling in order to provide the best Gaussian distribution for the dataset. Student’s t-test, non-parametric Mann-Whitney U test, Chi-square test, and Fisher’s exact test were performed in R to investigate prenatal clinical characteristics. Partial least squares discriminant analysis (PLS-DA) was performed using the MetaboAnalyst 5.0 package for R to screen for significant metabolites and identify metabolic profile differences between twin groups (http://www.metaboanalyst.ca). The linear logistic regression model was implicated in detecting differential metabolites between comparisons without the influence of confounding factor (e.g. gestation age and the way of conception) using the general linear model (glm) logistic regression package in R. Generalized estimating equation (GEE) modeling was utilized to examine metabolite correlations with birthweight discordance both within and between twin pairs. Metabolite pathway activity was calculated based on the KEGG database using Metaboanalyst 5.0. Metabolic networks of interest were reconstructed using the metascape package in Cytoscape (Version 3.9.1). Heatmaps, line graphs, and circus plots were made using the ggplot2 and GOplot R-packages [22, 23]. The intra- and inter-observer reliabilities of the CRL ultrasound measurement were determined using the Bland-Altman method [24] using Irr R package.

Results

Population characteristics

Characteristics of women pregnant with MC or DC twins included in this study are listed in Table 1. The postnatal outcomes of MC twins and DC twins are summarized in Table 2. No differences were observed between groups regarding pre-gestational maternal body mass index, weight gain during pregnancy, mode of delivery, maternal age, neonatal sex, Apgar score (1 and 5 min), and birth weight discrepancy with co-twins (comparison 1 and 2). In contrast, the gestational age at delivery, average birthweight, and mode of conception were significantly different between MC and DC groups (comparison 3). This was in part due to the fact that MC pregnancies were generally delivered earlier to avoid the higher risk of complications than DC pregnancies. In vitro fertilization (IVF) dominated the mode of conception in our DC pregnancies (79.17%, p < 0.001). Indicators of fetal development, including fetal abdominal circumference (p < 0.01), height (p < 0.05), fetal head circumference (p < 0.05), and birthweight (p < 0.01), were different between the larger twin and the smaller twin of DC groups.

Table 1 Comparison of the clinical characteristics of MC and DC groups
Table 2 Comparison of postnatal outcomes in the MC and DC twins

Analysis of the umbilical cord plasma metabolome profiles in MC and DC twins

The metabolites of twin umbilical cord plasma samples were identified using our in-house MCF mass spectral library and NIST library (https://www.nist.gov/nist-research-library) with the inter-assay coefficient of variation in QC samples ranging from 0.9 to 23.0% (See Supplementary Table 1). Representative metabolites are indicated on a GC-chromatogram displayed in Supplementary Fig. 2. Umbilical cord plasma of DC co-twins demonstrated better global metabolomic separation of larger and smaller twins compared to the MC co-twins (Fig. 2a and b). Lower plasma levels of EDTA and higher plasma levels of nicotinamide discriminated larger from smaller twins within DC pairs with a p-value less than 0.05 (C1, Fig. 3), while octadecane, hexadecane,2,6,11,15-tetramethyl, and 2-hydroxyglutaramic acid was significantly higher in smaller MC twins (p < 0.05) relative to their larger co-twin (C2, Fig. 3). PLS-DA score plots (Fig. 2) illustrate that the metabolic profiles of MC twins as a group were generally different from DC twins (Fig. 2C). The results of comparison 3 (C3) between DC and MC displayed a valid model performance (Accuracy = 0.93, R2 = 0.92, Q2 = 0.59); the major three latent variables accounted for 10.9%, 5.6%, and 4.4% of the variation in the metabolite levels. The most distinct separations and the most valid LOOCV were demonstrated when comparing DC and MC (Fig. 2, C). Adjusted logistic regression was performed for all comparisons to account for the potential influencing factors of gestational age at delivery and mode of conception. This revealed 23 significant plasma metabolites that contributed to the separation of the MC twins and DC twins with a p-value and q-value less than 0.05 and 0.05 respectively (C3; Fig. 3). Among them, two amino acids and one saturated fatty acid were found at higher levels in DC twins.

Fig. 2
figure 2

Partial least squares discriminant analysis (PLS-DA) of the umbilical cord plasma metabolome between the three twin comparisons, including a measure of prediction model performance (right bar graphs). The right bar graphics are leave-one-out cross-validations (LOOCV), where R2 indicates how well the model explains the data and Q2 indicates the reproducibility of the PLS-DA model. The red asterisk indicates the best classifier. List of abbreviations; MC = Monochorionic; DC = Dichorionic; L = Larger twin; S = Smaller twin

Fig. 3
figure 3

The heatmap illustrates the differences in the umbilical cord plasma metabolome between the three groups of DC twins and MC twins. The relative concentrations of umbilical cord plasma metabolites are shown via a log2(foldchange). Red color blocks represent higher metabolites levels in dividend groups than the divisor groups, whereas blue color blocks represent lower metabolites levels in dividend groups than the divisor groups. Only the metabolites with a p-value less than 0.05 (The logistic regression adjusted for gestational age and the way of conception) were displayed. The metabolites with a significant p-value less than 0.05 and a q-value less than 0.05 were marked via a black asterisk. List of abbreviations: DC = Dichorionic; MC = Monochorionic; L = Larger twin; S = Smaller twin

Within and between pair correlation of metabolites with birthweight discordance and birthweight

A GEE regression model was applied to measure the level of correlation between plasma metabolites and birth weight discordance within pairs (smaller vs. larger twin) and between twin pairs, accounting for the individual (twin) and shared (maternal) factors (Fig. 4). Octadecane and 2,6,11,15-tetramethyl-hexadecane were negatively associated with birthweight discordance within MC twin pairs only, while decamethyl-cyclopentasiloxane and 2-oxoadipic acid were positively associated with birthweight discrepancy in DC twins only. 2 − hydroxyglutaramic acid levels were negatively associated with birthweight within both MC and DC twin pairs. A total of 70 metabolites were found to be associated with the birthweight following twin pair comparisons. This included unsaturated fatty acids, TCA cycle intermediates, antioxidants, as well as most amino acids, amino acids derivatives, TCA cycle intermediates derivatives, and organic acids (Fig. 4).

Fig. 4
figure 4

Correlation of birth weight within (weight discordance within larger and smaller co-twin, left column) and between (average birth weight between twin pairs, right column) twin pairs of umbilical cord plasma metabolites detected from DC and MC twins, analyzed using a generalized estimating equation (GEE). The red lines represent the 95% confidence intervals for the correlation of metabolites with weight discordance in the DC umbilical cord plasma. The blue lines represent the 95% confidence intervals for correlating metabolites with weight discordance in MC umbilical cord plasma. The center dotted line in each column indicates 0 correlation; metabolites to the right of the dotted line are positively correlated with weight discordance, whereas metabolites to the left of the dotted line are negatively correlated with weight discordance. In addition, the distance between the metabolites and the dotted line represents the strength of the correlation. Metabolites are classified in accordance with their chemical properties, and only the metabolites significantly correlated with birth weight (p-value < 0.05) are plotted

Pathway enrichment analysis of MC twins and DC twins

The differences in metabolic pathway activities derived from the identified metabolites in umbilical cord plasma samples appear to be associated with chorionicity. The majority of metabolic pathways we observed were upregulated in comparison 3 (DC/MC) (Fig. 5). Lipid metabolism, one metabolism of the endocrine system, four pathways associated with amino acids metabolism, three metabolisms of the nervous system, and two energy metabolisms were found in lower levels in the MC compared to the DC twins. However, no significant metabolic pathway differences were detected in comparison 2 (MC-L/MC-S). Vitamin digestion and absorption were upregulated in larger DC twins (p = 0.069). The significant pathways were linked to their shared metabolites and reconstructed into a metabolic network and a chord plot based on the KEGG metabolic framework via Cytoscape (Fig. 6). Eight metabolites were significantly different between DC and MC twins including tryptophan, 2-oxobutanoate, cysteine, glutathione, dodecanoic acid, oxaloacetate, 2-oxoglutarate, and 4-methyl-2-oxopentanoate (Fig. 6a). Four significant metabolites were involved in the glutamate metabolism including cysteine, glutathione, oxaloacetate, and 2-oxoglutarate (Fig. 6b).

Fig. 5
figure 5

The predicted metabolic activity in the umbilical cord plasma associated with DC twins and MC twins was illustrated using log2(fold change). The black dotted line in each column indicates metabolic pathways in the divisor groups that were adjusted to 0. The red plots at positive values represent upregulated metabolic activity in dividend groups compared to the divisor groups, whereas the red plots at negative values represent downregulated metabolic activity in dividend groups compared to the divisor groups. The red dot sizes represent the enrichment ratio of pathway computed by metabolite hits. Only the metabolic pathways with a significant p-value less than 0.05 (Logistic regression) and a q-value less than 0.1 (false discovery rate) are plotted. The pathways with a significant p-value less than 0.05 and a q-value less than 0.05 were marked via red asterisks. List of abbreviations: DC = Dichorionic; MC = Monochorionic; L = Larger twin; S = Smaller twin

Fig. 6
figure 6

The metabolic networks represented in the umbilical cord plasma metabolome of uncomplicated MC twins. (a) The two-dimensional network was constructed using the metabolic pathways that were significantly different between the MC twins and the DC twins. The green circles are significantly different metabolites between the DC-L and DC-S abbreviated as C1. The orange circles are metabolites that were significantly different between the MC-L and MC-S abbreviated as C2. The yellow circles are metabolites that were significantly different between the DC and MC abbreviated as C3. The blue squares are the metabolic reactions associated with significant metabolites. All purple circles are unidentified metabolites that were directly linked to identified metabolites. The arrowheads indicate the direction of the metabolic reactions. (b) A chord plot displays how metabolites with p < 0.05 participate in different significant metabolic pathways

Inter- and intra-observer variations for ultrasound measurements

The determination of gestational age using CRL ultrasound indicated that inter-observer reliability agreement between observers 1 and 2 was 0.994 (95% CI 0.986–0.998), whilst the intra-observer agreement of CRL measurement was 0.997 (95% CI 0.992–0.999) and 0.996 (95% CI 0.989–0.998) for observer 1 and 2 correspondingly (Supplementary Fig. 3). These outcomes indicated good reproducibility of the CRL ultrasound measurement to determine gestational age in our hospital.

Discussion

In this study, we aimed to identify metabolomic differences in newborn twins associated with chorionicity (MC vs. DC) and to relate these to birthweight and birthweight discordance. The different metabolic profiles identified highlight the difference in the intrauterine growth environment experienced by twins in association with the mode of placentation. We observed more similarities in infant outcomes and metabolite profiles in MC co-twin pairs compared to DC co-twin pairs. Higher concentrations of most amino acids and organic acids were revealed in DC twins compared to MC twins associated with chorionicity. The accumulations of most nutritional metabolites and exogenous substances were positively correlated with the birthweight of both DC and MC twins. In contrast, 2 − hydroxyglutaramic acid (2-HG) was negatively correlated with birthweight discordance of both DC and MC twins.

The most interesting findings to emerge from the within-twin pair’s analysis were that higher plasma concentrations of 2-HG only tended to be associated with the metabolite differences identified in MC co-twins. In contrast, the metabolite nicotinamide contributed to differ in DC twin pairs only. GSH metabolism, tryptophan metabolism, and fatty acid metabolism were significantly downregulated in MC twins relative to DC twins in the between-pair analyses (Fig. 7). These metabolic alterations in the umbilical cord blood appear to serve as an essential reflection of the unequal blood distribution resulting in differentiation of chorionicity, thus influencing twins’ growth.

Fig. 7
figure 7

Changes in the metabolic phenotypes of twin cord blood between differential chorionicities. Uncomplicated monochorionic twin pregnancy shows three types of anastomoses: the veno-venous (VV) anastomosis is bi-directional occurred between the red veins (red star), large arterio-arterial (AA) anastomosis is bi-directional occurred between the blue arteries (blue star), and arterio-venous (AV) anastomosis is unidirectional from blue arteries to red veins (purple star). The umbilical cord plasma discrepancies within twin pairs and between twin pairs were exhibited. Within twin pairs: a higher 2-HG level was observed in the MC smaller twin compared to the larger one; a higher nicotinamide level was observed in the DC larger twin compared to the smaller one. Between MC and DC twins: the GSH metabolism, tryptophan metabolism, and linolenic acid metabolism were downregulated in the MC twins. List of abbreviations: DC = Dichorionic; MC = Monochorionic; L = Larger twin; S = Smaller twin; 2-HG, 2 − hydroxyglutaramic acid; α-KG, α-Ketoglutaric acid; HGDH, 2-Hydroxyglutarate dehydrogenase; GSH, glutathione

Differential metabolites within MC twin pairs

Growth in monochorionic twins is determined by the same genetic potential, as well as by the placental sharing and vascular anastomoses [25,26,27]. The present study demonstrated that less birthweight discrepancy was observed in uncomplicated MC twin pairs compared to DC twin pairs. After eliminating the potential influence of gestation age and IVF-ET, we found that the spatial images generated by PLS-DA displayed more similar metabolite profiles within MC twins than in DC twins (Fig. 2). Indeed, placental vascular anastomoses that connect the two circulations of monochorionic twin pregnancies account for the equalization of oxygen and nutritional supply between larger and smaller twins [25, 26]. Surprisingly, higher cord blood concentrations of 2-HG, octadecane, and 2,6,11,15-tetramethyl hexadecane were found in MC smaller babies than in their larger siblings (Fig. 3). A negative correlation between 2-HG levels and twin birth weight was also observed in this study (Fig. 4). Literature reported that a higher level of 2-HG has been identified as a result of abnormal metabolism under homeostasis [28]. The accumulation of 2-HG acted as a competitive inhibitor of alpha-ketoglutarate-dependent dioxygenases, causing profound metabolic and epigenetic dysregulation [29, 30]. A recent study showed that preventing the reduction of 2-HG levels after fertilization would impede the erasure of histone modifications such as H3K4me3 [31]. This highlights the potential role associated with epigenetic network remodelling of 2-HG in fetal development.

Differential metabolites within DC twin pairs

The second interesting observation of the current study is that nicotinamide was upregulated in DC larger compared to smaller DC twins, while no such disparity was observed within MC twin pairs (Fig. 3). Nicotinamide is an amide form of vitamin B3 as an indispensable nutrient often supplemented to pregnant women for embryonic development. It has been reported that nicotinamide exerts anti-inflammatory and antioxidative properties in the human placenta [32]. Inflammation plays a pivotal role during human pregnancy. F. Li et al. [33] demonstrated that dietary nicotinamide benefits both dams and pups with preeclampsia characterized by an excessive inflammatory response. A cross-sectional study design that included 626 mothers and their singleton offspring also reported that maternal vitamin B3 intake in pregnancy contributes to increasing birthweight [34]. It will be interesting to validate that adequate niacinamide supplementation may benefit the birthweight of twins.

The metabolic disparities between MC and DC twins related to vascular anastomoses

A hypothesis proposed by Sebire et al. [35] is that anastomoses between placental circulations in monochorionic twins develop randomly during early pregnancy. Subsequently, placental growth is influenced by spontaneous closure or disruption of anastomoses. Asymmetrical loss of anastomoses may trigger TTTS due to disturbed hemodynamics [35]. An imbalance of angiogenesis is more vulnerable to oxidative stress, which could increase the likelihood of intrauterine growth restriction, fetal death, spontaneous preterm labour, and preeclampsia. The findings from our umbilical plasma metabolome analysis illustrated that antioxidants (glutathione, GSH) were downregulated in MC twin pairs (Fig. 3). GSH is an intracellular tripeptide composed of glutamate, glycine, and cysteine, which protects the organism from oxidative stress by eliminating reactive oxygen species (ROS). Altered cysteine and glycine derivatives (Glycine, N-ethyl-N-(2-methoxyethoxycarbonyl)-, 2-methoxyethyl) were discovered in comparison 3, as illustrated in Fig. 3. The alterations of other essential metabolites involved in glutamate metabolism were also observed, including 2-oxoglutaric acid and oxaloacetic acid (Fig. 6b). In addition, our reconstructed metabolic network indicated that glutathione interacted with cysteine during pregnancy (Fig. 6a). Several studies supported that the regulation of antioxidant levels and their metabolites may be related to the maintenance of the redox state and blood vessel formation [17, 36, 37]. These findings support that the disproportionate vascular anastomoses in the MC twins are associated with oxidative stress, thus affecting the regulation of the antioxidants and their metabolites.

Interstitial vascular anastomose in twin pregnancy has also been known as a strong factor influencing the blood distribution and results in disturbed nutrient supply. A prospective observational study provided evidence that underlying hemodynamic changes also occur in uncomplicated MC twins [38]. Previous studies reported that the higher rates of abnormal placental vascular observed in MC twins compared to DC twins, even in the uncomplicated MC twin pregnancies, were associated with neonatal encephalopathy and cerebral palsy [39,40,41]. We observed that tryptophan metabolism was downregulated in MC twins (Fig. 5). Tryptophan is an essential amino acid, which must be obtained exclusively from a dietary source in humans. Disrupted tryptophan metabolism was associated with the decrease of neuronal protection for fetal growth and immunosuppressive against fetal rejection [42]. Consistent with the literature, this research found that the metabolic pathways related to the nervous system were impaired in MC twins (Fig. 5). Furthermore, we observed that three unsaturated fatty acids and two saturated fatty acids were downregulated in the MC twins (Fig. 3). Among them, linolenic acid and its derivatives were positively associated with the birthweight of twins (Fig. 4). It has been reported that prenatal administration of linolenic acid benefits neurogenesis and the cognitive abilities of mice with Down syndrome [43]. These results demonstrated that the placental vascular anastomoses in MC twins were related to the imbalanced nutrient distribution, thereby affecting their downstream metabolic pathway such as nervous system metabolism.

The accumulations of exogenous substances in MC and DC twins

Lastly, we found that a higher caffeine level was discovered in MC twins and the enrichment of caffeine in umbilical cord plasma was positively correlated with birth weight both in MC and DC twins (Figs. 3 and 4). Several studies suggested that caffeine was consistently elevated during pregnancy, and the elevation might be due to a slower caffeine metabolism in pregnant women rather than an increase in coffee intake [44, 45]. Together, the caffeine concentration of twins’ cord blood increases with gestational age and positively correlates with progressive fetal weight. Secondly, the xenobiotic compound decamethyl-cyclopentasiloxane (D5) was positively correlated with birth weight discordance between DC co-twins (Fig. 4). It is one of the cyclic siloxanes widely used in cosmetics and body care products. Although limited research has reported that the other cyclic siloxanes may impair the fertility and reproduction of female rats [46, 47]. The role of D5 in fetal development remains largely unknown. Therefore, our study suggested that the twins may suffer from a slower metabolism of some exogenous substances, like caffeine and D5. The mechanistic function of these accumulated metabolites in twins’ growth was clearly worth revealing.

Lastly, this study has presented several limitations. The metabolic profile of umbilical blood was collected at delivery, which is metabolically different from that at early or middle gestation period. Another concern was that the sample size was small (low power for cause-effect relationship) due to the high rate of complications in twin pregnancies. In the future, the sample size and different biospecimens should be expanded. Maternal blood composition, placental pathophysiological structure, and fetal metabolism should be further investigated to better understand the influence of chorionicity on maternal-fetal transfer and the growth discordance of twins.

Conclusions

To our knowledge, the present study is the first to investigate the metabolic profile of the umbilical cord plasma between uncomplicated MC and DC twin pregnancies. The main findings fill the gap in the scientific literature on the association between the metabolites in cord plasma and the development of uncomplicated twins. The discrepancy of umbilical cord plasma metabolome within and between MC and DC twins seems to have resulted from the alterations of metabolite allocation, such as 2-HG, nicotinamide, GSH, tryptophan, linolenic acid, and is related to discordant birthweight. The underlying mechanisms of vascular anastomosis and accumulated exogenous substances in the development of twins were still worthy of further investigation.

Data availability

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

Abbreviations

AC:

Abdominal circumference

AMDIS:

Automated Mass Spectral Deconvolution and Identification System

AUC:

Area Under Curve

CRL:

Crown rump length

DC:

Dichorionic

EFW:

Estimated fetal weight

FDR:

False discovery rate

GC-MS:

Gas chromatography-mass spectrometry

GEE:

Generalized estimating equation.

GLM:

General linear model

GSH:

Glutathione

HC:

Head circumference

HG:

Hydroxyglutaramic acid

ISUOG:

International Society of Ultrasound in Obstetrics and Gynecology

IUFD:

Intrauterine fetal death

IVF-ET:

In Vitro Fertilization & Embryo Transfer

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LOOCV:

Leave-one-out cross validation

MC:

Monochorionic

MCF:

Methyl chloroformate

PLS-DA:

Partial least squares discriminant analysis

QC:

Quality control

ROC:

Receiver operating characteristic

ROS:

Reactive oxygen species

sIUGR:

Selective intrauterine growth restriction

TAPS:

Twin anemia polycythemia sequence

TTTS:

Twin-to-twin transfusion syndrome

References

  1. National Collaborating Centre for Ws, Children’s H. : National Institute for Health and Clinical Excellence: Guidance. Multiple Pregnancy: The Management of Twin and Triplet Pregnancies in the Antenatal Period. edn. London:RCOG PressCopyright © 2011, National Collaborating Centre for Women’s and Children’s Health.; 2011.

  2. Gao W. Chinese twins register center. 2019.

  3. Deng C, Dai L, Yi L, Li X, Deng K, Mu Y, Wang K, Tao J, Li Q, Xu L. Temporal trends in the birth rates and perinatal mortality of twins: a population-based study in China. PLoS ONE. 2019;14(1):e0209962.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Google Scholar 

  5. DeJesus Allison SO, Javitt MC, Glanc P, Andreotti RF, Bennett GL, Brown DL, Dubinsky T, Harisinghani MG, Harris RD, Mitchell DG, et al. ACR appropriateness Criteria® multiple gestations. Ultrasound Q. 2012;28(2):149–55.

    Article  PubMed  Google Scholar 

  6. Peter C, Wenzlaff P, Kruempelmann J, Alzen G, Bueltmann E, Gruessner SE. Perinatal morbidity and early neonatal mortality in twin pregnancies. Open J Obstet Gynecol. 2013;03(01):78–89.

    Article  Google Scholar 

  7. Ortibus E, Lopriore E, Deprest J, Vandenbussche FP, Walther FJ, Diemert A, Hecher K, Lagae L, De Cock P, Lewi PJ, et al. The pregnancy and long-term neurodevelopmental outcome of monochorionic diamniotic twin gestations: a multicenter prospective cohort study from the first trimester onward. Am J Obstet Gynecol. 2009;200(5):494e491–498.

    Article  Google Scholar 

  8. Russo FM, Pozzi E, Pelizzoni F, Todyrenchuk L, Bernasconi DP, Cozzolino S, Vergani P. Stillbirths in singletons, dichorionic and monochorionic twins: a comparison of risks and causes. Eur J Obstet Gynecol Reprod Biol. 2013;170(1):131–6.

    Article  PubMed  Google Scholar 

  9. Hack KE, Derks JB, Elias SG, Franx A, Roos EJ, Voerman SK, Bode CL, Koopman-Esseboom C, Visser GH. Increased perinatal mortality and morbidity in monochorionic versus dichorionic twin pregnancies: clinical implications of a large dutch cohort study. BJOG. 2008;115(1):58–67.

    Article  CAS  PubMed  Google Scholar 

  10. Zhao D, Lipa M, Wielgos M, Cohen D, Middeldorp JM, Oepkes D, Lopriore E. Comparison between Monochorionic and Dichorionic Placentas with Special attention to vascular anastomoses and placental share. Twin Res Hum Genet. 2016;19(3):191–6.

    Article  PubMed  Google Scholar 

  11. Cheong-See F, Schuit E, Arroyo-Manzano D, Khalil A, Barrett J, Joseph KS, Asztalos E, Hack K, Lewi L, Lim A et al. Prospective risk of stillbirth and neonatal complications in twin pregnancies: systematic review and meta-analysis. BMJ. 2016:i4353.

  12. Aviram A, Lipworth H, Asztalos EV, Mei-Dan E, Melamed N, Cao X, Zaltz A, Hvidman L, Barrett JFR. Delivery of monochorionic twins: lessons learned from the Twin Birth Study. Am J Obstet Gynecol. 2020;223(6):916.e911–916.e919.

  13. Litwinska E, Syngelaki A, Cimpoca B, Sapantzoglou I, Nicolaides KH. Intertwin discordance in fetal size at 11–13 weeks’ gestation and pregnancy outcome. Ultrasound Obstet Gynecol. 2020;55(2):189–97.

    Article  CAS  PubMed  Google Scholar 

  14. Wishart DS. Metabolomics for investigating physiological and pathophysiological processes. Physiol Rev. 2019;99(4):1819–75.

    Article  CAS  PubMed  Google Scholar 

  15. Cosmi E, Visentin S, Favretto D, Tucci M, Ragazzi E, Viel G, Ferrara SD. Selective intrauterine growth restriction in monochorionic twin pregnancies: markers of endothelial damage and metabolomic profile. Twin Res Hum Genet. 2013;16(4):816–26.

    Article  PubMed  Google Scholar 

  16. Wang L, Han TL, Luo X, Li S, Young T, Chen C, Wen L, Xu P, Zheng Y, Saffery R, et al. Metabolic biomarkers of Monochorionic Twins Complicated with Selective Intrauterine Growth Restriction in Cord plasma and placental tissue. Sci Rep. 2018;8(1):15914.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Yang J, Wei Y, Qi H, Yin N, Yang Y, Li Z, Xu L, Wang X, Yuan P, Li L, et al. Neonatal hair profiling reveals a metabolic phenotype of monochorionic twins with selective intrauterine growth restriction and abnormal umbilical artery flow. Mol Med. 2020;26(1):37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Salomon LJ, Alfirevic Z, Bilardo CM, Chalouhi GE, Ghi T, Kagan KO, Lau TK, Papageorghiou AT, Raine-Fenning NJ, Stirnemann J, et al. ISUOG practice guidelines: performance of first-trimester fetal ultrasound scan. Ultrasound Obstet Gynecol. 2013;41(1):102–13.

    Article  CAS  PubMed  Google Scholar 

  19. Khalil A, Rodgers M, Baschat A, Bhide A, Gratacos E, Hecher K, Kilby MD, Lewi L, Nicolaides KH, Oepkes D, et al. ISUOG Practice Guidelines: role of ultrasound in twin pregnancy. Ultrasound Obstet Gynecol. 2016;47(2):247–63.

    Article  CAS  PubMed  Google Scholar 

  20. Smart KF, Aggio RB, Van Houtte JR, Villas-Bôas SG. Analytical platform for metabolome analysis of microbial cells using methyl chloroformate derivatization followed by gas chromatography-mass spectrometry. Nat Protoc. 2010;5(10):1709–29.

    Article  CAS  PubMed  Google Scholar 

  21. Han TL, Cannon RD, Gallo SM, Villas-Bôas SG. A metabolomic study of the effect of Candida albicans glutamate dehydrogenase deletion on growth and morphogenesis. NPJ Biofilms Microbiomes. 2019;5(1):13.

    Article  PubMed  PubMed Central  Google Scholar 

  22. H. W: ggplot2: elegant graphics for data analysis. New York: Springer 2009,NA.

  23. Walter W, Sánchez-Cabo F, Ricote M. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics. 2015;31(17):2912–4.

    Article  CAS  PubMed  Google Scholar 

  24. Souka AP, Pilalis A, Papastefanou I, Salamalekis G, Kassanos D. Reproducibility study of crown-rump length and biparietal diameter measurements in the first trimester. Prenat Diagn. 2012;32(12):1158–65.

    Article  PubMed  Google Scholar 

  25. Denbow ML, Cox P, Taylor M, Hammal DM, Fisk NM. Placental angioarchitecture in monochorionic twin pregnancies: relationship to fetal growth, fetofetal transfusion syndrome, and pregnancy outcome. Am J Obstet Gynecol. 2000;182(2):417–26.

    Article  CAS  PubMed  Google Scholar 

  26. Lewi L, Cannie M, Blickstein I, Jani J, Huber A, Hecher K, Dymarkowski S, Gratacós E, Lewi P, Deprest J. Placental sharing, birthweight discordance, and vascular anastomoses in monochorionic diamniotic twin placentas. Am J Obstet Gynecol. 2007;197(6):587e581–588.

    Article  Google Scholar 

  27. Fick AL, Feldstein VA, Norton ME, Wassel Fyr C, Caughey AB, Machin GA. Unequal placental sharing and birth weight discordance in monochorionic diamniotic twins. Am J Obstet Gynecol. 2006;195(1):178–83.

    Article  PubMed  Google Scholar 

  28. Engqvist MK, Eßer C, Maier A, Lercher MJ, Maurino VG. Mitochondrial 2-hydroxyglutarate metabolism. Mitochondrion. 2014;19 Pt B:275–281.

  29. Xu W, Yang H, Liu Y, Yang Y, Wang P, Kim SH, Ito S, Yang C, Wang P, Xiao MT, et al. Oncometabolite 2-hydroxyglutarate is a competitive inhibitor of α-ketoglutarate-dependent dioxygenases. Cancer Cell. 2011;19(1):17–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Chung C, Sweha SR, Pratt D, Tamrazi B, Panwalkar P, Banda A, Bayliss J, Hawes D, Yang F, Lee HJ, et al. Integrated metabolic and epigenomic reprograming by H3K27M mutations in diffuse intrinsic pontine gliomas. Cancer Cell. 2020;38(3):334–349e339.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Zhao J, Yao K, Yu H, Zhang L, Xu Y, Chen L, Sun Z, Zhu Y, Zhang C, Qian Y, et al. Metabolic remodelling during early mouse embryo development. Nat metabolism. 2021;3(10):1372–84.

    Article  CAS  Google Scholar 

  32. Lappas M, Permezel M. The anti-inflammatory and antioxidative effects of nicotinamide, a vitamin B(3) derivative, are elicited by FoxO3 in human gestational tissues: implications for preterm birth. J Nutr Biochem. 2011;22(12):1195–201.

    Article  CAS  PubMed  Google Scholar 

  33. Li F, Fushima T, Oyanagi G, Townley-Tilson H, Sato E, Nakada H, Oe Y, Hagaman J, Wilder J, Li M, et al. Nicotinamide benefits both mothers and pups in two contrasting mouse models of preeclampsia. Proc Natl Acad Sci USA. 2016;113(47):13450–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Souza R, Miranda C, Dos Santos LC. Maternal vitamin B(3) and C intake in pregnancy influence birth weight at term. Nutrition. 2021;91–92:111444.

  35. N.J. Sebire ea: Twin-to-Twin Transfusion Syndrome Results From Dynamic Asymmetrical Reduction in Placental Anastomoses:A Hypothesis. 2001.

  36. Schneider D, Hernández C, Farías M, Uauy R, Krause BJ, Casanello P. Oxidative stress as common trait of endothelial dysfunction in chorionic arteries from fetuses with IUGR and LGA. Placenta. 2015;36(5):552–8.

    Article  CAS  PubMed  Google Scholar 

  37. Herrera EA, Cifuentes-Zúñiga F, Figueroa E, Villanueva C, Hernández C, Alegría R, Arroyo-Jousse V, Peñaloza E, Farías M, Uauy R, et al. N-Acetylcysteine, a glutathione precursor, reverts vascular dysfunction and endothelial epigenetic programming in intrauterine growth restricted guinea pigs. J Physiol. 2017;595(4):1077–92.

    Article  CAS  PubMed  Google Scholar 

  38. Torres X, Bennasar M, García-Otero L, Martínez-Portilla RJ, Valenzuela-Alcaraz B, Crispi F, Goncé A, Gratacós E, Figueras F, Martínez JM. Uncomplicated Monochorionic Twins: Two Normal Hearts Sharing One Placenta. J Clin Med. 2020;9(11).

  39. Weiner E, Barber E, Feldstein O, Dekalo A, Schreiber L, Bar J, Kovo M. Placental histopathology differences and neonatal outcome in Dichorionic-Diamniotic as compared to Monochorionic-Diamniotic Twin Pregnancies. Reprod Sci. 2018;25(7):1067–72.

    Article  PubMed  Google Scholar 

  40. Redline RW. Severe fetal placental vascular lesions in term infants with neurologic impairment. Am J Obstet Gynecol. 2005;192(2):452–7.

    Article  PubMed  Google Scholar 

  41. Redline RW. Placental pathology and cerebral palsy. Clin Perinatol. 2006;33(2):503–16.

    Article  PubMed  Google Scholar 

  42. Badawy AA. Tryptophan metabolism, disposition and utilization in pregnancy. Biosci Rep. 2015;35(5).

  43. García-Cerro S, Rueda N, Vidal V, Puente A, Campa V, Lantigua S, Narcís O, Velasco A, Bartesaghi R, Martínez-Cué C. Prenatal administration of oleic acid or linolenic acid reduces neuromorphological and cognitive alterations in Ts65dn down syndrome mice. J Nutr. 2020;150(6):1631–43.

    Article  PubMed  Google Scholar 

  44. Liang L, Rasmussen MH, Piening B, Shen X, Chen S, Röst H, Snyder JK, Tibshirani R, Skotte L, Lee NC, et al. Metabolic Dynamics and Prediction of gestational age and time to delivery in pregnant women. Cell. 2020;181(7):1680–1692e1615.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Knutti R, Rothweiler H, Schlatter C. Effect of pregnancy on the pharmacokinetics of caffeine. Eur J Clin Pharmacol. 1981;21(2):121–6.

    Article  CAS  PubMed  Google Scholar 

  46. Quinn AL, Dalu A, Meeker LS, Jean PA, Meeks RG, Crissman JW, Gallavan RH Jr, Plotzke KP. Effects of octamethylcyclotetrasiloxane (D4) on the luteinizing hormone (LH) surge and levels of various reproductive hormones in female Sprague-Dawley rats. Reprod Toxicol. 2007;23(4):532–40.

    Article  CAS  PubMed  Google Scholar 

  47. Siddiqui WH, Stump DG, Plotzke KP, Holson JF, Meeks RG. A two-generation reproductive toxicity study of octamethylcyclotetrasiloxane (D4) in rats exposed by whole-body vapor inhalation. Reprod Toxicol. 2007;23(2):202–15.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by grants from Chongqing Municipal Education Commission (KJZD-K202100407), Chongqing Science & Technology Commission (cstc2021jcyj-msxmX0213), the National Natural Science Foundation of China (No. 81871185), Peking University New Engineering Crossover Youth Project (PKU2023XGK015), the Natural Science Foundation of Beijing (7212133), and the Kuanren Talents Program of the Second Affiliated Hospital of Chongqing Medical University.

Author information

Authors and Affiliations

Authors

Contributions

XL and JY performed the research; XL, YY, TLH, and FL interpreted and analyzed the data; XL, JY, TLH, RR, LH, XD, YW, RS, BN and HM reviewed and edited the manuscript; XL and TLH wrote the manuscript. YW and TLH were responsible for conception and design of the project. YW, LH and TLH were responsible for supervision and project administration.

Corresponding authors

Correspondence to Yuan Wei, Lina Hu or Ting-Li Han.

Ethics declarations

Ethics approval and consent to participate

This research was conducted in accordance with the Declaration of Helsinki and approved by the Ethical Committee of Peking University Third Hospital and Second Affiliated Hospital of Chongqing Medical University. Informed consent was signed by each pregnant participant included in this study.

Consent for publication

Informed consent was signed by each pregnant participant included in this study.

Competing interests

The authors declare no financial and non-financial competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1 Table 1:

Identified metabolites

Supplementary Material 2 Fig. 1:

Flowchart of the selection of study participant

Supplementary Material 3 Fig. 2:

Representative total ion chromatogram (TIC) of the umbilical cord plasma metabolome

Supplementary Material 4 Fig. 3:

Intra-observer 1 variability (a), intra-observer 2 variability (b), and inter-observer variability (c) for determining gestational age. The proportional difference from the average for the crown?rump lengths (CRL) of 20 larger twins were diagnosed by two independent registered sonographers (observer 1 and 2). The upper and lower red dot lines mean the 2.5th and 97.5th percentiles for limits of agreement. The middle black line is the mean

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

Liu, X., Yang, J., Ran, R. et al. Chorionicity-associated variation in metabolic phenotype of cord blood in twin. Nutr Metab (Lond) 20, 31 (2023). https://doi.org/10.1186/s12986-023-00744-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12986-023-00744-1

Keywords