Combined Metabolomic Analysis of Plasma and Urine Reveals AHBA, Tryptophan and Serotonin Metabolism as Potential Risk Factors in Gestational Diabetes Mellitus (GDM)

Gestational diabetes mellitus during pregnancy has severe implications for the health of the mother and the fetus. Therefore, early prediction and an understanding of the physiology are an important part of prenatal care. Metabolite profiling is a long established method for the analysis and prediction of metabolic diseases. Here, we applied untargeted and targeted metabolomic protocols to analyze plasma and urine samples of pregnant women with and without GDM. Univariate and multivariate statistical analyses of metabolomic profiles revealed markers such as 2-hydroxybutanoic acid (AHBA), 3-hydroxybutanoic acid (BHBA), amino acids valine and alanine, the glucose-alanine-cycle, but also plant-derived compounds like sitosterin as different between control and GDM patients. PLS-DA and VIP analysis revealed tryptophan as a strong variable separating control and GDM. As tryptophan is biotransformed to serotonin we hypothesized whether serotonin metabolism might also be altered in GDM. To test this hypothesis we applied a method for the analysis of serotonin, metabolic intermediates and dopamine in urine by stable isotope dilution direct infusion electrospray ionization mass spectrometry (SID-MS). Indeed, serotonin and related metabolites differ significantly between control and GDM patients confirming the involvement of serotonin metabolism in GDM. Clustered correlation coefficient visualization of metabolite correlation networks revealed the different metabolic signatures between control and GDM patients. Eventually, the combination of selected blood plasma and urine sample metabolites improved the AUC prediction accuracy to 0.99. The detected GDM candidate biomarkers and the related systemic metabolic signatures are discussed in their pathophysiological context. Further studies with larger cohorts are necessary to underpin these observations.

Introduction

Gestational diabetes mellitus (GDM) is defined as glucose intolerance with onset or new recognition in pregnancy (Kautzky-Willer et al., 2004, 2016b). The prevalence of GDM is rising worldwide, reaching up to 25.1% (Zhu and Zhang, 2016). Therefore, GDM is the most common form of metabolic complication in pregnancy (Erem et al., 2015) and is associated with fetal macrosomia (Witkop et al., 2009), hyperbilirubinemia, and in consequence shoulder dystocia (Xiong et al., 2001; HAPO Study Cooperative Research Group, 2008; Yogev et al., 2009) as well as maternal morbidity (hypertension, polyhydramnion, infection). Women who had GDM have an elevated risk to develop Diabetes mellitus Type 2 (T2DM) or cardiovascular disease (Tobias et al., 2011), as well as hyperlipidemia or obesity (Bartha et al., 2008; Clausen et al., 2009; Landon et al., 2009; Gillman et al., 2010; Harreiter et al., 2014) in later life, whereas their children have a higher risk for obesity or impaired fasting glucose (Silverman et al., 1995). Further, GDM affects the psychological health of the mother and the child: GDM increases a woman's risk of postpartum depression 4-fold (Hinkle et al., 2016), and postpartum depression is decreased with treatment of GDM (Crowther et al., 2005; Beucher et al., 2010). Intrauterine exposure to hyperglycemia is linked to an increased risk for neuropsychiatric and neurodevelopmental disorders of the offspring (Xiang et al., 2015; Nahum Sacks et al., 2016).

Early detection of GDM and treatment can reduce the risk for mother and child. The current gold standard of diagnosing GDM is an oral glucose tolerance test (oGTT) between 24 and 28 weeks of gestation. However, detection of women at risk even earlier during pregnancy would be important to enable early lifestyle modification or even drug treatment in order to improve perinatal outcomes of these women. In addition other markers than glucose could be useful in identification of women and neonates at greatest risk. Therefore, the identification of new reliable and easily accessible biomarkers for earlier diagnosis of women with metabolic alterations during pregnancy would be of great value.

When the fetus starts to grow the maternal metabolism changes. More and more energy is required to ensure the growth of the unborn child. For this reason the maternal metabolic state must be modified in different and multiple ways to afford the energy demand. Among the altered mechanisms, an adaptation of hormones appears, such as insulin, serotonin (5-hydroxytryptamin, 5HT), hepatocyte growth factor (HGF) and cortisol (Ernst et al., 2011). Moreover transcription factors and cell cycle regulators cause a change in the metabolism of the mother during pregnancy (Ernst et al., 2011). These modulations are essential physiological elements for the normal progress of pregnancy, but can implicate pathological diabetic condition if they are disrupted. Although the altered mechanisms in GDM and overt diabetes are similar and pathways can be deduced, the conditions and the pathophysiology differ, and more research is needed on the most common but heterogeneous form, namely GDM (International Association of Diabetes in Pregnancy Study Group Working Group on Outcome et al., 2015; Kautzky-Willer et al., 2016a; Simmons et al., 2016; Rosta et al., 2017).

Metabolomics has the capacity to detect early deregulations and disruptions in metabolism associated with diseases or disorders. To investigate physiological processes and to develop (early) diagnostics, metabolomics is one of the most promising technologies (Bain et al., 2009; Pinto et al., 2015; Allalou et al., 2016). In case of glucose disrupted states (like GDM), there are novel findings regarding hitherto inconspicious hormones like melatonin or serotonin which arouse interest (Ernst et al., 2011) and legitimate closer examinations of the metabolome and hormone levels in GDM. Therefore we used untargeted and targeted metabolomic technology for the combined investigation of the metabolome in 32 pregnant women with and without GDM in blood plasma and urine and performed a metabolite profiling. Differences in the profiles between the case and control group were detected pointing to potential biomarkers and physiological processes for early GDM. Complementary to these metabolite markers the body mass index (BMI) and week of pregnancy (see Supplementary Table S1) were monitored. Combined regression analysis of the latter, blood plasma and urine metabolites improved the AUC (Area under curve) prediction accuracy to 0.99. In this paper we discuss corresponding hypotheses and assumptions for detected alterations in metabolic pathways for a better understanding of the metabolomic changes occuring in GDM.

Materials and Methods

Study Population

Participants were recruited from the outpatient clinic at the Medical University of Vienna. We investigated 14 women with GDM and 18 women without GDM (nGDM). The participants were 21–41 years old and in 12^th−26^th week of pregnancy (see Supplementary Table S1 also for BMI). The study was performed in accordance with the ethical principles of the Declaration of Helsinki II and was approved by the local ethics committee (Ethics Committee of the Medical University of Vienna). All participants gave their written informed consent.

GDM Definition, BMI and Sample Preparation

GDM was diagnosed with a standard 2 h 75 g oral glucose tolerance test (oGTT), according to the International Association of Diabetes and Pregnancy Study Groups (IADPSG) criteria (International Association of Pregnancy Study Groups Consensus Panel et al., 2010; Colagiuri et al., 2014). Eighteen control and fourteen GDM cases were defined. Urine samples were taken before the oGTT and urine was stored at −20°C. Blood samples were taken at three time points of the oGTT (0, 60, and 120 min after glucose intake). Plasma samples were prepared by centrifugation. Coagulation was avoided with an anticoagulant, ethylenediaminetetraacetic acid (EDTA). Samples were stored at −20°C until measurement. BMI was calculated as weight (kg) divided by the square of height (m²).

Plasma Sample Extraction, Derivatization, GC-MS Analysis, Identification and Quantification

Extraction and derivatization were performed slightly modified according to Weckwerth et al. (2004b). Chemicals were purchased from Sigma-Aldrich (Austria) if not stated otherwise: Methanol (CHROMASOLV^®, HPLC grade), chloroform (anhydrous ≥99%, stabilized with 0.5–1% ethanol). Water was double-distilled and deionized (Milli-Q water^® Advantage A10, Austria).

Samples were extracted in batches of 12, each including a pooled quality control sample and a blank processed in the same way. The pooled sample was a mixture of randomly chosen blood aliquots of half of the women, containing most of the metabolites expected in plasma (Dunn et al., 2011). Extraction was done by adding 200 μl of plasma to 1.7 ml of pre-chilled (−20°C) extraction mixture: methanol/chloroform/water (2.5:1:0.5 (v/v/v)) (MCW) and internal standard was spiked (1 μmol of D-Sorbitol-¹³C₆, 98 atom % ¹³C, Campro Scientific, Germany). Sample tubes were agitated for 10 s and incubated for 8 min on ice. Samples were centrifuged for 4 min at 14,000 g at 4°C. The one-phase-supernatant was collected in a new Eppendorf tube, shortly vortexed and divided into two equal aliquots, corresponding to 100 μl plasma each. Samples were dried using a speed vac (Savant^TM, Thermo Scientific, Austria) and stored at −80°C until measurements. One aliquot of each sample was used for further GC-MS analyses and metabolites were derivatized before measurement in two steps. Dried extracts were acclimated at room temperature for 10 min, dissolved in 20 μl solution of 40 mg/ml methoxyamine hydrochloride (CH3ONH2*HCL) in pyridine and incubated for 90 min at 30°C on a thermo shaker (550 rpm). Subsequently, 80 μl of N-methyl-N-(trimethylsilyl) trifluoroacetamide (MSTFA) (Macherey Nagel, Germany) were added and incubated for 30 min at 37°C. After centrifugation (2 min at 14,000 g), the supernatant was transferred to glass vials with micro inserts, closed with crimp caps and measured by gas chromatography coupled to mass spectrometry (GC-MS). GC-MS analyses and data validation were performed according to previous publications (Mari et al., 2013; Prezelj et al., 2016) with slight modifications on a ThermoFisher Trace GC coupled to a Triple Quadrupole mass analyzer (Thermo Scientific TSQ Quantum GC™, Bremen, Germany). Each batch included randomly chosen samples, a blank, a pooled quality control sample and a pure non derivatized alkane standard mixture of even-numbered n-alkanes (C10-C40, each 50 mg/L in hexane) for retention index (RI) determination. One microliter of sample was injected in splitless mode at a constant injector temperature of 230°C using a deactivated stainless steel Siltek liner (Restek Corp., USA). GC separation was performed on a HP-5MS capillary column (30 m × 0.25 mm × 0.25 μm) (Agilent Technologies, CA) at a constant helium flow rate of 1 mL min⁻¹. Initial oven temperature was set to 70°C and held for 1 min, followed by a ramp to 76°C at 1°C min⁻¹ and a second ramp at 6°C min⁻¹ to 350°C held for 1 min. Transfer line temperature was set to 340°C and post run temperature to 325°C for 10 min. The quadrupole mass analyzer was used in full scan mode with a scan range of m/z 40–600 Th and a scan time of 250 msec. Electron impact (EI) ionization was performed at 70 eV with 50 μA emission current and ion source temperature was set to 250°C. Metabolite derivatives were identified by matching retention time as well as mass spectra with those of reference standards and by comparison of alkane based retention indices with an in house mass spectral library, as well as the GMD library (Kopka et al., 2005). Metabolites were considered as annotated with a spectral match factor higher than 850 (NIST MS Search 2.0 Program algorithm) and RI-deviation lower than 4%. Deconvolution and RI-deviation calculation was performed with AMDIS (Stein, 1999) and quantification with LC-Quan 2.6.0 (Thermo Fisher Scientific Inc.). Peak areas of a specific ion of a compound (quant m/z) were normalized to the ¹³C-sorbitol peak within each run. More detailed information and a list of quantified analytes can be found in Supplementary Table S1.

Stable Isotope Diluted Direct Infusion Electrospray Ionization Mass Spectrometry (SID-MS) Analysis of Serotonin Metabolism in Urine

If not stated otherwise solvents were purchased from Sigma-Aldrich (Austria) in high quality (CHROMASOLV^®, HPLC grade), formic acid (ROTIPURAN^®) was purchased from Carl Roth (Germany), acetic acid from Fisher Scientific (Austria) and water was double-distilled Milli-Q water as stated above.

Urine samples were purified prior analyses by solid phase extraction (SPE) using a slightly modified protocol according to Moriarty et al. (2011). Reaction tubes (15 ml, Greiner bio-one) were prepared by cutting a hole into the lid in the size of a SPE C18 cartridge (InertSep^TM 100 mg/1 ml, GL Sciences, Japan). SPE cartridges were inserted into the tubes and conditioned with 3 × 1 ml methanol, followed by 3 × 1 ml of acidified water to ph 3.5 with acetic acid. Before loading, urine samples were diluted by adding 500 μl sample to 800 μL acidified water (ph 3.5) in a 2 ml reaction tube including internal standard serotonin-d₄ (1 μMol) (98 atom% D, CDN Isotopes, Canada). After vortexing and centrifugation for 5 min at 14,000 g the supernatant was loaded on the SPE cartridge and washed with 1 ml 5% (v/v) methanol solution. Metabolites were eluted with 5 × 1 ml 0.1 M ammonium acetate in methanol. SPE cartridge tubes were centrifuged after each solvent or sample loading step at 1,000 g for 30 sec at 4°C. The collected eluate was dried for 2 h to absolute dryness under a gentle stream of nitrogen (N₂ 99.999%) and stored at −20°C.

Purified and dried urine samples were dissolved in 200 μL 0.1% formic acid in methanol, centrifuged at 21,000 g and supernatant were transferred to glass vials with micro-inserts which were crimped with pre-perforated lids. Samples were kept at 4°C during all preparation steps of extraction, purification and analyses.

Targeted SID-MS detection and quantification of serotonin-melatonin-tryptophan pathway metabolites were performed on an Orbitrap Elite mass spectrometer (Thermo Fisher Scientific, USA). Direct infusion of samples was performed using a nano UHPLC pump equipped with an autosampler (Dionex UltiMate 3000 RSLCnano UHPLC pump, Thermo Fisher Scientific) by an isocratic flow. To prevent clogging of the nano spray needle, a PicoChip^TM nano emitter system for Infusion (New objective Inc., USA) was used and ionization was conducted using a nano spray ionisation source (NSI). In the following, instrument descriptions and detailed parameters are given. Injection volume was 5 μl, isocratic flow was performed at a flow rate of 500 nL/min, with 60% mobile phase A: 0.1% formic acid and 40% mobile phase B (90% ACN, 10% H₂O, +0.1% formic acid) with a total run time of 15 min. Tip size of the PicoChip^TM nano emitter was 15 μm and NSI source parameters were as follows: source voltage 1.9 kV, source current 100 μA, sheath gas 0, aux gas 0 and capillary temperature 275°C.

Accurate mass analysis was performed using the Orbitrap FTMS mass analyzer using the lock mass option in MS and MS/MS mode. Ions of cyclomethicone N5 (m/z = 371.101230) were used for internal mass calibration. Data dependent MS² scan experiments of a target parent mass list was performed at a mass resolution of 120,000; scan event 1 was performed in full scan mode and a scan range of m/z 110–600 Th, MS₂ fragmentation of triggered precursor masses in scan event 2 was performed by collision induced dissociation (CID) with 50 eV normalized collision energy. Precursor ions of the target analytes serotonin, 5-hydroxyindolic acetic acid, N-acetylserotonin, 5-methoxytryptamin, melatonin, 6-hydroxymelatonin, L-tryptophan, 5-hydroxytryptophan and dopamine, as well as observed MS² fragments are summarized in Table 1. Mass calibration was performed once a week. For quantification average intensities of MS2 product ions were used and normalized to the stable isotope labeled internal standard serotonin-d₄.

TABLE 1

Table 1. Metabolite and mass spectral information of targeted SID-MS analysis.

Statistical Analysis

Log transformation and data normalization using internal standards to avoid interday tuning differences were applied. Two tailed unpaired t-test was used to compare the differences in metabolites between the control and the GDM-group. A p-value of less than 0.05 was regarded to be statistically significant. Further, the data were analysed with principal component analysis (PCA) and ANOVA integrated in the statistical toolbox COVAIN (Sun and Weckwerth, 2012) to find group separations and significant metabolite changes between control and GDM. Additionally, the data matrices are analyzed with partial least squares discriminant analysis (PLS-DA) using SIMCA-PR (Umetrics, Sweden). The PLS-DA model employed 7 cross validation groups, assigning every 7th observation to the same group, and grouping similar observations in the same groups. In Supplementary Figure S1 Q2(cum) and R2(cum) are given as indicators for goodness of fit and the predictive quality of all three time points. Variable importance selection was conducted with a VIP-analysis (variable importance in the projection).

Clustered correlation coefficient matrix visualization was performed with an in house written Matlab script (Matlab R2016b, version 9.1, ^® Natick, Massachusetts, United States). The script is available upon request from the corresponding author.

BMI, week of pregnancy and blood and urine sample metabolites were analyzed by LASSO regression (least absolute shrinkage and selection operator). For 21 patients (10 GDM, 11 normal) with serotonin and intermediates measurements, we built statistical models to predict GDM from a combination of metabolomics (131 metabolites) as well as serotonin data (8 compounds). First, metabolomics data were log-transformed. Second, all variables were standardized by the z-score transformation. Third, LASSO regression combined with a linear SVM (support vector machine) classifier, was applied to select the best subset of variables that achieve highest prediction accuracy. LASSO is a regularization regression method that penalizes the absolute size of the regression coefficients to avoid overfitting that is common in many regression problems. Here, LASSO was applied with a range of lambda values, or the regularization parameter. For each lambda value, a subset of variables are selected. Then a linear SVM model was built using these variables for the two-class (GDM = 1, nGDM = 0) classification problem. The concept of SVM is to maximize the margin that distinguish groups of data by fitting a “hyperplane.” It proves to work well with biological data where the sample size is small. The SVM model was evaluated by five-fold cross validation (CV). After trying all lambda values, the best subset of variables were selected from the model with minimal CV error. All the analysis was done by the Statistics and Machine Learning Toolbox in Matlab R2016b (version 9.1, ^® Natick, Massachusetts, United States).

Results

Plasma Metabolite Analysis of Control and GDM Patients during an Oral Glucose Tolerance Test (oGTT)

We applied an universal integrative protocol for extraction and analysis of metabolites from plasma samples (Weckwerth et al., 2004b). After extraction with MCW (methanol, chloroform, water in a one phase mixture) no further phase separation was performed and the complete MCW extract was injected into GC-MS. With the modified protocol 131 annotated and putative metabolites were detected (see Supplementary Table S1). The metabolite data were further analyzed by PLS-DA (see Figure 1). Healthy pregnant women are labeled by black squares (class 1) and diabetic individuals are labeled with red dots (class 2). PLS-DA revealed a separation between case and control groups in all three time points of the oGTT (Figures 1A–C). To detect the metabolites, which are most responsible for the separation in the PLS-scatter plots a VIP analysis has been performed. The pattern of the corresponding VIP analysis (Figure 1D) illustrates the metabolites importantance projection for the sample separation seen in the PLS-DA plot in Figure 1B. VIP values larger than one indicate important metabolites, which were grouped into classes. The VIP plot (Figure 1D), sorted by importance, showed the metabolites tryptophan and 3-hydroxybutanoic acid as most discriminatory between GDM and control. Further interesting metabolite markers discriminating GDM and control are discussed below.

FIGURE 1

Figure 1. PLS-DA model built for two classes, control and GDM for three different time points of the oral glucose tolerance test (oGTT). Healthy pregnant women are labeled by black squares (class 1) and GDM individuals are labeled with red dots (class 2). (A) 0 h, (B) 1 h, (C) 2 h (D) VIP projection of the variables according to the scores plot of B.

Amino Acids, Fatty Acids, Organic Acids, Sugars, and Steroids Discriminate Control and GDM

In Supplementary Figure S1 selected boxplots of statistically significant organic acids, branched-chain amino acids (BCAA), and further amino acids are shown. Valine, alanine and β-alanine differ significantly between patients with and without GDM. Several fatty acids, hydroxy acids and other organic acids as well as ketone bodies and sugars show differences. 2- and 3-hydroxybutanoic acid (α- and β-hydroxybutyric acid, AHBA and BHBA) differ significantly (see boxplots Supplementary Figure S1).

The intermediates of the tricarboxylic acid cycle (TCA) show significant differences between the GDM- group and the control-group. The potentially identified plant sterol β-sitosterol, derived from diet, and cholesterol were found significantly changed. β-Sitosterol had also high loadings in the VIP analysis (see Figure 1D).

Serotonin Metabolism Is Changed in GDM

The strong VIP loadings of tryptophan pointed us to serotonin metabolism. Because serotonin and intermediates are only weakly covered by GC-MS we applied a stable isotope dilution direct-infusion method (SIDE-MS assay, see section Materials and Methods). We applied this assay to the analysis of urine samples from the GDM and the control group. Urine analysis is a non-invasive technique complementing any other profiling method with lowest costs. In Figure 2A a principal component analysis demonstrates that serotonin metabolism is altered in GDM patients in comparison to the control-group. In Figure 2B boxplots are shown of selected intermediates with high loadings.

FIGURE 2

Figure 2. Measurement of extracellular serotonin metabolic intermediates and dopamine in urine of GDM patients. (A) Principal component analysis of serotonin metabolic intermediates. (B) Boxplots of serotonin metabolic intermediates in urine of GDM vs. control groups.

Integrative Analysis of Plasma and Urine Metabolomics

Integrative analysis of plasma and urine metabolite profiles was performed in two steps: (i) a systemic analysis using metabolite correlation networks of control and GDM cases, and (ii) using a linear regression method to determine the best set of candidate biomarker for accurate prediction of GDM cases.

In Figures 3A,B clustering of correlation coefficients between all metabolites from plasma and from urine samples is shown. The clustered correlation map is very different between control and GDM cases indicating a large reprogramming of metabolism in GDM patients. This pattern can be interpreted as a signature of GDM. Details are discussed below and in the discussion. In Figures 3C,D detailed clustered correlation coefficient maps of the serotonin-related metabolites in urine samples are depicted. Here, a remarkable difference is the clustering of serotonin, dopamine, 5-HIAA, N-acetyl-serotonin and tryptophan (cluster 1) and 5-methoxytryptamin, melatonin and 6-hydroxy-melatonin (cluster 2) in the control samples whereas cluster 1 is completely dissolved and cluster 2 is maintained in GDM cases (see Figures 3C,D).

FIGURE 3

Figure 3. Metabolite correlation network analysis visualized as clustered correlation coefficient matrix. (A) Clustered heat map shows a significant metabolic signature of control samples. (B) GDM case show a different pattern compared to the control samples indicating a dramatic reprogramming of metabolism in GDM disease. (C) Detailed view of serotonin/melatonin metabolites in urine samples of control. Cluster 1 is a highly correlated serotonin cluster and cluster 2 a highly correlated melatonin cluster in normal metabolic conditions. (D) Detailed view of serotonin/melatonin metabolites in urine samples of GDM case. Cluster 1 disappeared in GDM cases and cluster 2 is conserved from control to GDM.

To select for the set of most discriminatory variables from blood and urine sample metabolites a LASSO (least absolute shrinkage and selection operator) regression method was employed. For the 21 patients (10 GDM, 11 nGDM) with serotonin and metabolic intermediate measurements, we built statistical models to predict GDM from a combination of 131 plasma metabolites as well as 8 urine metabolites related to serotonin metabolism. LASSO regression wrapped with a linear SVM (support vector machine) classifier, was applied to select the best subset of variables that achieves the highest prediction accuracy. We found that the best subset of metabolites plus BMI achieves a good prediction performance with an AUC (Area under curve) value 0.94 of the ROC (Receiver operating characteristic) curve. These metabolites are glycolic acid, urea, methionine, erythronic acid, an unknown organic acid (potentially tartaric acid), glutamine, an unknown carbohydrate, unknown, alanin, serin and tryptophan. Next, we tested whether the prediction accuracy can be improved by combining the serotonin data from urine analysis. Following the same procedure, the best model adds five serotonin compounds, namely, serotonin, 5-HIAA, L-tryptophan, melatonin and 6-hydroxymelatonin, and achieves a higher AUC value of 0.99. The comparison between these two models is shown in Figure 4A.

FIGURE 4

Figure 4. Receiver operating characteristic (ROC) analysis of selected metabolites from plasma and urine samples of GDM vs. control patients. (A) Area under curve (AUC) for selected metabolites from plasma analysis and combined analysis of selected analysis of plasma and urine analysis. The final AUC is 0.99. Selection of variables was performed by LASSO regression (see section Materials and Methods). 11 Metabolites from blood samples and 5 metabolites from urine samples were selected (for further information see results and discussion). (B) Correlation analysis of serotonin, associated metabolic intermediates and dopamine measured in urine with BMI. The strongest correlation is found for dopamine.

Correlation between BMI and Metabolites of Serotonin Metabolism in Urine Samples

Pearson correlation coefficients (PCC) were calculated between 8 urine target compounds and BMI (Figure 4B). Significant associations were found between BMI and dopamine (PCC = 0.50, p-value = 0.022), and 5-HIAA (PCC = 0.49, p-value = 0.023), respectively.

Discussion

Changes in Serotonin Metabolism in GDM vs. Control Group

L-tryptophan and its metabolites serotonin (5HT) and melatonin have notable functions in regulating growth and development of the fetus (Glover, 2015; St-Pierre et al., 2016; Wu et al., 2016) and are involved in a magnifique number of physiological pathways and adaption processes during pregnancy (Sano et al., 2016). As L-trytophan (TRP) is highest ranked in the PLS-analysis (Figure 1D), we searched for differences regarding serotonin metabolism. Therefore, we applied a serotonin assay by using stable isotope dilution direct electrospray ionization mass spectrometry (SID-MS, see section Materials and Methods) to analyze urine samples of the same patients. In the measurements the control and the GDM patients showed significantly different concentrations of serotonin and intermediates of this pathway in urine samples (see Figure 2): serotonin levels were higher in women with GDM compared with women without GDM. Besides the well-known role of serotonin in mood and feeding behavior, recent results renewed the interest in the role of serotonin in metabolic diseases (Almaça et al., 2016). Conditions of metabolic demand lead to the production of serotonin in the pancreatic islet (Kim et al., 2010; Goyvaerts et al., 2015, 2016). Further, studies showed that serotonin via autocrine signaling increases beta- cell function and ß-cell mass during insulin resistant states (e.g., pregnancy). In mouse models evidence suggests that increased serotonin in islets during pregnancy drives ß-cell expansion (Kim et al., 2010). Additionally, serotonin levels were shown to rise in regard to BMI (Almaça et al., 2016). A recent study showed that ß-cell-derived serotonin inhibits glucagon secretion, a hormone which elevates the blood glucose levels. Based on this, our results could be explained by (1) a different lifestyle (nutrition: protein intake) of the patients with GDM and those without GDM (2) a higher BMI of women with GDM compared to the control group could be a reason either, or (3) higher serotonin levels in the case of disrupted metabolism, like GDM, could substantiate that serotonin is a factor in maintaining normoglycemia, the increase presenting a possible compensatory mechanism.

In summary, changes in L-tryptophan in the GDM group pointed to an altered serotonin metabolism which indeed is changed in GDM patients. Accordingly, alterations in serotonin metabolism may be fundamental in the pathogenesis of GDM. Nevertheless, such findings need to be confirmed by further studies analyzing larger cohorts.

Glucose-L-Alanine Cycle

Some of the significantly changed metabolites are found in the glucose-L-alanine cycle, which plays a role in glycolysis and gluconeogenesis. These findings appear to be relevant for the energy metabolism in GDM. The glucose-alanine cycle is a biochemical route between muscle and liver metabolism. In the muscle, protein is degraded resulting in glutamic acid. Subsequently, the alanine aminotransferase transfers the amino group from glutamate to pyruvate, which arises from glycolysis. L-alanine is formed and then transferred from muscle to the liver through the blood. Finally, the amino acid is used in the liver for gluconeogenesis. Several other metabolites show an alteration in concentrations from the GDM group (see Supplementary Figure S1). Furthermore, we analysed also the alanine-amino-transferase activity and identified a significant change between the case and control groups confirming the metabolomics data (data not shown). A recent study suggests an up-regulated glycolysis and higher alanine concentration in T2DM patients (Huang and Joseph, 2012). Like in our study another recent study discovered increased concentrations of alanine and lactate as potential gluconeogenic precursors (Galazis et al., 2012) as well as alanine and glutamate as significant marker in GDM (Bentley-Lewis et al., 2015).

Strong Markers of GDM Disease −2- and 3-Hydroxybutanoic Acid

Ketone bodies play a role in the fatty acid biosynthesis. 3-hydroxybutanoic acid (β-hydroxybutyrate; BHBA) is an organic acid, which is used for the biosynthesis of fatty acids. BHBA reached statistical significance between the case and control group. Furthermore, 2-hydroxybutanoic acid (α-hydroxybutyrate; AHBA) shows the lowest p-value of all metabolites between control and GDM and therefore is one of the strongest metabolic alterations (Supplementary Figure S1).

Our results confirm already well documented conclusions: AHBA is suggested as an early biomarker for insulin resistance and proposed to identify insulin resistance earlier than current diagnostic tests (Gall et al., 2010). Previous research has studied liver and plasma metabolome in mice to identify early alterations in development of insulin resistance and found significant changes in metabolites, of which AHBA is also significant in our study (Li et al., 2009). Altogether AHBA tends to be important for decreased insulin sensitivity in GDM and may be critical in the prevention and treatment of diabetes as potent biomarker.

Propanoate Metabolism

In consideration of four metabolites—lactic acid, AHBA, beta-alanin, and valine—of the propanoate pathway which show a significant change in case and control (see Supplementary Table S1 and Supplementary Figure S1) it seems reasonable to propose that the whole propanoate pathway is changed in the case group. On the basis of our data we suggest an adaption in the propanoate metabolism as described in an in-vivo study for DM by Huang in 2006 (Huang et al., 2006). Literature-derived evidence comprises a lipidomic analysis which found statistically significant pathways, among them the propanoate pathway (Zhao et al., 2013).

Degradation of Valine and Fatty Acids

Among the amino acids especially valine is one with the highest significance. Further, intermediates of the citric cycle are significantly changed in GDM (see Supplementary Figure S1). It seems very likely that the rate of degradation of amino acids is varying between case and control group. From protein breakdown L-valine is generated and finally succinyl-CoA appears to contribute to the citric acid cycle in muscle. In liver succinyl-CoA can be used for synthesis of glucose. Therefore, there is considerable evidence of an alteration in gluconeogenesis rates in the pathogenesis of GDM. 2,4-dihydroxybutyrate is described in the HMDB as usually absent in normal human urine extracts or present only in trace amounts in neonates (Wishart et al., 2013). However, in our study this metabolite was identified as being significantly changed between case and control group in blood samples.

A metabolic Signature of GDM by Combined Plasma and Urine Metabolomics Analysis

The analysis of metabolite correlation networks is a powerful method for the description of systemic biochemical regulation (Steuer et al., 2003; Weckwerth, 2003, 2011; Weckwerth et al., 2004a; Weckwerth and Morgenthal, 2005; Nägele et al., 2014). Recently, we demonstrated that differential metabolite correlation or covariance networks reflect biochemical regulation depending on the genotype or other factors (Weckwerth et al., 2004a; Sun and Weckwerth, 2012; Nägele et al., 2014). To reveal a picture of the metabolite correlation network in control vs. GDM case samples we used here a clustered correlation coefficient matrix visualization (Figure 3). Based on this visualization clearly distinguishable metabolite signatures differentiate control and GDM cases (Figures 3A,B). In a more detailed analysis of serotonin metabolites in urine samples a highly significant reprogramming of serotonin/melatonin metabolism is visible. Whereas serotonin and tryptophan and related metabolites form a highly correlated cluster (cluster 1 in Figure 3C) this cluster is compromised in GDM metabolism. In contrast, the melatonin cluster (cluster 2 in Figures 3C,D) is conserved. There are several reports that serotonin metabolism and especially serotonin transport is impaired in GDM (Li et al., 2014). Subsequently this would lead to changes in biochemical correlation networks as shown in our study. Further studies are necessary to underline the dynamics of these different metabolic signatures and especially reveal causal relationships between changed metabolite levels and enzymatic regulation (Sun and Weckwerth, 2012; Nägele et al., 2014).

To select for the best set of predictive candidate biomarker we used LASSO regression (see sections Materials and Methods and Results). The combined analysis of plasma metabolomics and targeted metabolomics of serotonin metabolism in urine revealed an improvement in prediction accuracy for GDM (see Figure 4A). Recently, a study of type 2 diabetes (T2DM) showed no improvement of prediction accuracy by integrating NMR metabolite profiling of blood and urine samples (Friedrich et al., 2015). In contrast, our targeted approach of profiling specific serotonin-related metabolites in urine indeed improved the prediction accuracy in GDM. Thereby the specificity of serotonin metabolism in GDM is supported.

Integrative regression analysis revealed further interesting metabolite markers. Glycolate was also recently detected as a T2DM metabolite marker in woman (Friedrich et al., 2015). In another study on T1D in rats urea was found to be a significant marker in plasma samples (Zhang et al., 2008). Several amino acids are part of the identified metabolic signature for GDM such as methionine, glutamine, alanine, serine, and tryptophan. Tryptophan is directly related to serotonin metabolism. Consequently, the extension of the metabolic signature by metabolites derived from the serotonin pathway improved the prediction accuracy (Figure 3A). Another very interesting metabolic marker is erythronic acid. This acid is a side-product from the degradation of so called advanced glycation end-products (AGE) such as fructosamine (Jakus and Rietbrock, 2004). AGE's play an important role in the pathophysiology of diabetes. We further correlated serotonin-related metabolites and others measured in the urine samples with the BMI of the cohort and found a strong correlation with dopamine. An intimate relationship between dopamine and obesity is postulated in the “reward deficiency syndrome” (Blum et al., 2014). We will investigate this relationship in more detail in future studies.

Conclusion

Plasma metabolomics is a comprehensive technology for the rapid and in depth analysis of all kinds of metabolic diseases. The amount of information by this approach is still not fully understood. Here, we applied this technology for the early diagnosis of GDM. Based on the physiological interpretation of identified metabolic markers we developed the hypothesis that serotonin is involved in GDM pathogenesis. We tested this hypothesis by targeted profiling of serotonin-derived metabolites in urine samples and the integration of plasma and urine metabolic markers improved the prediction accuracy of GDM in our study. This workflow reveals the power of metabolomics screening of metabolic diseases, especially in the context of physiological interpretation of identified significant changes in metabolites (McCabe and Perng, 2017). In future studies, we will test our hypothesis in more detail and will extend the study to larger cohorts to unambiguously validate or falsify the potential metabolic marker from our study. Furthermore, our study lays the ground to investigate GDM in more detail on a biochemical and physiological basis.

Author Contributions

WW and AK-W conceived the study. SD, LF, ML, NH, ST, and HD performed the experiments. SD, LF, NH, HD, XS, GB, and WW analyzed the data and performed statistics. WW, WJ, KL, and AK-W provided the reagents, materials and analytical tools. WW and ML wrote the manuscript. All authors revised the manuscript and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank the members of the Department Ecogenomics and Systems Biology at the University of Vienna for critical and helpful discussions. We acknowledge the strong support of the University of Vienna for the Vienna Metabolomics Center.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2017.00084/full#supplementary-material

References

Allalou, A., Nalla, A., Prentice, K. J., Liu, Y., Zhang, M., Dai, F. F., et al. (2016). A Predictive metabolic signature for the transition from gestational diabetes mellitus to Type 2 diabetes. Diabetes 65, 2529–2539. doi: 10.2337/db15-1720

PubMed Abstract | CrossRef Full Text | Google Scholar

Almaça, J., Molina, J., Menegaz, D., Pronin, A. N., Tamayo, A., Slepak, V., et al. (2016). Human beta cells produce and release serotonin to inhibit glucagon secretion from alpha cells. Cell Rep 17, 3281–3291. doi: 10.1016/j.celrep.2016.11.072

PubMed Abstract | CrossRef Full Text | Google Scholar

Bain, J. R., Stevens, R. D., Wenner, B. R., Ilkayeva, O., Muoio, D. M., and Newgard, C. B. (2009). Metabolomics applied to diabetes research: moving from information to knowledge. Diabetes 58, 2429–2443. doi: 10.2337/db09-0580

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartha, J. L., González-Bugatto, F., Fernández-Macías, R., González-González, N. L., Comino-Delgado, R., and Hervías-Vivancos, B. (2008). Metabolic syndrome in normal and complicated pregnancies. Eur. J. Obstet. Gynecol. Reprod. Biol. 137, 178–184. doi: 10.1016/j.ejogrb.2007.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bentley-Lewis, R., Huynh, J., Xiong, G., Lee, H., Wenger, J., Clish, C., et al. (2015). Metabolomic profiling in the prediction of gestational diabetes mellitus. Diabetologia 58, 1329–1332. doi: 10.1007/s00125-015-3553-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Beucher, G., Viaris De Lesegno, B., and Dreyfus, M. (2010). Maternal outcome of gestational diabetes mellitus. Diab. Metab. 36, 522–537. doi: 10.1016/j.diabet.2010.11.006

PubMed Abstract | CrossRef Full Text

Blum, K., Thanos, P. K., and Gold, M. S. (2014). Dopamine and glucose, obesity, and reward deficiency syndrome. Front. Psychol. 5:919. doi: 10.3389/fpsyg.2014.00919

PubMed Abstract | CrossRef Full Text | Google Scholar

Clausen, T. D., Mathiesen, E. R., Hansen, T., Pedersen, O., Jensen, D. M., Lauenborg, J., et al. (2009). Overweight and the metabolic syndrome in adult offspring of women with diet-treated gestational diabetes mellitus or type 1 diabetes. J. Clin. Endocrinol. Metab. 94, 2464–2470. doi: 10.1210/jc.2009-0305

PubMed Abstract | CrossRef Full Text | Google Scholar

Colagiuri, S., Falavigna, M., Agarwal, M. M., Boulvain, M., Coetzee, E., Hod, M., et al. (2014). Strategies for implementing the WHO diagnostic criteria and classification of hyperglycaemia first detected in pregnancy. Diab. Res. Clin. Pract. 103, 364–372. doi: 10.1016/j.diabres.2014.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Crowther, C. A., Hiller, J. E., Moss, J. R., McPhee, A. J., Jeffries, W. S., Robinson, J. S., et al. (2005). Effect of treatment of gestational diabetes mellitus on pregnancy outcomes. N. Engl. J. Med. 352, 2477–2486. doi: 10.1056/NEJMoa042973

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunn, W. B., Broadhurst, D., Begley, P., Zelena, E., Francis-McIntyre, S., Anderson, N., et al. (2011). Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat. Protocols 6, 1060–1083. doi: 10.1038/nprot.2011.335

PubMed Abstract | CrossRef Full Text | Google Scholar

Erem, C., Kuzu, U. B., Deger, O., and Can, G. (2015). Prevalence of gestational diabetes mellitus and associated risk factors in Turkish women: the Trabzon GDM Study. Arch. Med. Sci. 11, 724–735. doi: 10.5114/aoms.2015.53291

PubMed Abstract | CrossRef Full Text | Google Scholar

Ernst, S., Demirci, C., Valle, S., Velazquez-Garcia, S., and Garcia-Ocana, A. (2011). Mechanisms in the adaptation of maternal beta-cells during pregnancy. Diab. Manag. (Lond) 1, 239–248. doi: 10.2217/dmt.10.24

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedrich, N., Budde, K., Suhre, K., Volker, U., John, U., Felix, S. B., et al. (2015). Sex differences in urine metabolites related with risk of diabetes using NMR spectroscopy: results of the study of health in pomerania. Metabolomics 11, 1405–1415. doi: 10.1007/s11306-015-0795-6

CrossRef Full Text | Google Scholar

Galazis, N., Iacovou, C., Haoula, Z., and Atiomo, W. (2012). Metabolomic biomarkers of impaired glucose tolerance and type 2 diabetes mellitus with a potential for risk stratification in women with polycystic ovary syndrome. Euro. J. Obstet. Gynecol. Reprod. Biol. 160, 121–130. doi: 10.1016/j.ejogrb.2011.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Gall, W. E., Beebe, K., Lawton, K. A., Adam, K.-P., Mitchell, M. W., Nakhle, P. J., et al. (2010). α-hydroxybutyrate is an early biomarker of insulin resistance and glucose intolerance in a nondiabetic population. PLoS ONE 5:e10883. doi: 10.1371/journal.pone.0010883

PubMed Abstract | CrossRef Full Text | Google Scholar

Glover, V. (2015). Prenatal stress and its effects on the fetus and the child: possible underlying biological mechanisms. Adv. Neurobiol. 10, 269–283. doi: 10.1007/978-1-4939-1372-5_13

PubMed Abstract | CrossRef Full Text | Google Scholar

Gillman, M. W., Oakey, H., Baghurst, P. A., Volkmer, R. E., Robinson, J. S., and Crowther, C. A. (2010). Effect of treatment of gestational diabetes mellitus on obesity in the next generation. Diab. Care 33, 964–968. doi: 10.2337/dc09-1810

PubMed Abstract | CrossRef Full Text | Google Scholar

Goyvaerts, L., Lemaire, K., Arijs, I., Auffret, J., Granvik, M., Van Lommel, L., et al. (2015). Prolactin receptors and placental lactogen drive male mouse pancreatic islets to pregnancy-related mRNA changes. PLoS ONE 10:e0121868. doi: 10.1371/journal.pone.0121868

PubMed Abstract | CrossRef Full Text | Google Scholar

Goyvaerts, L., Schraenen, A., and Schuit, F. (2016). Serotonin competence of mouse beta cells during pregnancy. Diabetologia 59, 1356–1363. doi: 10.1007/s00125-016-3951-2

PubMed Abstract | CrossRef Full Text | Google Scholar

HAPO Study Cooperative Research Group (2008). Hyperglycemia and adverse pregnancy outcomes. N. Engl. J. Med. 358, 1991–2002. doi: 10.1056/NEJMoa0707943

PubMed Abstract | CrossRef Full Text

Harreiter, J., Dovjak, G., and Kautzky-Willer, A. (2014). Gestational diabetes mellitus and cardiovascular risk after pregnancy. Womens Health (Lond) 10, 91–108. doi: 10.2217/WHE.13.69

PubMed Abstract | CrossRef Full Text | Google Scholar

Hinkle, S. N., Buck Louis, G. M., Rawal, S., Zhu, Y., Albert, P. S., and Zhang, C. (2016). A longitudinal study of depression and gestational diabetes in pregnancy and the postpartum period. Diabetologia 59, 2594–2602. doi: 10.1007/s00125-016-4086-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Kim, Y., Caramori, M. L., Moore, J. H., Rich, S. S., Mychaleckyj, J. C., et al. (2006). Diabetic nephropathy is associated with gene expression levels of oxidative phosphorylation and related pathways. Diabetes 55, 1826–1831. doi: 10.2337/db05-1438

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, M., and Joseph, J. W. (2012). Metabolomic analysis of pancreatic beta-cell insulin release in response to glucose. Islets 4, 210–222. doi: 10.4161/isl.20141

PubMed Abstract | CrossRef Full Text | Google Scholar

International Association of Diabetes in Pregnancy Study Group Working Group on Outcome, D., Feig, D. S., Corcoy, R., Jensen, D. M., Kautzky-Willer, A., Nolan, C. J., et al. (2015). Diabetes in pregnancy outcomes: a systematic review and proposed codification of definitions. Diab. Metab. Res. Rev. 31, 680–690.

Google Scholar

International Association of Pregnancy Study Groups Consensus Panel, Metzger, B. E., Gabbe, S. G., Persson, B., and Buchanan, T. A. (2010). International association of diabetes and pregnancy study groups recommendations on the diagnosis and classification of hyperglycemia in pregnancy. Diabetes Care 33, 676–682. doi: 10.2337/dc09-1848

PubMed Abstract | CrossRef Full Text

Jakus, V., and Rietbrock, N. (2004). Advanced glycation end-products and the progress of diabetic vascular complications. Physiol. Res. 53, 131–142.

PubMed Abstract | Google Scholar

Kautzky-Willer, A., Bancher-Todesca, D., and Birnbacher, R. (2004). [Gestational diabetes mellitus]. Acta Med. Austriaca 31, 182–184.

PubMed Abstract | Google Scholar

Kautzky-Willer, A., Harreiter, J., Bancher-Todesca, D., Berger, A., Repa, A., Lechleitner, M., et al. (2016a). [Gestational diabetes mellitus]. Wien. Klin. Wochenschr. 128(Suppl. 2), S103–S112.

Kautzky-Willer, A., Harreiter, J., and Pacini, G. (2016b). Sex and gender differences in risk, pathophysiology and complications of Type 2 diabetes mellitus. Endocr. Rev. 37, 278–316.

PubMed Abstract | Google Scholar

Kim, H., Toyofuku, Y., Lynn, F. C., Chak, E., Uchida, T., Mizukami, H., et al. (2010). Serotonin regulates pancreatic beta cell mass during pregnancy. Nat. Med. 16, 804–808. doi: 10.1038/nm.2173

PubMed Abstract | CrossRef Full Text | Google Scholar

Kopka, J., Schauer, N., Krueger, S., Birkemeyer, C., Usadel, B., Bergmuller, E., et al. (2005). GMD@CSB.DB: the golm metabolome database. Bioinformatics 21, 1635–1638. doi: 10.1093/bioinformatics/bti236

PubMed Abstract | CrossRef Full Text | Google Scholar

Landon, M. B., Spong, C. Y., Thom, E., Carpenter, M. W., Ramin, S. M., Casey, B., et al. (2009). A multicenter, randomized trial of treatment for mild gestational diabetes. N. Engl. J. Med. 361, 1339–1348. doi: 10.1056/NEJMoa0902430

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Xu, Z., Lu, X., Yang, X., Yin, P., Kong, H., et al. (2009). Comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry for metabonomics: biomarker discovery for diabetes mellitus. Anal. Chim. Acta 633, 257–262. doi: 10.1016/j.aca.2008.11.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Hadden, C., Singh, P., Mercado, C. P., Murphy, P., Dajani, N. K., et al. (2014). GDM-associated insulin deficiency hinders the dissociation of SERT from ERp44 and down-regulates placental 5-HT uptake. Proc. Natl. Acad. Sci. U.S.A. 111, E5697–E5705. doi: 10.1073/pnas.1416675112

PubMed Abstract | CrossRef Full Text | Google Scholar

Mari, A., Lyon, D., Fragner, L., Montoro, P., Piacente, S., Wienkoop, S., et al. (2013). Phytochemical composition of Potentilla anserina L. analyzed by an integrative GC-MS and LC-MS metabolomics platform. Metabolomics 9, 599–607. doi: 10.1007/s11306-012-0473-x

PubMed Abstract | CrossRef Full Text | Google Scholar

McCabe, C. F., and Perng, W. (2017). Metabolomics of diabetes in pregnancy. Curr. Diab. Rep. 17:57. doi: 10.1007/s11892-017-0890-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Moriarty, M., Lee, A., O'Connell, B., Kelleher, A., Keeley, H., and Furey, A. (2011). Development of an LC-MS/MS method for the analysis of serotonin and related compounds in urine and the identification of a potential biomarker for attention deficit hyperactivity/hyperkinetic disorder. Anal. Bioanal. Chem. 401, 2481–2493. doi: 10.1007/s00216-011-5322-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Nägele, T., Mair, A., Sun, X., Fragner, L., Teige, M., and Weckwerth, W. (2014). Solving the differential biochemical Jacobian from metabolomics covariance data. PLoS ONE 9:e92299. doi: 10.1371/journal.pone.0092299

PubMed Abstract | CrossRef Full Text | Google Scholar

Nahum Sacks, K., Friger, M., Shoham-Vardi, I., Abokaf, H., Spiegel, E., Sergienko, R., et al. (2016). Prenatal exposure to gestational diabetes mellitus as an independent risk factor for long-term neuropsychiatric morbidity of the offspring. Am. J. Obstet. Gynecol. 215, 380.e1-7. doi: 10.1016/j.ajog.2016.03.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinto, J., Almeida, L. M., Martins, A. S., Duarte, D., Barros, A. S., Gahano, E., et al. (2015). Prediction of gestational diabetes through NMR metabolomics of maternal blood. J. Proteome Res. 14, 2696–2706. doi: 10.1021/acs.jproteome.5b00260

PubMed Abstract | CrossRef Full Text | Google Scholar

Prezelj, N., Covington, E., Roitsch, T., Gruden, K., Fragner, L., Weckwerth, W., et al. (2016). Metabolic Consequences of Infection of Grapevine (Vitis vinifera L.) cv. “Modra frankinja” with Flavescence Doree Phytoplasma. Front. Plant Sci. 7:711. doi: 10.3389/fpls.2016.00711

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosta, K., Al-Aissa, Z., Hadarits, O., Harreiter, J., Nadasdi, A., Kelemen, F., et al. (2017). Association Study with 77 SNPs Confirms the Robust Role for the rs10830963/G of MTNR1B Variant and Identifies Two Novel Associations in Gestational Diabetes Mellitus Development. PLoS ONE 12:e0169781. doi: 10.1371/journal.pone.0169781

PubMed Abstract | CrossRef Full Text | Google Scholar

Sano, M., Ferchaud-Roucher, V., Kaeffer, B., Poupeau, G., Castellano, B., and Darmaun, D. (2016). Maternal and fetal tryptophan metabolism in gestating rats: effects of intrauterine growth restriction. Amino Acids 48, 281–290. doi: 10.1007/s00726-015-2072-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Silverman, B. L., Metzger, B. E., Cho, N. H., and Loeb, C. A. (1995). Impaired glucose tolerance in adolescent offspring of diabetic mothers. Relationship to fetal hyperinsulinism. Diab. Care 18, 611–617. doi: 10.2337/diacare.18.5.611

PubMed Abstract | CrossRef Full Text | Google Scholar

Simmons, D., Devlieger, R., van Assche, A., Jans, G., Galjaard, S., Corcoy, R., et al. (2016). Effect of physical activity and/or healthy eating on GDM risk: the DALI Lifestyle Study. J. Clin. Endocrinol. Metab. 102, 903–913. doi: 10.1210/jc.2016-3455

PubMed Abstract | CrossRef Full Text | Google Scholar

Stein, S. E. (1999). An integrated method for spectrum extraction and compound identification from gas chromatography/mass spectrometry data. J. Am. Soc. Mass Spectrometry 10, 770–781. doi: 10.1016/S1044-0305(99)00047-1

CrossRef Full Text | Google Scholar

Steuer, R., Kurths, J., Fiehn, O., and Weckwerth, W. (2003). Observing and interpreting correlations in metabolomic networks. Bioinformatics 19, 1019–1026. doi: 10.1093/bioinformatics/btg120

PubMed Abstract | CrossRef Full Text | Google Scholar

St-Pierre, J., Laurent, L., King, S., and Vaillancourt, C. (2016). Effects of prenatal maternal stress on serotonin and fetal development. Placenta 48(Suppl. 1), S66–S71. doi: 10.1016/j.placenta.2015.11.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Sumner, L. W., Amberg, A., Barrett, D., Beale, M. H., Beger, R., Daykin, C. A., et al. (2007). Proposed minimum reporting standards for chemical analysis. Metabolomics 3, 211–221. doi: 10.1007/s11306-007-0082-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X., and Weckwerth, W. (2012). COVAIN: a toolbox for uni- and multivariate statistics, time-series and correlation network analysis and inverse estimation of the differential Jacobian from metabolomics covariance data. Metabolomics 8, 81–93. doi: 10.1007/s11306-012-0399-3

CrossRef Full Text | Google Scholar

Tobias, D. K., Hu, F. B., Forman, J. P., Chavarro, J., and Zhang, C. (2011). Increased risk of hypertension after gestational diabetes mellitus: findings from a large prospective cohort study. Diab. Care 34, 1582–1584. doi: 10.2337/dc11-0268

PubMed Abstract | CrossRef Full Text | Google Scholar

Weckwerth, W. (2003). Metabolomics in systems biology. Annu. Rev. Plant Biol. 54, 669–689. doi: 10.1146/annurev.arplant.54.031902.135014

PubMed Abstract | CrossRef Full Text | Google Scholar

Weckwerth, W. (2011). Unpredictability of metabolism–the key role of metabolomics science in combination with next-generation genome sequencing. Anal. Bioanal. Chem. 400, 1967–1978. doi: 10.1007/s00216-011-4948-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Weckwerth, W., and Morgenthal, K. (2005). Metabolomics: from pattern recognition to biological interpretation. Drug Discov. Tod. 10, 1551–1558. doi: 10.1016/S1359-6446(05)03609-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Weckwerth, W., Loureiro, M. E., Wenzel, K., and Fiehn, O. (2004a). Differential metabolic networks unravel the effects of silent plant phenotypes. Proc. Natl. Acad. Sci. U.S.A. 101, 7809–7814. doi: 10.1073/pnas.0303415101

PubMed Abstract | CrossRef Full Text | Google Scholar

Weckwerth, W., Wenzel, K., and Fiehn, O. (2004b). Process for the integrated extraction, identification and quantification of metabolites, proteins and RNA to reveal their co-regulation in biochemical networks. Proteomics 4, 78–83. doi: 10.1002/pmic.200200500

PubMed Abstract | CrossRef Full Text | Google Scholar

Wishart, D. S., Jewison, T., Guo, A. C., Wilson, M., Knox, C., Liu, Y., et al. (2013). HMDB 3.0–The human metabolome database in 2013. Nucleic Acids Res. 41, D801–D807. doi: 10.1093/nar/gks1065

PubMed Abstract | CrossRef Full Text | Google Scholar

Witkop, C. T., Neale, D., Wilson, L. M., Bass, E. B., and Nicholson, W. K. (2009). Active compared with expectant delivery management in women with gestational diabetes: a systematic review. Obstet. Gynecol. 113, 206–217. doi: 10.1097/AOG.0b013e31818db36f

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, H. H., Choi, S., and Levitt, P. (2016). Differential patterning of genes involved in serotonin metabolism and transport in extra-embryonic tissues of the mouse. Placenta 42, 74–83. doi: 10.1016/j.placenta.2016.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, A. H., Wang, X., Martinez, M. P., Walthall, J. C., Curry, E. S., Page, K., et al. (2015). Association of maternal diabetes with autism in offspring. JAMA 313, 1425–1434. doi: 10.1001/jama.2015.2707

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong, X., Saunders, L. D., Wang, F. L., and Demianczuk, N. N. (2001). Gestational diabetes mellitus: prevalence, risk factors, maternal and infant outcomes. Int. J. Gynaecol. Obstet. 75, 221–228. doi: 10.1016/S0020-7292(01)00496-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Yogev, Y., Metzger, B. E., and Hod, M. (2009). Establishing diagnosis of gestational diabetes mellitus: impact of the hyperglycemia and adverse pregnancy outcome study. Sem. Fetal Neonatal Med. 14, 94–100. doi: 10.1016/j.siny.2009.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Nagana Gowda, G. A., Asiago, V., Shanaiah, N., Barbas, C., and Raftery, D. (2008). Correlative and quantitative 1H NMR-based metabolomics reveals specific metabolic pathway disturbances in diabetic rats. Anal. Biochem. 383, 76–84. doi: 10.1016/j.ab.2008.07.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Mao, J., Ai, J., Shenwu, M., Shi, T., Zhang, D., et al. (2013). Integrated lipidomics and transcriptomic analysis of peripheral blood reveals significantly enriched pathways in type 2 diabetes mellitus. BMC Med. Genomics 6(Suppl. 1):S12. doi: 10.1186/1755-8794-6-S1-S12

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y., and Zhang, C. (2016). Prevalence of gestational diabetes and risk of progression to Type 2 diabetes: a global perspective. Curr. Diab. Rep. 16:7. doi: 10.1007/s11892-015-0699-x

PubMed Abstract | CrossRef Full Text | Google Scholar