- Open Access
Discovery of metabolic biomarkers for gestational diabetes mellitus in a Chinese population
Nutrition & Metabolism volume 18, Article number: 79 (2021)
Gestational diabetes mellitus (GDM), one of the most common pregnancy complications, can lead to morbidity and mortality in both the mother and the infant. Metabolomics has provided new insights into the pathology of GDM and systemic analysis of GDM with metabolites is required for providing more clues for GDM diagnosis and mechanism research. This study aims to reveal metabolic differences between normal pregnant women and GDM patients in the second- and third-trimester stages and to confirm the clinical relevance of these new findings.
Metabolites were quantitated with the serum samples of 200 healthy pregnant women and 200 GDM women in the second trimester, 199 normal controls, and 199 GDM patients in the third trimester. Both function and pathway analyses were applied to explore biological roles involved in the two sets of metabolites. Then the trimester stage-specific GDM metabolite biomarkers were identified by combining machine learning approaches, and the logistic regression models were constructed to evaluate predictive efficiency. Finally, the weighted gene co-expression network analysis method was used to further capture the associations between metabolite modules with biomarkers and clinical indices.
This study revealed that 57 differentially expressed metabolites (DEMs) were discovered in the second-trimester group, among which the most significant one was 3-methyl-2-oxovaleric acid. Similarly, 72 DEMs were found in the third-trimester group, and the most significant metabolites were ketoleucine and alpha-ketoisovaleric acid. These DEMs were mainly involved in the metabolism pathway of amino acids, fatty acids and bile acids. The logistic regression models for selected metabolite biomarkers achieved the area under the curve values of 0.807 and 0.81 for the second- and third-trimester groups. Furthermore, significant associations were found between DEMs/biomarkers and GDM-related indices.
Metabolic differences between healthy pregnant women and GDM patients were found. Associations between biomarkers and clinical indices were also investigated, which may provide insights into pathology of GDM.
Gestational diabetes mellitus (GDM), defined as diabetes diagnosed during pregnancy, affects approximately 15% of pregnant women globally . Several risk factors are highly correlated with the development of GDM, including maternal obesity, advanced age, family history of diabetes mellitus, and history of abnormal glucose metabolism [2, 3]. Studies have found that GDM can lead to several short-term and long-term complications, including obesity, impaired glucose metabolism, and cardiovascular diseases, for both mothers and infants . It was reported that in women with previous GDM the risk of developing diabetes after delivery was more than seven times that in women with normal glucose tolerance (NGT). For the offspring, they are more than two to eight times as likely to develop obesity, the metabolic syndrome, type 2 diabetes, and impaired insulin sensitivity and secretion . Although the pathogenesis of GDM is still not well understood, metabolomics has introduced new insights into the pathology of GDM and has yielded potential biomarkers related to GDM .
Metabolomics is the comprehensive analysis of low molecular weight compounds, known as metabolites, in biological systems. As the end products of metabolic processes, metabolites can reflect the internal physiological status of the organism that changes in response to environmental factors. Metabolic profiling is a useful approach for qualitative and quantitative metabolite studies of cell, bio-fluids, and tissue . In recent decades, metabolomics has been applied in various aspects, including disease diagnosis and treatment, metabolic pathway elucidation, biomarker discovery, and drug safety evaluation [8, 9]. For example, an untargeted metabolomics study has revealed a large number of pregnancy-related metabolic pathways and metabolites, five of which were able to predict gestational age in high accordance with ultrasound .
Many omics-based studies have been conducted for identifying biomarkers of GDM and for exploring the underlying mechanisms of its development [11, 12]. Previous proteomics studies have confirmed the downregulation of adiponectin among GDM women [13, 14]. Adiponectin could mediate antidiabetic metabolic effects through phosphorylation and activation of the 5'-AMP-activated protein kinase and acetyl coenzyme A carboxylase, thereby increasing fatty acid oxidation and glucose uptake in vivo. Conversely, downregulation of adiponectin may induce insulin resistance as well as GDM [15, 16]. Like proteomics, metabolomics could provide a deeper insight in the pathogenesis of GDM. Akturk et al. observed that asymmetric dimethylarginine was elevated in women with GDM during late pregnancy . Cetin et al. found that in GDM pregnancies, valine, methionine, phenylalanine, isoleucine, leucine, ornithine, glutamate, proline, and alanine were increased while glutamine was significantly decreased . Chen et al. revealed that 2-aminobutyric acid was associated with an increased likelihood of GDM in China . Moreover, 26 serum metabolites investigated by Liu et al. contributed to GDM, including 1-methyladenosine, homovanillic acid sulfate, and glucosamine compared with healthy pregnant women. These identified biomarkers are involved in some metabolic pathways that mainly participate in lipid, carbohydrate, and amino acid metabolisms . Although these investigations of metabolite profiles have identified branched-chain amino acids (BCAAs), aromatic amino acids, sulfur-containing amino acids, and other metabolites, their findings still lack consistency. Therefore, systemic analysis of the metabolites of GDM is necessary for obtaining more clues for GDM diagnosis and mechanism research.
In this study, we aimed to explore the metabolic difference between GDM women and normal pregnant women in the second and third trimester, using ultra-performance liquid chromatography coupled to tandem mass spectrometry (UPLC-MS/MS) system. Additionally, we expected to reveal the association of clinical indices with differentially expressed metabolites (DEMs) found in such metabolomics analysis.
Materials and methods
Study participant recruitment
All the samples were obtained from 200 healthy pregnant women and 200 GDM women in the second trimester, 199 normal controls, and 199 GDM patients in the third trimester at the University of Hong Kong-Shenzhen Hospital from 2016 to 2018. GDM patients were matched 1:1 with normal pregnant women. Matching was based on maternal age (± 3 years), pregestational BMI (± 3 kg/m2), and gestational week (± 3 weeks). Participants with cancer, kidney disease, heart disease, hepatic disease, other metabolic diseases or patients using medications that might affect glucolipid metabolism were excluded. The diagnostic criteria for GDM women were based on the standards recommended by the International Association of the Diabetes and Pregnancy Study Group. Glucose (75 g) was used for 2 h for conducting oral glucose tolerance test (OGTT). GDM was defined when fasting plasma glucose (FPG) was ≥ 5.1 mmol/L or 1-h plasma glucose was ≥ 10.0 mmol/L or 2-h plasma glucose was ≥ 8.5 mmol/L. All values for the OGTT less than the thresholds were considered normal. The study was approved by the ethics committee of the University of Hong Kong-Shenzhen Hospital (13). It was conducted according to relevant regulations, and informed consent was signed by every participant.
Sample collection and serum metabolomics
Age, height, and weight were recorded for every participant. General background information, including family history of diabetes, reproductive history, and medical history, were collected. Information about pregnancy outcomes was acquired after delivery, from the maternal and infant medical records in the hospital. Body mass index (BMI) was calculated as body weight (in kg)/height squared (in m2). Serum samples obtained from 24 to the end of 27 gestational weeks were defined as the second-trimester group, while those after 28 gestational weeks were defined as the third-trimester group. Blood samples were drawn in the morning after an overnight fast through the antecubital vein. Glucose levels were measured using the hexokinase method on a Roche Cobas 701 analyzer (Roche, Ltd, Basel, Switzerland). HbA1c values were tested using an Arkray HA-8160 analyzer (Arkray, Ltd, Kyoto, Japan). Serum lipids, including total triglyceride (TG), total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C) levels were measured with a Siemens ADVIA2400 fully automated chemistry analyzer (Siemens AG, Munich, Germany).
Metabolomic profiling was performed using Metabo-Profile (Shanghai, China). Serum samples were stored at − 80 °C until analysis. All of the standards of targeted metabolites were accurately weighed and prepared in water or methanol to obtain individual stock solution with a concentration of 5 mg/mL. Appropriate amount of each stock solution was mixed to create stock calibration solutions. A mixture of stable isotope labeled internal standards were prepared in methanol at a concentration of 50 µM/L. 25 µL of serum was added to a 96-well plate and the plate was transferred to the Biomek 4000 automation workstation (Biomek 4000, Beckman Coulter, Inc., California, USA). Approximately 120 µL of ice-cold methanol with partial internal standards was automatically added to each sample for extracting the metabolites. After vortexing for 5 min, the mixture was centrifuged for 30 min at 4000 g; 30 µL of supernatant was transferred to a new 96-well plate, and 20 µL of freshly prepared derivative reagents (3-Nitrophenylhydrazine) was added to each well. After derivatization for 60 min at 30 °C, 330 µL of ice-cold 50% methanol solution was added for dilution. The samples were stored at − 20 °C for 20 min and were centrifuged at 4000 g for 30 min at 4 °C. Approximately 135 µL of supernatant was mixed with 10 µL of internal standards in each well of a new 96-well plate. Serial dilutions of derivatized stock standards were added to the left wells, and the plate was ready for analysis.
Chromatographic separation was performed on an ACQUITY UPLC BEH C18 VanGuard pre-column (2.1×5 mm, 1.7 µm) and an ACQUITY UPLC BEH C18 analytical column (2.1×100 mm, 1.7 µm). The temperature of the column and the sample manager was set at 40 °C and 10 °C, respectively. The mobile phase A was water with 0.1% formic acid, while B was a mixture of acetonitrile and isopropanol (70:30). Gradient conditions were 0–1 min, 5% B; 1–11 min, 5%–78% B; 11–13.5 min, 78%–95% B; 13.5–14 min, 95%–100% B; 14–16 min, 100% B; 16–16.1 min, 100%–5% B; 16.1–18 min, 5% B. The flow rate was 0.4 mL/min with a 5 µL injection volume.
The mass spectrometer was operated in positive electrospray ionization (ESI +) mode with a capillary voltage of 1.5 kV as well as the negative electrospray ionization (ESI −) mode with a capillary voltage of 2 kV. The temperature of the ion source and desolvation was 150 °C and 550 °C, respectively. The desolvation gas flow was set at 1000 L/h. Raw data files generated by UPLC-MS/MS were processed using Masslynx software (v4.1, Waters, Milford, MA, USA) for performing peak integration, calibration, and quantitation for each metabolite. An overview of the research design and analysis workflow is presented in Fig. 1.
Bioinformatics analysis and metabolic biomarker recognition
For the analysis of clinical characteristics, participants were divided into four groups: the second-trimester GDM group (group 1), second-trimester normal group (group 2), third-trimester GDM group (group 3), and third-trimester normal group (group 4). Clinical characteristics between the same gestational-stage groups (group1 vs. group 2 and group3 vs. group 4) were compared using the Student’s paired t-test or signed-rank test for continuous variables and chi-square tests for categorical variables with SAS 9.3 (SAS Institute, Cary, NC, USA). After normality testing for continuous variables, the data was presented as mean ± standard deviation (SD) or median (interquartile range). The level of significance was set at p < 0.05.
To observe changes in metabolites between the same gestational-stage groups, the univariate and multivariate analyses were conducted, including differential expression analysis using t-test, fold-change with volcano plot, principal component analysis (PCA), partial least square discriminant analysis (PLS-DA), and orthogonal partial least square discriminant analysis (OPLS-DA). Data processing was performed using the iMAP platform (v1.0, Metabo-Profile, Shanghai, China). p < 0.05 and |log2FC|> 0 in the univariate analysis and variable importance in the projection (VIP) > 1 in the multivariate analysis were the standards for screening DEMs. To investigate the biological meaning of DEMs, enrichment analysis was performed with the small molecule pathway database (SMPDB) with p < 0.05 regarded as the significance level, and pathway analysis was also performed using p values to account for significance and false discovery. A threshold of p < 0.05 was used for considering the function or pathway to be impactful. Metabolites were further analyzed using random forest (RF) for identifying potential biomarkers based on which logistic regression (LR) models were constructed for evaluating predictive efficiency. Area under the curve (AUC) for receiver operating characteristic (ROC) curves were calculated for assessing the performance of potential biomarkers with LR for GDM. The determination of sensitivity and specificity at the optimal cut-off point was defined by the minimum distance to the top-left corner. Longitudinal analysis for catching continue-changed metabolic biomarkers with the disease was performed by extracting the shared candidate metabolic biomarkers from above classification-model evaluated biomarkers in two pregnancy stages. The weighted gene co-expression network analysis (WGCNA)  was used for inferring the association between metabolite modules and clinical indices.
Clinical characteristic of study participants
The clinical characteristics of all participants are presented in Table 1. Compared with normal pregnant women, pregnant women who developed GDM had higher FPG, 1-h, and 2-h blood glucose levels after an OGTT and smaller changes in BMI during pregnancy, in both the second and third trimesters. HbA1c, TG, TC, and LDL were statistically different between the case and control groups in the second trimester, which were consistent with previous research conclusions [22, 23]. TC, delivery gestational week, neonatal weight, and neonatal length were statistically different between GDM patients and healthy pregnant women in the third trimester. Women with GDM were more likely to have a family history of diabetes. Age, BMI, and gestational week were similar between the two groups in each trimester.
Metabolomics profiling of study participants
As shown in the metabolomics profiles (Fig. 2a, b), the case group (GDM) did not separate from the control group (normal) in the second and third trimester with PCA; moreover, PLS-DA showed a relatively clear discrimination for both the trimesters (Fig. 2c, d). Finally, the results of OPLS-DA (Fig. 2e, f) indicated the possibility of evaluating the differences between GDM patients and normal controls with metabolite abundance. The results showed that R2 and Q2 of the OPLS-DA model in the second trimester were 0.347 and 0.165, while in the third trimester were 0.324 and 0.201, respectively. Permutation tests (n = 200) were employed for validating the predictive ability of the built OPLS-DA models (Additional file 1: Figure S1). The R2 and Q2 values derived from the permuted data were lower than the original values, which demonstrated that the OPLS-DA model did not overfit. The VIP values of the OPLS-DA model and criteria, including |log2FC|> 0 with p < 0.05, were further employed for determining the DEMs. A volcano plot provided a quick way to display the differences in metabolite expression levels between normal pregnant women and GDM patients with statistical significance (Fig. 3). As shown in the volcano plot, of 200 metabolites found in this study, 57 metabolites in the second-trimester group (Additional file 4: Table S1) and 72 metabolites in the third-trimester group (Additional file 5: Table S2) were considered DEMs. The top five DEMs in the second-trimester group were 3-methyl-2-oxovaleric acid, 3-hydroxybutyric acid, palmitic acid, alpha-hydroxyisobutyric acid, and acetic acid. As for the third-trimester group, they were ketoleucine, alpha-ketoisovaleric acid, pyruvic acid, L-tryptophan, and succinic acid. Most DEMs were down-regulated during the second trimester (Fig. 3a), while many DEMs were positively-regulated during the third trimester (Fig. 3b), suggesting different metabolomic dysfunction at different trimester stages. Amino acids, fatty acids, and organic acids account for around 70% of DEMs found in the second- and third-trimester stages (Fig. 3e, f).
Metabolic enrichment of biological function and pathway relevant to GDM
The DEMs for each comparison group were evaluated using enrichment analysis with SMPDB (Fig. 4a, b). In the second-trimester group, the alpha linolenic acid and linoleic acid metabolism pathways had the highest fold enrichment, lowest p value (p < 0.001), and FDR of < 0.1. Other significant functions included beta oxidation of very long chain fatty acids and valine-leucine-isoleucine degradation (Fig. 4a and Table 2, Additional file 6: Table S3). In the third-trimester group, functions such as urea cycle, ammonia recycling, glycine and serine metabolism, valine-leucine-isoleucine degradation, arginine and proline metabolism, alanine metabolism, glutamate metabolism, aspartate metabolism, glucose-alanine cycle, phenylalanine and tyrosine metabolism, and carnitine synthesis were significantly associated with the corresponding DEMs (Fig. 4b and Table 3, Additional file 7: Table S4).
Pathway analysis was also performed for investigating the function of DEMs. 32 pathways were observed, nine of which were significantly enriched in the second-trimester group (Fig. 4c, Additional file 8: Table S5), among which two pathways, valine-leucine-isoleucine biosynthesis and valine-leucine-isoleucine degradation, played key roles in reflecting the changes in metabolites. For the third-trimester group, 48 pathways were found, of which 21 were significantly enriched with DEMs (Fig. 4d, Additional file 9: Table S6).
Furthermore, remarkable differences also exist between different trimester stages on function and pathway levels (Fig. 4 and Additional file 2: Fig. S2), suggesting that stage-specific biomarkers and diagnostic models should be considered.
Selection of potential metabolic biomarkers for GDM
After observing metabolomics differences between the groups and reliable functional enrichment analysis, it was necessary to establish a diagnostic model for predicting the presence of GDM in pregnant women and for selecting potential metabolite characteristics with the importance determined using machine learning algorithms [24, 25].
For the second-trimester group, samples were first divided into training data with 70% samples and test data with 30% samples. The RF model was then learned on the training data for obtaining the importance score for each metabolite, based on which the candidate metabolite biomarkers were selected with top large importance scores and with overlaps to DEMs as much as possible. Next, the LR model was constructed on the training data with such metabolite biomarkers (Additional file 10: Table S7), with a training AUC of 0.808. Finally, the second-trimester group-specific LR model was validated on the test data and achieved a testing AUC of 0.807 (Fig. 5a). Similarly, for the third-trimester group, the metabolite biomarkers (Additional file 10: Table S7) were selected using RF, and an LR model was built, which had learning performance as AUC 0.819 and validation performance as AUC 0.810 on the test data (Fig. 5b). In practice, patients diagnosed as positive by our new diagnostic method have about 30% probability of actual disease when the prevalence of GDM is 11.91% . The positive predictive value in the second-trimester group is 0.328, and the negative predictive value is 0.964 (when sensitivity and specificity are all 0.783, the best diagnostic bounds), while the positive predictive value in the third-trimester group is 0.269, and the negative predictive value is 0.958 (when sensitivity is 0.767 and specificity is 0.717). In addition, almost all identified biomarker candidates were DEMs (Table 4). Of note, there are three common potential metabolic biomarkers (3-hydroxybutyric acid, isobutyric acid, and isovaleric acid) discovered both in second- and third-trimester groups, whose trend changes longitudinally across pregnancy were displayed in Additional file 3: Figure S3. As the pregnancy progressed, these three metabolic biomarkers significantly increased among normal pregnant women, while in GDM patients only 3-hydroxybutyric acid was increased. We also assess the goodness of fit between the predicted and real probabilities with calibration plot (Fig. 5c, d), evaluating how similar the predicted probability values are to the actual probabilities.
Clinical relevance of metabolic biomarkers for GDM
Finally, to confirm the clinical relevance of metabolic biomarkers associated with GDM, WGCNA was performed for inferring the association between metabolite modules and clinical indices. As shown in Fig. 6a, four modules were detected for the second-trimester group. Module turquoise was significantly associated with GDM (i.e., group index) and many other important clinical indices, including pre-pregnancy BMI, OGTT, TC, TG, and LDL. Additionally, this module significantly included 11 DEMs (P = 9.99e−04), compared to module blue containing three DEMs (P = 0.666), module brown containing two DEMs (P = 0.580), and module grey containing 41 DEMs (P = 0.868) (Fig. 6c). Similarly, five modules were found for the third-trimester group (Fig. 6b), where the module turquoise contains 21 DEMs (P = 6.737e−08), module blue contains five DEMs (P = 0.944), module brown contains six DEMs (P = 0.090), module grey contains 37 DEMs (P = 0.992), and module yellow contains three DEMs (P = 0.606) (Fig. 6d). At this time, module turquoise was associated with group, OGTT, TC, TG, and LDL; and module yellow was positively associated with the group, pre-pregnancy BMI, and OGTT. It should be noted that the biomarker candidates are almost DEMs (Table 4); thus, they have similar contribution to particular modules and their associations with GDM.
Using UPLC-MS/MS for metabolomics analysis, detailed information was obtained on the metabolic changes in normal pregnant women and GDM patients in our study. The changes in serum metabolites were further investigated using univariate and multivariate statistical analyses. Thirteen candidate metabolite biomarkers found in the second-trimester group and thirteen found in the third-trimester group contributed to GDM when compared with healthy pregnant women. According to metabolic enrichment and pathway analyses, valine-leucine-isoleucine degradation in the second-trimester group was consistently found in both analyses. In the third trimester, valine-leucine-isoleucine degradation and glycine, serine, arginine, proline, alanine, glutamate, aspartate, and phenylalanine metabolisms were consistent in both the analyses. Metabolic biomarkers have been found by RF, and LR models based on which showed high predictive efficiency. Furthermore, these biomarkers demonstrate remarkable relationship with clinical indices.
BCAAs, which consist of valine, leucine, and isoleucine, were higher in the GDM group than in the NGT group in our study. Although some studies have found that the levels of BCAAs did not differ significantly between GDM and NGT [27, 28], several studies have shown that elevated BCAAs in GDM patients might serve as biomarkers for GDM [29,30,31]. Phosphatidylinositol-3-kinase (PI3K)/ protein kinase B (Akt)/mammalian target of rapamycin (mTOR) signaling is involved in functions such as protein and lipid biosynthesis . Mammalian target of rapamycin complex 1 (mTORC1) is controlled by the small GTPase Ras homolog enriched in brain (Rheb) and tuberous sclerosis complex (TSC) complex (TSC1 and TSC2) . When insulin exists, it activates insulin receptor and insulin receptor substrate (IRS), which further promotes the activation of PI3K and Akt. Akt inhibits the TSC complex, allowing Rheb to stimulate mTORC1 to mediate cellular functions . Ribosomal S6 protein kinase 1 (S6K1), one of the mTORC1 substrate, participates in the Rheb/mTOR/S6K pathway. However, constitutive activation of the Rheb/mTOR/S6K pathway can induce a negative feedback and thus cause insulin resistance . BCAAs might be associated with insulin resistance in GDM women [36, 37]. One potential mechanism is that elevated BCAAs levels lead to activation of the mTOR/S6K1 pathway and serine phosphorylation of IRS1, contributing to inhibition of insulin-induced PI3K activation and insulin resistance . Also, BCAAs dysmetabolism could result in the accumulation of toxic BCAA metabolites that cause β-cell mitochondrial dysfunction and highly susceptibility to insulin resistance .
In our study, glycine and serine metabolisms were lower in the GDM group, which is consistent with findings of previous studies. Takashina et al. observed that fasting glucose and 2-h plasma glucose levels or the homeostasis model assessment of insulin resistance negatively correlated with glycine levels, and the homeostasis model assessment for the β-cell function index negatively correlated with glycine and serine levels . Moreover, oral glycine has been reported to increase insulin secretion without affecting insulin sensitivity .
Additionally, arginine and proline levels were also lower in the GDM group. It has been reported that arginine and its metabolites promote insulin secretion  and improve insulin resistance in humans . Arginine plays multiple beneficial roles against metabolic abnormalities, but it might also induce oxidative stress . Proline is absorbed and metabolized into glutamine, which may enter the tricarboxylic cycle and ultimately be converted into glucose. One study demonstrated that the ingestion of proline with glucose attenuated the glucose area response without affecting insulin response and decreased glucagon levels compared to glucose alone .
Alpha linolenic acid is a precursor of polyunsaturated fatty acids, which mainly contain omega-6 and omega-3 fatty acids. A meta-analysis by Zhong et al. investigated the efficacy of omega-3 fatty acid for GDM, revealing that omega-3 fatty acids supplementation in GDM patients could reduce FPG and HOMA-IR score . Omega-3 fatty acids were involved in the mechanism of increasing β-oxidation of fatty acids, improving antioxidant functions and insulin action, and reducing lipogenesis .
Short-chain fatty acids (SCFAs) refer to those fatty acids that have one to six carbons. Previous study has indicated that increased expression of free fatty acid receptor-2 and alteration of its endogenous ligands SCFAs contributed to glucose homeostasis by improving insulin secretion throughout gestation . SCFAs also diminished late gestational androgen excess through suppression of adenosine deaminase /xanthine oxidase pathway, which protected against glucose dysmetabolism and poor fetal outcome . SCFAs can not only impact metabolism, but also influence cardiovascular diseases. SCFAs are recognized by G protein-coupled receptors, one of which is G protein-coupled receptor 41(Gpr41). Another SCFA receptor olfactory receptor 78, together with Gpr41, participated in the modulation of blood pressure. While Gpr41 lowered baseline blood pressure [50, 51], olfactory receptor 78 induced the hypertensive effect of SCFAs . Acetate, one of major SCFAs, decreases among GDM patients in our study. It was believed to enhance glucose incorporation and lipid metabolism through activating AMP-activated protein kinase [52, 53]. Additionally, supplementation of acetate significantly increased the number of activated Treg cells , which could attenuate cardiac hypertrophy and fibrosis and improve electric remodeling in hypertension/Angiotensin II–induced cardiac damage .
Of note, 3-hydroxybutyric acid was selected as a potential metabolic biomarker in both the trimester groups. As a classic ketone body, the levels of 3-hydroxybutyric acid increases because of the oxidation of free fatty acids and excess acetyl-CoA. ATP production from fatty acids and carbohydrate oxidation happens out of control, resulting in increased acetyl-CoA levels. A study in diabetic rats showed that inefficient utilization and mobilization of glucose may contribute to the elevation of 3-hydroxybutyric acid levels .
In the WGCNA, carnitine was included in the module grey for both the trimesters. Previous studies have shown that using 2 g/day of L-carnitine resulted in a reduction of TC and LDL, and its mechanism may be related to the phenomena of insulin resistance and lipotoxicity . However, a higher dose of L-carnitine had different effects, since it contributed to the elevation of TG, apolipoprotein-A1, and apolipoprotein-B100 levels . While BCAAs are in the module grey in the second-trimester group, isoleucine belongs to the module grey, leucine and valine are in the module turquoise for the third-trimester group. It is widely accepted that BCAAs transaminase helps in the conversion of isoleucine and valine into branched-chain α-ketoacids, which are further transformed into propinonyl-CoA by the branched-chain α-ketoacid dehydrogenase complex. Propinonyl-CoA can become methylmalonyl-CoA with relevant carboxylase. Methylmalonyl-CoA mutase (MUT) is an enzyme that catalyzes the conversion of methylmalonyl-CoA to succinyl-CoA . Based on experiments with mice, decreased Mut expression led to higher body weight, hyperinsulinemia, elevated fasting glucose and increased triglyceride .
Actually, this study has some limitations. First, in a cross-sectional study design, metabolites were only measured at one point; thus, further prospective cohort studies are needed for establishing the dynamic association of these metabolites with GDM. Second, though metabolite biomarkers in the first trimester allow early diagnosis and timely intervention, we mainly focused on the metabolic changes during the second and third trimesters due to lack of enough participants recruited in the first trimester stage. Third, while the second trimester starts at week 14 of pregnancy and lasts through the end of week 27, participants were recruited from 24 weeks to the end of 27 gestational weeks for the second-trimester group; thus, important information may have missed out. Fourth, the precise molecular mechanisms underlying the development of GDM remain unclear and mechanistic studies need to be conducted for clarifying the exact roles of these discovered metabolites in GDM.
In conclusion, with the analysis of serum samples, our study suggested that specific metabolomic profile existed among GDM patients. Several key metabolites, such as glycine, serine, proline, and 3-hydroxybutyric acid, are associated with GDM in the second or third trimester, and potential biomarkers for GDM have been identified. These metabolites mainly participated in fatty acid and amino acid metabolism, which may shed light on the pathology of GDM. However, further research is needed to confirm our findings and explore the underlying molecular mechanisms.
Availability of data and materials
The datasets used or analyzed during the current study are available from the corresponding author (Cheng Hu) on reasonable request.
Protein kinase B
Area under the curve
Branched-chain amino acids
Body mass index
Differentially expressed metabolites
Fasting plasma glucose
Gestational diabetes mellitus
G protein-coupled receptor 41
High-density lipoprotein cholesterol
Insulin receptor substrate
Low-density lipoprotein cholesterol
Mammalian Target of Rapamycin
Mammalian Target of Rapamycin Complex 1
Normal glucose tolerance
Oral glucose tolerance test
Orthogonal partial least square discriminant analysis
Principal component analysis
Partial least square discriminant analysis
Ras homolog enriched in brain
Receiver operating characteristic
Ribosomal S6 protein kinase 1
Short-chain fatty acids
Small molecule pathway database
Tuberous sclerosis complex
Ultra-performance liquid chromatography coupled to tandem mass spectrometry
Variable Important for the Projection
Weighted gene co-expression network analysis
Zhu Y, Zhang C. Prevalence of gestational diabetes and risk of progression to type 2 diabetes: a global perspective. Curr Diab Rep. 2016;16:7.
Jovanovic L, Pettitt DJ. Gestational diabetes mellitus. JAMA. 2001;286:2516–8.
The L. Gestational diabetes in England: cause for concern. Lancet. 2019;393:1262.
McIntyre HD, Catalano P, Zhang C, Desoye G, Mathiesen ER, Damm P. Gestational diabetes mellitus. Nat Rev Dis Primers. 2019;5:47.
Damm P, Houshmand-Oeregaard A, Kelstrup L, Lauenborg J, Mathiesen ER, Clausen TD. Gestational diabetes mellitus and long-term consequences for mother and offspring: a view from Denmark. Diabetologia. 2016;59:1396–9.
Huynh J, Xiong G, Bentley-Lewis R. A systematic review of metabolite profiling in gestational diabetes mellitus. Diabetologia. 2014;57:2453–64.
Nicholson JK, Lindon JC. Systems biology: metabonomics. Nature. 2008;455:1054–6.
Vogeser M, Seger C. A decade of HPLC-MS/MS in the routine clinical laboratory–goals for further developments. Clin Biochem. 2008;41:649–62.
Hsieh Y. HPLC-MS/MS in drug metabolism and pharmacokinetic screening. Expert Opin Drug Metab Toxicol. 2008;4:93–101.
Liang L, Rasmussen MH, Piening B, Shen X, Chen S, Rost H, Snyder JK, Tibshirani R, Skotte L, Lee NC, et al. Metabolic dynamics and prediction of gestational age and time to delivery in pregnant women. Cell. 2020;181:1680–92.
Lorenzo-Almoros A, Hang T, Peiro C, Soriano-Guillen L, Egido J, Tunon J, Lorenzo O. Predictive and diagnostic biomarkers for gestational diabetes and its associated metabolic and cardiovascular diseases. Cardiovasc Diabetol. 2019;18:140.
McIntyre HD, Kapur A, Divakar H, Hod M. Gestational diabetes mellitus-innovative approach to prediction, diagnosis, management, and prevention of future NCD-mother and offspring. Front Endocrinol (Lausanne) 2020, 11:614533.
Ravnsborg T, Svaneklink S, Andersen LLT, Larsen MR, Jensen DM, Overgaard M: First-trimester proteomic profiling identifies novel predictors of gestational diabetes mellitus. PLoS One 2019, 14:e0214457.
Ravnsborg T, Andersen LL, Trabjerg ND, Rasmussen LM, Jensen DM, Overgaard M. First-trimester multimarker prediction of gestational diabetes mellitus using targeted mass spectrometry. Diabetologia. 2016;59:970–9.
Yamauchi T, Kamon J, Minokoshi Y, Ito Y, Waki H, Uchida S, Yamashita S, Noda M, Kita S, Ueki K, et al. Adiponectin stimulates glucose utilization and fatty-acid oxidation by activating AMP-activated protein kinase. Nat Med. 2002;8:1288–95.
Yamauchi T, Kamon J, Ito Y, Tsuchida A, Yokomizo T, Kita S, Sugiyama T, Miyagishi M, Hara K, Tsunoda M, et al. Cloning of adiponectin receptors that mediate antidiabetic metabolic effects. Nature. 2003;423:762–9.
Akturk M, Altinova A, Mert I, Dincel A, Sargin A, Buyukkagnici U, Arslan M, Danisman N. Asymmetric dimethylarginine concentrations are elevated in women with gestational diabetes. Endocrine. 2010;38:134–41.
Cetin I, de Santis MS, Taricco E, Radaelli T, Teng C, Ronzoni S, Spada E, Milani S, Pardi G. Maternal and fetal amino acid concentrations in normal pregnancies and in pregnancies with gestational diabetes mellitus. Am J Obstet Gynecol. 2005;192:610–7.
Chen X, de Seymour JV, Han TL, Xia Y, Chen C, Zhang T, Zhang H, Baker PN. Metabolomic biomarkers and novel dietary factors associated with gestational diabetes in China. Metabolomics. 2018;14:149.
Liu T, Li J, Xu F, Wang M, Ding S, Xu H, Dong F. Comprehensive analysis of serum metabolites in gestational diabetes mellitus by UPLC/Q-TOF-MS. Anal Bioanal Chem. 2016;408:1125–35.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559.
Kwon SS, Kwon JY, Park YW, Kim YH, Lim JB. HbA1c for diagnosis and prognosis of gestational diabetes mellitus. Diabetes Res Clin Pract. 2015;110:38–43.
Wang J, Li Z, Lin L. Maternal lipid profiles in women with and without gestational diabetes mellitus. Medicine (Baltimore) 2019, 98:e15320.
Yu XT, Wang ZJ, Zeng T. Essential gene expression pattern of head and neck squamous cell carcinoma revealed by tumor-specific expression rule based on single-cell RNA sequencing. Biochimica Et Biophys. Acta-Mol. Basis Dis. 2020, 1866.
Yu XT, Gan ZJ, Xu Y, Wan SB, Li M, Ding SJ, Zeng T. Identifying essential methylation patterns and genes associated with stroke. IEEE Access. 2020;8:96669–76.
Nguyen CL, Pham NM, Binns CW, Duong DV, Lee AH. Prevalence of gestational diabetes mellitus in Eastern and Southeastern Asia: a systematic review and meta-analysis. J Diabetes Res. 2018;2018:6536974.
Bentley-Lewis R, Huynh J, Xiong G, Lee H, Wenger J, Clish C, Nathan D, Thadhani R, Gerszten R. Metabolomic profiling in the prediction of gestational diabetes mellitus. Diabetologia. 2015;58:1329–32.
Pappa KI, Vlachos G, Theodora M, Roubelaki M, Angelidou K, Antsaklis A. Intermediate metabolism in association with the amino acid profile during the third trimester of normal pregnancy and diet-controlled gestational diabetes. Am J Obstet Gynecol. 2007;196(65):e61-65.
Park S, Park JY, Lee JH, Kim SH. Plasma levels of lysine, tyrosine, and valine during pregnancy are independent risk factors of insulin resistance and gestational diabetes. Metab Syndr Relat Disord. 2015;13:64–70.
Andersson-Hall U, Gustavsson C, Pedersen A, Malmodin D, Joelsson L, Holmang A. Higher concentrations of BCAAs and 3-HIB are associated with insulin resistance in the transition from gestational diabetes to type 2 diabetes. J Diabetes Res. 2018;2018:4207067.
Fitch WL, King JC. Plasma amino acid, glucose, and insulin responses to moderate-protein and high-protein test meals in pregnant, nonpregnant, and gestational diabetic women. Am J Clin Nutr. 1987;46:243–9.
Bertacchini J, Heidari N, Mediani L, Capitani S, Shahjahani M, Ahmadzadeh A, Saki N. Targeting PI3K/AKT/mTOR network for treatment of leukemia. Cell Mol Life Sci. 2015;72:2337–47.
Dibble CC, Cantley LC. Regulation of mTORC1 by PI3K signaling. Trends Cell Biol. 2015;25:545–55.
Yu Y, Yoon SO, Poulogiannis G, Yang Q, Ma XM, Villen J, Kubica N, Hoffman GR, Cantley LC, Gygi SP, Blenis J. Phosphoproteomic analysis identifies Grb10 as an mTORC1 substrate that negatively regulates insulin signaling. Science. 2011;332:1322–6.
Shah OJ, Wang Z, Hunter T. Inappropriate activation of the TSC/Rheb/mTOR/S6K cassette induces IRS1/2 depletion, insulin resistance, and cell survival deficiencies. Curr Biol. 2004;14:1650–6.
Kimball SR, Jefferson LS. Signaling pathways and molecular mechanisms through which branched-chain amino acids mediate translational control of protein synthesis. J Nutr. 2006;136:227S-231S.
O’Connell TM. The complex role of branched chain amino acids in diabetes and cancer. Metabolites. 2013;3:931–45.
Um SH, D’Alessio D, Thomas G. Nutrient overload, insulin resistance, and ribosomal protein S6 kinase 1, S6K1. Cell Metab. 2006;3:393–402.
Lynch CJ, Adams SH. Branched-chain amino acids in metabolic signalling and insulin resistance. Nat Rev Endocrinol. 2014;10:723–36.
Takashina C, Tsujino I, Watanabe T, Sakaue S, Ikeda D, Yamada A, Sato T, Ohira H, Otsuka Y, Oyama-Manabe N, et al. Associations among the plasma amino acid profile, obesity, and glucose metabolism in Japanese adults with normal glucose tolerance. Nutr Metab (Lond). 2016;13:5.
Gonzalez-Ortiz M, Medina-Santillan R, Martinez-Abundis E, von Drateln CR. Effect of glycine on insulin secretion and action in healthy first-degree relatives of type 2 diabetes mellitus patients. Horm Metab Res. 2001;33:358–60.
Schmidt HH, Warner TD, Ishii K, Sheng H, Murad F. Insulin secretion from pancreatic B cells caused by L-arginine-derived nitrogen oxides. Science. 1992;255:721–3.
Monti LD, Galluccio E, Villa V, Fontana B, Spadoni S, Piatti PM. Decreased diabetes risk over 9 year after 18-month oral L-arginine treatment in middle-aged subjects with impaired glucose tolerance and metabolic syndrome (extension evaluation of L-arginine study). Eur J Nutr. 2018;57:2805–17.
Carvalho AM, Oliveira AA, Loureiro AP, Gattas GJ, Fisberg RM, Marchioni DM. Arginine intake is associated with oxidative stress in a general population. Nutrition. 2017;33:211–5.
Nuttall FQ, Gannon MC, Jordan K. The metabolic response to ingestion of proline with and without glucose. Metabolism. 2004;53:241–6.
Zhong N, Wang J. The efficacy of omega-3 fatty acid for gestational diabetes: a meta-analysis of randomized controlled trials. Gynecol Endocrinol. 2019;35:4–9.
Backes J, Anzalone D, Hilleman D, Catini J. The clinical relevance of omega-3 fatty acids in the management of hypertriglyceridemia. Lipids Health Dis. 2016;15:118.
Fuller M, Priyadarshini M, Gibbons SM, Angueira AR, Brodsky M, Hayes MG, Kovatcheva-Datchary P, Backhed F, Gilbert JA, Lowe WL Jr, Layden BT. The short-chain fatty acid receptor, FFA2, contributes to gestational glucose homeostasis. Am J Physiol Endocrinol Metab. 2015;309:E840-851.
Usman TO, Areola ED, Badmus OO, Kim I, Olatunji LA. Sodium acetate and androgen receptor blockade improve gestational androgen excess-induced deteriorated glucose homeostasis and antioxidant defenses in rats: roles of adenosine deaminase and xanthine oxidase activities. J Nutr Biochem. 2018;62:65–75.
Natarajan N, Hori D, Flavahan S, Steppan J, Flavahan NA, Berkowitz DE, Pluznick JL. Microbial short chain fatty acid metabolites lower blood pressure via endothelial G protein-coupled receptor 41. Physiol Genomics. 2016;48:826–34.
Pluznick JL, Protzko RJ, Gevorgyan H, Peterlin Z, Sipos A, Han J, Brunet I, Wan LX, Rey F, Wang T, et al. Olfactory receptor responding to gut microbiota-derived signals plays a role in renin secretion and blood pressure regulation. Proc Natl Acad Sci USA. 2013;110:4410–5.
Maruta H, Yoshimura Y, Araki A, Kimoto M, Takahashi Y, Yamashita H. Activation of AMP-activated protein kinase and stimulation of energy metabolism by acetic acid in L6 myotube cells. PLoS One 2016, 11:e0158055.
Sakakibara S, Yamauchi T, Oshima Y, Tsukamoto Y, Kadowaki T. Acetic acid activates hepatic AMPK and reduces hyperglycemia in diabetic KK-A(y) mice. Biochem Biophys Res Commun. 2006;344:597–604.
Kaye DM, Shihata WA, Jama HA, Tsyganov K, Ziemann M, Kiriazis H, Horlock D, Vijay A, Giam B, Vinh A, et al. Deficiency of prebiotic fiber and insufficient signaling through gut metabolite-sensing receptors leads to cardiovascular disease. Circulation. 2020;141:1393–403.
Kvakan H, Kleinewietfeld M, Qadri F, Park JK, Fischer R, Schwarz I, Rahn HP, Plehm R, Wellner M, Elitok S, et al. Regulatory T cells ameliorate angiotensin II-induced cardiac damage. Circulation. 2009;119:2904–12.
Oshida Y, Iwao N, Ohsawa I, Sato J, Nakao T, Sato Y. Effect of insulin on intramuscular 3-hydroxybutyrate levels in diabetic rats. Horm Metab Res. 1998;30:70–1.
Vidal-Casariego A, Burgos-Pelaez R, Martinez-Faedo C, Calvo-Gracia F, Valero-Zanuy MA, Luengo-Perez LM, Cuerda-Compes C. Metabolic effects of L-carnitine on type 2 diabetes mellitus: systematic review and meta-analysis. Exp Clin Endocrinol Diabetes. 2013;121:234–8.
Rahbar AR, Shakerhosseini R, Saadat N, Taleban F, Pordal A, Gollestan B. Effect of L-carnitine on plasma glycemic and lipidemic profile in patients with type II diabetes mellitus. Eur J Clin Nutr. 2005;59:592–6.
Sivanand S, Vander Heiden MG. Emerging roles for branched-chain amino acid metabolism in cancer. Cancer Cell. 2020;37:147–56.
Lerin C, Goldfine AB, Boes T, Liu M, Kasif S, Dreyfuss JM, De Sousa-Coelho AL, Daher G, Manoli I, Sysol JR, et al. Defects in muscle branched-chain amino acid oxidation contribute to impaired lipid metabolism. Mol Metab. 2016;5:926–36.
We thank Metabo-profile Biotechnology (Shanghai, China) for this metabolomics profiling and metabolomic analysis.
This study was supported by the National Natural Science Foundation of China (61803360 and 81974118), the Shanghai Municipal Science and Technology Major Project (2017SHZDZX01), the Postdoctoral Science Foundation of China (2019M660088 and 2020T13408), and the Shanghai Science and Technology Innovation Action Plan (20XD1433300).
Ethics approval and consent to participate
The study was approved by the ethics committee of the University of Hong Kong-Shenzhen Hospital. It was conducted according to the Declaration of Helsinki, and informed consent was signed by every participant.
Consent for publication
The authors have no conflicts of interest to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Figure S1: Permutation test in the second trimester group (a) and the third trimester group (b).
. Figure S2: Overlap statistic of different findings in two trimester groups.
. Figure S3: Longitudinal changes of three potential metabolic biomarkers. *P < 0.05.
. Table S1: List of DEMs in the second trimester group.
. Table S2: List of DEMs in the third trimester group.
. Table S3. Function enrichments of DEMs in the second trimester group.
. Table S4: Function enrichments of DEMs in the third trimester group.
. Table S5: Pathway enrichments of DEMs in the second trimester group.
. Table S6: Pathway enrichments of DEMs in the third trimester group.
. Table S7: Candidate metabolite biomarkers in the second- and third-trimester groups.
About this article
Cite this article
Lu, W., Luo, M., Fang, X. et al. Discovery of metabolic biomarkers for gestational diabetes mellitus in a Chinese population. Nutr Metab (Lond) 18, 79 (2021). https://doi.org/10.1186/s12986-021-00606-8
- Gestational diabetes mellitus
- Metabolic biomarkers