Metabolomic discovery of biomarkers of type 2 diabetes (T2D) risk may reveal etiological pathways and help to identify individuals at risk for disease. We prospectively investigated the association between serum metabolites measured by targeted metabolomics and risk of T2D in the European Prospective Investigation into Cancer and Nutrition (EPIC)-Potsdam (27,548 adults) among all incident cases of T2D (n = 800, mean follow-up 7 years) and a randomly drawn subcohort (n = 2,282). Flow injection analysis tandem mass spectrometry was used to quantify 163 metabolites, including acylcarnitines, amino acids, hexose, and phospholipids, in baseline serum samples. Serum hexose; phenylalanine; and diacyl-phosphatidylcholines C32:1, C36:1, C38:3, and C40:5 were independently associated with increased risk of T2D and serum glycine; sphingomyelin C16:1; acyl-alkyl-phosphatidylcholines C34:3, C40:6, C42:5, C44:4, and C44:5; and lysophosphatidylcholine C18:2 with decreased risk. Variance of the metabolites was largely explained by two metabolite factors with opposing risk associations (factor 1 relative risk in extreme quintiles 0.31 [95% CI 0.21–0.44], factor 2 3.82 [2.64–5.52]). The metabolites significantly improved T2D prediction compared with established risk factors. They were further linked to insulin sensitivity and secretion in the Tübingen Family study and were partly replicated in the independent KORA (Cooperative Health Research in the Region of Augsburg) cohort. The data indicate that metabolic alterations, including sugar metabolites, amino acids, and choline-containing phospholipids, are associated early on with a higher risk of T2D.
Type 2 diabetes (T2D) is characterized by impaired insulin sensitivity of several tissues and inadequate insulin secretion from β-cells (1). A detailed understanding of the pathophysiology of T2D is a prerequisite for the development of preventive strategies. In particular, the identification of early metabolic alterations is promising in the study of etiological pathways and may further help to identify high-risk individuals. A number of biomarkers have been proposed as indicators for the estimation of T2D risk, such as fasting plasma glucose and glycated hemoglobin A1c (HbA1c) (2), triglycerides (3), HDL cholesterol (4), inflammatory markers (5), adiponectin (5,6), liver enzymes (7), and fetuin-A (8). However, most biomarkers fail to grasp the complexity of T2D etiology (3). Design and advancement of high-throughput analytical techniques determined the emergence of metabolomics, which is the simultaneous study of numerous low-molecular weight compounds, namely metabolites. Metabolites represent intermediates and end products of metabolic pathways that reflect more rapidly physiological dysfunctions than current biomarkers and, thus, may mirror earlier stages of T2D (9). Cross-sectional studies have linked alterations in metabolic profiles with obesity (10), glucose tolerance (11), and prevalent diabetes (12–14). The most prominent metabolic shifts involved blood acylcarnitines and branched-chain amino acids (BCAAs). Recently, observations from a prospective study found that a set of five amino acids was predictive for T2D, representing pioneering work in the emerging field of systems epidemiology (15).
In the current study, we investigated whether a targeted metabolomic approach involving a broader spectrum of metabolites and a larger number of study participants may help to identify metabolites associated with the risk of T2D and the mechanisms involved. Therefore, we profiled 163 serum metabolites in originally healthy individuals who were consecutively followed up for incident T2D in two large-scale prospective cohort studies in Germany and studied cross-sectional relationships of the identified metabolites with insulin sensitivity and secretion in precisely phenotyped participants. In addition, we evaluated the usefulness of the metabolites for T2D risk prediction compared with the German Diabetes Risk Score (DRS) (16) and established biomarkers.
RESEARCH DESIGN AND METHODS
European Prospective Investigation into Cancer and Nutrition-Potsdam study.
The European Prospective Investigation into Cancer and Nutrition (EPIC)-Potsdam is part of the ongoing multicenter EPIC study and comprises 27,548 participants from the general population in the area of Potsdam in eastern Germany who were mainly 35–65 years of age at time of recruitment between 1994 and 1998 (17). At baseline, participants underwent an examination by qualified staff, including medical history, blood pressure measurement, and anthropometry (18). Participants also completed sociodemographic and lifestyle questionnaires and a validated food frequency questionnaire. In addition, 30 mL blood were drawn (random sampling) and immediately processed (19). Only participants with morning appointments were asked to fast overnight. Blood was fractionated into serum, plasma, buffy coat, and erythrocytes; aliquotted into straws of 0.5 mL each; and stored in tanks of liquid nitrogen at −196°C until analysis. Besides metabolomic profiling, other biomarkers have been measured in baseline blood samples as described previously (7,20,21). Every 2–3 years, follow-up questionnaires were sent to participants to identify incident cases of T2D, with response rates of ∼95% (22,23). Once a participant was identified as a potential case, disease status was further verified with medical records, including the correct diagnosis (International Classification of Diseases, 10th revision, E11, non–insulin-dependent diabetes), the date of the diagnosis, and the means of diagnosis confirmation. This verification was achieved by sending a standard inquiry form to the treating physician. Consent was obtained from all study participants a priori, and the study was approved by the ethics committee of the Medical Society of the State of Brandenburg.
We constructed a case-cohort study within EPIC-Potsdam, including all incident cases of T2D of the full cohort identified up to 31 August 2005 (n = 849, mean follow-up 7 years), and a subcohort (n = 2,500) randomly drawn from the EPIC-Potsdam study population. Because the subcohort was representative of the full cohort, it included 2,415 noncases and 85 cases of incident T2D (i.e., internal cases). The remaining 764 cases did not belong to the subcohort (external cases). By randomly selecting the subcohort and using the appropriate statistic for this study design, the biomarkers only needed to be measured in the case-cohort sample; however, the results are expected to be generalizable to the full cohort (7). The case-cohort design was previously chosen based on its advantages, including a reduced chance of selection bias for the control group (24).
For the present analysis, we further excluded participants with prevalent T2D at baseline (n = 110), with missing or nonverified data on incident or prevalent T2D (n = 13), with missing blood samples or biomarker measurements (n = 64), and with missing covariate information (n = 80). Thus, the analytical sample included 2,282 individuals of the subcohort and 800 individuals with incident T2D.
Cooperative Health Research in the Region of Augsburg study.
The Cooperative Health Research in the Region of Augsburg (KORA) study consists of population-based surveys and follow-up periods in the area of Augsburg in southern Germany. A total of 4,261 individuals between 25 and 74 years of age participated in the S4 survey between 1999 and 2001 (25). In the KORA cohort, blood was drawn into serum gel tubes after a fasting period of at least 8 h. Blood samples were rested for coagulation for 30 min at room temperature; serum was obtained by centrifugation at 2,750g at 15°C for 10 min and stored in a freezer at −80°C. A total of 3,080 individuals took part in the F4 follow-up survey during the years 2006–2008 (26). The identification of incident T2D was based on an oral glucose tolerance test (OGTT) or a validated physician diagnosis (27). All KORA participants gave written informed consent, and the KORA study was approved by the ethics committee of the Bavarian medical association.
A subcohort of 876 S4 participants 55–74 years of age without T2D at baseline and with metabolomics data available was included in the current study. Of them, 91 developed incident T2D during the 7-year follow-up.
Tübingen Family study for T2D.
The Tübingen Family (TüF) study is in an ongoing investigation of the pathophysiology of T2D in southern Germany (28). Individuals meeting at least one of the following criteria were included in the study: a family history of T2D, a BMI >27 kg/m2, and previous impaired glucose tolerance or gestational diabetes mellitus. They were considered healthy according to a physical examination and routine laboratory tests. Written informed consent was obtained from all participants, and the medical ethics committee of the University of Tübingen approved the protocol.
All individuals underwent a 75-g OGTT. Venous plasma samples were drawn at 0, 30, 60, 90, and 120 min for plasma glucose, insulin, C-peptide, and metabolomic analyses (minute 0). Insulin sensitivity was calculated from the OGTT (29). The plasma glucose and C-peptide areas under the curve (AUCs) during the OGTT were calculated by applying the trapezoid method. Insulin secretion was calculated from AUCC-peptide/AUCglucose. The present analysis included 76 Caucasians from the TüF study who had measurements of insulin sensitivity and insulin secretion as well as metabolomics data available.
Serum metabolite concentrations.
Serum concentrations of metabolites were determined with the AbsoluteIDQ p150 and p180 Kits (Biocrates Life Sciences AG, Innsbruck, Austria) using the flow injection analysis tandem mass spectrometry (FIA-MS/MS) technique (30). The metabolomic method simultaneously quantified 163 metabolites, including 41 acylcarnitines (Cx:y), 14 amino acids, 1 hexose (sum of six-carbon monosaccharides without distinction of isomers), 92 glycerophospholipids (lyso-, diacyl-, and acyl-alkyl-phosphatidylcholines), and 15 sphingomyelins. To ensure valid measurements, metabolites below the limit of detection (n = 30) and those with very high analytical variance (n = 6) in our samples were excluded, leaving 127 metabolites for the present analysis.
Metabolomic measurements were performed in the Genome Analysis Center at the Helmholtz Zentum München. Sample preparation was done according to the manufacturer’s protocol (Biocrates user’s manual UM-P150) and has been described previously (30). In brief, after centrifugation, 10 μL serum were inserted into a filter on a 96-well sandwich plate, which already contained stable isotope-labeled internal standards. Amino acids were derivated with 5% phenylisothiocyanate reagent. Metabolites and internal standards were extracted with 5 mmol/L ammonium acetate in methanol. The solution was then centrifuged through a filter membrane and diluted with mass spectrometry running solvent. Final extracts were analyzed by FIA-MS/MS, and metabolites were quantified in µmol/L by appropriate internal standards. The method has been validated, and analytical specifications were provided in the Biocrates manual AS-P150. The manufacturer selected the metabolites based on the robustness of their measurements. The uncertainty of the measurements was <10% for most of the metabolites. Regarding accuracy, all included metabolites were found in the range of 80–115% of their theoretical values. The median analytical variance of EPIC-Potsdam samples was a 7.3% within-plate coefficient of variation and a 11.3% between-plates coefficient of variation (31). To account for run-order effects, serum samples were randomly analyzed together, regardless of the case status. We have shown previously that most of the metabolites had moderate to high intraclass correlation coefficients measured in participants over a 4-month period, indicating reasonable reliability of the measurements (31).
Step 1: Identification of metabolites associated with T2D risk in EPIC-Potsdam.
Cox proportional hazards regression with weighting as suggested by Prentice (32) and robust sandwich covariance estimates to account for the case-cohort design were used to calculate multivariable-adjusted hazard ratios as a measure of relative risk (RR) and 95% CI, with age as the underlying time scale from recruitment to study exit (T2D diagnosis or censoring) of each participant. We considered z score–standardized metabolite concentrations (mean 0 [SD 1]) as the exposure variable and calculated a multivariable-adjusted model to select metabolites associated with T2D risk. This model was adjusted for age, sex, alcohol intake from beverages (nonconsumers; women >0–6, 6–12, and >12 g/day; and men >0–12, 12–24, and >24 g/day), smoking (never, former, current ≤20 cigarettes/day, current >20 cigarettes/day), cycling and sports (h/week), education (no degree/vocational training, trade/technical school, university degree), coffee intake (cups/day), red meat intake (g/day), whole-grain bread intake (g/day), prevalent hypertension (yes/no), BMI (kg/m2), and waist circumference (cm). Because the metabolomic approach is exploratory, the P values from Cox regression were corrected for multiple testing (n = 127) using the Bonferroni-Holm procedure (33), and a corrected P < 0.05 (two-sided testing) was considered significant to select metabolites.
We next calculated a model that included the covariates and all the identified metabolites and used stepwise Cox regression to select the independent predictors. To also account for the intercorrelation of some of the metabolites, we conducted a principal component analysis (PCA). In brief, the PCA aggregates the individual metabolites based on their degree of correlation with one another to a smaller number of metabolite factors (principal components). These metabolite factors are extracted in a way that explains the major fraction of the variance of individual metabolites. We included all metabolites associated with T2D risk in the PCA, based the PCA on the correlation matrix of metabolites, and used an orthogonal rotation procedure with the varimax method. We retained two metabolite factors according to the scree test and because they accounted for most of the observed variance. Thus, the proportion of explained variance of factor 1 and factor 2 were 34.2 and 16.1%, respectively. To investigate the association between metabolite factors and T2D risk, we estimated RR and 95% CI across quintiles of metabolite factors, particularly choosing quintiles of metabolite factors to facilitate the interpretation. We also investigated a possible effect modification of sex or fasting status on the association between metabolite factors and T2D risk by including multiplicative interaction terms into the models. We then repeated the PCA to include only fasting blood samples in order to evaluate whether the metabolite factors were different from those obtained from random blood samples. Finally, we calculated hazard rates of T2D during different periods of follow-up and tested whether they were different with a test of heterogeneity (34).
Step 2: Additional analyses, risk prediction, and replication in KORA.
For significant metabolites found in step 1, we calculated multivariable-adjusted models with additional adjustment for blood glucose, HbA1c, HDL cholesterol, and triglycerides. We also adjusted for the amino acids phenylalanine, tyrosine, and isoleucine, which have been found to be associated with T2D risk (35). Using data of the EPIC-Potsdam subcohort, we calculated Spearman partial correlation coefficients between identified metabolites and established T2D biomarkers. Data of the TüF study were used to calculate Spearman partial correlation coefficients between identified metabolites and measures of insulin sensitivity and secretion. We calculated measures of discrimination and calibration in different multivariable-adjusted logistic regression models using the DRS (16) as the reference model and adding established T2D biomarkers and the identified metabolites. Receiver operating characteristic (ROC) AUCs were compared using the DeLong test (36).
The results were replicated in the prospective KORA study, and metabolite factors were recalculated in KORA using the linear factor equations retrieved from the PCA in the EPIC-Potsdam sample. The risk estimates from EPIC-Potsdam and KORA were combined using a meta-analytical approach (37). Power calculation suggested a detectable RR per SD of 1.26 (38). Additionally, multivariable-adjusted RRs of T2D were calculated for the amino acids that were recently identified by Wang et al. (35) and that were also measured in the EPIC-Potsdam study. The statistical analyses were conducted with SAS version 9.2 (SAS Institute, Inc, Cary, NC) and in the R statistical environment (www.r-project.org).
Baseline characteristics of the EPIC-Potsdam sample are presented in Table 1. Of all the metabolites measured using a targeted metabolomic approach, 34 were significantly associated with T2D risk in the EPIC-Potsdam study after correction for multiple testing (Supplementary Table 1). These relations were also independent of relevant dietary and lifestyle factors as well as anthropometry and hypertension. Among these 34 metabolites, 14 were identified to be significantly associated with T2D risk independent of the others (Table 2). Specifically, hexose, phenylalanine, and diacyl-phosphatidylcholines C32:1, C36:1, C38:3, and C40:5 were significantly positively associated with T2D risk, whereas glycine, sphingomyelin C16:1, lysophosphatidylcholine C18:2, and acyl-alkyl-phosphatidylcholines C34:3, C40:6, C42:5, C44:4, and C44:5 were significantly inversely related to T2D risk. Using PCA, we identified two metabolite factors that included multiple metabolites and explained most of their variation (Fig. 1). These metabolite factors showed significant and opposing associations with T2D risk. When comparing extreme quintiles, metabolite factor 1, which mainly contains acyl-alkyl-phosphatidylcholines, sphingomyelins, and lysophosphatidylcholines, was associated with a significant 69% reduced risk of T2D (RR 0.31 [95% CI 0.21–0.44]) (Table 3), whereas metabolite factor 2, consistent of diacyl-phosphatidylcholines, BCAA and aromatic amino acids, propionylcarnitine, and hexose, was associated with a significant 3.82-fold increased risk of T2D (3.82 [2.64–5.52]). When we restricted the PCA to fasting samples (n = 429), very similar metabolite factors could be generated with only minor differences in factor loadings (Supplementary Table 2). We also estimated the joint effects of both metabolite factors by summing them (factor 1 received a negative sign because it was inversely associated with T2D risk) and calculating RR of T2D across quintiles of combined factors. The RR (95% CI) of T2D from quintile 1 to quintile 5 of summed metabolite factors was as follows: 1.0, 1.47 (0.96–2.24), 2.29 (1.54–3.38), 2.67 (1.79–4.0), and 6.69 (4.50–9.96) (P for trend <0.0001).
Adjustment for established T2D biomarkers only marginally affected the magnitude of risk association for most of the metabolites (Table 2). An exception was that the inverse association between acyl-alkyl-phosphatidylcholines and sphingomyelin C16:1 and T2D risk was attenuated and no longer significant after adjustment for blood glucose, HbA1c, HDL cholesterol, and triglycerides. The positive association between hexose and T2D risk was considerably weakened but remained significant after adjustment for plasma glucose. When we also adjusted for phenylalanine, tyrosine, and isoleucine, the associations for metabolite factor 1 were unchanged. The associations for metabolite factor 2, which included these three amino acids, were weakened but remained significant (data not shown).
We further observed that metabolites were linked to established T2D biomarkers (Supplementary Table 3). Specifically, metabolite factor 1, which was inversely associated with T2D risk, was negatively correlated to plasma glucose, HbA1c, and triglycerides and positively related to HDL cholesterol and adiponectin. Metabolite factor 2, which was positively associated with T2D risk, was positively correlated with triglycerides and liver enzymes. Data from the TüF study revealed that acyl-alkyl-phosphatidylcholines, lysophosphatidylcholine C18:2, and glycine were positively associated with insulin sensitivity, whereas hexose and diacyl-phosphatidylcholines were inversely related to insulin sensitivity (Table 4). Furthermore, phenylalanine was positively associated with insulin secretion, whereas hexose, sphingomyelin C16:1, and acyl-alkyl-phosphatidylcholines were inversely related to insulin secretion. The potential of identified metabolites to discriminate between T2D cases and noncases was comparable to that of the DRS (16) (ROC AUC 0.849 and 0.847, respectively, P for difference = 0.838) (Fig. 2). When the metabolites were added to established risk prediction models of T2D, discrimination was slightly but significantly improved up to a ROC AUC of 0.912, and these models were well calibrated (Fig. 2 and Supplementary Table 4). Replication in KORA revealed significant associations with T2D risk for metabolite factor 2 and hexose (Table 3 and Supplementary Table 5). In KORA, similar trends as in EPIC-Potsdam were seen for metabolite factor 1, acyl-alkyl-phosphatidylcholines, glycine, lysophosphatidylcholine C18:2, and sphingomyelin C16:1, with borderline significance. Although the risk estimates for diacyl-phosphatidylcholines were considerably lower in KORA than in EPIC-Potsdam, there was no significant heterogeneity between studies. We also calculated the RR of T2D for the BCAA and aromatic amino acids, which have recently been reported to be associated with T2D risk in the Framingham Offspring cohort and the Malmö Diet and Cancer study, to facilitate the comparison (35). In EPIC-Potsdam, isoleucine, valine, tyrosine, and phenylalanine were positively associated with T2D risk (RR per SD 1.30 [95% CI 1.17–1.43], 1.27 [1.16–1.40], 1.31 [1.18–1.45], 1.35 [1.22–1.49], respectively); leucine was not measured in EPIC-Potsdam. When combining isoleucine, tyrosine, and phenylalanine, the RR of T2D from lowest to highest quartile was 1.0, 1.13 (0.82–1.54), 1.45 (1.07–1.98), and 2.18 (1.62–2.95), respectively (P for trend < 0.0001).
We conducted several sensitivity analyses. In EPIC-Potsdam, a small proportion of the participants (14.3%) had fasted. The proportion of incident T2D cases was equally distributed among fasting and nonfasting participants (26.8% and 26.7%, respectively). Additional adjustment for fasting status did not change the results. Further, we did not observe an effect modification of fasting status on the association between metabolite factors 1 and 2 and T2D risk (P for interaction = 0.115 and 0.688, respectively). We observed no interaction with sex (P = 0.407 and 0.441, respectively), and in both men and women, the risk associations were similar. Hazard rates of T2D in different periods of follow-up were not different (factor 1 P = 0.126, factor 2 P = 0.994), indicating that follow-up time did not affect the association between metabolite factors and risk of T2D. To ensure that the metabolite changes preceded the onset of T2D and were not attributed to prediabetic conditions, we repeated the analysis, excluding all cases of T2D that occurred shortly after the baseline examination during the first 2 years of follow-up (n = 208). The risk associations were slightly lower, but not markedly different.
In this prospective investigation using a targeted metabolomic approach at population level, we found increased concentrations of hexose; phenylalanine; and diacyl-phosphatidylcholines C32:1, C36:1, C38:3, and C40:5 and reduced concentrations of glycine; sphingomyelin C16:1; acyl-alkyl-phosphatidylcholines C34:3, C40:6, C42:5, C44:4, and C44:5; and lysophosphatidylcholine C18:2 to be independently predictive of T2D in EPIC-Potsdam. The results agree with data from cross-sectional studies showing that patients with T2D had increased concentrations of sugar metabolites (13), acylcarnitines (14), and BCAA (13) and reduced concentrations of glycine (12). We were able to further replicate the results of Wang et al. (35), who recently reported that BCAA and aromatic amino acids predicted T2D in the prospective Framingham Offspring cohort and the Malmö Diet and Cancer study. In agreement with Wang et al. (35), we found higher concentrations of phenylalanine, isoleucine, tyrosine, and valine to be associated with increased risk of T2D and glycine to be associated with reduced risk of T2D. However, in the present study, BCAA and aromatic amino acids were linked to each other, and only phenylalanine was independently associated with T2D risk when accounting for the other metabolites. BCAAs may serve as substrates for the glucose-alanine cycle in skeletal muscle. Through alanine aminotransferase–catalyzed transamination reactions, this may result in increased substrate availability for hepatic gluconeogenesis, thereby increasing hepatic glucose production (39). Conversely, glycine is a gluconeogenic amino acid; therefore, reduced serum glycine may also reflect increased gluconeogenesis. Alternative theories suggest that glycine depletion may reflect glutathione consumption driven by oxidative stress (40) or abundance of incompletely oxidized fuels that are excreted as urinary acylglycine conjugates (41–43). The frequently observed increase of BCAA in subjects with insulin resistance is also believed to be the result of reduced activities of key BCAA catabolic enzymes in liver and adipose tissue (44). Furthermore, amino acids may directly cause muscular insulin resistance by disrupting insulin signaling (45).
The positive association between hexose and T2D risk remained significant after adjustment for glucose. This observation could be an artifact from the different methods used to measure hexose and glucose. However, it has to be noted that hexose represented not only glucose but also the sum of all six-carbon monosaccharides. Previous studies have shown that in addition to glucose, fructose levels were elevated in individuals with T2D (12) and that intake of fructose was positively associated with risk of insulin resistance and T2D (46). In insulin-resistant conditions, the body aims to compensate for decreased glucose uptake of peripheral tissues through increased pancreatic insulin secretion (1). However, at the stage of overt insulin resistance, this system will eventually be exhausted, as caused by β-cell dysfunction and, subsequently, insulin secretion decreases (1). Phenylalanine, which was positively correlated to insulin secretion in the present study, may be involved in pathways to compensate early stages of insulin resistance through stimulation of insulin secretion. In contrast, increased hexose concentrations may indicate manifest insulin resistance and defect of β-cells.
It is noteworthy that we observed significant associations between choline-containing phospholipids (i.e., diacyl-, acyl-alkyl-, and lysophosphatidylcholines and sphingomyelins) and T2D risk. Diacyl-phosphatidylcholines consist of glycerol linked to phosphocholine and two fatty acid residues, and removal of one fatty acid produces lysophosphatidylcholines. The corresponding acyl-alkyl-phosphatidylcholines comprise an ether linkage to one alkyl chain and one polyunsaturated fatty acid (47). Sphingomyelins are built of a ceramide core linked to one fatty acid and a phosphocholine or phosphoethanolamine (Fig. 3). Together, these phospholipids make up the main constituent of cellular membranes and may be involved in cellular signal transduction (48). In addition, they represent a major fraction of the human plasma lipidome because they are most abundant in all lipoproteins (49). Diacyl-phosphatidylcholines are particularly essential for hepatic secretion of triglyceride-rich VLDL particles and HDL (48), whereas acyl-alkyl-phosphatidylcholines may act as serum antioxidants to prevent lipoprotein oxidation (50). Their hepatic synthesis requires dietary choline (48). It was previously shown that choline-deficient mice on a high-fat diet showed reduced phosphatidylcholine biosynthesis and accumulated hepatic fat, but at the same time, they had reduced fasting insulin and improved glucose tolerance (51). In addition, impaired hepatic phosphatidylcholine biosynthesis led to reduced levels of plasma triglycerides and HDL cholesterol in vivo (52). Accordingly, phosphatidylcholines and sphingomyelins were positively related to plasma HDL cholesterol in the present study. Furthermore, acyl-alkyl-phosphatidylcholines were inversely correlated to plasma triglycerides, opposite to diacyl-phosphatidylcholines, and higher levels of acyl-alkyl-phosphatidylcholines but not diacyl-phosphatidylcholines were linked to improved insulin sensitivity and reduced insulin secretion. Previous studies reported that acyl-alkyl-phosphatidylcholine levels were lower in obese subjects and subjects with insulin resistance (50,53). These mechanisms may contribute to the antithetical association between two phosphatidylcholine subclasses and T2D risk found in the present study and may indicate a key role of the type of linkage between phospholipid core and fatty acid residue. Furthermore, those phosphatidylcholines containing fatty acids with a lower number of carbons and double bonds were positively associated with T2D risk, contrary to those with a higher number of carbons and double bonds. Similar observations have recently been reported for fatty acid compositions of erythrocyte membrane phospholipids (54) and triglycerides (55), suggesting that lipids with a shorter chain length and saturated fatty acid residues may trigger development of T2D, whereas those containing longer chains and unsaturated fatty acids may offer protection.
In summary, the present data suggest that the identified metabolites could be part of different pathways involved in the early genesis of T2D. Therefore, these novel candidates could be useful in clinical practice to identify high-risk individuals earlier in order to delay or prevent disease onset. Furthermore, serum metabolites in this study predicted risk of T2D in a similar manner to a combination of classic risk factors; thus, their measurement may be a useful approach to predict T2D risk for individuals in the future. The metabolites could also serve as markers for specific metabolic pathways that are deranged and, thereby, allow the implementation of individualized preventive and therapeutic strategies. However, future investigations are warranted to calculate in detail the individual risks and to better understand the metabolic effects of these biomarkers and their biological mechanisms.
The primary strength of this study is that, to our knowledge, we were among the first to adopt a targeted metabolomic approach at population level and included a large sample from three independent, well-described study populations. Furthermore, our targeted metabolomic platform covered a wide variety of metabolites with known identity and quantitative measurements. Because we used a prospective design with consecutive follow-up, we were able to investigate time-dependent exposure-disease associations.
The study, however, had several limitations. First, because we used independent study populations, the conditions of biosample collection, storage, and preparation were not necessarily the same, which may be a source of variation. In the KORA and TüF studies, fasting blood samples were collected from all participants, whereas in EPIC-Potsdam only a small proportion of participants provided fasting blood samples. Nevertheless, we did not observe an effect modification of fasting status on the association between metabolite factors and T2D risk, and we could reproduce very similar metabolite factors comparing fasting to nonfasting samples. Furthermore, the metabolomic analyses were based on serum samples in the EPIC-Potsdam and KORA studies but on plasma samples in the TüF study. As previously reported (56), the correlation between these serum and plasma metabolites was high; however, the absolute metabolite concentrations were higher in serum, which could lead to systematic changes. Second, the analytical method detected most of the metabolites with high specificity; however, it may not have detected all possible interferences among metabolites. Of the metabolites that we identified, diacyl-phosphatidylcholine C38:3 and sphingomyelin C16:1 may be interfering compounds for sphingomyelin C24:1 and diacyl-phosphatidylcholine C30:2, respectively. Third, we only had a limited number of incident T2D cases available from KORA. We may not have had sufficient statistical power, and the replication results have to be interpreted with caution. Fourth, there is a chance that reverse causation may explain the results, implying that overt diabetic conditions that were undiagnosed may have caused these metabolite changes. When we accounted for this issue, the results remained robust. Last, because this was an observational study, we cannot prove causality but only show associations. However, the identified metabolites were also correlated to established T2D biomarkers as well as to measures of insulin sensitivity and secretion in a different population, which underlines the biological plausibility of the results.
In conclusion, this prospective investigation using metabolomics data of independent study populations identified sugar metabolites, amino acids, and choline-containing phospholipids to be independently associated with risk of T2D. Beyond the classic pathways, these candidates point toward a novel role of phospholipid and lipoprotein metabolism in T2D pathophysiology. Future studies should further elucidate the biological mechanisms.
See accompanying commentary, p. 349.
The current study was supported by a grant from the Federal Ministry of Education and Research, Germany (Bundesministerium für Bildung und Forschung), to the German Center for Diabetes Research (Förderkennzeichen 01GI0922 to D.Z.D.). N.S. is supported by a Heisenberg professorship from the Deutsche Forschungsgemeinschaft.
No potential conflicts of interest relevant to this article were reported.
A.Fl., H.B., and T.P. designed the study. A.Fl., N.S., Z.Y., and K.M. analyzed the data. A.Fl., N.S., Z.Y., K.M., D.D., H.-G.J., A.Fr., H.-U.H., M.H.A., A.P., M.R., C.P., R.W.-S., T.I., M.B.S., J.A., H.B., and T.P. discussed the data and interpreted the results. A.Fl. wrote the first draft of the manuscript. N.S., Z.Y., K.M., D.D., H.-G.J., A.Fr., H.-U.H., M.H.A., A.P., M.R., C.P., R.W.-S., T.I., M.B.S., J.A., H.B., and T.P. provided their expertise and contributed to the writing of the manuscript. A.Fr., A.P., T.I., and H.B. collected the data. C.P. and J.A. conducted and supervised the metabolomic measurements. All authors take full responsibility for the contents of the manuscript. A.Fl. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
The authors thank all the study participants of EPIC-Potsdam, KORA, and TüF. Special thanks go to Sven Knüppel and Wolfgang Bernigau (Department of Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke) for statistical advice and to Martin Floegel (Max Born Institute for Non-Linear Optics, Berlin) for his support in figure formatting.