Evidence That Multiple Defects in Lipid Regulation Occur before Hyperglycemia during the Prodrome of Type-2 Diabetes

Background Blood-vessel dysfunction arises before overt hyperglycemia in type-2 diabetes (T2DM). We hypothesised that a metabolomic approach might identify metabolites/pathways perturbed in this pre-hyperglycemic phase. To test this hypothesis and for specific metabolite hypothesis generation, serum metabolic profiling was performed in young women at increased, intermediate and low risk of subsequent T2DM. Methods Participants were stratified by glucose tolerance during a previous index pregnancy into three risk-groups: overt gestational diabetes (GDM; n = 18); those with glucose values in the upper quartile but below GDM levels (UQ group; n = 45); and controls (n = 43, below the median glucose values). Follow-up serum samples were collected at a mean 22 months postnatally. Samples were analysed in a random order using Ultra Performance Liquid Chromatography coupled to an electrospray hybrid LTQ-Orbitrap mass spectrometer. Statistical analysis included principal component (PCA) and multivariate methods. Findings Significant between-group differences were observed at follow-up in waist circumference (86, 95%CI (79–91) vs 80 (76–84) cm for GDM vs controls, p<0.05), adiponectin (about 33% lower in GDM group, p = 0.004), fasting glucose, post-prandial glucose and HbA1c, but the latter 3 all remained within the ‘normal’ range. Substantial differences in metabolite profiles were apparent between the 2 ‘at-risk’ groups and controls, particularly in concentrations of phospholipids (4 metabolites with p≤0.01), acylcarnitines (3 with p≤0.02), short- and long-chain fatty acids (3 with p< = 0.03), and diglycerides (4 with p≤0.05). Interpretation Defects in adipocyte function from excess energy storage as relatively hypoxic visceral and hepatic fat, and impaired mitochondrial fatty acid oxidation may initiate the observed perturbations in lipid metabolism. Together with evidence from the failure of glucose-directed treatments to improve cardiovascular outcomes, these data and those of others indicate that a new, quite different definition of type-2 diabetes is required. This definition would incorporate disturbed lipid metabolism prior to hyperglycemia.


Introduction
The metabolic basis of type 2 diabetes mellitus (T2DM) has traditionally had hyperglycemia as its sine qua non, despite generally being accompanied by a long prior history of (central) obesity together with relative physical inactivity. Evidence suggests that blood vessel dysfunction, either overt or inducible, is detectable prior to rises in blood glucose [1][2][3], as occurs in the disease itself [4]. Debate over whether glucose is the direct cause of the blood vessel damage has not yet been resolved. Many lines of evidence suggest that hyperglycemia may not be the earliest metabolic change in the complications of T2DM. One, based on current treatment results in clinical trials, is that complications are not prevented by glycemic control, intensive or not [5][6][7], confirmed by the latest very large trials of dipeptidyl peptidase-4 (DPP-4) inhibitors [8]. Earlier evidence suggested that microvascular components were delayed more by lowered blood pressure [6,[9][10][11] than by tight blood glucose control. The ACCORD (Action to Control Cardiovascular Risk in Diabetes) trial reported no overall difference in microvascular outcomes in diabetic subjects at risk of cardiovascular events, and intensive glycemic treatment was associated with higher mortality [12]. HMG CoA reductase inhibitors with its anti-flammatory and anti-thrombotic effects [13,14], have been used to target successfully total and LDL-cholesterol [15,16]. Despite such reductions, intriguingly statin treatment may marginally increase glycemia [17][18][19]. A second line of evidence is that the impaired blood vessel responsiveness is in both large arteries in vivo [1] and smaller arteries in tissue biopsy studies, which occurs even when blood glucose is normal [20]. A third set of arguments, coupled to an extensive literature reviewed elsewhere, is that iron and copper dysregulation are implicated in diseases that manifest in changes in both lipid and carbohydrate metabolism (and their attendant comorbidities) [21,22] Metabolomics is a systems biology strategy for exploring the low molecular weight metabolites present in the metabolome of an organism [23]. It portrays a dynamic interaction of a phenotype with the environment, across genomic and post-transcriptional regulation [23] and has been applied to study cardiovascular diseases [24][25][26][27][28][29] including heart failure [30], myocardial ischemia [31,32], myocardial infarction [33] and preeclampsia [34]. Its application in the investigation of glucose intolerance [35][36][37][38] has led to the identification of new metabolic biomarkers and has highlighted the influence of drugs on the metabolic profile of subjects diagnosed with glucoregulatory disorders [39,40]. Animal studies using targeted metabolomic approaches have confirmed that mitochondrial overload and incomplete fatty acid oxidation in skeletal muscle occur in both major types of diabetes [41].
Gestational glycemic status including overt gestational diabetes mellitus (GDM) increases susceptibility to subsequent development of the T2DM 'phenotype' postnatally [42], although a confounding factor is obesity [43].
Here, we examined the early metabolic natural history of 'prediabetes' by comparing the serum metabolic profiles of women from three backgrounds, systematically determined in the third trimester of pregnancy. However, here, we chose a data-driven approach free of specific hypotheses [44] to determine which metabolite classes might be so changed on a number of pathways. All these women were followed for some two years postnatally when serum samples for metabolomic analysis were taken. Our main hypothesis was that the metabolome at follow-up would differ significantly between those women at high risk of T2DM (having had previous GDM) compared to those who remained normoglycemic throughout pregnancy and a third group who were normoglycemic during pregnancy but in the upper quartile of the glycemic distribution. Samples of these women were included in the vascular sub-study [3].

Study Population
The Hyperglycemia and Pregnancy Outcome (HAPO) study was a multi-centre study investigating the impact of glycemia below (but not including) overt diabetes in singleton pregnancies of women not taking anti-hypertensive drugs nor any other chronic therapies. Inclusion criteria were that women were at least 3 months pregnant, were to deliver at our local maternity hospital, and had completed a 75 g oral glucose tolerance test (GTT) at 24-32 weeks gestation.
To establish our sampling frame (Figure 1), we used the glycemic distribution from the first 957 participants recruited at the Manchester site of the Hyperglycemia and Pregnancy Outcome (HAPO) study [45]. From the group of 250 of these women who were initially followed up, we selected 100 women, including all 18 with previous GDM and 82 additional participants by computer-generated random sampling, who were stratified into the upper quartile of the original glycaemic distribution (UQ group) or below that distribution's median (control group). They were then matched for confounding factors of age, BMI and ethnicity in that order. There was no prior nor current use of statins/other cardiometabolic medications in these young women. The three final study groups were as follows: i) the 18 women who fulfilled the WHO definition of overt GDM at their HAPO GTT (GDM group); ii) 39 women with an index gestational fasting plasma glucose (FPG) value $4.8 but ,5.5 mmol?L 21 , and/or a 2 hr glucose value of $6.8 but ,7.8, mmol?L 21 (i.e. these were the upper quartile ('UQ') cut-off values for the whole distribution in the 957 original women, below GDM, forming the UQ group); and iii) 43 women whose gestational FPG had been #4.5, and 2hr plasma glucose #5.8, mmol?L 21 (from the lower half of the original GTT distribution) (Control Group) -see Figure 1.
Follow-up was performed at a mean of 22 months after the index pregnancy when fasting blood serum samples for metabolic profiling and, if possible, 2-h GTTs were, repeated.

Anthropometric measurements
All anthropometric measurements were taken by trained staff following WHO guidelines [46]. Total body-fat estimation was via a widely employed bioimpedance method (Bodystat 1500, Bodystat Ltd, UK).

Biochemical measurements
Blood samples were centrifuged, and serum and plasma aliquotted, immediately frozen and maintained at 280uC for later analysis of lipids and hormones. Blood samples were analyzed for glucose, triglyceride (TG), total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), NEFA, insulin, adiponectin, and leptin. We measured glucose by the glucose oxidase method on a Beckman Synchron LX system. HbA 1c was assayed by reversed phase cation exchange chromatography (Menarini Diagnostic, UK). Serum  adiponectin and leptin were measured using ELISA (R&D Systems, Minneapolis, MN), and insulin with monoclonal-based ELISAs (Mercodia, Sweden). Serum TC and TG were measured by the CHOD/PAP and GPO/PAP methods respectively on a Cobas Mira S analyzer (ABX Diagnostics, Shefford, UK); all reagents were obtained from the same source. HDL-C was measured by a second-generation homogenous method using PEG-modified enzymes (Roche Diagnostics, Lewes, UK). LDL-C was calculated using the Friedewald formula. A calculated LDL-C value of ,0.1 mmol.L 21 was set as the detection limit for cholesterol. Finally, we measured non-esterified free fatty acids (NEFA) in plasma using an enzymatic endpoint assay (WAKO Chemicals, Richmond VA) with a detection limit of 0.01 mmol.L 21 .

Preparation of serum samples for metabolomic analysis
Fasting serum samples taken at follow-up, with group of origin blinded to the analyst, were thawed on ice and prepared as previously described [47,48]. Samples were deproteinised by mixing 200-mL plasma with 600-mL methanol followed by vortexmixing (15 s) and centrifugation (15 min, 13,865 g). 370-mL aliquots of each supernatant were transferred to two Eppendorf tubes and lyophilised (HETO VR MAXI vacuum centrifuge attached to a Thermo Svart RVT 4104 refrigerated vapour trap; Thermo Life Sciences, Basingstoke, UK). The aliquots were for separate positive-and negative-ion Ultra Performance Liquid Chromatography-Mass Spectrometry (UPLC-MS) analyses, respectively. A pooled quality control (QC) sample [47] was also prepared by mixing 50-mL aliquots of serum from each of 100 subject samples followed by vortex mixing for one minute. 200-mL aliquots of the pooled QC sample were deproteinised and lyophilised as described above.

UPLC-MS analysis
Each sample was reconstituted in 100-mL water and vortex mixed (15 s) and then centrifuged (15 min, 13,865 g). Samples were analysed in a random order using Ultra Performance Liquid Chromatography (UPLC; Waters, Elstree, UK) coupled to an electrospray hybrid LTQ-Orbitrap mass spectrometer (Thermo-Fisher Scientific, Bremen, Germany). Each sample was analysed twice, once in negative-ion mode and once in positive-ion mode.
The analytical conditions [47,49] and application of QC samples [47] were applied as described previously.

Raw data processing and data analysis
Raw data files (.RAW) were converted to the NetCDF format using the File converter program in XCalibur (ThermoFisher Scientific, Bremen, Germany). Deconvolution of data was performed using XCMS as described previously [48], as were signal correction and quality assurance procedures [34].

Statistical analysis
Statistical analyses were carried out using STATA version 12 (Stata Corporation, College Station, Texas) or programs written in the MatlabH scripting language (version 7.8; http://www. mathworks.com/). Summary statistics of non-normally distributed continuous variables are presented as geometric means -derived from log-transformed data. Univariate analysis was performed using the Mann-Whitney U test, a non-parametric method for assessing whether two independent samples come from the same distribution. We used maximum-likelihood multinomial logit models to assess the relationship between levels of adiponectin, leptin, indices of adiposity (BMI), smoking status, triglyceride, nonesterified fatty acids (NEFA) as well as cholesterol and the likelihood of having GDM or the UQ of glycemia compared to the control group. Missing values were ignored.

Annotation of putative metabolites matched to features
Metabolic features characterized by measuring both the accurate m/z and retention time, and corresponding putative molecular annotations were assigned by standard methods as described [50]. One or more molecular formulae within available databases were assigned to each feature with mass accuracy of 63 ppm. These were subsequently searched against The Manchester Metabolomics Database, which has been constructed with information from the Human Metabolome Database (http:// www.hmdb.ca/, v2.0) and Lipidmaps (http://www.lipidmaps. org/). This is a level 2 annotation according to the proposed reporting standards of the Metabolomics Standards Initiative [51]. In these types of raw metabolomic data, a single metabolic feature can be assigned to one or more metabolites due to uncertainty caused by possible isomerism, resulting in a non-specific annota-   [51]. Where more than one putative structure can be assigned to any analytical feature corresponding to a particular molecular mass (that is, more than one molecule of the particular mass could occur in physiology), each possible annotation has been listed with 'AND/OR' as the conjunction. To minimise the influence of false discovery we  PCA showed no clustering related to class or sub-clustering of subjects from one or multiple classes therefore these data were not included in the manuscript. Similarly, Partial Least Squares -Discriminant analyses (PLS-DA) was also performed but no validated models were constructed and therefore these data were not reported.

Subject group characteristics
Standard anthropometric and metabolic parameters were measured in all participants, and stratified by glycemic status (Table 1). No significant between-group differences were present in age, ethnicity, BMI nor smoking status at follow-up. Small but significant differences in fasting and two-hour serum glucose concentrations occurred during pregnancy as expected. Significant between-group differences were observed at follow-up in waist circumference, adiponectin, fasting glucose, post-prandial glucose and HbA 1c , with means and all ranges are still within the 'normal' range (as defined by WHO), between control and both UQ and GDM women.
In maximum-likelihood multinomial logit models, increasing adiponectin concentrations (60% reduction in risk per mg.

Between-group differences in metabolite concentrations
3,552 metabolomic features were judged suitable for univariate analysis after raw metabolite data and related quality assurance processes had been performed. Levels of numerous metabolites differed significantly between groups. Data are presented (Tables 2 to 4), according to a metabolite classification system in which each molecule is listed as a member either of a structural class (e.g. 'short-chain fatty acids and related metabolites') or a functional class according to its participation in a defined metabolic process (e.g. ,participating in. 'tetrahydrofolate metabolism'). Each metabolite has been listed only once as a member of a single class. If a metabolite was detected more than once, the feature with the lowest p value was reported. Within each class, data have been separated into those with higher and lower ratios and are then presented in order from lowest to highest p-value. Figure 2 shows in ascending order of fold difference, the top 32 metabolites for all three comparisons.
By comparing the control and UQ groups, 173 of 3552 metabolic features were statistically different (p,0.05). Of these, 43 unique metabolic features were annotated (Table 2). 35  (Figure 4), were present at lower concentrations in the UQ than in controls, as were levels of certain vitamin D metabolites and the anorectic pentapeptide, enterostatin [52].
In comparisons between the Control and GDM groups, 392 of 3,552 metabolic features differed significantly. Of these 392 metabolites, 69 unique metabolic features were annotated (Table 3). Here the picture was more evenly balanced, with about one half of the metabolites in each class higher (35 of 69 metabolites) and the remainder lower in the UQ than the control group.
For the UQ versus GDM comparison, 401 of 3552 metabolic features differed significantly, of these, 72 unique metabolic features were annotated, ( Table 4). Many of the metabolites annotated in Table 3 recurred in Table 4, but notably the metabolite 2-Oxoglutaramate, an important biomarker of hepatic encephalopathy and other hyperammonemic diseases [53,54] was twice as abundant in previous GDM compared those in the UQ subgroup (2.11 (0.79, 3.96; p = 0.015)). A notable fold change in the metabolite N-(aminomethyl)urea was observed between the groups (Table 3 and Table 4).

Discussion
The pathophysiological metabolic changes in the very early stages of type 2 diabetes, before measurable hyperglycemia, remain comparatively little known or understood. Our current results provide compelling evidence for the occurrence of significant metabolic defects that antedate the onset of hyperglycemia, even if marginal differences in glycemia well within the normal range were present. These metabolic defects may exert effects that can lead to or cause subsequent glucoregulatory decompensation deteriorating to 'hyperglycemia', which currently defines the disease.
The particular metabolic pathways suggested by this study are defects in those regulating systemic lipid metabolism [55] and hormone secretion/responsiveness [56]; they appear to antedate and could therefore 'cause' or lead to overt hyperglycemia. Hormones currently implicated in the development of T2DM include the beta-cell hormones insulin and amylin [55,56], and the adipocyte hormones leptin [57] and adiponectin [58]. Early damage to blood vessels [1,59] and pancreatic islet beta-cells [60], for example, provide evidence for metabolic defects that antedate diabetes. Copper homeostasis and iron status are also related to GDM [61][62][63]. For example, high body iron stores, leading to unliganded iron, cause hydroxyl radical formation via Fenton chemistry and are significantly associated with a greater risk of T2DM [22,[64][65][66][67]. Here, 29 of the women were included in a vascular sub-study where there was a gradation of declining endothelial function of resistance blood vessels ex-vivo, poorest in the 12 of the 18 women with prior GDM studied here and less marked in those with UQ, compared with controls defined the same way [3]. Those vascular findings parallel the metabolic changes reported here.
To address questions of what metabolic markers identify the pathogenic pathways to T2DM and from them potential new strategies for disease prevention, we compared the 2 at-risk groups with controls to quantify specific metabolic differences between groups. The data suggest that some pathogenic processes may have begun by the time women reached the UQ state, with others underway when they further deteriorate, previously indicated by being GDM. Several distinct if overlapping molecular processes may underpin these successive degrees of regulatory impairment represented by the two increased-risk states. Dividing the complex time-dependent process into stages produces artificial categories but enables identification of earlier-and later-onset pathways.
Twenty-two months after their index pregnancy, when originally profiled by their glucose tolerance, the women had this status re-assessed by fasting plasma glucose and hemoglobin A 1c values. In contrast to their within-pregnancy glucose tolerance, glycaemic indices at re-testing were not different between the UQ and GDM groups, although both were marginally defective compared to controls, yet still within the usual, 'currently normal' glycaemic range. Pair-wise between-group comparisons pinpointed relatively circumscribed subsets of defined metabolite classes related to elevated diabetes risk. Those metabolite classes perturbed in the UQ compared with control women included: phospholipid subclasses, in particular phosphatidylcholines; LCFA; LCFAcarnitines; SCFA and SCFA-metabolites. Other perturbed classes included diglycerides; bile acids; steroids; prostanoids; and amino Prominent differences in phosphatidylcholines were identified in both the control/UQ and UQ/GDM contrasts. Diacyl-phosphatidylcholines has been shown to be independently associated with increased risk of type 2 diabetes in a prospective study of type 2 diabetes patients in the European Prospective Investigation into Cancer and Nutrition (EPIC)-Potsdam cohort [27]. Phospholipids are highly insoluble in aqueous media so these molecules will have originated in membranous structures in plasma, namely lipoproteins; this suggests that differences in phosphatidylcholine composition are related and could contribute to glucoregulatory transitions preceding hyperglycemia. Alterations in additional lipid classes including those of steroids/bile acids, and diglycerides are also probably related to changes in lipoprotein metabolism. Consistent with these findings, diabetes itself is associated with prominent changes in plasma lipoprotein content [68,69]. This disturbance in phospholipid metabolism cannot be localised or characterised further here since the observed changes could reflect alterations in any or all of the HDL, LDL, or VLDL fractions. Prominent alterations in LDL-particle composition have previously been identified in diabetes pathogenesis [70], lipoproteinbound phospholipids are reportedly targets of glycoxidationmediated damage [71], and oxidized phospholipids can become pathogenic [72,73]. Such direct effects of lipid alterations on blood vessels possibly underlie the major benefits of statin treatment in T2DM, although statins are also thought to be anti-inflammatory [22,74]. Whether, and through what pathways, statins may lead to increases in glycemia [17,18] remain unanswered questions relevant to this early pathogenesis. Altered lipoproteins are also implicated in the mechanisms that lead to or cause beta-cell dysfunction in diabetes [75]. Follow-up proteomic and metabolomic studies of purified lipoprotein fractions from different classes of at-risk patients would now help identify the specific molecules more clearly and may in time be useful in improving the performance of classification models based on standard factors [76]. The data here clearly point to early alteration in lipoprotein metabolism in the chain of events that culminate in diabetes and its complications.
A lysophospholipid-related signal may also be present, particularly in the UQ/GDM and control/GDM comparisons, indicating the onset of pro-inflammatory stress, which contributes to tissue damage. Plasma lyosphosphospholipid content is another potential biomarker for monitoring oxidative damage caused, for example [77] by perturbed regulation of catalytically-active copper metabolism before and in diabetes [78]. Lysophospholipid measurements could help monitor progression of tissue damage in people at risk of developing diabetes, and perhaps the response to preventive/therapeutic interventions.
Another significantly perturbed lipid-related signal here was for LCFA and LCFA-carnitines. Both classes tended to be lower in the UQ than in controls. Gall et al reported that medium-chain acylcarnitines such as decanoylcarnitine decreased in concentration with increasing insulin resistance and dysglycemia [35]. In the population-based Cooperative Health Research in the Region of Augsburg (KORA) cohort, three metabolites, namely glycine, lysophosphatidylcholine (LPC) (18:2) and acetylcarnitine had significantly altered levels in IGT individuals as compared to those with normal glucose tolerance [79]. Acylcarnitines are biosynthesized solely in mitochondria, where they transport fatty acids into the organelle for beta-oxidation, so decreases in their plasma levels might reflect increased mitochondrial utilisation [80]. Here, serum levels of both LCFA and LCFA-carnitines were lower in UQ compared to control women, consistent with increased rates of tissue fatty acid utilisation in the UQ group. Such changes can occur in the glucose-sparing fuel economy that emerges in diabetes [81]. Preferential fatty-acid utilisation may contribute to systemic hyperglycemia as recognised long ago [82]. Our data indicate that such utilisation begins much earlier in the pathogenic process than hitherto recognised. The lowering of LCFA and LCFA-carnitines coincided with a small increase in fasting plasma glucose in the UQ group, consistent with substitution of LCFA for glucose in mitochondrial oxidation. Perturbations in LCFA metabolism have been implicated in the pathogenesis of beta-cell damage in diabetes [83]; the early onset of altered LCFA here may lead to or cause beta-cell dysfunction/ damage [57]. Acyl carnitine levels were elevated in pregnant women who went on to develop pre-eclampsia [34]. By contrast, this pattern is no longer evident in the UQ/GDM comparison, where LCFA tended to be higher, probably consistent with their impaired mitochondrial oxidation, typical of insulin resistance in the former (and fatter) GDM group.
Another complex metabolic alteration change more prominent in the UQ/GDM comparison is a tendency to increased numbers of SCFA and SCFA-metabolites. Elevated SCFA and SCFAmetabolites suggest their defective utilisation, as in diabetes [81], again occurring earlier than hitherto realised. Shikimate 3phosphate an obligatory intermediate in the anabolic pathway for biosynthesis of the essential aromatic amino acids, is potentially a microbial metabolite not produced in human cellular metabolism [84]. Some SCFA-metabolites identified may originate from microbial biosynthesis. The identification of microbial metabolites in human plasma with possible links to defective glucoregulation could point to between-group differences in their production by gut microflora and/or uptake from the gut. Other identified metabolic features, as in terpenoid/quinones and teasterone/ typhasterol may be of plant origin, consistent with possible differences in dietary intake and/or uptake from the gut.
We also found ( Table 1) that significantly lower circulating adiponectin levels occurred before measurable alterations in insulin or leptin levels. Adiponectin deficiency occurs from infancy, as found in the children of this cohort [85] and may influence GDM [28] and T2DM [86][87][88]. It is associated with defective glycosylation and functionality, such as impaired ability to stimulate hepatic or muscle mitochondrial fatty acid oxidation via AMP kinase [58,89]. Adiponectin deficiency could provide a central link between perturbed phosphatidylcholine metabolism and mitochondrial lipid utilisation here. However, whether changes in production/secretion and/or signalling of known hormones including adiponectin really antedate or rather result from the described metabolic changes remains uncertain. It is certainly known that adiponectin deficiency can cause these changes but together with the exact nature and origin of the adiponectin deficiency observed here, requires further longitudinal study.
In summary, we identified here a rather consistent pattern of metabolic perturbations in groups of women whose diabetes risk was stratified a priori by differences in their degree of glucoregulatory impairment during a previous pregnancy. The data point to a time-line in the molecular pathologies ultimately leading to type 2 diabetes; the changes found in the control/UQ comparison likely precede those in the UQ/GDM comparison (e.g. perturbed plasma phospholipids and altered lipoprotein metabolism). A second early alteration was the relative fall in plasma LCFA and LCFA-carnitines, along with minor increases in fasting plasma glucose and HbA 1c levels. Those are consistent with a glucose-sparing mitochondrial fuel economy, related to the increased abdominal circumference in the UQ and GDM groups.
Many changes occurred in clusters of metabolite classes, for example phospholipids, lysophospholipids, LCFA, LCFA-carnitines, and SCFA/SCFA-metabolites, pointing to mechanisms that affect large subsets of these metabolite classes (e.g. transcription factors), long before the emergence of overt disease. Differences in relative timings of activation in different potential pathways to the onset/progression of T2DM pathogenesis were also observed. Modified lysophospholipid metabolism possibly implies elevated pro-inflammatory stress; lowered LCFA/LCFA-carnitine levels are consistent with early metabolic fuel substitution leading to preferential mitochondrial oxidation of LCFA as opposed to glucose, providing an early hyperglycemic stimulus; a widespread increase in SCFA/SCFA-metabolites suggest potential early defects in their generation and/or defective mitochondrial utilisation.
Finally, we found early adiponectin deficiency which may initiate or contribute to several of the metabolic disturbances, The results point to a probable defect in adipose tissue regulation contributing to the initiation of T2DM pathogenesis; further characterisation of the early changes in adiponectin synthesis and post-translational modifications and its causes will be useful. Our current conclusions are reminiscent in several respects of those from a recent study of the antecedents of type 1 diabetes wherein dysregulation of lipid and amino acid metabolism preceded islet autoimmunity in children who later progressed to overt disease [37].
Our study paves the way for targeted investigation of the pathogenic biochemical pathways that lead to or cause type 2 diabetes and more effective prevention and therapy [90], notably of blood vessel damage. Further longitudinal studies of diabetes development as we are doing here will be needed for assessing those at risk in general populations. Our study highlights the important role of metabolic profiling in discovery studies related to diabetes. Although metabolite identifications are not definitive they provide mechanistic information to guide further targeted studies. The major perturbations in this hypothesis-generating stage affected large subsets of metabolite classes showing covariation between metabolites. Therefore, no corrections for multiple comparisons were applied. Finally, whether these patterns of metabolic derangements after prior GDM may lead to or cause the T2DM in general populations needs testing.