Serum metabolomic signatures of gestational diabetes in South Asian and white European women.

INTRODUCTION: This study aimed to identify serum metabolomic signatures associated with gestational diabetes mellitus (GDM), and to examine if ethnic-specific differences exist between South Asian and white European women. RESEARCH DESIGN AND METHODS: Prospective cohort study with a nested case-control analysis of 600 pregnant women from two Canadian birth cohorts; using an untargeted approach, 63 fasting serum metabolites were measured and analyzed using multisegment injection-capillary electrophoresis-mass spectrometry. Multivariate logistic regression modeling was conducted overall and by cohort.
RESULTS: The proportion of women with GDM was higher in South Asians (27.1%) compared with white Europeans (17.9%). Several amino acid, carbohydrate, and lipid pathways related to GDM were common to South Asian and white European women. Elevated circulating concentrations of glutamic acid, propionylcarnitine, tryptophan, arginine, 2-hydroxybutyric acid, 3-hydroxybutyric acid, and 3-methyl-2-oxovaleric acid were associated with higher odds of GDM, while higher glutamine, ornithine, oxoproline, cystine, glycine with lower odds of GDM. Per SD increase in glucose concentration, the odds of GDM increased (OR=2.07, 95% CI 1.58 to 2.71), similarly for metabolite ratios: glucose to glutamine (OR=2.15, 95% CI 1.65 to 2.80), glucose to creatinine (OR=1.79, 95% CI 1.39 to 2.32), and glutamic acid to glutamine (OR=1.46, 95% CI 1.16 to 1.83). South Asians had higher circulating ratios of glucose to glutamine, glucose to creatinine, arginine to ornithine, and citrulline to ornithine, compared with white Europeans.
CONCLUSIONS: We identified a panel of serum metabolites implicated in GDM pathophysiology, consistent in South Asian and white European women. The metabolic alterations leading to larger ratios of glucose to glutamine, glucose to creatinine, arginine to ornithine, and citrulline to ornithine in South Asians likely reflect the greater burden of GDM among South Asians compared with white Europeans. © Author(s) (or their employer(s)) 2022. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

Emerging epidemiological studies in pregnant women have identified metabolites associated with gestational diabetes mellitus (GDM); however, few have examined the metabolome of South Asian women with GDM.

We identified a panel of serum metabolites involved in amino acid, carbohydrate, nitrogen, and lipid pathways that were altered in those with GDM, generally consistent in two ethnically distinct cohorts of South Asian and white European women.

The larger ratios of glucose to glutamine, glucose to creatinine, arginine to ornithine, and citrulline to ornithine in South Asian women may reflect early metabolic remodeling that contributes to the greater burden of GDM among South Asians compared with white Europeans.

Identification of metabolites implicated in GDM pathophysiology may help with targeted prevention strategies and to identify populations at high risk of incident metabolic complications (eg, type 2 diabetes).

Introduction

Gestational diabetes mellitus (GDM) is defined as hyperglycemia first recognized during pregnancy that is not attributable to previous diabetes.1 2 GDM can lead to several metabolic changes and clinical complications for pregnant women and their offspring, including stillbirth, preterm birth, neonatal macrosomia, and higher adiposity in offspring.1–3 Evidence is mounting of lasting intergenerational metabolic sequelae of GDM, where women and offspring have higher rates of incident type 2 diabetes, obesity, hypertension, dyslipidemia, and cardiovascular disease later in life.1 2 4 5 Globally, GDM rates are highest among South Asian (ie, people who originate from the Indian subcontinent), Latin American, Middle Eastern, and Indigenous populations.6 7 Specifically, the proportion of women with GDM is almost twofold higher among South Asians living in Canada compared with white Europeans.3 The reason for the increased risk of GDM among South Asian women is not well understood, but is hypothesized to be related to higher genetic risk, greater adiposity, and differences in dietary intake in pregnancy.3 6 8–10

The field of metabolomics is contributing to rapid improvements in our understanding of diseases of complex etiology at critical stages of development, such as GDM. This is because the human metabolome reflects the interplay of many of the complex factors that contribute to the disease phenotype, including gene expression, environmental exposures, and physiological responses.11 12 To date, 61 metabolomic studies have investigated the pathophysiology of GDM, but most were small studies conducted primarily in either white European (42% of studies) or East Asian (30% of studies) populations (online supplemental tables S1 and S2). Despite being a high-risk population, only seven studies to date8 13–18 have examined the metabolome of South Asian women with GDM, of which most included fewer than 20 GDM cases (online supplemental table S1). Non-targeted metabolomic profiling in a large population-based sample may reveal new insights into why South Asian women are disproportionately impacted by GDM as compared with other populations.8 Thus, we sought to identify serum metabolites associated with GDM, to characterize metabolic pathways related to GDM, and to determine if these metabolites differ between South Asian and white European pregnant women.

Research design and methods

Study participants

We compared GDM cases to non-cases from two prospective birth cohorts of pregnant women in Canada: the South Asian Birth Cohort (START) study and the Family Atherosclerosis Monitoring in Early Life (FAMILY) study.9 19 Each cohort enrolled pregnant women from Ontario, Canada, in their second or third trimester and followed the women and infants prospectively. Recruitment for the START study began in July 2011, and enrolled 1012 women with 1002 newborns.9 Recruitment for the FAMILY study occurred between 2004 and 2009, and enrolled 857 women and 901 newborns.19 A sample of 600 women (300 from START and 300 from FAMILY), selected for a nutrition substudy to be balanced across low, intermediate, and high diet quality scores within each cohort, and who had prior non-targeted metabolomics data, were eligible for this analysis.12

GDM classification

All women underwent a 75 g oral glucose tolerance test (OGTT) between 24 and 28 weeks of gestation to diagnose GDM. In START and FAMILY, GDM (ie, ‘case’ status) was diagnosed if any one of the following criteria was met: (1) The OGTT returned a result consistent with GDM based on the International Association of the Diabetes and Pregnancy Study Groups (IADPSG) definition: fasting glucose ≥5.1 mmol/L, 1-hour glucose ≥10.0 mmol/L, and/or 2-hour glucose ≥8.5 mmol/L, (2) a woman self-reported GDM at baseline or post partum, (3) a review of a woman’s birth chart indicated positive GDM status, or (4) a woman reported taking insulin during pregnancy, without prior diabetes diagnosis. In START, we additionally assessed GDM using the South Asian-specific Born in Bradford (BiB) criteria: fasting glucose level ≥5.2 mmol/L, or a 2-hour postload level of ≥7.2 mmol/L. The BiB GDM criteria were developed and validated in 4821 South Asian women living in the UK and showed that lower glucose thresholds predicted adverse perinatal outcomes, including infant adiposity and high infant birth weight.20

Serum metabolomics analysis by multisegment injection-capillary electrophoresis-mass spectrometry

Overnight fasted serum samples were collected from women and stored in liquid nitrogen at the Hamilton Clinical Research Laboratory. The workflow for the metabolomics analysis in START and FAMILY was the same and is described in detail elsewhere.21 In brief, multisegment injection-capillary electrophoresis-mass spectrometry (MSI-CE-MS) was used for the identification and quantification of polar ionic metabolites measured consistently in serum filtrate samples with rigorous quality control (QC).21 Most serum metabolites were unambiguously identified by spiking authentic standards, and subsequently quantified (μM) using a calibration curve, where ion responses were normalized to a single internal standard (ie, relative peak area (RPA)).21 Unknown serum metabolites were annotated based on their most likely molecular formula, relative migration time, and mode of detection with most compounds identified following acquisition of high-resolution tandem mass spectrometry (MS/MS) spectra at different collision energies (10, 20, 40 V).21 In START and FAMILY, 63 and 62 serum metabolites, respectively, met the selection criteria using MSI-CE-MS, namely they were frequently measured (>75% of serum samples) with adequate technical precision (coefficient for variation for QC <30%). Maternal serum metabolome coverage was similar in both cohorts, except for an unknown cation (m/z: 334.689, identified as a doubly charged peptide) only detected in START (online supplemental table S3). Ratios of serum metabolites previously identified to be related to glucose homeostasis, type 2 diabetes, or GDM were also investigated, including the ratios of glutamic acid to glutamine, isoleucine to glutamine, leucine to glutamine, valine to glutamine, glucose to creatinine, glucose to glutamine, arginine to ornithine, citrulline to ornithine, tyrosine to methionine, alanine to glycine, and Fischer’s ratio.22–27 Fischer’s ratio is a ratio of branched-chain amino acids (BCAAs) leucine, valine, isoleucine to aromatic amino acids (AAs) phenylalanine and tyrosine.27 Ratios for other serum metabolites with better discriminative ability for GDM as compared with glucose were also reported in our work.

Statistical analysis

Metabolites were autoscaled and natural logarithm transformed to correct for skewness. Student’s two-tailed t-tests were performed to compare means and Mann-Whitney U tests for medians using QC-based batch-corrected metabolite concentrations (or RPA for metabolites with no authentic standards commercially available) between pregnant women with and without GDM, and between South Asian and white European women. Each natural log-transformed metabolite was tested independently for its association with GDM, and those with univariate p<0.10 were candidates for subsequent multivariable logistic regression. There was no adjustment for multiple testing during this exploratory phase of analysis. Unknown serum metabolites significant at p<0.10 were eligible for additional MS/MS experiments. In multivariable logistic regression, the association between metabolite and GDM was expressed per 1 SD increase of log-transformed metabolite RPA, adjusted for ethnicity, maternal age, maternal height, pre-pregnancy body mass index (BMI), gestational age at recruitment, family history of diabetes, smoking history, and a diet quality score.12 These covariates were selected as they were found to be associated with GDM and various metabolites in previous analyses.2 8 9 12 15 16 28 29 Standardized protocols were used to measure anthropometrics (eg, weight, height). Participants self-reported ethnicity, with options of Hispanic, European, South Asian, Arab, East/Southeast Asia, African, First Nation, Mixed, and Other. Only participants of European and South Asian ethnicities were included for this study. Maternal age was calculated by subtracting the date of birth from the date of baseline visit. Pre-pregnancy BMI was calculated at baseline, using measured height and self-reported or measured weight. Gestational age at recruitment was defined as the number of weeks the woman has been pregnant, as reported from her expected delivery date. Participants had a family history of diabetes if at least one of their immediate family members (parents and full-blooded siblings) was diagnosed with diabetes. Smoking history was self-reported, with options of never smoking, quitting before pregnancy, quitting during pregnancy, or currently smoking. The diet quality score was calculated as the sum of daily servings of ‘healthy’ foods (eg, fermented dairy, vegetables, fruits) less the sum of daily servings of ‘unhealthy foods’ (eg, processed meats, refined grains), as previously described.12 To assess the robustness of regression models, we conducted a k-fold cross-validation (k=10). The associations between serum metabolites and GDM status were calculated overall, as well as stratified by cohort to examine any differences between South Asian and white European women. Multiplicative interactions between serum metabolites and ethnicity were tested. Analyses were conducted in SAS V.9.4 and R V.4.0.2.

Pathway analysis

Pathway analysis was performed using freely available online software, MetaboAnalyst V.5.0.30 Metabolites were generalized logarithm transformed, autoscaled, and then mapped into metabolic pathways using the Kyoto Encyclopedia of Genes and Genomes database. Pathway enrichment analysis compares the number of significant metabolites with an expected value within a specific pathway to determine the significant pathways (p<0.05) with false discovery rate correction. A pathway importance value was calculated for each metabolite based on its position relative to other molecules within the same pathway.30 Topology analysis calculates the pathway impact as the sum of the importance measures of the matched metabolites normalized by the sum of the importance measures of all the metabolites of each pathway. Pathways were retained if p<0.05 and impact was >0.10, thresholds consistent with prior MetaboAnalyst pathway analyses.31–33

Sensitivity analyses

The primary analysis was a complete-case analysis that assessed the associations of metabolites with GDM status, adjusted for ethnicity, maternal age, maternal height, pre-pregnancy BMI, gestational age at recruitment, family history of diabetes, smoking history, and diet quality. Sensitivity analyses compared this model with: (1) results of unadjusted models; (2) models using a different method of GDM diagnosis in START (IADPSG vs BiB); (3) using multiple imputation (m=10) because 7.7% of women (n=46) had missing data on at least one covariate; and (4) models additionally adjusting for total plant protein intake (% of total energy). To assess the consistency of our results, we also conducted partial least-squares discriminant analysis (PLS-DA), a supervised multivariate data analysis method for ranking serum metabolites of significance based on their variable importance in projection (VIP >1.5) score.

Results

Participant characteristics

Women with GDM were older, shorter, had higher pre-pregnancy BMI, and more likely to have a family history of diabetes than women without GDM (table 1). Despite lower pre-pregnancy BMI, South Asian women were more likely to have GDM (27.1%) than white European women (17.9%) (p=0.008; table 1). South Asian women were almost twice as likely to report a family history of diabetes (58.0%), compared with white European women (32.7%) with GDM (p<0.001). While both cohorts had similar total protein intakes, South Asian women with GDM consumed more plant protein (7.2% vs 4.8% of total energy, p<0.001) and less animal protein (8.7% vs 11.4% of total energy, p<0.001), respectively, compared with white European women with GDM.

Table 1

Participant characteristics for nested case–control metabolite analysis, IADPSG criteria for GDM diagnosis

	Overall, START and FAMILY (n=590)		START, South Asian mothers (n=299)		FAMILY, white European mothers (n=291)
	GDM (n=133, 22.5%)	No GDM (n=457)	GDM (n=81,27.1%)	No GDM (n=218)	GDM (n=52, 17.9%)	No GDM (n=239)
Age, mean (SD)	32.2 (4.9)	30.9 (4.4)	31.1 (4.2)	29.6 (3.5)	33.9 (5.4)	32.0 (4.8)
Gestational age at recruitment, mean (SD)	27.6 (3.1)	28.2 (3.3)	26.4 (2.1)	26.7 (1.6)	29.4 (3.6)	29.6 (3.8)
Pre-pregnancy BMI, mean (SD)	27.3 (6.6)	24.7 (5.1)	25.3 (5.2)	24.5 (3.8)	30.5 (7.3)	25.8 (5.8)
Mean height, cm (SD)	162.1 (6.2)	164.0 (6.7)	161.2 (6.0)	162.3 (6.5)	163.6 (6.4)	165.6 (6.5)
Gravida, mean (SD)	2.6 (1.6)	2.3 (1.3)	2.5 (1.3)	2.3 (1.1)	2.7 (1.9)	2.3 (1.5)
Family history of diabetes, n (%)	64 (48.1)	134 (29.3)	47 (58.0)	90 (41.3)	17 (32.7)	44 (18.4)
Unknown	7 (5.3)	20 (4.4)	2 (2.5%)	2 (0.9)	5 (9.6)	18 (7.5)
Parity, n (%)
Primiparous	49 (36.8)	196 (42.9)	27 (33.3)	79 (36.2)	22 (42.3)	117 (49.0)
Multiparous	84 (63.2)	260 (56.9)	54 (66.7)	139 (63.8)	30 (57.7)	121 (50.6)
Unknown	0 (0)	1 (0.2)	0 (0)	0 (0)	0 (0)	1 (0.4)
Glucose, mean (SD), mmol/L
Fasting, plasma	5.0 (0.8)	4.3 (0.3)	5.0 (0.8)	4.3 (0.3)	5.1 (0.7)	4.3 (0.3)
1-hour GTT	10.0 (1.9)	7.2 (1.4)	10.2 (2.0)	7.3 (1.4)	9.7 (1.8)	7.1 (1.4)
2-hour GTT	8.0 (2.0)	5.9 (1.2)	8.5 (2.1)	6.0 (1.2)	7.3 (1.6)	5.9 (1.2)
Maternal ethnicity, n (%)
South Asian	81 (60.9)	218 (47.7)	81 (100)	218 (100)	0 (0)	0 (0)
White European	52 (39.1)	239 (52.3)	0 (0)	0 (0)	52 (100)	239 (100)
Physical activity during pregnancy, n (%)
Sedentary	24 (18.0)	79 (17.3)	17 (21.0)	49 (22.5)	7 (13.5)	30 (12.6)
Mild	79 (59.4)	266 (58.2)	46 (56.8)	127 (58.3)	33 (63.5)	139 (58.2)
Moderate	29 (21.8)	108 (23.6)	18 (22.2)	41 (18.8)	11 (21.2)	67 (28.0)
Strenuous	1 (0.8)	4 (0.9)	0 (0)	1 (0.5)	1 (1.9)	3 (1.3)
Social disadvantage index, n (%)*
High	11 (8.3)	36 (7.9)	8 (9.9)	31 (14.2)	3 (5.8)	5 (2.1)
Moderate	34 (25.6)	89 (19.5)	24 (29.6)	65 (29.8)	10 (19.2)	24 (10.0)
Low	68 (51.1)	270 (59.1)	33 (40.7)	96 (44.0)	35 (67.3)	174 (72.8)
Unknown	20 (15.0)	62 (13.6)	16 (19.8)	26 (11.9)	4 (7.7)	36 (15.1)
Smoking history, n (%)
Never smoked	107 (80.5)	376 (82.3)	81 (100)	218 (100)	26 (50.0)	158 (66.1)
Quit before pregnancy	11 (8.3)	35 (7.7)	0 (0.0)	0 (0.0)	11 (21.2)	35 (14.6)
Quit during pregnancy	11 (8.3)	32 (7.0)	0 (0.0)	0 (0.0)	11 (21.2)	32 (13.4)
Current smoker	3 (2.3)	10 (2.2)	0 (0.0)	0 (0.0)	3 (5.8)	10 (4.2)
Unknown	1 (0.8)	4 (0.9)	0 (0.0)	0 (0.0)	1 (1.9)	4 (1.7)
Diet quality index, mean (SD)	5.4 (7.2)	4.1 (8.0)	7.8 (6.9)	6.9 (8.5)	1.7 (6.2)	1.5 (6.6)
Protein, mean (SD), % of total energy	16.4 (2.8)	16.1 (2.6)	16.1 (2.7)	15.7 (2.4)	16.9 (2.8)	16.5 (2.6)
Plant protein, mean (SD), % of energy	6.3 (1.8)	5.9 (1.7)	7.2 (1.5)	7.1 (1.4)	4.8 (1.2)	4.9 (1.2)
Animal protein, mean (SD), % of energy	9.7 (3.3)	9.8 (3.1)	8.7 (3.2)	8.4 (2.9)	11.4 (2.8)	11.1 (2.7)

*The social disadvantage index was constructed based on the sum of scores from household income, marital status, and employment status, and was previously validated.41 For income: 2 points were given if household income was <$30 000, 1 point if between $30 000 and $50 000, and 0 point if >$50 000. If married or common law 0 point was given, while 1 point was given if widowed, separated, divorced, or never married. If unemployed or retired 2 points were given, if employed 0 point was given. If the sum of points from household income, marital status, and employment status was 0–1 point participants were considered to have low disadvantage index, if 2–3 points moderate disadvantage index, and 4–5 points as high disadvantage index.

BMI, body mass index; FAMILY, Family Atherosclerosis Monitoring in Early Life; GDM, gestational diabetes mellitus; GTT, glucose tolerance test; IADPSG, International Association of the Diabetes and Pregnancy Study Groups; START, South Asian Birth Cohort.

Association between serum metabolites and GDM

Twenty-eight candidate serum metabolites were associated with GDM at the initial statistical threshold of p<0.10 (online supplemental tables S4 and S5; figure S1). Nine serum metabolites were associated with GDM in multivariable logistic regression models (figure 1). Higher levels of glucose, glutamic acid, propionylcarnitine, pyruvic acid, 2-hydroxybutyric acid, 3-hydroxybutyric acid, and 3-methyl-2-oxovaleric acid were associated with higher odds of GDM (figure 1). Greater ratios of glucose to glutamine, glucose to creatinine, glutamic acid to glutamine, and arginine to ornithine also corresponded with higher odds of GDM (figure 1). In contrast, higher levels of unknown m/z 129.066 (C5H8N2O2) were associated with lower odds of GDM (figure 1). With additional MS/MS experiments, this cation was tentatively classified as an AA analog, which was correlated with serum creatinine (r=0.44) and glutamine (r=0.32)31 (online supplemental figure S2). Our 10-fold cross-validation of multivariable logistic regression models demonstrated good predictive performance (area under the curve ranging from 0.69 to 0.76, sensitivity 0.94–0.97, online supplemental table S6), but the PLS-DA model was not robust (average Q2=0.05, R2=0.18).

Figure 1

Associations between serum metabolites and odds of gestational diabetes mellitus (GDM) (IADPSG diagnosed), per 1 SD increase in metabolite concentration, START and FAMILY (n=590 women). *Reported as most likely molecular formula (m/z). Adjusted for ethnicity, maternal age, maternal height, pre-pregnancy BMI, gestational age at recruitment, family history of diabetes, smoking history, and diet quality score. †Identified as a cation amino acid analog.21 ‡Identified as doubly charge peptide, containing alanine, glutamic acid, and histidine. BMI, body mass index; FAMILY, Family Atherosclerosis Monitoring in Early Life; IADPSG, International Association of the Diabetes and Pregnancy Study Groups; RPA, relative peak area; START, South Asian Birth Cohort.

Seventeen pathways were associated with GDM, including several AA metabolism pathways, as well as carbohydrate metabolism, nitrogen metabolism, and lipid metabolism pathways (table 2; figure 2; online supplemental figure S3).

Table 2

Metabolic pathways significantly impacted by GDM, with p<0.05 and impact >0.10, sorted by impact, KEGG database*

Pathway name†	Significant metabolites (univariatep<0.10)	Measured metabolites	Total metabolites	Metabolites involved	P value	FDR	Impact
Valine, isoleucine, and leucine (BCAAs) biosynthesis	1	6	8	Valine, leucine, 3-methyl-2-oxovaleric acid, α-ketoisovaleric acid, isoleucine, threonine	0.040	0.073	1.00
Alanine, aspartic acid, and glutamic acid metabolism	4	7	28	Aspartic acid, asparagine, alanine, glutamic acid, glutamine, citric acid, pyruvic acid	0.006	0.021	0.53
Glutamine and glutamic acid metabolism	2	2	6	Glutamic acid, glutamine	<0.001	0.004	0.50
Arginine biosynthesis	4	6	14	Glutamic acid, arginine, citrulline, aspartic acid, ornithine, glutamine	0.001	0.005	0.48
Arginine and proline metabolism	4	7	38	Arginine, creatine, glutamic acid, ornithine, pyruvic acid, proline, guanidinoacetic acid	0.001	0.005	0.37
Nitrogen metabolism	2	2	6	Glutamic acid, glutamine	<0.001	0.004	0.25
Histidine metabolism	1	4	16	Glutamic acid, histidine, methyl-l-histidine, aspartic acid	0.026	0.059	0.22
Pyruvic acid metabolism	2	2	22	Pyruvic acid, lactic acid	0.039	0.073	0.21
Glyoxylic acid and dicarboxylic acid metabolism	4	6	32	Citric acid, serine, glycine, glutamic acid, pyruvic acid, glutamine	0.004	0.016	0.18
Tryptophan metabolism	1	1	41	Tryptophan	0.027	0.059	0.14
Synthesis and degradation of ketone bodies	1	1	5	3-Hydroxybutyric acid	0.002	0.01	0.14
Tyrosine metabolism	1	2	42	Tyrosine, pyruvic acid	0.026	0.059	0.14
Citric acid cycle (tricarboxylic acid cycle)	2	2	20	Citric acid, pyruvic acid	0.028	0.059	0.14
Butyric acid metabolism	2	2	15	3-Hydroxybutyric acid, glutamic acid	<0.001	0.003	0.13
Cysteine and methionine metabolism	2	4	33	Serine, methionine, cystine, pyruvic acid	0.01	0.029	0.13
Glutathione metabolism	4	4	28	Glycine, glutamic acid, oxoproline, ornithine	0.001	0.006	0.12
Glycolysis/ gluconeogenesis	3	3	26	Glucose, pyruvic acid, lactic acid	<0.001	<0.001	0.10

*Adapted from MetaboAnalyst results (https://www.metaboanalyst.ca/).

†There were four additional pathways with p<0.05, where pathway impact was <0.10: valine, leucine and isoleucine degradation, porphyrin and chlorophyll metabolism, propanoate metabolism, and pyrimidine metabolism.

BCAA, branched-chain amino acid; FDR, false discovery rate; GDM, gestational diabetes mellitus; KEGG, Kyoto Encyclopedia of Genes and Genomes.

Figure 2

Pathway impact by significance, metabolites significantly associated with gestational diabetes mellitus (GDM), South Asian Birth Cohort (START) and Family Atherosclerosis Monitoring in Early Life (FAMILY) (n=590 women). The size of pathway symbols represent the significance level, and the color is based on impact factor. (Adapted from MetaboAnalyst30 results (https://www.metaboanalyst.ca/)).

Ethnic-specific differences

Comparing the two groups, the ratios of glucose to glutamine, glucose to creatinine, citrulline to ornithine, and arginine to ornithine were 0.4–1.0 SD higher in South Asians compared with white European women (online supplemental figure S4). Among women with GDM, there were significant differences in the adjusted metabolite concentrations between South Asians (n=81) and white Europeans (n=52) (online supplemental figures S5 and S6). Concentrations of several AAs, carboxylic acids, and hydroxy acids were higher among white European women with GDM, except for 2-aminooctanoic acid, tryptophan betaine, arginine, glucose to creatinine, and arginine to ornithine ratios (online supplemental figure S6). Despite ethnic differences in the distribution of many serum metabolites, ORs associating metabolites with GDM were of similar magnitude and direction in both groups (table 2; online supplemental tables S7 and S8).

Sensitivity analysis

Multiple imputation, compared with complete case, or further adjusting for plant protein intake, did not meaningfully change the results (online supplemental table S9). For most serum metabolites, adjusted models, compared with unadjusted, slightly attenuated the relationship between metabolites and GDM (online supplemental table S9). Multivariate-adjusted models using BiB-defined GDM in START and IADPSG in FAMILY resulted in the same nine serum metabolites associated with IADPSG-defined GDM, plus an additional five metabolites: tryptophan, ornithine, alanine to glycine ratio, and citrulline to ornithine ratio (online supplemental table S9). Results from a PLS-DA model with VIP score ranking were also consistent with the main analysis results (online supplemental figure S7).

Discussion

This metabolomics analysis in two ethnically distinct cohorts of women in pregnancy highlighted several pathways that were altered with GDM and common to both South Asian and white European women. These GDM-altered pathways included AA metabolism (eg, BCAA biosynthesis, tryptophan metabolism), carbohydrate metabolism, and fatty acid metabolism (eg, ketone body production, butyric acid metabolism). Further, several pathways involving glutamic acid and glutamine were consistently associated with GDM, such as glutathione metabolism and nitrogen metabolism. Within these pathways, higher concentrations of glutamic acid, propionylcarnitine, tryptophan, arginine, 2-hydroxybutyric acid, 3-hydroxybutyric acid, pyruvic acid, and 3-methyl-2-oxovaleric acid were associated with higher odds of GDM. Higher circulating glutamine, ornithine, oxoproline, cystine, glycine, and unknown m/z 129.066 (C5H8N2O2, identified as a cation AA derivative) were inversely related to GDM. This is the first study to consider serum metabolite ratios as predictors of GDM, discovering increased odds for GDM with increasing ratios of glucose to glutamine, glucose to creatinine, glutamic acid to glutamine, and arginine to ornithine. Increasing evidence suggests that GDM reflects an early stage of type 2 diabetes disease progression that manifests due to the metabolic stress caused by pregnancy.2 4 5 Overall, our study results support this hypothesis, as we identified GDM-altered pathways consistent with previous investigations in participants with GDM or type 2 diabetes2 15 16 34 35 (online supplemental table S1).

We observed a positive association between the glutamic acid to glutamine ratio and GDM. This ratio is an established type 2 diabetes risk marker, but ours is the first study to show that this ratio was also associated with GDM.22 34 Glutamic acid and glutamine also play important roles in seven of the 17 pathways we identified to be associated with GDM, including AA metabolism (alanine, aspartic acid, histidine, arginine, proline), nitrogen metabolism, and glyoxylic and dicarboxylic acid metabolism. Glutamine and glutamic acid provide carbon for glucose production in the kidney and liver, where glutamic acid promotes insulin secretion from pancreatic β-cells.22 Abnormal glutamic acid homeostasis is associated with increased pancreatic β-cell damage, insulin secretion, and gluconeogenesis, which can be mitigated, in part, by increased glutamine.22 36 In type 2 diabetes studies, even after adjusting for insulin resistance and BCAA levels, the glutamic acid to glutamine ratio remained a strong predictor of incident type 2 diabetes.22 In a systematic review of 19 clinical trials and animal studies, glutamine improved fasting plasma glucose, postprandial blood glucose, and increased insulin production in patients and animals with type 2 diabetes.36 Therefore, glutamine may be an important therapeutic target for GDM, as well as a marker of future metabolic risk; however, additional studies are warranted to investigate this hypothesis.36

While several studies identify elevated BCAAs as important risk markers for insulin resistance, pre-diabetes, and type 2 diabetes, the relationship between circulating BCAAs and GDM is inconsistent34 35 (online supplemental table S1). Seven previous studies have identified positive relationships between BCAAs and GDM, while 12 studies found no difference (online supplemental table S1). We did not identify significant associations between either BCAAs or Fischer’s ratio and GDM. We did, however, find that the BCAA biosynthesis pathway was related to GDM, along with significant positive associations of 3-methyl-2-oxovaleric acid and propionylcarnitine with GDM, two downstream leucine metabolism by-products, both previously linked to impaired fasting glucose and type 2 diabetes mellitus.34 35 We also identified that tryptophan metabolism was associated with GDM, consistent with prior investigations.15 16 35 Elevated tryptophan can result in higher production of the diabetogenic xanthurenic acid via the tryptophan-kynurenine pathway, which may explain why it was elevated in women with GDM.2

We observed that higher levels of 2-hydroxybutyric acid and 3-hydroxybutyric acid were associated with increased odds of GDM, consistent with previous studies (online supplemental table S1). It is postulated that the increased fatty acid catabolism, oxidative stress, and enhanced glutathione synthesis observed in GDM are related to higher circulating concentrations of 2-hydroxybutyric acid.37 This shift to fatty acid oxidation is further supported by increased 3-hydroxybutyric acid, which is a major ketone body in circulation in a fasted state.37 Elevated levels of 2-hydroxybutyric acid and 3-hydroxybutyric acid have been implicated as predictive biomarkers several years prior to the onset of type 2 diabetes, highlighting the potential clinical utility of these compounds for early detection of glucose intolerance in women during pregnancy prior to onset of GDM.37

Consistent with the BiB study of 4072 white European and 4702 South Asian women living in the UK, among all pregnant women, our study demonstrated that South Asian women had higher levels of glucose, while white Europeans had higher levels of histidine, creatinine, and pyruvic acid.8 Contrary to the BiB study, we did not observe higher levels of glutamine, alanine, BCAAs, or aromatic AAs in South Asian women.8 We did, however, observe higher levels of serine, arginine, tryptophan betaine, 3-methyl-2-oxovaleric acid, and acetylcarnitine. Variability between studies may reflect our lower sample size, or dietary pattern and socioeconomic status differences. In our study, the higher ratios of glucose to glutamine, glucose to creatinine, arginine to ornithine, and citrulline to ornithine in South Asians compared with white Europeans are important ethnic-specific findings. These four ratios are strongly related to GDM and incident type 2 diabetes and could reflect early metabolic perturbations that contribute to the larger burden of GDM, incident type 2 diabetes, and other endocrine disorders in South Asians compared with white Europeans.22–27 34 35 Elevated circulating arginine and larger ratios of arginine to ornithine, and citrulline to ornithine related to increased risk for GDM in South Asian women are signs of abnormal urea cycle metabolism.23 Increased arginine relative to ornithine levels can also indicate an impairment of nitric oxide synthesis and vasodilation properties, which has previously been documented in people with type 2 diabetes.23 Our study suggests that AAs in urea cycle may be candidate markers for predicting GDM, particularly among South Asian women.

Metabolite differences by ethnicity were robust after adjustments for age, height, pre-pregnancy BMI, gestational age at recruitment, family history of diabetes, smoking, and diet quality. However, differences in specific food groups may not have been fully captured in the diet quality score.12 For example, the higher amount of tryptophan betaine in South Asian women may be explained by higher legume intake, while higher 3-methylhistidine in white Europeans by their larger protein and meat intake.12 38 However, ethnic-specific metabolite concentrations were robust to further adjustments for total intake of plant protein (online supplemental figure S6). There is a growing body of literature that physical activity and sedentary behavior are associated with various metabolite concentrations, including BCAAs, alanine, proline, and lactic acid.39 Nevertheless, further adjustment for self-reported physical activity (data not shown) did not significantly change ethnic-specific metabolite concentrations in our study. Other lifestyle and genetic variations not explored herein may also contribute to observed differences.

Strengths and limitations

This study investigated the metabolic phenotype differences associated with GDM risk among South Asian women, an understudied population, as compared with white European women; an important strength compared with previous metabolomic studies (online supplemental tables S1 and S2). Both cohorts provided fasting measurements, used a validated and identical technical platform for metabolomic analyses with stringent QC and batch correction adjustment.21 We conducted internal cross-validation, and sensitivity analyses to demonstrate robustness and similar outcomes using complementary univariate and multivariate statistical models. We adjusted our analyses for known important confounders; however, we cannot exclude the possibility of residual confounding from those that were unmeasured. Our study was limited as only polar ionic metabolites in serum filtrates were analyzed, but not fatty acids or lipids.40 Also, blood samples were collected at the same time as the OGTT was measured, and thus some degree of reverse causation is possible. Our analysis included 62 serum metabolites measured in two birth cohorts; however, larger scale investigations may be needed to validate metabolic differences implicated in GDM risk in different ethnic populations.

Conclusion

In summary, we identified a panel of serum metabolites involved in AA, carbohydrate, nitrogen, and lipid pathways that were altered with GDM. Although there was a large overlap in metabolites associated with GDM in South Asian and white European women, the larger ratios of glucose to glutamine, glucose to creatinine, arginine to ornithine, and citrulline to ornithine in South Asians may reflect early metabolic remodeling that contributes to the greater burden of GDM among South Asians compared with white Europeans. Further work is necessary to understand if interventions to alter these metabolite ratios can help alleviate the burden of GDM. GDM research, prevention, diagnosis, and management predominately focus on glycemic control, but this may be an oversimplified framework in the era of precision medicine.