Metabolomic and genetic architecture of gestational diabetes subtypes

Aims/hypothesis: Physiological gestational diabetes mellitus (GDM) subtypes that may confer different risks for adverse pregnancy outcomes have been defined. The aim of this study was to characterise the metabolome and genetic architecture of GDM subtypes to address the hypothesis that they differ between GDM subtypes. Methods: This was a cross-sectional study of participants in the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) study who underwent an OGTT at approximately 28 weeks’ gestation. GDM was defined retrospectively using International Association of Diabetes and Pregnancy Study Groups/WHO criteria, and classified as insulin-deficient GDM (insulin secretion <25th percentile with preserved insulin sensitivity) or insulin-resistant GDM (insulin sensitivity <25th percentile with preserved insulin secretion). Metabolomic analyses were performed on fasting and 1 h serum samples in 3463 individuals (576 with GDM). Genome-wide genotype data were obtained for 8067 individuals (1323 with GDM). Results: Regression analyses demonstrated striking differences between the metabolomes for insulin-deficient or insulin-resistant GDM compared to those with normal glucose tolerance. After adjustment for covariates, 33 fasting metabolites, including 22 medium- and long-chain acylcarnitines, were uniquely associated with insulin-deficient GDM; 23 metabolites, including the branched-chain amino acids and their metabolites, were uniquely associated with insulin-resistant GDM; two metabolites (glycerol and 2-hydroxybutyrate) were associated with the same direction of association with both subtypes. Subtype differences were also observed 1 h after a glucose load. In genome-wide association studies, variants within MTNR1B (rs10830963, p=3.43×10⁻¹⁸, OR 1.55) and GCKR (rs1260326, p=5.17×10⁻¹³, OR 1.43) were associated with GDM. Variants in GCKR (rs1260326, p=1.36×10⁻¹³, OR 1.60) and MTNR1B (rs10830963, p=1.22×10⁻⁹, OR 1.49) demonstrated genome-wide significant association with insulin-resistant GDM; there were no significant associations with insulin-deficient GDM. The lead SNP in GCKR, rs1260326, was associated with the levels of eight of the 25 fasting metabolites that were associated with insulin-resistant GDM and ten of 41 1 h metabolites that were associated with insulin-resistant GDM. Conclusions/interpretation: This study demonstrates that physiological GDM subtypes differ in their metabolome and genetic architecture. These findings require replication in additional cohorts, but suggest that these differences may contribute to subtype-related adverse pregnancy outcomes. Graphical Abstract: (Figure presented.)