Skip to main content

Metabolomic profiling reveals early biomarkers of gestational diabetes mellitus and associated hepatic steatosis

Abstract

Background

This study aims to identify early metabolomic biomarkers of gestational diabetes mellitus (GDM) and evaluate their association with hepatic steatosis.

Methods

We compared maternal serum metabolomic profiles between women who developed GDM (n = 118) and matched controls (n = 118) during the first (10–14 gestational weeks) and second (24–28 gestational weeks) trimesters using ultra-performance liquid chromatography coupled with mass spectrometry. Mediation analysis was performed to evaluate the mediating role of metabolic dysfunction-associated steatotic liver disease (MASLD) in the relationship between metabolites and subsequent development of GDM. A refined prediction model was developed to predict GDM using established clinical factors and selected metabolites.

Results

Significant alterations in circulating metabolites, including amino acids, bile acids, and phospholipids, were observed in the GDM group compared to controls during early pregnancy. Mediation analysis revealed that several metabolites, including glycocholic acid (proportion mediated (PM) = 31.9%), butanoyl carnitine (PM = 25.7%), and uric acid (PM = 22.4%), had significant indirect effects on GDM incidence mediated by hepatic steatosis. The refined prediction model composed of clinical factors and selected metabolites in the first trimester demonstrated higher performance in predicting GDM development than the established prediction model composed solely of clinical factors (AUC, 0.85 vs. 0.63, p < 0.001).

Conclusions

Women who developed GDM exhibited altered metabolomic profiles from early pregnancy, which showed a significant correlation with GDM, with MASLD as a mediator. Selected metabolomic biomarkers may serve as predictive markers and potential targets for early risk assessment and intervention in GDM.

Graphical abstract

Research insights

What is currently known about this topic?

Gestational diabetes mellitus (GDM) is a common pregnancy complication with significant health risks. Early identification of women at high risk for GDM is crucial for timely intervention and improved outcomes.

What is the key research question?

What alterations in circulating metabolites during early pregnancy are associated with subsequent GDM development?

Does metabolic dysfunction-associated steatotic liver disease (MASLD) mediate the association between specific metabolites and GDM risk?

What is new?

Significant alterations in bile acids, amino acids, phosphatidylethanolamines, and phosphatidylinositols were observed in early pregnancy sera of women who later developed GDM. MASLD significantly mediated the effects of several metabolites on GDM risk, with mediation proportions ranging from 9.7 to 31.9%. A refined prediction model composed of clinical factors and metabolites significantly improved the performance in predicting GDM development.

How might this study influence clinical practice?

These results provide new insights into early metabolic alterations associated with GDM development and highlight the potential mediating role of MASLD. This comprehensive metabolomic approach may contribute to the development of improved risk prediction models and targeted interventions for GDM prevention.

Background

Gestational diabetes mellitus (GDM) is a common pregnancy complication, affecting approximately 14% of pregnancies worldwide [1]. GDM is associated with adverse maternal and fetal outcomes during pregnancy and increased long-term health risks for both mothers and children [2,3,4,5,6,7,8]. Early identification of women at high risk for GDM is crucial for timely intervention and improved outcomes [9,10,11].

Current screening strategies for GDM, typically conducted between 24 and 28 weeks of gestation, may miss opportunities for early intervention [12, 13]. Recent studies have explored the potential of metabolomics to identify early biomarkers for GDM risk. Among metabolites in early pregnancy, selected metabolites have been used to develop prediction models for GDM [14, 15]. Additionally, emerging evidence suggests a link between metabolic dysfunction-associated steatotic liver disease (MASLD), formerly known as non-alcoholic fatty liver disease, and the development of GDM [12, 16,17,18].

Metabolomics offers a comprehensive approach to understanding the metabolic alterations that precede GDM diagnosis. Previous studies have demonstrated associations between various metabolites, including amino acids, lipids, and microbiota-related compounds, and subsequent GDM development [15, 19]. However, many of these studies have been limited by their focus on specific metabolite classes or use of single analytical techniques.

In this study, we aimed to conduct a comprehensive metabolomic analysis of serum samples from early pregnancy to identify potential biomarkers of GDM risk. We utilized a multi-platform approach, incorporating bile acid (BA) analysis, polar metabolite profiling, and lipidomics. Furthermore, we investigated the potential mediating role of MASLD in the relationship between circulating metabolites and GDM development. Finally, a refined prediction model was developed to predict GDM using clinical factors and selected predictive metabolites.

Methods

Study subjects

This study used data from the ongoing “Fatty Liver in Pregnancy” registry (NCT02276144), a prospective cohort study initiated in 2014 at Incheon Seoul Women’s Hospital and Seoul Metropolitan Government Seoul National University Boramae Medical Center in Korea. The study was approved by the Institutional Review Board of Seoul Metropolitan Government Seoul National University Boramae Medical Center (no. 26-2014-18) and the Public Institutional Review Board designated by the Ministry of Health and Welfare (no. P01-201404-BM-03). All participants provided informed consent.

Between 2014 and 2020, 1,936 women were enrolled in the cohort. Singleton pregnant women who visited the antenatal clinic in the first trimester were invited to participate (Fig. 1). Enrolled women were assessed for hepatic steatosis using liver ultrasound at 10–14 weeks of gestation and were followed until delivery. GDM was diagnosed using the two-step approach, following the recommendations of the American College of Obstetricians and Gynecologists [20]. For the 50 g oral glucose challenge test (GCT), plasma glucose was measured one hour after the glucose load, with ≥ 140 mg/dL considered positive. Women with a positive GCT underwent a 100 g oral GCT, and GDM was diagnosed if at least two values exceeded the thresholds (> 95 mg/dL fasting, > 180 mg/dL at 1 h, > 155 mg/dL at 2 h, and > 140 mg/dL at 3 h).

Fig. 1
figure 1

Schematic overviews of the study

After applying exclusion criteria (viral hepatitis, previable birth, consent withdrawal, and incomplete follow-up), 1756 subjects remained eligible. From this population, we selected 118 women who developed GDM and 118 matched controls based on maternal age and parity.

Clinical parameters and assessments

Participants were enrolled before 14 weeks of gestation. Baseline characteristics, including clinical and demographic information, were collected at enrollment via self-report questionnaire. At 10–14 weeks of gestation, participants underwent ultrasound examination to assess for the presence of hepatic steatosis, and fasting blood samples were collected. The presence of hepatic steatosis was evaluated based on the detection of bright echogenic patterns of the liver by ultrasonography. A semi-quantitative grading system (grade 0–3) was used to classify the severity of hepatic steatosis [21, 22].

The second blood sample was obtained at 24–28 weeks of gestation during the 50 g glucose challenge test. In this case-control study, we performed lipidomic analysis, polar metabolite analysis, and BA analysis on serum samples collected at both 10–14 weeks and 24–28 weeks of gestation. Results were analyzed in relation to subsequent GDM development.

Serum BA analysis by UPLC/TQ MS

For serum BA analysis, 20 µL of serum samples were dispensed into a 96-well plate. All the samples were pooled in equal volumes to generate quality control (QC) samples, which were analyzed prior to sample acquisition and subsequently after every 11 samples. For protein precipitation, 80 µL of methanol was added to each well and vortexed for 10 min. After vortexing, the mixtures were incubated at 4 °C for 2 h, followed by centrifugation at 3300 × g for 20 min at 10 °C. 50 µL of supernatant were transferred into a new 96-well plate, and 50 µL of methanol: water (20:80, v/v) with internal standards was added to each well. BA Standard Mix 1– Unconjugated and BA Standard Mix 2– Conjugated (Cambridge Isotope Laboratories, Inc., USA) were used as internal standards. Serum extracts were analyzed by an ACQUITY UPLC system (Waters, Milford, MA, USA) coupled with a XEVO TQ-XS Triple Quadrupole MS system (Waters, Milford, MA, USA) in multiple reaction monitoring (MRM) mode. Chromatographic separation of BAs in the serum was performed using an ACQUITY UPLC BEH C18 column (2.1 × 50 mm, 1.7 μm; Waters) at 60 °C and a flow rate of 0.5 mL/min. The binary mobile phases comprised 5 mM ammonium acetate in water (solvent A) and methanol: acetonitrile (50:50 v/v, solvent B). The gradient was as follows: 35–40% B from 0 to 3.5 min, 40–45% B from 3.5 to 4.0 min, 45–50% B from 4.0 to 4.8 min, 50–55% B from 4.8 to 5.5 min, 55–80% B from 5.5 to 6.2 min, 80–85% B from 6.2 to 6.4 min, 85% B from 6.4 to 7.0 min, 85–35% B from 7.0 to 7.01 min, and 35% B from 7.01 to 9.0 min to equilibrate. Retention times and MRM transitions of each BA are shown in Supplementary Table S1. Measured BA concentrations below the limit of quantification (LOQ) were imputed as the LOQ of the respective BA. In the case of BA concentrations in which 30% or more are below the LOQ, those BAs, such as LCA, GLCA, TLCA and TUDCA, were excluded from further analysis.

Serum semi-target analysis by UPLC/TQ MS

For serum polar metabolite analysis, 20 µL of serum samples were dispensed into a 96-well plate. All the samples were pooled in equal volumes to generate quality control (QC) samples, which were analyzed prior to sample acquisition and after every 11 samples. For protein precipitation, 80 µL of methanol was added to each well and vortexed for 10 min. After vortexing, the mixtures were incubated at 4 °C for 2 h, followed by centrifugation at 3300 × g for 20 min at 10 °C. Then, 50 µL of supernatant was transferred into a new 96-well plate and 50 µL of methanol: water (20:80 v/v) with internal standards was added to each well. A mixed solution of betaine D11, glutamate 13C, leucine 13C6, octanoyl carnitine 13C, phenylalanine 13C6, succinate 13C2, taurine 13C, and uridine 13C was used as the internal standards. Serum semi-target analysis was conducted by 1290 Infinity LC (Agilent Technologies, USA) coupled to a 6490 Triple Quadrupole MS system (Agilent Technologies, USA) in the dynamic multiple reaction monitoring (dMRM) mode. Chromatographic separation of polar metabolites was performed using a Scherzo SM-C18 column (2.0 × 100 mm, 3 μm; Imtakt) at 25 °C and a flow rate of 0.2 mL/min. The binary mobile phases comprised 0.1% (v/v) formic acid in water (solvent A) and 0.1% (v/v) formic acid in methanol (solvent B). The gradient was as follows: 5% B from 0 to 3.0 min, 5–100% B from 3.0 to 14.0 min, 100% B from 14.0 to 16.0 min, 100–5% B from 16.0 to 17.0 min, and 5% B from 17.0 to 20.0 min to equilibrate.

Information on selected metabolites and MRM transitions was described previously [23]. The MassHunter Workstation (Ver. B.06.00; Agilent Technologies) was used for data acquisition and analysis. After normalization by internal standards, three metabolites (5-methyltetrahydrofolate, adenosine, and riboflavin) with coefficients of variation > 30 in quality control samples were excluded. Then, total 62 metabolites were used for further statistical analysis.

Serum lipidomic profiling by UPLC/QTOF MS

For serum lipidomic analysis, 20 µL of serum samples were dispensed into a 96-well plate. All the samples were pooled in equal volumes to generate quality control (QC) samples, which were analyzed prior to sample acquisition and after every 11 samples. For protein precipitation, 30 µL of water and 200 µL of isopropanol with internal standards were added to each well and vortexed for 10 min. After vortexing, the mixtures were incubated at 4 °C for 2 h followed by centrifugation at 3300 × g for 20 min at 10 °C. Then, 150 µL of supernatant were transferred into a new 96-well plate, and 20 µL of samples were transferred into new two conical plates (one each for positive and negative mode). 180 µL of isopropanol: water (80:20 v/v) was added to all wells. SPLASH Lipidomix (Avanti Polar Lipids, USA) was used as the internal standard. Lipidomic profiling analysis was conducted using an ACQUITY UPLC system (Waters, Milford, MA, USA) coupled with a XEVO G2-XS QTOF mass spectrometer equipped with an electrospray ionization (ESI) source (Waters, Milford, MA, USA). Chromatographic separation of lipids in the serum was performed using an ACQUITY UPLC CSH C18 column (2.1 × 100 mm, 1.7 μm; Waters) at 55 °C and a flow rate of 0.4 mL/min. The binary mobile phases comprised 10 mM ammonium formate and 0.1% (v/v) formic acid in acetonitrile: water (60:40 v/v, solvent A) and 0.1% (v/v) formic acid in isopropanol: acetonitrile (90:10 v/v, solvent B) for positive mode and 10 mM ammonium acetate in water: acetonitrile (60:40 v/v, solvent A) and isopropanol: acetonitrile (90:10 v/v, solvent B) for negative mode. The gradient was as follows: 40–43% B from 0 to 2 min, 43–50% B from 2 to 2.1 min, 50–54% B from 2.1 to 12 min, 54–70% B from 12 to 12.1 min, 70–99% B from 12.1 to 18 min, 99–40% B from 18 to 18.1 min, and 40% B from 18.1 to 20 min to equilibrate. The mass spectrometer was operated in positive and negative ion mode, and the mass range was set at 80–1500 m/z for analysis of the lipid extracts.

MS data produced from the sum of the total intensity of all detected peaks in a UPLC/QTOF-MS was processed with Progenesis QC (Waters) and normalized to QC samples. Features with coefficients of variation below 30% in pooled QC samples were selected to determine the precision of detection. Lipid metabolites were identified using several authentic standard compounds and online databases, such as the Human Metabolome Database (HMDB) (www.hmdb.ca), METLIN (https://metlin.scripps.edu) and LIPID MAPS (www.lipidmaps.org).

Statistical analysis

Statistical analyses were performed using SPSS Statistics version 28.0 (IBM Corporation, Armonk, NY, USA), MedCalc version 23.0.5 (MedCalc Software, Ostend, Belgium), and R version 3.6.0 (http://www.r-project.org). To assess the normality of continuous variables, the Kolmogorov-Smirnov test was conducted. Continuous variables following a normal distribution are expressed as means ± standard deviations and compared using independent t-tests, while those not following a normal distribution are presented as medians (interquartile ranges) and compared using Mann-Whitney U-tests. Categorical variables are reported as numbers (percentages) and compared using chi-square tests. For metabolite data, all values are expressed as medians (interquartile ranges) and compared using ranked analysis of covariance (ANCOVA). P-values were obtained from ranked analysis of covariance (ANCOVA), with adjustments for maternal age, pre-pregnancy BMI, nulliparity, and family history of diabetes. For multiple testing, p values were adjusted with the Benjamini–Hochberg method. Spearman’s correlation was used to assess relationships between variables. Multiple logistic regression was performed to evaluate the association between circulating metabolites and GDM, adjusting for maternal age, pre-pregnancy body mass index (BMI), nulliparity, and family history of diabetes. In addition, causal mediation analysis was performed to evaluate the direct and indirect effects of circulating metabolites on GDM, with MASLD as a mediator. The R-package “mediation” [24] was used, adjusting for confounders such as maternal age, pre-pregnancy BMI, nulliparity, and family history of diabetes in both outcome and mediator models. A p-value < 0.05 was considered statistically significant for both direct and indirect effects. Receiver operating characteristic (ROC) curves were used to evaluate the predictive performance of clinical factors and selected metabolites in the first trimester. To establish the prediction models for GDM development, the risk score using clinical factors was calculated based on baseline clinical factors in the first trimester described previously [25]. Baseline clinical factors included overweight status, family history of DM, age, history of pregnancy complications, history of polycystic ovary syndrome (PCOS), issues with insulin or blood glucose regulation, and a history of high blood pressure, high cholesterol, and/or heart disease. The areas under the ROC curves (AUCs) were compared using the DeLong test [26]. Flow chart of the statistical analysis was presented in Supplementary Fig. 1.

Results

Baseline characteristics of study participants

Of the 1,936 pregnant women enrolled in the cohort, 236 women (118 women who developed GDM and 118 matched controls) participated in this study (Fig. 1). Table 1 presents the baseline characteristics of the study population. As expected, due to matching, maternal age and parity were similar between GDM cases and controls. Women who developed GDM were more likely to have a family history of type 2 diabetes mellitus and had higher pre-pregnancy BMI. At enrollment, GDM cases also demonstrated higher BMI and blood pressure compared to controls. Ultrasound examinations at 10–14 weeks of gestation revealed a higher frequency of steatotic liver disease in women who subsequently developed GDM (35.2% vs. 10.7%, p < 0.001). Laboratory tests showed that GDM cases had higher median values of aminotransferases, cholesterol, low-density lipoprotein cholesterol, triglycerides, and gamma-glutamyl transferase (GGT).

Table 1 Baseline clinical, biochemical, and metabolic features of study participants in the first trimester

Association between circulating metabolites in the first trimester and GDM

BA, polar semi-targeted, and lipidomic analyses were performed on a total of 236 sera collected from 236 study participants in the first trimester before GDM diagnosis using UPLC/TQ-MS and UPLC/QTOF-MS. The quantified 11 BAs and 62 polar metabolites are summarized in Supplementary Tables S2 and S3. Lipidomic profiles encompassed all lipid classes: free fatty acids (FFAs), ceramides (Cers), sphingomyelins (SMs), lysophosphatidylcholines (LPCs), lysophosphatidylethanolamines (LPEs), phosphatidylcholines (PCs), phosphatidylethanolamines (PEs), phosphatidylinositols (PIs), diacylglycerols (DAGs), and triacylglycerols (TAGs). All annotated lipid metabolites are summarized in Supplementary Table S4.

Comprehensive metabolomic analysis revealed significant alterations in diverse circulating metabolites in early pregnancy serum samples from women who later developed GDM (Fig. 2). After adjusting for potential confounders, the most distinctive changes were observed in BAs, amino acids, PEs, and PIs. BAs, particularly conjugated primary BAs such as glycocholate (GCA), glycochenodeoxycholate (GCDCA), taurocholate (TCA), and taurochenodeoxycholate (TCDCA), showed characteristic changes in the GDM group (Fig. 2A). Among polar metabolites, levels of several amino acids, including branched-chain amino acids (leucine, isoleucine, valine), aromatic amino acids (tyrosine), and others (proline, arginine, alanine, threonine), were elevated in the GDM group (Fig. 2B). Lipidomic analysis revealed significant alterations in FFAs, LysoPEs, PEs, PIs, and TAGs between the GDM and control groups (Fig. 2C). Notably, several PE species (PE 34:1, PE 36:1, PE 36:3, PE 36:4, PE 40:6) showed marked increases in the GDM group.

Fig. 2
figure 2

Alterations of circulating metabolites in serum samples from pregnant women during the first trimester. (A) bile acid (BA) analysis, (B) polar semi-targeted analysis, and (C) lipidomic analysis. Fold changes were calculated by dividing the median value of GDM group by the median value of non-GDM group. Adjusted p-values (padj) were calculated using the Benjamini–Hochberg method from p-values obtained from ranked analysis of covariance (ANCOVA) test, with adjustments for maternal age, pre-pregnancy BMI, nulliparity, and family history of diabetes; *padj<0.05, **padj <0.01, and ***padj <0.001

Association between circulating metabolites in the first trimester and GDM risk factors

Correlation analysis revealed significant associations between altered metabolites and established GDM risk factors (Fig. 3). Amino acids, particularly branched-chain amino acids (BCAAs) and tyrosine, showed strong correlations with waist circumference, triglycerides (TGs), and GGT. BA levels primarily correlated with TGs and fasting blood glucose, while PE and PI levels showed strong associations with TGs and total cholesterol. Consequently, metabolites exhibiting significant changes associated with GDM development commonly demonstrated strong correlations with TGs.

Fig. 3
figure 3

Correlations between circulating metabolite levels and laboratory results in serum samples from pregnant women during the first trimester. (A) bile acid (BA) analysis, (B) polar semi-targeted analysis, and (C) lipidomic analysis. Statistical significance in Spearman’s correlation heat maps is indicated with asterisks (*p < 0.05, **p < 0.01, and ***p < 0.001)

Association between circulating metabolites in the first trimester and GDM based on MASLD status

Multiple logistic regression analysis further elucidated the relationships between specific metabolites and the risk of GDM and MASLD (Fig. 4). Several conjugated BAs (GCA, GCDCA, TCA, TDCA, TCDCA) were associated with increased GDM risk, while CA, GCA, and TCA correlated with MASLD risk (Fig. 4A). Among polar metabolites, various amino acids (leucine, valine, tyrosine, hydroxyproline, arginine, alanine, kynurenine) were associated with increased risks of both GDM and MASLD (Fig. 4B). Interestingly, carnitine and some short-chain acylcarnitines showed a distinct relationship with MASLD risk. Additionally, certain lipid metabolites, particularly from FFA, PE, and PI classes, demonstrated stronger associations with GDM risk compared to MASLD risk (Fig. 4C). Mediation analysis revealed that MASLD significantly mediated the effects of several metabolites on GDM risk (Table 2). A total of 56 metabolites demonstrated a significant total effect on the relationship between metabolites and GDM (Supplementary Table S5). Conjugated primary BAs (GCA and TCA), acylcarnitine (butanoyl carnitine), purine metabolites (inosine and uric acid), and various amino acids showed significant indirect effects on GDM incidence mediated by MASLD, with mediation proportions ranging from 9.7 to 31.9%. Additionally, several lipids, such as PE 40:6, PI 38:3, and SM d40:1, showed significant mediation effects through MASLD. Pathway analysis from MetaboAnalyst (https://www.metaboanalyst.ca/) confirmed that these significant metabolites are linked to biological processes such as valine, leucine, and isoleucine biosynthesis, as well as primary BA metabolism (Supplementary Fig. S2).

Fig. 4
figure 4

Associations between circulating metabolites and GDM / MASLD during the first trimester. (A) bile acid (BA) analysis, (B) polar semi-targeted analysis, and (C) lipidomic analysis. Multivariable–adjusted odds ratio, 95% CI and p-value were calculated from multiple logistic regression with adjustments for maternal age, pre-pregnancy BMI, nulliparity and family history of diabetes; *p < 0.05, **p < 0.01, and ***p < 0.001 for GDM and #p < 0.05, ##p < 0.01, and ###p < 0.001 for MASLD

Table 2 Mediation analysis on the association between circulating metabolites in the first trimester and GDM, with MASLD as a mediator

Association between circulating metabolites in the second trimester and GDM

To assess metabolic changes after GDM diagnosis, we also conducted BA, polar semi-targeted, and lipidomic analyses on sera collected from identical 236 study subjects in the second trimester. In the second trimester, women who developed GDM had higher values of BMI, AST, and TGs, with lower values of HDL-C than those who did not (Supplementary Table S6). The association between circulating metabolites and GDM risk differed based on MASLD status and gestational age (Fig. 5). Many metabolites showed further increases in GDM cases with concurrent MASLD at both the first and second trimesters.

Fig. 5
figure 5

Alterations in circulating metabolites according to MASLD status in serum samples from pregnant women during the first and second trimesters. P-values were obtained from the ranked analysis of covariance (ANCOVA) test, with adjustments for maternal age, pre-pregnancy BMI, nulliparity and family history of diabetes; *p < 0.05, **p < 0.01, and ***p < 0.001 for GDM and #p < 0.05, ##p < 0.01, and ###p < 0.001 for MASLD

Prediction models for GDM development

A refined prediction model was developed incorporating both clinical factors and selected metabolites from the first trimester. Metabolites were selected based on their significant effects on GDM risk mediated by MASLD, including the following 16 metabolites: GCA, TCA, butanoyl carnitine, inosine, uric acid, leucine, valine, tyrosine, glycine, hydroxyl-proline, arginine, alanine, kynurenine, PE 40:6, PI 38:3, and SM d40:1. This model demonstrated significantly improved performance in predicting GDM development compared to the established prediction model based solely on clinical factors (AUC 0.85 vs. 0.63, p < 0.001, Fig. 6). The inclusion of metabolomic biomarkers significantly enhanced the model’s ability to identify women at high risk for GDM early in pregnancy.

Fig. 6
figure 6

Prediction model evaluation using receiver operating characteristic (ROC) curves

Discussion

Main findings of the current study

This comprehensive metabolomic study identified prospective alterations in various metabolite classes during early pregnancy that were associated with subsequent GDM development. We observed that specific BAs, amino acids, and lipids were linked to both GDM and MASLD risk. Furthermore, mediation analysis suggested that certain circulating metabolites may influence GDM development through MASLD-related pathways.

Metabolic profiling and the risk of GDM

Our findings aligned with previous studies demonstrating associations between amino acids, lipids, and purine metabolites in early pregnancy and GDM incidence [19, 27]. However, our multi-platform approach provided a more comprehensive metabolomic profile, revealing additional potential biomarkers and metabolic pathways involved in GDM pathogenesis.

The observed alterations in amino acid levels, particularly branched-chain amino acids (BCAAs), were consistent with their established role in insulin resistance and pancreatic β-cell dysfunction [28]. Elevated BCAA levels have been consistently associated with insulin resistance across various populations, including normoglycemic adults, obese individuals, and those with prediabetes or diabetes [28]. Our results suggested that increased BCAA levels in early pregnancy may indicate underlying insulin resistance and impaired β-cell function in women who later developed GDM.

Conversely, decreased levels of glutamine, glycine, and serine were also linked to an increased risk of GDM. While the literature has reported inconsistent findings regarding the relationship between amino acid levels and GDM, these discrepancies might be due to variations in genetic background, dietary intake, sample collection time, and environmental factors [19, 27].

Interestingly, we found that several amino acids, including BCAAs, exhibited significant mediation effects through MASLD. This finding aligns with previous research demonstrating altered amino acid profiles in MASLD patients, independent of obesity or insulin resistance [29]. The greater changes observed in GDM cases with concurrent MASLD suggested that pregnant women with MASLD may be at particularly high risk for GDM development.

Our study also identified novel associations between BA metabolism and GDM risk. While previous research on BAs in GDM has been limited and inconsistent [30, 31], our findings suggested that elevated levels of conjugated primary BAs, particularly GCA and TCA, in early pregnancy may contribute to disturbances in hepatic BA homeostasis and increase GDM risk. This is supported by recent evidence linking intrahepatic cholestasis of pregnancy, characterized by aberrant BA metabolism, to an increased incidence of GDM [32, 33]. Furthermore, our results indicated that chemically distinct BAs may exert varying effects on hyperglycemia across different species.

In the current study, elevated levels of uric acid, the end product of purine metabolism, were associated with increased risks of both GDM and MASLD, consistent with previous studies [34, 35]. Uric acid has been linked to insulin resistance, inflammation, and oxidative stress [36, 37], all of which may contribute to the pathogenesis of both GDM and MASLD.

Our findings regarding altered lipid profiles, particularly PEs and PIs, provided new insights into the lipid metabolism disturbances associated with GDM. Elevated levels of specific PE and PI species in early pregnancy may reflect underlying insulin resistance and altered lipid metabolism [38, 39]. These observations are consistent with previous studies reporting dysregulated lipid profiles in women who later developed GDM [40].

Altered metabolic profiling for Predicting GDM Risk in the context of hepatic steatosis

The development of a prediction model incorporating both clinical factors and selected metabolites represented a significant advancement in early GDM risk assessment. Our refined model demonstrated superior performance compared to established models based solely on clinical factors [25], with an AUC of 0.85 versus 0.63 (p < 0.001). This finding aligns with recent studies highlighting the potential of integrating metabolomic biomarkers with traditional risk factors to improve GDM prediction [41, 42]. Although external validation is necessary to confirm the model’s generalizability, identifying a comparable external cohort is particularly challenging due to the unique characteristics of our study population. Nevertheless, our results demonstrate that incorporating metabolic biomarkers enhances predictive performance. Future studies with independent cohorts will be essential to further validate and strengthen our findings.

MASLD, GDM, and metabolites

Previous studies have consistently reported an association between hepatic steatosis and GDM [12, 16,17,18]. However, the underlying mechanisms linking hepatic steatosis to GDM risk remain unclear. One possibility is that these two metabolic conditions share a common underlying dysfunction, such as insulin resistance [43, 44]. Additionally, proinflammatory cytokines released from intra-abdominal visceral tissue might be transferred to the liver, promoting hepatic fat accumulation and ultimately leading to hepatic steatosis [45, 46]. This process, along with the release of adipokines and increased lipolytic activity, may make visceral adipose tissue more metabolically detrimental, contributing to both hepatic steatosis and GDM [47, 48].

The significant mediation effects of MASLD on the relationship between several metabolites and GDM risk highlight the complex interplay between hepatic steatosis and glucose metabolism during pregnancy. This finding supports the emerging concept of MASLD as a potential driver of metabolic dysfunction in pregnancy [17, 18], rather than merely a consequence. We observed stronger associations between certain metabolites and GDM risk in women with concurrent MASLD, suggesting that this subgroup may benefit from more intensive monitoring and intervention strategies, such as lifestyle modifications.

Strengths and limitations of the current study

Our study has several strengths, including its prospective design, comprehensive metabolomic profiling approach, and inclusion of MASLD assessment. The use of multiple analytical platforms allowed for a more holistic view of the metabolic alterations preceding GDM development. Additionally, the mediation analysis provided novel insights into potential mechanistic pathways linking metabolic disturbances to GDM risk.

However, some limitations should be acknowledged. First, the relatively small sample size, particularly for MASLD cases, may limit the generalizability of our findings. Second, while we adjusted for several potential confounders, residual confounding cannot be ruled out. Third, the use of non-fasting samples during the oral glucose loading test in the second trimester may have influenced some metabolite levels, although this limitation was mitigated by consistent sampling conditions across all participants. Lastly, we opted for a 1:1 matching strategy to ensure a well-balanced comparison between cases and controls. However, a higher control-to-case ratio would have provided additional statistical power.

Future research should focus on external validation of our findings in larger, more diverse cohorts and explore longitudinal changes in metabolomic profiles throughout pregnancy. Additionally, mechanistic studies are needed to elucidate the biological pathways underlying the observed associations between specific metabolites, MASLD, and GDM development.

Conclusions

In conclusion, our comprehensive metabolomic analysis identified novel early pregnancy biomarkers associated with subsequent GDM development and highlighted the potential mediating role of MASLD. These findings enhance our understanding of the complex metabolic perturbations preceding GDM and may inform the development of improved risk prediction models and targeted interventions. The integration of metabolomic biomarkers with clinical risk factors shows promise for enhancing early GDM risk assessment, potentially enabling more timely and personalized preventive strategies.

Data availability

No datasets were generated or analysed during the current study.

Abbreviations

ADMA:

Asymmetric dimethylarginine

ALT:

Alanine aminotransferase

AST:

Aspartate aminotransferase

BA:

Bile acid

BMI:

Body mass index

CA:

Cholate

CDCA:

Chenodeoxycholate

Cer:

Ceramide

95% CI:

95% confidence intervals

DAG:

Diacylglycerol

DCA:

Deoxycholate

DMG:

Dimethylglycine

FC:

Fold change

FFA:

Free fatty acid

GABA:

γ-aminobutyrate

GCA:

Glycocholate

GCDCA:

Glycochenodeoxycholate

GDCA:

Glycodeoxycholate

GDM:

Gestational diabetes mellitus

GGT:

Gamma-glutamyl transferase

GLCA:

Glycolithocholate

GPC:

Sn-glycerol-3-phosphocholine

GUDCA:

Glycoursodeoxycholate

HDL-C:

High-density lipoprotein-cholesterol

LCA:

Litocholate

LDL-C:

Low-density lipoprotein-cholesterol

LysoPC:

Lysophosphatidylcholine

LysoPE:

Lysophosphatidylethanolamine

MASLD:

Metabolic dysfunction-associated steatotic liver disease

MRM:

Multiple reaction monitoring

OR:

Odd ratio

PC:

Phosphatidylcholine

PE:

Phosphatidylethanolamine

PI:

Phosphatidylinositol

SAH:

S-adenosylhomocysteine

SM:

Sphingomyelin

TAG:

Triacylglycerol

TCA:

Taurocholate

TCDCA:

Taurochenodeoxycholate

TDCA:

Taurodeoxycholate

TLCA:

Taurolithocholate

TGs:

Triglycerides

TMAO:

Trimethylamine-N-oxide

UDCA:

Ursodeoxycholate

References

  1. Sweeting A, Hannah W, Backman H, Catalano P, Feghali M, Herman WH, Hivert MF, Immanuel J, Meek C, Oppermann ML, et al. Epidemiology and management of gestational diabetes. Lancet. 2024;404(10448):175–92.

    Article  PubMed  Google Scholar 

  2. Nakshine VS, Jogdand SD. A comprehensive review of gestational diabetes mellitus: impacts on maternal health, fetal development, childhood outcomes, and long-term treatment strategies. Cureus. 2023;15(10):e47500.

    PubMed  PubMed Central  Google Scholar 

  3. Wicklow B, Retnakaran R. Gestational diabetes mellitus and its implications across the life span. Diabetes Metab J. 2023;47(3):333–44.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Sweeting A, Wong J, Murphy HR, Ross GP. A clinical update on gestational diabetes mellitus. Endocr Rev. 2022;43(5):763–93.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Kramer CK, Campbell S, Retnakaran R. Gestational diabetes and the risk of cardiovascular disease in women: a systematic review and meta-analysis. Diabetologia. 2019;62(6):905–14.

    Article  PubMed  Google Scholar 

  6. Lee SM, Shivakumar M, Park JW, Jung YM, Choe EK, Kwak SH, Oh S, Park JS, Jun JK, Kim D, et al. Long-term cardiovascular outcomes of gestational diabetes mellitus: a prospective UK Biobank study. Cardiovasc Diabetol. 2022;21(1):221.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Plows JF, Stanley JL, Baker PN, Reynolds CM, Vickers MH. The pathophysiology of gestational diabetes Mellitus. Int J Mol Sci 2018, 19(11).

  8. Hod M, Kapur A, Sacks DA, Hadar E, Agarwal M, Di Renzo GC, Cabero Roura L, McIntyre HD, Morris JL, Divakar H. The international federation of gynecology and obstetrics (FIGO) initiative on gestational diabetes mellitus: a pragmatic guide for diagnosis, management, and care. Int J Gynaecol Obstet. 2015;131(Suppl 3):S173–211.

    Article  PubMed  Google Scholar 

  9. Wang J, Huang P, Hou F, Hao D, Li W, Jin H. Predicting gestational diabetes mellitus risk at 11–13 weeks’ gestation: the role of extrachromosomal circular DNA. Cardiovasc Diabetol. 2024;23(1):289.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Frontino G, Succurro E, Corcoy R, Scialabba F, Poloniato A, Scavini M. Editorial: current and future trends in gestational diabetes diagnosis, care and neonatal outcomes. Front Endocrinol. 2023;14:1270472.

    Article  Google Scholar 

  11. Lin H, Chen X, Wang L, Zhu T, Feng X, Liu X, Chen H, Pan S. Unraveling the role of microRNAs: potential biomarkers for gestational diabetes mellitus revealed through RNA sequencing analysis. Arch Gynecol Obstet. 2024;310(2):1255–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Lee SM, Hwangbo S, Norwitz ER, Koo JN, Oh IH, Choi ES, Jung YM, Kim SM, Kim BJ, Kim SY, et al. Nonalcoholic fatty liver disease and early prediction of gestational diabetes mellitus using machine learning methods. Clin Mol Hepatol. 2022;28(1):105–16.

    Article  PubMed  Google Scholar 

  13. Avalos GE, Owens LA, Dunne F. Applying current screening tools for gestational diabetes mellitus to a European population: is it time for change? Diabetes Care. 2013;36(10):3040–4.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Razo-Azamar M, Nambo-Venegas R, Meraz-Cruz N, Guevara-Cruz M, Ibarra-González I, Vela-Amieva M, Delgadillo-Velázquez J, Santiago XC, Escobar RF, Vadillo-Ortega F, et al. An early prediction model for gestational diabetes mellitus based on metabolomic biomarkers. Diabetol Metab Syndr. 2023;15(1):116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhan Y, Wang J, He X, Huang M, Yang X, He L, Qiu Y, Lou Y. Plasma metabolites, especially lipid metabolites, are altered in pregnant women with gestational diabetes mellitus. Clin Chim Acta. 2021;517:139–48.

    Article  CAS  PubMed  Google Scholar 

  16. Lee SM, Kwak SH, Koo JN, Oh IH, Kwon JE, Kim BJ, Kim SM, Kim SY, Kim GM, Joo SK, et al. Non-alcoholic fatty liver disease in the first trimester and subsequent development of gestational diabetes mellitus. Diabetologia. 2019;62(2):238–48.

    Article  CAS  PubMed  Google Scholar 

  17. Lee SM, Jung YM, Choi ES, Kwak SH, Koo JN, Oh IH, Kim BJ, Kim SM, Kim SY, Kim GM, et al. Metabolic dysfunction-associated fatty liver disease and subsequent development of adverse pregnancy outcomes. Clin Gastroenterol Hepatol. 2022;20(11):2542–e25502548.

    Article  PubMed  Google Scholar 

  18. Lee SM, Cho GJ, Wi WY, Norwitz ER, Koo BK, Lee J, Jung YM, Kwak SH, Park CW, Jun JK, et al. Metabolic dysfunction-associated fatty liver disease as a risk factor for adverse outcomes in subsequent pregnancy: a nationwide cohort study. Hepatol Int. 2023;17(2):367–76.

    Article  PubMed  Google Scholar 

  19. Chen Q, Francis E, Hu G, Chen L. Metabolomic profiling of women with gestational diabetes mellitus and their offspring: review of metabolomics studies. J Diabetes Complications. 2018;32(5):512–23.

    Article  PubMed  Google Scholar 

  20. Practice Bulletin No. Gestational diabetes mellitus. Obstet Gynecol. 2013;137(2 Pt 1):406–16.

    Google Scholar 

  21. Saadeh S, Younossi ZM, Remer EM, Gramlich T, Ong JP, Hurley M, Mullen KD, Cooper JN, Sheridan MJ. The utility of radiological imaging in nonalcoholic fatty liver disease. Gastroenterology. 2002;123(3):745–50.

    Article  PubMed  Google Scholar 

  22. Taylor KJ, Riely CA, Hammers L, Flax S, Weltin G, Garcia-Tsao G, Conn HO, Kuc R, Barwick KW. Quantitative US attenuation in normal liver and in patients with diffuse liver disease: importance of fat. Radiology. 1986;160(1):65–71.

    Article  CAS  PubMed  Google Scholar 

  23. Jung Y, Ahn E, Park T, Hwang GS. Longitudinal study investigating serum metabolites and their association with type 2 diabetes risk in a Korean population. Diabetes Obes Metab. 2023;25(8):2120–30.

    Article  CAS  PubMed  Google Scholar 

  24. Tingley D, Yamamoto T, Hirose K, Keele L, Imai K. Mediation: R package for causal mediation analysis. J Stat Softw. 2014;59(5):1–38.

    Article  Google Scholar 

  25. Artzi NS, Shilo S, Hadar E, Rossman H, Barbash-Hazan S, Ben-Haroush A, Balicer RD, Feldman B, Wiznitzer A, Segal E. Prediction of gestational diabetes based on nationwide electronic health records. Nat Med. 2020;26(1):71–6.

    Article  CAS  PubMed  Google Scholar 

  26. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44(3):837–45.

    Article  CAS  PubMed  Google Scholar 

  27. Xie J, Li L, Xing H. Metabolomics in gestational diabetes mellitus: a review. Clin Chim Acta. 2023;539:134–43.

    Article  CAS  PubMed  Google Scholar 

  28. Lynch CJ, Adams SH. Branched-chain amino acids in metabolic signalling and insulin resistance. Nat Rev Endocrinol. 2014;10(12):723–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Gaggini M, Carli F, Rosso C, Buzzigoli E, Marietti M, Della Latta V, Ciociaro D, Abate ML, Gambino R, Cassader M, et al. Altered amino acid concentrations in NAFLD: impact of obesity and insulin resistance. Hepatology. 2018;67(1):145–58.

    Article  CAS  PubMed  Google Scholar 

  30. Zhou Q, Wang Y, Gu Y, Li J, Wang H, Leng J, Li W, Yu Z, Hu G, Ma RCW et al. Genetic variants associated with beta-cell function and insulin sensitivity potentially influence bile acid metabolites and gestational diabetes mellitus in a Chinese population. BMJ Open Diabetes Res Care 2021, 9(1).

  31. Li J, Huo X, Cao YF, Li SN, Du Z, Shao P, Leng J, Zhang C, Sun XY, Ma RCW, et al. Bile acid metabolites in early pregnancy and risk of gestational diabetes in Chinese women: a nested case-control study. EBioMedicine. 2018;35:317–24.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Martineau MG, Raker C, Dixon PH, Chambers J, Machirori M, King NM, Hooks ML, Manoharan R, Chen K, Powrie R, et al. The metabolic profile of intrahepatic cholestasis of pregnancy is associated with impaired glucose tolerance, dyslipidemia, and increased fetal growth. Diabetes Care. 2015;38(2):243–8.

    Article  PubMed  Google Scholar 

  33. Martineau M, Raker C, Powrie R, Williamson C. Intrahepatic cholestasis of pregnancy is associated with an increased risk of gestational diabetes. Eur J Obstet Gynecol Reprod Biol. 2014;176:80–5.

    Article  PubMed  Google Scholar 

  34. Yue C, Ying C, Li X. Elevated serum uric acid is associated with gestational diabetes mellitus: an observational cohort study. J Clin Endocrinol Metab. 2023;108(7):e480–6.

    Article  PubMed  Google Scholar 

  35. Brennan P, Clare K, George J, Dillon JF. Determining the role for uric acid in non-alcoholic steatohepatitis development and the utility of urate metabolites in diagnosis: an opinion review. World J Gastroenterol. 2020;26(15):1683–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Copur S, Demiray A, Kanbay M. Uric acid in metabolic syndrome: does uric acid have a definitive role? Eur J Intern Med. 2022;103:4–12.

    Article  CAS  PubMed  Google Scholar 

  37. Law KP, Zhang H. The pathogenesis and pathophysiology of gestational diabetes mellitus: deductions from a three-part longitudinal metabolomics study in China. Clin Chim Acta. 2017;468:60–70.

    Article  CAS  PubMed  Google Scholar 

  38. Zhang Z, Zhou Z, Li H. The role of lipid dysregulation in gestational diabetes mellitus: early prediction and postpartum prognosis. J Diabetes Investig. 2024;15(1):15–25.

    Article  PubMed  Google Scholar 

  39. Chang W, Hatch GM, Wang Y, Yu F, Wang M. The relationship between phospholipids and insulin resistance: from clinical to experimental studies. J Cell Mol Med. 2018;23:702–10.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Wang Y, Huang Y, Wu P, Ye Y, Sun F, Yang X, Lu Q, Yuan J, Liu Y, Zeng H, et al. Plasma lipidomics in early pregnancy and risk of gestational diabetes mellitus: a prospective nested case-control study in Chinese women. Am J Clin Nutr. 2021;114(5):1763–73.

    Article  CAS  PubMed  Google Scholar 

  41. Zhu Y, Barupal DK, Ngo AL, Quesenberry CP, Feng J, Fiehn O, Ferrara A. Predictive metabolomic markers in early to mid-pregnancy for gestational diabetes mellitus: a prospective test and validation study. Diabetes. 2022;71(8):1807–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Lu L, Koulman A, Petry CJ, Jenkins B, Matthews L, Hughes IA, Acerini CL, Ong KK, Dunger DB. An unbiased lipidomics approach identifies early second trimester lipids predictive of maternal glycemic traits and gestational diabetes mellitus. Diabetes Care. 2016;39(12):2232–9.

    Article  CAS  PubMed  Google Scholar 

  43. De Souza LR, Berger H, Retnakaran R, Vlachou PA, Maguire JL, Nathens AB, Connelly PW, Ray JG. Non-alcoholic fatty liver disease in early pregnancy predicts dysglycemia in mid-pregnancy: prospective study. Am J Gastroenterol. 2016;111(5):665–70.

    Article  PubMed  Google Scholar 

  44. Hagström H, Höijer J, Ludvigsson JF, Bottai M, Ekbom A, Hultcrantz R, Stephansson O, Stokkeland K. Adverse outcomes of pregnancy in women with non-alcoholic fatty liver disease. Liver Int. 2016;36(2):268–74.

    Article  PubMed  Google Scholar 

  45. Porepa L, Ray JG, Sanchez-Romeu P, Booth GL. Newly diagnosed diabetes mellitus as a risk factor for serious liver disease. CMAJ. 2010;182(11):E526–531.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Després JP. Abdominal obesity as important component of insulin-resistance syndrome. Nutrition. 1993;9(5):452–9.

  47. Meng K, Lee CH, Saremi F. Metabolic syndrome and ectopic fat deposition: what can CT and MR provide? Acad Radiol. 2010;17(10):1302–12.

    Article  PubMed  Google Scholar 

  48. Larter CZ, Chitturi S, Heydet D, Farrell GC. A fresh look at NASH pathogenesis. Part 1: the metabolic movers. J Gastroenterol Hepatol. 2010;25(4):672–90.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

This study was supported by Korea Basic Science Institute (No. C523100 and C523400), the National Research Foundation of Korea (NRF) grants funded by the Korean government (MIST) (RS-2022-NR068419; RS-2024-00345936; 2021R1A2C2005820; RS-2021-NR056442; RS-2022-NR067269; RS-2023-00223831), the Research Program funded by the Korea National Institute of Health (KNIH) research project (2022ER090202), and Patient-Centered Clinical Research Coordinating Center (PACEN) funded by the Ministry of Health & Welfare, Republic of Korea (grant number: RS-2021-KH120178).

Author information

Authors and Affiliations

Authors

Contributions

Y. J., S.M.L., G.S.H., J.S.P., and W.K. conceived the idea and designed the study. S.M.L., J.N.K., I.H.O., K.H.K., B.J.K., S.M.K., S.Y.K., G.M.K., S.K.J., and D.H.L were responsible for the Incheon and Boramae cohort, and Y.J., J.L., Y.B., and G.S.H. acquired the metabolomics data and performed the analysis. Y.J., J.L., Y.K., and W. L. conducted the statistical analysis. Y.J. and S.M.L. wrote the manuscript, and G.S.H., J.S.P., and W.K. supervised the study. J.H.M., B.K.K., S.S., E.R.N., G.S.H., J.S.P., and W.K. critically reviewed and edited the manuscript. All authors contributed to the discussion, reviewed the manuscript, and approved the final version of the manuscript.

Corresponding authors

Correspondence to Geum-Sook Hwang, Joong Shin Park or Won Kim.

Ethics declarations

Competing interests

The authors declare no 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.

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jung, Y., Lee, S.M., Lee, J. et al. Metabolomic profiling reveals early biomarkers of gestational diabetes mellitus and associated hepatic steatosis. Cardiovasc Diabetol 24, 125 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12933-025-02645-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12933-025-02645-4

Keywords