Changing Metabolic Signatures of Amino Acids and Lipids During the Prediabetic Period in a Pig Model With Impaired Incretin Function and Reduced β-Cell Mass

GIPR^dn transgenic pigs have been generated as described previously (5). GIPR^dn transgenic pigs and age-matched controls were studied at age 2.5 or 5 months. All animals were housed in single pens under controlled conditions and had free access to water. Animals were fed a commercial diet (Supplementary Table 1) once daily and were fasted for 18 h prior to the glucose tolerance tests (GTTs). Physiological data of GIPR^dn transgenic pigs and controls are summarized in Table 1.

Intravenous GTTs (IVGTTs) were accomplished as described previously (5). Blood was collected at −10, 0, 1, 3, 5, 10, 20, 30, 40, 60, and 90 min relative to the glucose load. Blood samples were placed on ice immediately after collection, centrifuged at 1,500g for 15 min, and plasma was stored at −80°C. Plasma glucose levels were determined using an AU 400 autoanalyzer (Olympus, Hamburg, Germany). Plasma insulin levels were measured using a porcine insulin radioimmunoassay kit (Millipore, Billerica, MA). The homeostasis model assessment (HOMA) of β-cell function index (HOMA-β) and HOMA of insulin resistance index (HOMA-IR) were calculated as previously described (12). The insulinogenic index was calculated as the ratio of the increment of plasma insulin (µU/mL) to glucose concentration (mg/dL) 30 min relative to the oral glucose load (7).

Metabolites were quantified at the Helmholtz Center Munich using the AbsoluteIDQ p150 Kit (Biocrates Life Sciences AG, Innsbruck, Austria) covering amino acids (14), hexoses (H), free carnitine (C0), acylcarnitines (40), glycerophosphocholines (92), and sphingolipids (15). Leu and Ile are summarized as xLeu. Absolute quantification of amino acids, acylcarnitines, and hexoses in the biological sample was achieved by reference to internal standards labeled with stable isotopes. The semiquantitative assessment of lipids was achieved by several nonlabeled internal standards in every group (13). Concentrations of metabolites are reported in µmol/L unless indicated otherwise.

Total RNA was isolated from shock-frozen liver tissue samples from 5-month-old male GIPR^dn transgenic and control animals (n = 3 per group) using a tissue homogenizer (Silent Crusher M, Heidolph Instruments, Schwabach, Germany) and the Trizol method. Hybridization probes were prepared using the WT Expression (Ambion, Austin, TX) and Terminal Labeling Kits (Affymetrix, Santa Clara, CA) according to the manufacturers’ instructions. Labeled probes were hybridized to Affymetrix custom Gene ST arrays (SNOWBALLs520824F), which were scanned using an Affymetrix GeneChip Scanner 3000. Raw data were normalized using robust multiarray analysis (14), and the gene list was ranked by differential expression and screened for overrepresentation of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways using gene set enrichment analysis (GSEA) (15).

Metabolite data were statistically evaluated by ANOVA (Linear Mixed Models; PROC MIXED; SAS 8.2), taking the fixed effects of genotype (GT; transgenic, control), age (2.5 months, 5 months), time (relative to glucose application), and the interaction GT × Age, as well as the random effect of animal into account (16). To account for multiple testing, only P values ≤0.0001 were considered significant in the initial screening for significant effects. If specific hypotheses were tested, P values <0.05 were considered significant if not indicated otherwise. Correlation analyses were performed using PROC CORR (SAS 8.2).

Of the 14 plasma amino acids quantified, 6 (His, Trp, Orn, Ser, Phe, and xLeu) were significantly influenced by age and 7 (Phe, Orn, Val, xLeu, His, Arg, and Tyr) by the interaction GT × Age (i.e., the effect of GT was different in younger vs. older animals). In addition, the plasma concentrations of 6 amino acids (Arg, Orn, Phe, Ser, Tyr, and xLeu) were significantly affected by time (Supplementary Table 2) (i.e., their concentrations changed during the course of the IVGTT) (Supplementary Fig. 1). To account for this fact, least squares means (LSMs) and SEs were calculated for GT × Age (Fig. 1). Consistently, the plasma concentrations of Phe, Orn, Val, xLeu, His, Arg, and Tyr were increased in 2.5-month-old but decreased in 5-month-old GIPR^dn transgenic pigs versus age-matched controls. The concentration of Trp was not different between the younger transgenic and control pigs but was significantly reduced in the older transgenic animals (Fig. 1).

To evaluate previously proposed hypotheses (17) of an increased turnover of glucogenic amino acids during prediabetic stages, the sum of the main glucogenic amino acids Ser + Gly was tested. While in 2.5-month-old animals there was no difference between the transgenic and the control group, 5-month-old transgenic pigs displayed a significantly reduced level of these glucogenic amino acids (LSM: 737 vs. 841 μmol/L in controls; GT × Age: P = 0.0043).

Amino acid profiling can potentially provide surrogate information on metabolic aberrations that may contribute to vascular endothelial dysfunction. Increase of reactive oxygen species, as observed in various diabetic conditions and models (18), may lead to consumption of tetrahydrobiopterin (BH4), the cofactor for the enzymes phenylalanine hydroxylase (EC 1.14.16.1) and nitric oxide synthase (EC 1.14.13.39). In addition, increased arginase activity has been reported in diabetic conditions and may contribute to vascular endothelial dysfunction by decreasing l-arginine availability to endothelial nitric oxide synthase (19). Potential aberrations in these enzymes were addressed by calculating Phe/Tyr and Orn/Arg ratios as indirect markers for phenylalanine hydroxylase and arginase activities. While both ratios were significantly affected by GT × Age (P < 0.0001), they did not reveal any sign of endothelial dysfunction in our model (not shown).

Different sphingolipid metabolites (e.g., ceramide) have been suggested to potentially be important regulators of β-cell survival, proliferation, and function (20,21). Since disturbed cellular sphingomyelin (SM) biosynthesis and metabolism might partially also be reflected in extracellular fluids, we included sphingomyelins in our screening panel. Sums of all SM (Fig. 2A) and hydroxysphingomyelin (SM-OH) (Fig. 2B) concentrations that have been measured (Supplementary Table 3) were calculated. In controls, circulating SM concentrations increased significantly with age (P = 0.0002); this increment was only marginal in GIPR^dn transgenic pigs (GT × Age: P = 0.0393). This was reflected by reduced concentrations of SM C24:0 (−17%), SM C24:1 (−9%), SM C16:0 (−15%), and SM C16:1 (−19%) in 5-month-old GIPR^dn transgenic pigs versus age-matched controls (Supplementary Table 3).

SM-OH concentrations were one order of magnitude lower than SM concentrations and were higher in 2.5-month-old but lower in 5-month-old GIPR^dn transgenic pigs compared with age-matched controls (GT × Age: P = 0.0004). The concentrations of SM (OH) C22:1 (−18%) and SM (OH) C22:2 (−10%) were most prominently reduced in 5-month-old GIPR^dn transgenic pigs.

Aberrant diacylglycerol levels have been reported in the context of prediabetes (22,23). Although the plasma level of PC aa C42:0 was increased in 2.5-month-old GIPR^dn transgenic pigs, decreased concentrations of PC aa C42:0, PC aa C42:4, and PC aa C42:5 were observed in 5-month-old transgenic animals when compared with age-matched controls (GT × Age: P ≤ 0.0001) (Table 2). In addition, the plasma concentrations of a number of diacylglycerols were age dependent, with higher concentrations in 2.5- versus 5-month-old animals (Supplementary Table 4).

The plasma concentrations of several ether phospholipids were increased in 2.5-month-old GIPR^dn transgenic pigs, whereas decreased concentrations of specific ether phospholipids were measured in 5-month-old transgenic pigs compared with age-matched controls (GT × Age: P ≤ 0.0001) (Table 2). In addition, a number of ether phospholipid concentrations were affected by age (Supplementary Table 5).

Specific plasma lysophosphatidylcholine concentrations were found to be decreased in 2.5-month-old (lysoPC a C16:1 and lysoPC a C20:3) or 5-month-old (lysoPC a C18:0) GIPR^dn transgenic pigs versus age-matched controls (GT × Age: P ≤ 0.0001) (Table 2). A significant effect of age was observed for plasma concentrations of several lysophosphatidylcholines (Supplementary Table 6). The plasma concentration of lysoPC a C18:0 was significantly (P < 0.0001) affected by the intravenous glucose load (Supplementary Fig. 2A).

Accumulation of plasma acylcarnitines was reported to be associated with absolute or relative insulin deficiency (24,25). Irrespective of age, the plasma concentration of glutaconyl-l-carnitine (C5:1-DC) was significantly increased in GIPR^dn transgenic pigs compared with controls (GT: P < 0.0001). In addition, a trend of increased levels of hydroxypropionyl-l-carnitine (C3-OH; GT: P = 0.0003), valeryl-l-carnitine (C5; GT: P = 0.0036), hexenoyl-l-carnitine (C6:1; GT: P = 0.0006), and octenoyl-l-carnitine (C8:1; GT: P = 0.0018) was detected in 2.5- and 5-month-old GIPR^dn transgenic pigs compared with age-matched controls (Supplementary Table 7). The plasma concentration of butyryl-l-carnitine (C4) was reduced in 2.5-month-old but increased in 5-month-old GIPR^dn transgenic pigs compared with controls (GT × Age: P < 0.0001). The plasma concentration of propionyl-l-carnitine (C3) was significantly influenced by age (P < 0.0001), with higher levels in 2.5- compared with 5-month-old animals. In addition, there was a trend of higher C3 levels in 2.5-month-old GIPR^dn transgenic pigs versus age-matched controls, while this difference was not present in the 5-month age-group (GT × Age: P = 0.0010). The plasma concentration of decanoyl-l-carnitine (C10) was significantly influenced by the glucose infusion during the IVGTT (Supplementary Fig. 2B).

The differential metabolite profiles observed in our study could arise as the result of acute glucose challenge (intravenous glucose infusion) or might represent chronic alterations in glucose metabolism within the prediabetic state. To differentiate between these possibilities, the data analysis was restricted to a subset of metabolites and time points. Only the parameters that were significantly influenced by the interaction GT × Age and only the values from fasted plasma samples (10 min and 0 min before glucose infusion) were taken into account.

It is important that a significant effect of the interaction GT × Age was confirmed for 16 out of 26 metabolites, including all metabolite classes except acylcarnitines. This observation suggests that a significant proportion of the metabolite changes observed in our study are indicative of a chronic alteration of glucose metabolism, rendering them promising candidate biomarkers to monitor the prediabetic period. These metabolites are summarized in Table 3.

Correlations were calculated by using data from the 5-month-old animals, since GIPR^dn transgenic pigs exhibited a 35% reduction in total β-cell volume at this stage. Metabolic parameters listed in Table 3 were included in the analysis, and all of those depicted in Fig. 3 correlated positively with total β-cell volume. Furthermore, Arg and the majority of plasma lipids selected correlated positively with parameters of insulin release. A negative correlation was revealed between lysoPC a C18:0 and, less clearly, PC ae C40:3, PC aa C42:4, and Val and area under the curve (AUC) glucose, in particular with AUC glucose (10–60).

To evaluate if some of the plasma metabolite alterations were accompanied by expression changes of targets in corresponding metabolic pathways, we performed a holistic transcriptome study of liver samples from the three male GIPR^dn transgenic pigs and the three control boars in the 5-month age-group. GSEA was used as a tool, enabling sensitive detection of gene expression changes within specific pathways (26). In accordance with several plasma metabolite alterations, GSEA indeed revealed liver transcript changes to be significantly enriched in a number of pathways already predicted by metabolite alterations. Those include Gly, Ser, and Thr metabolism, sphingolipid metabolism, and biosynthesis of unsaturated fatty acids (Table 4). A representative GSEA plot for genes involved in sphingolipid metabolism and the corresponding pathway with the observed transcriptome changes are shown in Fig. 4.

A number of targeted and nontargeted metabolomic studies have been performed to address metabolic correlates of clinically manifest diabetes (reviewed in Lucio et al. and Suhre et al.) (27,28). In addition, several studies performed metabolomic screens addressing classical risk factors, such as obesity, or pathogenetic mechanisms, such as insulin resistance. For instance, a monozygotic twin study screened for metabolomic correlates of acquired obesity and detected primarily increased levels of lysophosphatidylcholines, while ether phospholipids were decreased (29). Impaired insulin sensitivity in nondiabetic human subjects was found to be associated with altered metabolite profiles in lipid metabolic pathways, steroid hormone biosynthesis, and bile acid metabolism (27). Metabolomic profiling of obese versus lean humans revealed that serum concentrations of eight amino acids (Ala, Val, Leu/Ile, Phe, Tyr, Glu/Gln, Asp/Asn, and Arg) were markedly increased, while Gly was decreased, and defined a branched-chain amino acid (BCAA)-related metabolomic signature in obese subjects (30). A recent prospective study found the plasma concentrations of three BCAAs (Leu, Ile, and Val) and two aromatic amino acids (Phe and Tyr) to be significantly associated with the future risk of diabetes (31). However, the patients included in the latter study were characterized by an already increased BMI and HOMA-IR. Human subjects with even earlier stages of prediabetes characterized by β-cell dysfunction are only rarely diagnosed and, thus, difficult to recruit. Our genetically engineered pig model provides a unique time window to investigate the course of gradual β-cell mass reduction in relation to changes of intermediary metabolism. In contrast to human patients, the data on β-cell mass in pigs can also be obtained as quantitative information. A characteristic biochemical feature of this animal model is the progressive reduction of the insulinogenic index with increasing life span. This parameter is correlated with impaired first-phase insulin release (7), which occurs in very early stages of β-cell dysfunction and is an important factor in the development of postprandial hyperglycemia (32).

Changes in plasma amino acid concentrations have been originally described in insulin-resistant obese patients (reviewed in Felig et al. and Muoio and Newgard) (33,34). In particular, elevated levels of BCAAs have been confirmed in a number of subsequent studies of obese and/or insulin-resistant patient collectives (30,35,36) and in clinically manifest diabetes (28). Of interest, plasma levels of Val, Leu, Ile, Tyr, and Phe—known to be associated with diabetes risk in human subjects with increased BMI and reduced insulin sensitivity (31)—and those of Orn and Arg were elevated in 2.5-month-old GIPR^dn transgenic pigs, which show only a reduced incretin effect but are neither obese nor insulin resistant according to HOMA-IR (Table 1). This observation is remarkable because it demonstrates that a mild chronic disturbance of glucose homeostasis caused by reduced GIP function might be sufficient to result in a plasma amino acid profile suggestive for an increased risk to develop diabetes even in the absence of obesity. The elevated levels of seven amino acids in 2.5-month-old GIPR^dn transgenic pigs could be related to the distinct pattern of insulin secretion in this model. After an oral glucose load, initial insulin secretion is delayed, although the overall amount of insulin secreted during the 120-min period of the test reaches levels comparable with controls (5). In this respect, our model reflects the condition of impaired glucose tolerance preceding type 2 diabetes in human subjects where plasma insulin levels were reduced at 30 min after glucose ingestion but normal or increased at 120 min (37). The bell-shaped curve of insulin secretion in GIPR^dn transgenic pigs implicates a longer phase of insulin above a level that could inhibit gluconeogenesis and, thus, the use of glucogenic amino acids. A further possible explanation is that in comparison with the steep and fast increase of insulin secretion after oral glucose load, the slower increase in plasma insulin levels seen in GIPR^dn transgenic pigs could be a weaker anabolic signal, leading to less efficient amino acid uptake in peripheral tissues and clearance from the circulation (38,39).

It is surprising that initially elevated plasma levels of these amino acids (at 2.5 months) were reduced in 5-month-old GIPR^dn transgenic pigs, which show a more pronounced prediabetic phenotype (Table 1). The reduction in pancreatic β-cell mass and impaired insulin secretion after oral glucose load (5) and probably also after food intake might explain this effect. Of interest, the plasma concentrations of Arg, His, Phe, Val, and xLeu were correlated with total β-cell mass in 5-month-old pigs. During evolution toward diabetes, a continual mitigation of the inhibiting effect of insulin on hepatic gluconeogenesis would be expected to result in a higher demand and consumption of glucogenic amino acids. This hypothesis is indeed supported by the fact that the sum of the major glucogenic amino acids (Gly + Ser) was significantly reduced in plasma of 5-month-old GIPR^dn transgenic pigs. In addition, the liver transcriptome changes in 5-month-old pigs, showing enrichment of differentially expressed genes in corresponding pathways of amino acid metabolism, might also contribute to the observed alterations in plasma amino acid concentrations.

Specific amino acid ratios (Tyr to Phe and Orn to Arg) did not reveal signs of endothelial dysfunction as expected in view of our model. The GIPR^dn transgenic pig mimics different stages of the prediabetic period but is not yet manifesting with overt clinical signs.

Sphingolipids are essential components of plasma membranes and endosomes and are critically involved in cell surface protection, protein and lipid transport and sorting, and cellular signaling cascades (40). Reduced plasma SM concentrations were recently observed in cynomolgus monkeys spontaneously developing hyperglycemia and hyperlipidemia (41).

We thus questioned if early stages of altered glucose homeostasis and β-cell dysfunction are associated with altered plasma concentrations of SM. It is interesting that GIPR^dn transgenic pigs did not show a comparable age-dependent increase in saturated SM species as observed in controls. In the more progressed prediabetic condition of 5-month-old GIPR^dn transgenic pigs, we found significant reductions in plasma levels of nearly all of the SM metabolites quantified (Supplementary Table 3). This observation is consistent with a recent study exploring specific sphingolipids as potential biomarkers in two different mouse models of insulin-dependent diabetes, streptozotocin-administered mice and Ins2^Akita mutant mice (42). In this study, plasma SM-C24:1 levels were consistently reduced in both models, and there was a trend of reduced SM C24:0 levels in the streptozotocin-administered animals. While both mouse models are characterized by a massive loss of β-cells, 5-month-old GIPR^dn transgenic pigs show only a modest disturbance of glucose homeostasis and a 35% reduction in β-cell mass. The concordance of SM changes in prediabetic GIPR^dn transgenic pigs with two mouse models with clinically manifest diabetes suggests that reductions in specific SMs could signify a proximal event in the cascade toward diabetes rather than being a consequence of manifesting disease. The fact that reduction of SM synthesis increases insulin sensitivity (43) raises the possibility that decreased SM levels in 5-month-old GIPR^dn transgenic pigs represent an adaptive response to reduced insulin secretion. It is interesting that there was a trend of a smaller HOMA-IR in 5-month-old GIPR^dn transgenic versus control pigs (see Table 1). The observed liver transcriptome changes related to SM metabolism (see Table 4 and Fig. 4) would support this concept. However, SM synthase 1 deficient mice show impaired insulin secretion (44). Thus, reduced SM levels, if present in β-cells, may also contribute to the impaired insulin secretion observed in 5-month-old GIPR^dn transgenic pigs.

Since plasma lipid profiles were found to be altered in obesity (29), an important risk factor for diabetes, we quantified the plasma concentrations of 38 diacylglycerols, 39 ether phospholipids, and 15 lysophosphatidylcholines in GIPR^dn transgenic and control pigs. While ether phospholipids can serve as endogenous antioxidants (45,46), increased levels of lysophosphatidylcholines have been described in the context of proinflammatory and proatherogenic conditions (reviewed in Pietiläinen et al.) (29). Acquired obesity was found to be associated with increases in lysophosphatidylcholines but decreases in ether phospholipids (29). Moreover, the levels of plasmalogens (ether lipids where the first position of glycerol binds a vinyl residue [from a vinyl alcohol] with the double bond next to the ether bond) were found to be decreased by ∼50% in LDLs from type 2 diabetic patients (47).

Of interest is the finding that in 2.5-month-old GIPR^dn transgenic pigs, the concentrations of PC aa C42:0 and of several ether phospholipids were increased compared with age-matched controls. In contrast, decreased plasma concentrations of specific diacylglycerols and ether phospholipids were observed in 5-month-old GIPR^dn transgenic pigs. The latter may be interpreted as an early stage of loss of metabolic control associated with reduced antioxidant potential by ether phospholipids. Furthermore, as observed for the plasma amino acid concentrations, the reversed plasma concentrations of some glycerophospholipids in 2.5- and 5-month-old GIPR^dn transgenic pigs are apparently indicative of progression within the prediabetic state.

Lysophosphatidylcholines were not elevated in GIPR^dn transgenic pigs. This is consistent with the fact that the metabolic alterations observed in 2.5- and 5-month-old GIPR^dn transgenic pigs are relatively mild and do not cause increased lipid peroxidation or other changes that might result in increased plasma lysophosphatidylcholine levels as observed in obesity and clinical diabetes (29).

Carnitine and its esters procure fatty acids for β-oxidation and therefore prevent fatty acid accumulation. Absolute or relative insulin deficiency results in increased and incomplete β-oxidation and accumulation of plasma acylcarnitines (24,25). In a recent study exploring metabolomic differences between obese and lean humans and their contribution to insulin resistance (30), 4 out of 37 acylcarnitine species measured were found to be increased in the plasma of obese subjects: propionyl-l-carnitine (C3), valeryl-l-carnitine (C5), hexenoyl-l-carnitine (C6), and octenoyl-l-carnitine (C8:1). Increased levels of short-chain acylcarnitines (C2, C6, and C8) also were observed in the diabetic state, and a positive correlation with plasma HbA_1c suggested a link with the severity of the disease (24). It is interesting that some short-chain acylcarnitines were also increased in GIPR^dn transgenic pigs, even at the age of 2.5 months. This observation indicates that in the absence of obesity, insulin resistance, or other severe disturbances of glucose homeostasis, short-chain acylcarnitines are sensitive parameters whose elevation suggests early stages and progression of prediabetes.

Metabolite profiles of 2.5- and 5-month-old GIPR^dn transgenic pigs representing an early prediabetic state characterized by a reduced incretin effect and progressive reduction of glucose tolerance, insulin secretion, and pancreatic β-cell mass in part show similar changes as previously detected in obese, insulin-resistant, or diabetic patients. Thus, metabolite concentration alterations detected even in young GIPR^dn transgenic pigs appear to be sensitive indicators of very early stages of prediabetes. It is important that >50% of the metabolites altered in GIPR^dn transgenic pigs versus controls in an age- and disease progression–specific manner during an IVGTT were also detected if only fasting plasma samples were analyzed. Therefore, these parameters appear to be associated with long-term glucose control and do not arise only as the consequence of an intravenous glucose bolus. This observation also is important in practical terms because time-consuming GTTs are not essential. The age-dependent inverse alteration of specific plasma amino acid and lipid concentrations (elevated in 2.5-month-old and reduced in 5-month-old animals) going along with the progression of phenotypic changes in GIPR^dn transgenic pigs strongly argues for repeated metabolomic profiling in regular intervals, which is easily achievable in a large animal model. The fact that metabolic signatures observed in fasting plasma samples of 5-month-old GIPR^dn transgenic pigs were 1) significantly correlated with total β-cell mass and 2) associated with liver transcriptome changes in corresponding metabolic pathways support their potential usage as biomarkers.

In summary, our findings document the potential of a genetically designed large animal model resembling early stages of prediabetes in the absence of confounding obesity or insulin resistance for the discovery of candidate biomarkers that may serve as indicators of early stages of prediabetes.

This study was supported by the Deutsche Forschungsgemeinschaft (GRK 1029), the Federal Ministry of Education and Research (Leading-Edge Cluster m⁴–Personalized Medicine and Targeted Therapies), the Bayerische Forschungsstiftung (492/02), the Diabetes Hilfs- und Forschungsfonds Deutschland, a grant from the German Federal Ministry of Education and Research (BMBF) to the German Center Diabetes Research (DZD e.V.), BMBF Grant 03IS2061B (project Gani_Med), and BMBF Grant 0315494A (project SysMBo).

No potential conflicts of interest relevant to this article were reported.

S.R. researched data, contributed to discussion, and wrote, reviewed, and edited the manuscript. W.R.-M. and C.P. researched data, contributed to discussion, and reviewed and edited the manuscript. S.K. and H.B. performed the microarray study, including the bioinformatics analysis. J.A., B.G., and K.S. contributed to discussion and reviewed and edited the manuscript. A.A.R. and E.W. contributed to discussion and wrote, reviewed, and edited the manuscript. S.R. and E.W. are the guarantors of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

The authors thank Julia Scarpa (Institute of Experimental Genetics, Genome Analysis Center, Helmholtz Zentrum München), Katharina Sckell (Institute of Experimental Genetics, Genome Analysis Center, Helmholtz Zentrum München), and Arsin Sabunchi (Institute of Bioinformatics and Systems Biology, Helmholtz Zentrum München) for metabolomics measurements performed at the Helmholtz Zentrum München, Genome Analysis Center, Metabolomics Core Facility.