Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Postpartum plasma metabolomic profile among women with preeclampsia and preterm delivery: implications for long-term health

BMC Medicine volume 18, Article number: 277 (2020) Cite this article

Abstract

Background

Preeclampsia and preterm delivery (PTD) are believed to affect women’s long-term health including cardiovascular disease (CVD), but the biological underpinnings are largely unknown. We aimed to test whether maternal postpartum metabolomic profiles, especially CVD-related metabolites, varied according to PTD subtypes with and without preeclampsia, in a US urban, low-income multi-ethnic population.

Methods

This study, from the Boston Birth Cohort, included 980 women with term delivery, 79 with medically indicated PTD (mPTD) and preeclampsia, 52 with mPTD only, and 219 with spontaneous PTD (sPTD). Metabolomic profiling in postpartum plasma was conducted by liquid chromatography-mass spectrometry. Linear regression models were used to assess the associations of each metabolite with mPTD with preeclampsia, mPTD only, and sPTD, respectively, adjusting for pertinent covariates. Weighted gene coexpression network analysis was applied to investigate interconnected metabolites associated with the PTD/preeclampsia subgroups. Bonferroni correction was applied to account for multiple testing.

Results

A total of 380 known metabolites were analyzed. Compared to term controls, women with mPTD and preeclampsia showed a significant increase in 36 metabolites, mainly representing acylcarnitines and multiple classes of lipids (diacylglycerols, triacylglycerols, phosphocholines, and lysophosphocholines), as well as a decrease in 11 metabolites including nucleotides, steroids, and cholesteryl esters (CEs) (P < 1.3 × 10−4). Alterations of diacylglycerols, triacylglycerols, and CEs in women with mPTD and preeclampsia remained significant when compared to women with mPTD only. In contrast, the metabolite differences between women with mPTD only and term controls were only seen in phosphatidylethanolamine class. Women with sPTD had significantly different levels of 16 metabolites mainly in amino acid, nucleotide, and steroid classes compared to term controls, of which, anthranilic acid, bilirubin, and steroids also had shared associations in women with mPTD and preeclampsia.

Conclusion

In this sample of US high-risk women, PTD/preeclampsia subgroups each showed some unique and shared associations with maternal postpartum plasma metabolites, including those known to be predictors of future CVD. These findings, if validated, may provide new insight into metabolomic alterations underlying clinically observed PTD/preeclampsia subgroups and implications for women’s future cardiometabolic health.

Peer Review reports

Background

Pregnancy is believed to be a vulnerable time period for women’s long-term health. While ample evidence suggests that preterm delivery (PTD) and other pregnancy complications such as preeclampsia are important risk factors for the health of the offspring, much less appreciated is that these same traits may also affect the long-term health of the mother. Pregnancy complications and PTD are associated with risk of having the same conditions in subsequent pregnancies, but also lead to an increased long-term risk of adverse health outcomes, including cardiovascular diseases (CVDs) [1,2,3,4], the leading causes of death for women in the USA [5].

All-cause PTDs are associated with a 1.5- to 3-fold higher risk of CVD [1,2,3,4, 6]. Previous studies have shown that much of this association is accounted for by metabolic disorders during pregnancy, especially preeclampsia (which explains about 24-26% of the association) [1], or by metabolic disorders postpartum (which may explain about 13-15% of the association) [3]; however, a substantial portion of the risk remains unaccounted for. The PTD-CVD associations persist even in pregnancies uncomplicated by preeclampsia or hypertensive disorders [3, 6, 7], indicating the need to explore other pathways via which PTD and CVD are linked. PTD is heterogeneous, with a variety of different PTD subtypes (i.e., spontaneous vs. medically indicated PTD), with and without preeclampsia, which may affect CVD risk through different pathways. However, such data are quite limited.

We hypothesized that women who experienced different PTD subtypes with and without preeclampsia may exhibit CVD risk factors or biomarkers early in the postpartum period, long before the development of CVD. The identification of such early biomarkers may provide new insight into the biological underpinning of the link between PTD/preeclampsia and women’s future cardiometabolic health and provide targets for prevention and intervention in high-risk women. Metabolomic profiling, via the high-throughput assessment of a large number of circulating small molecule metabolites across multiple pathways, may offer advantages in elucidating biological pathways and biomarker discovery. Recent metabolomic studies have successfully identified multiple metabolites of different classes as early biomarkers for CVD, including branched-chain and aromatic amino acids (i.e., leucine, isoleucine, valine, phenylalanine, and tyrosine) [8,9,10], short- and medium-chain acylcarnitines [11,12,13], and metabolites of lipid classes [14,15,16,17]. Little is known about whether such metabolites vary in women with PTD vs. term delivery with or without preeclampsia.

In this study, by applying state-of-the-art metabolomic approaches, we quantified maternal metabolomic profiles in postpartum plasma (collected within 1–3 days after delivery) of 1411 women enrolled from the Boston Birth Cohort (BBC), an inner-city multi-ethnic longitudinal study cohort. We aimed to identify both shared and divergent metabolomic patterns (especially in those metabolites associated with CVD risk) for different PTD subtypes with and without preeclampsia, and further to test whether such PTD/preeclampsia-associated metabolomic patterns may vary by maternal characteristics such as ethnicity and parity.

Methods

Study population

The study population was a subset of the BBC, an ongoing longitudinal cohort study begun in 1998 at the Boston Medical Center in Massachusetts. Detailed information about the BBC has been reported elsewhere [18, 19]. Briefly, mothers who delivered singleton live births were eligible for the study and were invited to participate within 1 to 3 days after delivery. We excluded pregnancies that were a result of in vitro fertilization, multiple gestations, those with fetal chromosomal abnormalities or major birth defects, and preterm deliveries due to non-obstetric factors (such as trauma). The BBC is a low-income patient population, with a relatively high proportion of PTD and preeclampsia. After obtaining written informed consent, each mother was interviewed using a standardized questionnaire to gather dietary and epidemiologic data, and their electronic medical records (EMRs) were abstracted. Maternal blood samples were obtained within 24–72 h after delivery. The study protocol was approved by the Institutional Review Boards of Boston University Medical Center and the Johns Hopkins Bloomberg School of Public Health.

Definition of preterm delivery and preeclampsia

Gestational age was assessed by early prenatal ultrasound (< 20 weeks) or based on the first day of the last menstrual period as recorded in maternal EMRs if early prenatal ultrasound was not available [19]. Medically indicated PTD (mPTD) was defined as a delivery by medical induction or Cesarean section (c-section) at < 37 weeks without uterine contractions or rupture of membranes. The indications for medical induction or c-section included gestational complication (mainly preeclampsia), placenta abruption, placenta previa, oligohydramnios, previous c-section, intrauterine growth restriction, absent end diastolic flow, not reassuring fetal heartbeat, fetal distress, signs of reduced amniotic fluid index, or persistent fetal tachycardia. Spontaneous PTD (sPTD) was defined as a delivery occurring secondary to documented active preterm labor (uterine contractions with cervical effacement and dilation at < 37 weeks) or premature rupture of membranes at < 37 weeks without uterine contractions or both. We further divided sPTD as early sPTD (spontaneous delivery at < 33 weeks) or late sPTD (spontaneous delivery at 33–366/7 weeks). Physician-diagnosed preeclampsia, eclampsia, and hemolysis, elevated liver enzymes, and low platelet syndromes (HELLP) were extracted from the maternal EMRs and were further confirmed by review of all relevant medical records to meet the definition of preeclampsia according to the recent American College of Obstetricians and Gynecologist (ACOG) criteria [20].

We classified all participating women into one of the following PTD/preeclampsia subgroups: (1) women of term delivery without preeclampsia (including physician-diagnosed preeclampsia, eclampsia, and HELLP), as term controls; (2) women with both mPTD and preeclampsia; (3) women with mPTD but no preeclampsia (referred to as “mPTD only” hereafter); and (4) women with sPTD but no preeclampsia (referred to as “sPTD” hereafter, and further subdivided by early or late sPTD in supplemental analyses). Among the 1411 women with available metabolomic profiles, 9 were removed due to missing preeclampsia status and 72 (61 preeclamptic women with term delivery and 11 preeclamptic women with sPTD) were removed as they did not meet the criteria for either the control or the PTD subgroups as described above.

Epidemiological and clinical data

Using a standard maternal questionnaire interview [18], maternal epidemiological factors were collected including race/ethnicity, maternal age at delivery, highest education level, parity, smoking during pregnancy, alcohol drinking during pregnancy, illicit drug use during pregnancy, lifetime stress, previous history of PTD, and maternal birthplace (US-born/non-US-born). Maternal pregestational body mass index (BMI), calculated as self-reported weight (kg) divided by height squared (m2), was categorized into four groups: normal weight (< 25.0 kg/m2), overweight (25–29.9 kg/m2), obesity (≥ 30 kg/m2), and unknown. Clinical complication before pregnancy, including chronic hypertension and pregestational diabetes, was defined based on archived EMRs.

Metabolomic profiles and data processing steps

Maternal plasma samples, randomly distributed in each plate, were shipped to the Broad Institute of MIT and Harvard (Cambridge, MA, USA) for metabolomic profiling. Details of this technique can be found elsewhere [21, 22]. Briefly, two liquid chromatography tandem mass spectrometry techniques were applied, including hydrophilic interaction liquid chromatography (HILIC) analyses of water-soluble metabolites in the positive ionization mode and C8 chromatography with positive ion mode analyses of polar and non-polar plasma lipids. Metabolites in different lipid classes, including cholesterol esters (CEs), diacylglycerols (DAGs), triacylglycerols (TAGs), lysophosphatidylcholines (LPCs), lysophosphatidylethanolamines (LPEs), phosphatidylcholines (PCs), and phosphatidylethanolamines (PEs), were further classified based on the number of total acyl chain carbon atoms and double bond contents and were annotated as “C[number of total acyl chain carbon atoms]:[number of double bonds in fatty acid meoeties] [lipid class],” accordingly.

A pooled plasma study reference sample, comprised from all of the studied plasma samples, was evenly distributed throughout the study samples (every ~ 20 study samples) as a quality control step. Using these reference samples, a coefficient of variance (CV) was calculated for each metabolite. Internal standard peak areas were monitored for quality control and to ensure system performance throughout the analyses. A total of 432 known metabolites were successfully quantified, and 40 metabolites with CV > 20%, 10 metabolites unquantifiable in > 10% participants, and 2 xenobiotics were removed from downstream analyses. The remaining 380 metabolites, with any non-detectable values being imputed as one half of the minimal value, were inverse normally transformed to render the distributions approximately Gaussian and to remove the impact of outliers.

Data analyses

Maternal characteristics were compared between each of the three PTD/preeclampsia subgroups (as defined above) and term controls, using ANOVA tests for continuous variables and χ2 tests (or Fisher’s exact tests, when necessary) for categorical variables. To identify whether metabolites varied significantly in each PTD/preeclampsia subgroup compared to term controls, we fit a linear regression model with each metabolite as the outcome and the four-categorical index variable (term controls/mPTD and preeclampsia/mPTD only/sPTD) as the exposure, adjusting for conventional and clinical covariates including maternal age at delivery, maternal ethnicity/race, maternal birthplace (US-born versus non-US-born), pregestational BMI category, pregestational diabetes, chronic hypertension, marital status, highest education level, parity, maternal smoking during pregnancy, illicit drug use, lifetime stress, and fetal sex. Bonferroni correction was used to account for multiple testing.

We then performed weighted gene coexpression network analysis (WGCNA) using the WGCNA package in R [23], to reduce the number of comparisons made and minimize issues of metabolite collinearity. This method utilizes an unsupervised network-based approach to group metabolites into “modules” based on their correlation patterns. After hierarchical clustering, highly interconnected metabolites were assigned to the same module and was annotated as a unique color. The remaining metabolites (n = 51) that were not assigned to any of the established modules were removed from the subsequent analyses. For each module, the top hub metabolite was identified, which was the metabolite having the highest connectivity with other metabolites in the same module. Each participant was then assigned a “score” for each module, which was calculated as the first principal component of the metabolites within the module. The score for each module was then analyzed as the outcome in the subsequent analyses, using the linear regression models with the adjustment of covariates as described above. All analyses were performed using R Software, version 3.6.1.

Results

Population characteristics of the study participants

This study included 980 women with term delivery (or term controls) and 350 women with any of the following three PTD/preeclampsia subgroups: (1) mPTD and preeclampsia (n = 79), (2) mPTD only (n = 52), and (3) sPTD (n = 219). The population characteristics across the four study groups are summarized in Table 1. Compared to term controls, women with both mPTD and preeclampsia were older and were more likely to be obese before pregnancy, be stressed, and have chronic hypertension and pregestational diabetes, while women with sPTD were more likely to be unmarried, to be US-born, to smoke during pregnancy, and to have illicit drug use than term controls (all P < 0.05). As expected, c-section was significantly more common in women with mPTD, either with (65.8%) or without preeclampsia (80.8%) than in term controls (31.8%).

Table 1 Population characteristics of the 1330 women enrolled in the Boston Birth Cohort, stratified by preterm delivery subgroups
Full size table

Altered plasma metabolites in relation to different PTD/preeclampsia subgroups

Figure 1 and Table 2 present the metabolites that were significantly different from term controls in at least one PTD/preeclampsia subgroup. After Bonferroni correction for multiple testing (P < 1.3 × 10−4) and with term controls as the reference, women with mPTD and preeclampsia had significantly different patterns in 47 metabolites, including higher levels of 36 metabolites that were mainly amino acids (n = 7), total carnitine or acylcarnitines (C2, C3, C5-DC, C8, n = 4), nucleotides (n = 4), DAGs (C34:1, C34:1 DAG/TAG fragment, C34:2, C36:2, C36:3, C38:4, and C38:5, n = 7), TAGs (C50:2, C52:2, C52:3, n = 3), LPCs (C22:5 and C22:6, n = 2), and PCs (C38:6, C40:6, and C40:9, n = 3) and lower levels of 11 metabolites including CEs (C14:0, C16:0, C18:0, C18:2, and C18:3, n = 5), anthranilic acid, bilirubin, cytosine, 1-methylguanine, and 2 steroids (Table 2 and Additional file 1: Table S1). However, metabolomic differences were modest in women with mPTD only, who had lower levels of 4 PEs (C32:1, C34:2, C36:2, and C36:4) than term controls. In comparison, women with sPTD had significantly lower plasma levels of 13 metabolites (3 amino acids, 3 nucleotides, 2 steroids, C9 carnitine, anthranilic acid, bilirubin, 4-acetamidobutanoate, and trimethylamine-N-oxide) and higher levels of 3 metabolites (pyroglutamic acid, niacinamide, and phosphocholine), compared to the term controls. No significant associations were observed between sPTD and lipid metabolites. The identified associations remained largely unchanged when we further adjusted for the timing of maternal blood collection (postpartum day 1, 2, or 3 and above) in the regression models (data not shown).

Fig. 1
figure 1

Manhattan plot for the metabolomic differences in women with each preterm delivery (PTD) subgroup compared to women with term delivery. The upper panel presents metabolites that were higher, and the lower panel presents metabolites that were lower in women with medically indicated PTD and preeclampsia (a), in women with medically indicated PTD only (b), and in women with spontaneous PTD (c), respectively. The analyses were adjusted for maternal age at delivery, maternal ethnicity/race, maternal birthplace, pregestational BMI category, pregestational diabetes, chronic hypertension, marital status, highest education level, parity, smoking during pregnancy, illicit drug use, lifetime stress, and fetal sex. X-axis represents the metabolite classes. AA amino acid, AC acylcarnitine, FA fatty acid and conjugate, QA quaternary amine, CE cholesterol ester, DAG diacylglycerol, LPC lysophosphatidylcholine, LPC_P lysophosphatidylcholine plasmalogen, LPE lysophosphatidylethanolamine, MAG monoacylglycerol, PC phosphatidylcholine, PC_P phosphatidylcholine plasmalogen, PE phosphatidylethanolamine, PE_P phosphatidylethanolamine plasmalogen, PI phosphatidylinositol, SM sphingomyelin, TAG triacylglycerol

Full size image
Table 2 Maternal metabolites that were significantly different in women of each preterm delivery (PTD) subgroups, compared to women with term delivery
Full size table

We further examined the 47 metabolites that were associated with the group of mPTD and preeclampsia, compared to women with mPTD only. At P < 0.001, women with both mPTD and preeclampsia had higher levels of four DAGs (C34:1, C34:2, C38:4, and C52:2) and two TAGs (C50:2 and C52:2) and lower levels of four CEs (C16:0, C18:0, C18:2, and C18:3) than women with mPTD only, suggesting specific associations with preeclampsia.

Metabolite patterns related to PTD/preeclampsia subgroups by metabolomic network analyses

With WGCNA network analysis, 329 metabolites were grouped into 8 metabolic network modules, annotated as different colors. Metabolites in each module are presented in Additional file 1: Table S2. Figure 2 presents the correlations between each PTD/preeclampsia subgroup and each module score. With the adjustment of covariates and with P < 0.006 (= 0.05/8 modules) as the significance cutoff, we observed that women with mPTD and preeclampsia had significantly higher scores for the “yellow” (C16:0 CE as the top hub metabolite that was inversely correlated with the module score, P = 7.8 × 10−7), “black” (C14:1 carnitine as the top hub metabolite, P = 1.6 × 10−5), and “green” modules (pseudouridine as the top hub metabolite, P = 0.0007) than term controls (Table 3). Women with mPTD only had a significantly lower score of the “blue” module (P = 0.004), while women with sPTD had a significantly lower score of the “green” module (P = 0.003, Table 3), compared to term controls.

Fig. 2
figure 2

Correlations and p-values (in parentheses) of each coexpressed metabolite module with different preterm delivery (PTD) subgroups. The metabolite modules were built via the WGCNA analysis on 1330 women using the “unsigned” network

Full size image
Table 3 Associations of the eight metabolite modules with each preterm delivery (PTD) subgroup, compared to women with term delivery
Full size table

Stratified analyses were further performed by maternal ethnicity and maternal parity. As shown in Additional file 1: Table S3 and Table S4, the associations between the PTD/preeclampsia subgroups and each metabolite module were largely comparable between Black and non-Black mothers and between nulliparity and multiparity mothers.

Exploratory analyses of early and late sPTD

When women with sPTD were split into early (< 33 weeks of gestation, n = 67) and late sPTD (33–366/7 weeks of gestation, n = 152), and compared to term controls, we observed that the identified sPTD-associated differences in 16 individual metabolites were predominantly observed in women with early sPTD rather than in women with late sPTD (Additional file 1: Table S5). In addition, another 20 metabolites were significantly different in women with early sPTD, including 6 lipids (4 CEs and 2 LPCs) that were higher and 4 lipids (3 PEs and C36:2 phosphatidylinositol) that were lower in mothers with early sPTD (Additional file 1: Table S5 and Additional file 2: Fig. S1). Of note, the identified two LPCs (C22:5 and C22:6) were also significantly higher in women with mPTD and preeclampsia, and two PEs (C34:2 and C36:2) that were lower in women with early sPTD were also significantly lower in women with mPTD only. In the WGCNA analyses, women with early sPTD had significantly lower scores of the “yellow” (C16:0 CE as the top hub metabolite) and “green” (pseudouridine as the top hub metabolite) modules than term controls (Additional file 1: Table S6).

Discussion

To our knowledge, this study represents the largest to date to explore postpartum metabolomic profiles significantly associated with PTD/preeclampsia subgroups from a high-risk US multi-ethnic population. A total of 47, 4, and 16 metabolites were significantly different in women with mPTD and preeclampsia, in women with mPTD only, and in women with sPTD, respectively, compared with term controls. The broad-scale metabolite differences in women with mPTD and preeclampsia predominantly involve acylcarnitines and lipids. The neutral lipids, including DAGs, TAGs, and CEs, were significantly higher in women with mPTD and preeclampsia than in women with mPTD only, suggesting potential specific associations with preeclampsia. This finding is in line with clinically observed significant and independent adverse impact of preeclampsia-related PTD on maternal cardiometabolic health [1, 3]. We also observed significant lipid metabolite alteration in mothers with early sPTD, including those lipids (LPCs and PEs) that were associated with medically indicated PTD, indicating shared pathways among different PTD subtypes in affecting maternal metabolic health. While a comprehensive discussion on the differences of all the metabolites identified in this study is beyond the word limit of this report, below we just highlight the major findings and relevant literature.

Differences in acylcarnitine associated with mPTD and preeclampsia

Significantly higher levels of multiple acylcarnitines (metabolites related to cellular energy production) were mainly observed in women with mPTD and preeclampsia. Acylcarnitines are responsible for the transport of long-chain fatty acids from the cytoplasm into the mitochondria for β-oxidation. Fatty acids play an important role during pregnancy as metabolic fuel for the placenta [24]. Elevated acylcarnitines have been reported as potential biomarkers for preeclampsia [25,26,27]. Using the network analyses via WGCNA, acylcarnitines together with metabolites in purine metabolism (6,8-dihydroxypurine, xanthine, and xanthosine), which all were higher in women with mPTD and preeclampsia, were classified into the “black” module, suggesting coexpression of these metabolites.

Alterations in lipid metabolites among PTD/preeclampsia subgroups

Neutral lipids, including DAGs, TAGs, and CEs, constitute the most abundant group of lipids and serve as energy and carbon storage. These lipids, or the score for the “yellow” module, were significantly altered in women with mPTD and preeclampsia compared to women with mPTD only, suggesting preeclampsia-specific associations. Accumulation of DAGs has been reported to play a role in insulin resistance in cell [28]. DAGs are precursors to TAGs which play an important role in the storage of energy and fatty acids. In our study, the major TAGs that were higher in women with mPTD and preeclampsia included TAG C50:2, C52:2, and C52:3. These TAGs are of short- or medium-carbon chains which are likely hydrolyzed to release saturated (i.e., C16:0 palmitic acid or C18:0 stearic acid) and monounsaturated (i.e., C18:1 oleic acid) fatty acids. Oxidation of saturated and monounsaturated fatty acids provides more energy than oxidation of polyunsaturated fatty acids of the same length [29]. Accumulation of such short- or medium-carbon chain TAGs may indicate dysfunctional energy metabolism and may induce pro-inflammatory signaling [30], which has been proposed as a feature of multiple PTD-related diseases, including diabetes [22, 31] and CVD [14]. Another type of neutral lipid, CEs, was inversely correlated with DAGs and TAGs, which may partly explain the lower levels of CEs (C14:0, C16:0, C18:0, C18:2, and C18:3) that we observed in women with mPTD and preeclampsia. Consistently, lower CEs have been linked to unhealthy lifestyles including a lack of physical activity (including C16:0 and C18:2) [32] and an insulinemic dietary pattern (including C14:0 and C18:3) [33]. However, the relationships between CEs and CVD risk have been inconclusive, including both positive [14] and inverse associations [34, 35]. In comparison, we observed that multiple CEs (C18:0, C20:3, C20:4, and C22:4) were positively associated with early sPTD and that the score for the "yellow" module was inversely associated with early sPTD in this study. Future studies are needed to explore whether and how differences in levels of such neutral lipids may explain the adverse impact of preeclampsia-related PTD and early sPTD on CVD risk.

Altered metabolites of amino acids and nucleotides by different PTD/preeclampsia subgroups

Different levels of amino acids and nucleotides were observed in women with mPTD and preeclampsia and in women with sPTD (especially those with early sPTD), including both shared and unique metabolites. For example, anthranilic acid and bilirubin were lower across the different PTD subtypes. Anthranilic acid is the metabolite of tryptophan. Tryptophan level during pregnancy was lower in mothers with PTD [36]. Anthranilic acid derivatives have been reported to have anti-inflammatory function [37, 38]. Collectively, this evidence may help to explain the inverse relationship between PTD and anthranilic acid. Bilirubin is a bile pigment with potent antioxidant properties that inhibit lipid oxidation and atherosclerotic development [39, 40]. A number of studies have reported an inverse association between bilirubin levels and risk of hypertension, diabetes, and CVD [41, 42]. Taken together, these findings suggest the hypothesis that inflammation and elevated oxidative stress are potential mechanisms via which both preeclampsia-related mPTD and early sPTD may influence future cardiometabolic health.

The main strengths of our study included a large sample size, a well-validated metabolomic platform and detailed information on covariates, and more importantly, our ability to define PTD/preeclampsia subgroups, which are clinically and pathogenically important to gain new insight. Some limitations of our study should also be acknowledged. First, the BBC was particularly designed to study preterm delivery in an inner-city, predominantly minority population; therefore, caution is needed when generalizing our findings to other populations. Second, there is lack of a replication cohort in this study. We cannot exclude the possibility that some of these findings may not be replicated in another cohort. Third, maternal metabolites were measured in plasma collected within 1–3 days after delivery, which was at a time of significant flux with respect to molecules in the circulation. Future studies are needed to investigate whether the identified metabolite differences (especially those associated with mPTD and preeclampsia) are persistent across different time periods postpartum and whether they are associated with future CVD risk. Fourth, we could not distinguish early-onset versus late-onset preeclampsia in this study. Given the possibility that prolonged exposure to preeclampsia may be associated with an increased risk of CVD, further studies are needed to test whether early-onset preeclampsia and late-onset preeclampsia have different associations with maternal postpartum metabolomic profiles. Fifth, pregestational cardiovascular risk factors are important risk factors of gestational complications including preeclampsia [43,44,45] and they may also lead to metabolomic differences postpartum. To minimize potential confounding effects, we adjusted for pregestational overweight/obesity, hypertension, and diabetes in our analyses. Pregestational dyslipidemia is another potential risk factor of preeclampsia [44, 45]. However, information on pregestational dyslipidemia was not available in this study. Future studies are needed to assess potential confounding by pregestational dyslipidemia. Finally, we were only able to further stratify preterm delivery into early vs. late for the sPTD group, but not for other mPTD/preeclampsia subgroups due to limited sample size.

Conclusion

We demonstrated individual and clusters of maternal postpartum metabolites that were significantly and specifically associated with clinically observed subgroups: mPTD and preeclampsia, mPTD only, and sPTD, as well as some metabolites that were shared among the subgroups. Some of the identified lipid metabolites and/or related pathways, especially those associated with the mPTD and preeclampsia subgroup, have been reported as early predictors of future risk of CVD. Our study represents an early attempt to better understand metabolomic underpinnings of PTD/preeclampsia. Additional studies are warranted to confirm our findings and further investigate the link of these identified metabolites in mediating or modifying a woman’s future cardiometabolic risk. A better understanding of the pathways from PTD/preeclampsia, postpartum metabolomic profiles, to maternal future cardiometabolic risk may lead to new risk assessment and intervention strategies and help to identify at-risk women early in the disease process when intervention may delay or prevent the onset of cardiometabolic disease.

Availability of data and materials

The datasets used and/or analyzed for the current study are available from the corresponding author on reasonable request, and after the Institutional Review Board review and approval.

Abbreviations

BBC:

Boston Birth Cohort

BMI:

Body mass index

CE:

Cholesteryl ester

CV:

Coefficient of variance

CVD:

Cardiovascular disease

DAG:

Diacylglycerol

EMR:

Electronic medical record

LPC:

Lysophosphatidylcholine

LPE:

Lysophosphatidylethanolamine

mPTD:

Medically indicated PTD

PC:

Phosphatidylcholine

PE:

Phosphatidylethanolamine

PI:

Phosphatidylinositol

PTD:

Preterm delivery

sPTD:

Spontaneous PTD

TAG:

Triacylglycerol

WGCNA:

Weighted gene coexpression network analysis

References

  1. Auger N, Potter BJ, He S, Healy-Profitos J, Schnitzer ME, Paradis G. Maternal cardiovascular disease 3 decades after preterm birth: longitudinal cohort study of pregnancy vascular disorders. Hypertension. 2020;75(3):788–795.

  2. Bonamy AK, Parikh NI, Cnattingius S, Ludvigsson JF, Ingelsson E. Birth characteristics and subsequent risks of maternal cardiovascular disease: effects of gestational age and fetal growth. Circulation. 2011;124(25):2839–46.

    Article  Google Scholar 

  3. Tanz LJ, Stuart JJ, Williams PL, Rimm EB, Missmer SA, Rexrode KM, Mukamal KJ, Rich-Edwards JW. Preterm delivery and maternal cardiovascular disease in young and middle-aged adult women. Circulation. 2017;135(6):578–89.

    Article  Google Scholar 

  4. Wu P, Gulati M, Kwok CS, Wong CW, Narain A, O'Brien S, Chew-Graham CA, Verma G, Kadam UT, Mamas MA. Preterm delivery and future risk of maternal cardiovascular disease: a systematic review and meta-analysis. J Am Heart Assoc. 2018;7(2):e007809.

  5. Murphy SL, Xu J, Kochanek KD, Curtin SC, Arias E. Deaths: final data for 2015. Natl Vital Stat Rep. 2017;66(6):1–75.

    PubMed  Google Scholar 

  6. Minissian MB, Kilpatrick S, Eastwood JA, Robbins WA, Accortt EE, Wei J, Shufelt CL, Doering LV, Merz CNB. Association of spontaneous preterm delivery and future maternal cardiovascular disease. Circulation. 2018;137(8):865–71.

    Article  Google Scholar 

  7. Catov JM, Wu CS, Olsen J, Sutton-Tyrrell K, Li J, Nohr EA. Early or recurrent preterm birth and maternal cardiovascular disease risk. Ann Epidemiol. 2010;20(8):604–9.

    Article  Google Scholar 

  8. Zheng Y, Yu B, Alexander D, Steffen LM, Nettleton JA, Boerwinkle E. Metabolomic patterns and alcohol consumption in African Americans in the Atherosclerosis Risk in Communities Study. Am J Clin Nutr. 2014;99(6):1470–8.

    Article  CAS  Google Scholar 

  9. Magnusson M, Lewis GD, Ericson U, Orho-Melander M, Hedblad B, Engstrom G, Ostling G, Clish C, Wang TJ, Gerszten RE, et al. A diabetes-predictive amino acid score and future cardiovascular disease. Eur Heart J. 2013;34(26):1982–9.

    Article  CAS  Google Scholar 

  10. Cheng S, Rhee EP, Larson MG, Lewis GD, McCabe EL, Shen D, Palma MJ, Roberts LD, Dejam A, Souza AL, et al. Metabolite profiling identifies pathways associated with metabolic risk in humans. Circulation. 2012;125(18):2222–31.

    Article  CAS  Google Scholar 

  11. Shah SH, Sun JL, Stevens RD, Bain JR, Muehlbauer MJ, Pieper KS, Haynes C, Hauser ER, Kraus WE, Granger CB, et al. Baseline metabolomic profiles predict cardiovascular events in patients at risk for coronary artery disease. Am Heart J. 2012;163(5):844–50 e841.

    Article  CAS  Google Scholar 

  12. Kalim S, Clish CB, Wenger J, Elmariah S, Yeh RW, Deferio JJ, Pierce K, Deik A, Gerszten RE, Thadhani R, et al. A plasma long-chain acylcarnitine predicts cardiovascular mortality in incident dialysis patients. J Am Heart Assoc. 2013;2(6):e000542.

    Article  Google Scholar 

  13. Rizza S, Copetti M, Rossi C, Cianfarani MA, Zucchelli M, Luzi A, Pecchioli C, Porzio O, Di Cola G, Urbani A, et al. Metabolomics signature improves the prediction of cardiovascular events in elderly subjects. Atherosclerosis. 2014;232(2):260–4.

    Article  CAS  Google Scholar 

  14. Stegemann C, Pechlaner R, Willeit P, Langley SR, Mangino M, Mayr U, Menni C, Moayyeri A, Santer P, Rungger G, et al. Lipidomics profiling and risk of cardiovascular disease in the prospective population-based Bruneck study. Circulation. 2014;129(18):1821–31.

    Article  CAS  Google Scholar 

  15. Alshehry ZH, Mundra PA, Barlow CK, Mellett NA, Wong G, McConville MJ, Simes J, Tonkin AM, Sullivan DR, Barnes EH, et al. Plasma lipidomic profiles improve on traditional risk factors for the prediction of cardiovascular events in type 2 diabetes mellitus. Circulation. 2016;134(21):1637–50.

    Article  CAS  Google Scholar 

  16. Ganna A, Salihovic S, Sundstrom J, Broeckling CD, Hedman AK, Magnusson PK, Pedersen NL, Larsson A, Siegbahn A, Zilmer M, et al. Large-scale metabolomic profiling identifies novel biomarkers for incident coronary heart disease. PLoS Genet. 2014;10(12):e1004801.

    Article  Google Scholar 

  17. Cavus E, Karakas M, Ojeda FM, Kontto J, Veronesi G, Ferrario MM, Linneberg A, Jorgensen T, Meisinger C, Thorand B, et al. Association of circulating metabolites with risk of coronary heart disease in a European population: results from the Biomarkers for Cardiovascular Risk Assessment in Europe (BiomarCaRE) consortium. JAMA Cardiol. 2019;4(12):1270–9.

    Article  Google Scholar 

  18. Wang G, Divall S, Radovick S, Paige D, Ning Y, Chen Z, Ji Y, Hong X, Walker SO, Caruso D, et al. Preterm birth and random plasma insulin levels at birth and in early childhood. JAMA. 2014;311(6):587–96.

  19. Wang X, Zuckerman B, Pearson C, Kaufman G, Chen C, Wang G, Niu T, Wise PH, Bauchner H, Xu X. Maternal cigarette smoking, metabolic gene polymorphism, and infant birth weight. JAMA. 2002;287(2):195–202.

  20. ACOG Practice Bulletin No. 202. Gestational hypertension and preeclampsia. Obstet Gynecol 2019, 133(1):e1-e25.

  21. Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, Lewis GD, Fox CS, Jacques PF, Fernandez C, et al. Metabolite profiles and the risk of developing diabetes. Nat Med. 2011;17(4):448–53.

    Article  Google Scholar 

  22. Rhee EP, Cheng S, Larson MG, Walford GA, Lewis GD, McCabe E, Yang E, Farrell L, Fox CS, O'Donnell CJ, et al. Lipid profiling identifies a triacylglycerol signature of insulin resistance and improves diabetes prediction in humans. J Clin Invest. 2011;121(4):1402–11.

    Article  CAS  Google Scholar 

  23. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

    Article  Google Scholar 

  24. Rakheja D, Bennett MJ, Foster BM, Domiati-Saad R, Rogers BB. Evidence for fatty acid oxidation in human placenta, and the relationship of fatty acid oxidation enzyme activities with gestational age. Placenta. 2002;23(5):447–50.

    Article  CAS  Google Scholar 

  25. Kenny LC, Broadhurst DI, Dunn W, Brown M, North RA, McCowan L, Roberts C, Cooper GJ, Kell DB, Baker PN, et al. Robust early pregnancy prediction of later preeclampsia using metabolomic biomarkers. Hypertension. 2010;56(4):741–9.

    Article  CAS  Google Scholar 

  26. Koster MP, Vreeken RJ, Harms AC, Dane AD, Kuc S, Schielen PC, Hankemeier T, Berger R, Visser GH, Pennings JL. First-trimester serum acylcarnitine levels to predict preeclampsia: a metabolomics approach. Dis Markers. 2015;2015:857108.

    Article  Google Scholar 

  27. Thiele IG, Niezen-Koning KE, van Gennip AH, Aarnoudse JG. Increased plasma carnitine concentrations in preeclampsia. Obstet Gynecol. 2004;103(5 Pt 1):876–80.

    Article  CAS  Google Scholar 

  28. Erion DM, Shulman GI. Diacylglycerol-mediated insulin resistance. Nat Med. 2010;16(4):400–2.

    Article  CAS  Google Scholar 

  29. DeLany JP, Windhauser MM, Champagne CM, Bray GA. Differential oxidation of individual dietary fatty acids in humans. Am J Clin Nutr. 2000;72(4):905–11.

    Article  CAS  Google Scholar 

  30. Milanski M, Degasperi G, Coope A, Morari J, Denis R, Cintra DE, Tsukumo DM, Anhe G, Amaral ME, Takahashi HK, et al. Saturated fatty acids produce an inflammatory response predominantly through the activation of TLR4 signaling in hypothalamus: implications for the pathogenesis of obesity. J Neurosci. 2009;29(2):359–70.

    Article  CAS  Google Scholar 

  31. Razquin C, Toledo E, Clish CB, Ruiz-Canela M, Dennis C, Corella D, Papandreou C, Ros E, Estruch R, Guasch-Ferre M, et al. Plasma lipidomic profiling and risk of type 2 diabetes in the PREDIMED trial. Diabetes Care. 2018;41(12):2617–24.

    Article  CAS  Google Scholar 

  32. Ding M, Zeleznik OA, Guasch-Ferre M, Hu J, Lasky-Su J, Lee IM, Jackson RD, Shadyab AH, LaMonte MJ, Clish C, et al. Metabolome-wide association study of the relationship between habitual physical activity and plasma metabolite levels. Am J Epidemiol. 2019;188(11):1932–43.

    Article  Google Scholar 

  33. Tabung FK, Balasubramanian R, Liang L, Clinton SK, Cespedes Feliciano EM, Manson JE, Van Horn L, Wactawski-Wende J, Clish CB, Giovannucci EL, et al. Identifying metabolomic profiles of insulinemic dietary patterns. Metabolites. 2019;9(6):120.

  34. Razquin C, Liang L, Toledo E, Clish CB, Ruiz-Canela M, Zheng Y, Wang DD, Corella D, Castaner O, Ros E, et al. Plasma lipidome patterns associated with cardiovascular risk in the PREDIMED trial: a case-cohort study. Int J Cardiol. 2018;253:126–32.

    Article  Google Scholar 

  35. Toledo E, Wang DD, Ruiz-Canela M, Clish CB, Razquin C, Zheng Y, Guasch-Ferre M, Hruby A, Corella D, Gomez-Gracia E, et al. Plasma lipidomic profiles and cardiovascular events in a randomized intervention trial with the Mediterranean diet. Am J Clin Nutr. 2017;106(4):973–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Virgiliou C, Gika HG, Witting M, Bletsou AA, Athanasiadis A, Zafrakas M, Thomaidis NS, Raikos N, Makrydimas G, Theodoridis GA. Amniotic fluid and maternal serum metabolic signatures in the second trimester associated with preterm delivery. J Proteome Res. 2017;16(2):898–910.

    Article  CAS  Google Scholar 

  37. Sharma S, Srivastava VK, Kumar A. Newer N-substituted anthranilic acid derivatives as potent anti-inflammatory agents. Eur J Med Chem. 2002;37(8):689–97.

    Article  CAS  Google Scholar 

  38. Duffy JC, Dearden JC, Rostron C. A QSAR study of anti-inflammatory N-arylanthranilic acids. J Pharm Pharmacol. 1996;48(9):883–6.

    Article  CAS  Google Scholar 

  39. Stocker R, Yamamoto Y, McDonagh AF, Glazer AN, Ames BN. Bilirubin is an antioxidant of possible physiological importance. Science. 1987;235(4792):1043–6.

    Article  CAS  Google Scholar 

  40. Stocker R, McDonagh AF, Glazer AN, Ames BN. Antioxidant activities of bile pigments: biliverdin and bilirubin. Methods Enzymol. 1990;186:301–9.

    Article  CAS  Google Scholar 

  41. Abbasi A, Deetman PE, Corpeleijn E, Gansevoort RT, Gans RO, Hillege HL, van der Harst P, Stolk RP, Navis G, Alizadeh BZ, et al. Bilirubin as a potential causal factor in type 2 diabetes risk: a Mendelian randomization study. Diabetes. 2015;64(4):1459–69.

    Article  CAS  Google Scholar 

  42. Wang L, Bautista LE. Serum bilirubin and the risk of hypertension. Int J Epidemiol. 2015;44(1):142–52.

    Article  Google Scholar 

  43. Magnussen EB, Vatten LJ, Lund-Nilsen TI, Salvesen KA, Davey Smith G, Romundstad PR. Prepregnancy cardiovascular risk factors as predictors of pre-eclampsia: population based cohort study. BMJ. 2007;335(7627):978.

    Article  CAS  Google Scholar 

  44. Harville EW, Viikari JS, Raitakari OT. Preconception cardiovascular risk factors and pregnancy outcome. Epidemiology. 2011;22(5):724–30.

    Article  Google Scholar 

  45. Egeland GM, Klungsoyr K, Oyen N, Tell GS, Naess O, Skjaerven R. Preconception cardiovascular risk factor differences between gestational hypertension and preeclampsia: cohort Norway study. Hypertension. 2016;67(6):1173–80.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors wish to thank the study participants in the BBC, the nursing staff at labor and delivery of the Boston Medical Center, and the field team for their contributions to the Boston Birth Cohort. Linda Rosen in the Clinical Data Warehouse assisted in obtaining relevant clinical data from electronic medical records. The Clinical Data Warehouse service is supported by Boston University’s Clinical and Translational Institute and the National Institutes of Health Clinical and Translational Science awards (grant U54-TR001012). We thank Tami Bartell for her scientific editing of the manuscript text. We thank Dr. Meir Stampfer for his careful review and very helpful edits and comments.

Funding

The Boston Birth Cohort (the parent study, PI, Wang X) was supported in part by the March of Dimes PERI grants (20-FY02-56, #21-FY07-605), the Health Resources and Services Administration (HRSA) of the U.S. Department of Health and Human Services (HHS) under grant number R40MC27443, Autism Field-initiated Innovative Research Studies Program; and grant number UJ2MC31074, Autism Single Investigator Innovation Program; and the National Institutes of Health (NIH) grants (R21ES011666, R01HD041702, R21HD066471, U01AI090727, R21AI079872, R01HD086013, R01HD098232, R21AI154233 and R01ES031272). Dr. Hong is also supported by NICHD (PI: Hong X, R03HD096136) and Hopkins Population Center (NICHD R24HD042854). This information or content and conclusions are those of the author and should not be construed as the official position or policy of, nor should any endorsements be inferred by HRSA, HHS, or the US Government.

Author information

Author notes

  1. Xiumei Hong and Boyang Zhang contributed equally to this work.

Authors and Affiliations

  1. Center on the Early Life Origins of Disease, Department of Population, Family and Reproductive Health, Johns Hopkins University Bloomberg School of Public Health, 615 N. Wolfe St, Baltimore, MD, 21205-2179, USA

    Xiumei Hong, Yuelong Ji, Guoying Wang & Xiaobin Wang

  2. Department of Biostatistics, Johns Hopkins University Bloomberg School of Public Health, 615 N. Wolfe St, Baltimore, MD, 21205, USA

    Boyang Zhang, Yan Zhang & Hongkai Ji

  3. Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA

    Liming Liang & Frank B. Hu

  4. Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA

    Liming Liang

  5. Broad Institute of MIT and Harvard University, Cambridge, MA, USA

    Clary B. Clish

  6. Integrated Research Center for Fetal Medicine, Department of Gynecology and Obstetrics, Johns Hopkins University School of Medicine, Baltimore, MD, 21287, USA

    Irina Burd

  7. Department of Pediatrics, Boston University School of Medicine and Boston Medical Center, Boston, MA, 02118, USA

    Colleen Pearson & Barry Zuckerman

  8. Department of Nutrition, Harvard T.H. Chan School of Public Health, Boston, MA, USA

    Frank B. Hu

  9. Channing Division of Network Medicine, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA, USA

    Frank B. Hu

  10. Division of General Pediatrics & Adolescent Medicine, Department of Pediatrics, Johns Hopkins University School of Medicine, Baltimore, MD, 21205, USA

    Xiaobin Wang

Authors
  1. Xiumei Hong
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Boyang Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Liming Liang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yuelong Ji
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Guoying Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Hongkai Ji
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Clary B. Clish
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Irina Burd
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Colleen Pearson
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Barry Zuckerman
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Frank B. Hu
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Xiaobin Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

XH, FH, and XW were responsible for the study concept and design. XH, LL, and XW were responsible for drafting the manuscript. CC supervised laboratory metabolomic assays and quality control. XH, BYZ, YJ, and LL assumed major responsibilities for statistical analyses and data presentations. GW and YJ prepared blood samples for the lab assay. XW was responsible for overseeing the acquisition of the epidemiological and clinical data as well as biospecimens. All the authors were responsible for the critical review and revision of the manuscript and contributed to data interpretations. XH and XW have full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiumei Hong.

Ethics declarations

Ethics approval and consent to participate

The study has been approved by the Institutional Review Boards of Boston Medical Center and Johns Hopkins Bloomberg School of Public Health. Written informed consent was obtained from each participant.

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.

Other maternal metabolites that were significantly different in women of each preterm delivery (PTD) subgroup, compared to women of term delivery. Table S2. Individual metabolites included in each of the 8 metabolite modules. Table S3. Associations of metabolite modules with different preterm delivery (PTD) subgroups, stratified by maternal parity. Table S4. Associations of metabolite modules with different preterm delivery (PTD) subgroups, stratified by maternal race / ethnicity. Table S5. Maternal metabolites that were significantly different in women with late sPTD and in women with early sPTD, separately, compared to women with term delivery. Table S6. Associations of metabolite modules with early and late sPTD, separately, compared to women with term delivery.

Additional file 2: Figure S1.

Manhattan plot for the metabolomic differences in women with early and late spontaneous PTD, separately, compared to women with term delivery.

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

Hong, X., Zhang, B., Liang, L. et al. Postpartum plasma metabolomic profile among women with preeclampsia and preterm delivery: implications for long-term health. BMC Med 18, 277 (2020). https://doi.org/10.1186/s12916-020-01741-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-020-01741-4

Keywords

BMC Medicine

ISSN: 1741-7015

Contact us