Figures
Abstract
Background
Although obesity and cardiometabolic traits commonly overlap, underlying pathways remain incompletely defined. The association of metabolite profiles across multiple cardiometabolic traits may lend insights into the interaction of obesity and metabolic health. We sought to investigate metabolic signatures of obesity and related cardiometabolic traits in the community using broad-based metabolomic profiling.
Methods and Results
We evaluated the association of 217 assayed metabolites and cross-sectional as well as longitudinal changes in cardiometabolic traits among 2,383 Framingham Offspring cohort participants. Body mass index (BMI) was associated with 69 of 217 metabolites (P<0.00023 for all), including aromatic (tyrosine, phenylalanine) and branched chain amino acids (valine, isoleucine, leucine). Additional metabolic pathways associated with BMI included the citric acid cycle (isocitrate, alpha-ketoglutarate, aconitate), the tryptophan pathway (kynurenine, kynurenic acid), and the urea cycle. There was considerable overlap in metabolite profiles between BMI, abdominal adiposity, insulin resistance [IR] and dyslipidemia, modest overlap of metabolite profiles between BMI and hyperglycemia, and little overlap with fasting glucose or elevated blood pressure. Metabolite profiles were associated with longitudinal changes in fasting glucose, but the involved metabolites (ornithine, 5-HIAA, aminoadipic acid, isoleucine, cotinine) were distinct from those associated with baseline glucose or other traits. Obesity status appeared to “modify” the association of 9 metabolites with IR. For example, bile acid metabolites were strongly associated with IR among obese but not lean individuals, whereas isoleucine had a stronger association with IR in lean individuals.
Conclusions
In this large-scale metabolite profiling study, body mass index was associated with a broad range of metabolic alterations. Metabolite profiling highlighted considerable overlap with abdominal adiposity, insulin resistance, and dyslipidemia, but not with fasting glucose or blood pressure traits.
Citation: Ho JE, Larson MG, Ghorbani A, Cheng S, Chen M-H, Keyes M, et al. (2016) Metabolomic Profiles of Body Mass Index in the Framingham Heart Study Reveal Distinct Cardiometabolic Phenotypes. PLoS ONE 11(2): e0148361. https://doi.org/10.1371/journal.pone.0148361
Editor: David Meyre, McMaster University, CANADA
Received: November 3, 2015; Accepted: January 19, 2016; Published: February 10, 2016
Copyright: © 2016 Ho et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Data are available on dbGap under study accession phs000007.v22.p8, and access can be requested via the dbGaP website (http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000007.v22.p8).
Funding: This work was partially supported by the National Heart, Lung and Blood Institute’s Framingham Heart Study (Contract No. N01-HC-25195 and HHSN268201500001I), National Institutes of Health grants K23-HL116780 (JEH), 2R01-DK081572 (MGL, REG, TJW), a Boston University School of Medicine Department of Medicine Career Investment Award (JEH), and the Ellison Foundation (SC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: Drs. Gerszten, Vasan, Larson, and Wang are named as co-inventors on a pending patent application relating to metabolite predictors of diabetes, entitled ‘Multiplexed Biomarkers of Insulin Resistance’, application number PCT/US2009/049831, publication number WO2010005982. No monetary benefit has been received from this application by any coauthors to date. There are no other declarations relevant to employment, consultancy, patents, products in development or modified products. Further, the authors’ would like to confirm that this does not alter their adherence to PLOS ONE policies on sharing data and materials.
Introduction
Currently one quarter of adults in the United States are obese, and it is predicted that more than half of the population will be obese by 2030 [1]. With the growing obesity epidemic, the incidence of diabetes mellitus has doubled over the last 30 years [2]. The increased risk of diabetes is most pronounced in those with a high body mass index (BMI) [2], underscoring the close link between obesity and metabolic disease. Recent advances in high-throughput technology have allowed for the systematic assessment of metabolic profiles, and have provided insights into metabolic pathways that appear dysregulated in cardiometabolic disease [3–5]. Prior studies have focused on cross-sectional associations of metabolite profiles and metabolic traits [3–5], with some recent studies relating selected metabolites to the development of future clinical disease [6–8].
Although obesity and metabolic disease traits commonly occur together, it has been increasingly recognized that a subset of obese individuals can be classified as 'metabolically healthy,’ and, in turn, lean individuals can be 'metabolically unhealthy'. This observation motivates the use of tools such as metabolomic profiling to provide a better understanding of the heterogeneity in metabolic risk among obese and lean individuals. Furthermore, we hypothesized that metabolite profiles could presage changes in metabolic traits over time.
The Framingham Heart study provides a unique community-based setting in which cardiometabolic traits have been systematically and longitudinally assessed in individuals across the spectrum of cardiometabolic risk. We investigated the comprehensive profile of non-lipid and lipid metabolites in a sample of over 2,300 community-based adults, in order to elucidate metabolic signatures of obesity and related cardiometabolic traits. We sought to determine metabolite signatures associated with cross-sectional cardiometabolic traits, as well as with their longitudinal trajectories. Using genetic markers, we explored the directionality of association between metabolites and BMI.
Materials and Methods
Study Sample
Participants of the Framingham Heart Study Offspring Study—a prospective, observational, community-based cohort—were examined [9]. Of 3,799 participants who attended the baseline examination cycle (1991–1995), profiling of polar positive-charge metabolites was performed on 2,526 participants with available blood samples, of whom 2,383 attended at least one subsequent examination among exam 6 (1995–1998), 7 (1998–2001), and 8 (2005–2008). Of these participants, 1,969 also underwent profiling of polar negative-charge metabolites, and 1,962 underwent lipid profiling.
Blood pressure was defined as the average of two seated resting measurements performed by a physician. BMI was calculated as weight divided by height2 (kg/m2), and waist circumference was measured at the level of the umbilicus at mid-respiration. Diabetes mellitus was defined as a fasting glucose ≥126 mg/dL, non-fasting blood glucose ≥200 mg/dL, or the use of antidiabetic medications. The homeostatic model assessment of insulin resistance was calculated as HOMA-IR [= fasting insulin (ulU/ml) x fasting glucose (mmol/ml) / 22.5] [10], and log-transformed due to a skewed distribution. Insulin resistance was estimated as outlined previously [11], using the top quartile of HOMA-IR from the whole FHS Offspring sample free of diabetes at baseline. Current smoking was defined as smoking regularly during the year prior to the examination. Total cholesterol, high-density lipoprotein (HDL) cholesterol, and triglycerides were measured after an overnight fast. All participants provided written informed consent and the study protocol was approved by the Boston University Medical Center Institutional Review Board.
Metabolite profiling
Baseline plasma samples were collected in EDTA after an overnight fast, immediately processed, and stored at -80°C until assayed. Targeted metabolite profiling was performed using liquid chromatography with tandem mass spectrometry. Positively charged polar metabolites (including amino acids and derivatives, urea cycle intermediates, and nucleotides), lipid metabolites (lysophosphatidylcholines, LPCs; lysophosphatidylethanolamines, LPEs; phosphatidylcholines, PCs; sphingomyelins, SMs; cholesterol esters, CEs; diacylglycerols, DAGs; and triacylglycerols, TAGs), and negatively charged polar metabolites (organic acids, sugars, and bile acids) were profiled as previously described [12,13]. Nomenclature for lipid analytes entailed an initial number, denoting the total number of carbons in the lipid acyl chain, and a second number following the colon, denoting the total number of double bonds. Details are provided in S1 File and a list of 217 metabolites measured using this platform, and descriptions are detailed in Table A in S1 File.
Statistical Analysis
Baseline characteristics were summarized for the sample, and correlations of metabolic traits (BMI, waist circumference, HOMA-IR, HDL cholesterol, log-triglycerides, fasting glucose, systolic and diastolic blood pressure, pulse pressure) at baseline and 3 subsequent examinations were evaluated. Metabolite concentrations were log-transformed for analysis due to right-skewed distributions.
Our primary analysis focused on metabolite associations with BMI; secondary analyses examined related metabolic traits (waist circumference, HOMA-IR, HDL, log-triglycerides, fasting glucose, systolic and diastolic blood pressure). For analyses examining HOMA-IR and fasting glucose, participants taking antidiabetic medications were excluded. Cross-sectional data for metabolic traits (response variables) and metabolites (predictors) were assessed using linear regression, adjusted for age, sex, BMI, log-triglycerides, and a multi-level indicator variable for batch number. Associations between baseline metabolites and longitudinal changes in metabolic traits were assessed using PROC GLIMMIX to account for repeated measures in a random effects model, with a class variable to indicate subsequent exam cycles. Models were adjusted for baseline metabolic trait, age, sex, BMI, log-triglycerides, and batch. In order to account for multiple testing, a P-value threshold of 0.05 / (number of metabolites tested) = 0.05 / 217 = 0.00023 was deemed statistically significant.
In exploratory analyses, we examined the association of metabolites with HOMA-IR in non-diabetic individuals stratified by obesity status. Analyses were performed using SAS, version 9.2 (SAS Institute, Cary, NC).
Results
A total of 2,383 participants (mean age 55 ± 10 years, 53% women) underwent profiling of 217 metabolites. The mean BMI was 27.5 ± 4.9 kg/m2, with median BMI of 26.9 kg/m2 (25th, 75th percentile 24.0, 30.0), and 25% were obese. The range of BMI is presented in S1 Fig Baseline median HOMA-IR was 2.7 mg*microU/dL*mL (25th, 75th percentile 0.9, 2.9), and 28% were insulin resistant, defined as the upper quartile among non-diabetic individuals (Table 1).
Body mass index (BMI) is associated with broad alterations in multiple biochemical pathways
BMI was significantly associated with 69 of 217 metabolites using a conservative Bonferroni-corrected threshold (P<0.00023, Tables 2 and 3). An additional 11 metabolites displayed “suggestive” associations (P<0.001, full results in S1 Table. The top metabolite groups associated with BMI included aromatic amino acids (tyrosine, P = 3.2x10-34; phenylalanine P = 5.7x10-25), branched chain amino acids (valine, P = 2.6x10-32; isoleucine, P = 2.9x10-25; leucine, P = 2.8x10-22), and other amino acids (alanine, P = 1.3x10-11; glycine, P = 4.2x10-10; proline, P = 1.7x10-4), all of which were positively associated with BMI.
In addition, BMI was positively associated with multiple metabolites in the citric acid cycle (isocitrate, P = 8.8x10-20; alpha-ketoglutarate, P = 6.3x10-11; aconitate, P = 5.2x10-6) and the tryptophan pathway (kynurenine, P = 3.3x10-19; kynurenic acid, P = 2.9x10-11). Other pathways associated with BMI included the urea cycle (citrulline, P = 8.6x10-5; ornithine, P = 7.2x10-5), nucleic acid metabolism (xanthosine, P = 3.8x10-7; uric acid, P = 3.3x10-6), and multiple creatine-related metabolites (carnitine, P = 1.6x10-7; choline, P = 1.4x10-8; glycero-phosphocholine, P = 3.6x10-7).
Lipid subspecies reveal a distinct pattern in obesity
Of 104 lipid subspecies profiled, 38 were significantly associated with BMI, including glycerol, an integral component of TAGs (P = 2.1x10-11), choline, and glycero-phosphocholine, both important components of PCs and SMs (P = 1.4x10-8; P = 3.6x10-7, respectively). In addition, BMI was associated with LPCs, LPEs (inversely associated with BMI), SMs, CEs (positively associated with BMI), PCs and TAGs (both positive and negative associations, P<0.00023 for all). A distinct pattern emerged within TAG subspecies; shorter carbon length and a lower number of double bonds were associated with higher BMI, whereas longer carbon length and higher number of double bonds were associated with a lower BMI (Fig 1A). The converse was observed with LPEs, where shorter carbon length and lower number of double bonds were associated with lower BMI (full results in S1 Table). No such pattern was discernible for PCs or LPCs.
X-axis represents number of double bonds, Y-axis is the beta-coefficient for the change in metabolic trait per 1 standard deviation increase in metabolite. Darker circles are shorter length TAGs, and lighter circles represent greater carbon length TAGs. Panel (A) shows body-mass index, (B) waist circumference, (C) HOMA-IR, (D) fasting glucose, (E) HDL cholesterol, and (F) systolic blood pressure.
Metabolite profiles of obesity are linked to insulin resistance and dyslipidemia
The majority of metabolites associated with BMI were concomitantly associated with abdominal adiposity, as assessed by waist circumference. There was also overlap in metabolic profiles associated with BMI and other metabolic traits, including insulin resistance as estimated by HOMA-IR, HDL cholesterol, log-triglyceride concentrations, and fasting glucose, as displayed in the heatmap (Figs 2 and 3, Tables 2 and 3). Of 31 non-lipid metabolites associated with BMI, 27 demonstrated at least suggestive associations (at P<0.001) with at least one other metabolic trait other than waist circumference, after adjusting for BMI (HOMA-IR, n = 16; HDL, n = 13; log-triglycerides, n = 23; fasting glucose, n = 8). Of these, 12 metabolites were associated with at least 4 metabolic traits: isoleucine, leucine, proline, valine, and aminoadipic acid were associated with BMI, HOMA-I77R, HDL, and triglycerides; alpha-hydroxybutyric acid was associated with BMI, fasting glucose, HDL, and triglycerides; fructose/glucose/galactose, kynurine, isocitrate, and lactate were associated with BMI, fasting glucose, HOMA-IR, and triglycerides; glycine, and quinolinic acid were associated with all 5 metabolic traits. By contrast, there was limited overlap of metabolites associated with BMI and blood pressure traits, with only 5 of the 31 non-lipid BMI metabolites found in association with either systolic and/or diastolic blood pressure.
TAG patterns associated with BMI were similar for waist circumference and HOMA-IR, with shorter carbon length and lower number of double bonds (e.g. more saturation) associated with worse metabolic profiles. By contrast, there was no distinct association between TAG characteristics and fasting glucose. Lastly, TAG patterns observed in association with HDL and blood pressure traits appeared J-shaped in distribution—with low carbon length, saturated TAGs associated with higher HDL and blood pressure on the one hand, and intermediate length and saturation TAGs associated with lower HDL and lower blood pressure on the other hand (Fig 1). Full results of metabolic traits and metabolic profile analyses are shown in S1 Table.
Baseline metabolite profiles associated with longitudinal changes in metabolic traits are distinct from cross-sectional findings
The correlation of baseline and longitudinal metabolic traits over approximately 14 years is presented in Table B in S1 File. There was a strong correlation of baseline BMI, waist circumference, and HDL cholesterol with subsequent measures (Pearson correlation coefficient r > 0.70 for all), but lower correlation of baseline fasting glucose, triglycerides, systolic and diastolic blood pressure with follow-up measures (r = 0.36–0.56 for baseline versus last exam).
Metabolites associated with longitudinal changes in fasting glucose were distinct from metabolites associated with fasting glucose cross-sectionally. Specifically, isoleucine (P = 1.1x10-4), ornithine (P = 2.2x10-5), 5-HIAA (P = 1.1x10-4), aminoadipic acid (P = 5.8x10-5), and cotinine (P = 3.5x10-5), all predicted future changes in fasting glucose, despite not being associated with baseline fasting glucose values (Table 4). By contrast, there were no metabolites associated with longitudinal changes in the remaining metabolic traits (BMI, waist circumference, HDL, triglycerides, blood pressure traits) that were not already associated in cross-sectional analyses.
Differential effect of metabolites on insulin resistance in obese and non-obese individuals
In exploratory analyses, we compared metabolite profiles associated with insulin resistance, as estimated by HOMA-IR, separately in obese (n = 511) and non-obese (n = 1680), non-diabetic individuals. Many metabolites associated with insulin resistance appear to be consistent across BMI strata (Tables A and B in S1 File). However, obesity status appeared to “modify” the association for 9 of the metabolites, with suggestive P-values for interaction (P≤0.01 for all, Table 5). This included amino acids and related metabolites (3-hydroxyphenylacetic acid, tryptophan, tyrosine, isoleucine), bile acids (glycocholate, taurocholate), and others (carnitine, creatine, inosine). All metabolites appeared more strongly associated with HOMA-IR among the non-obese, with exception of isoleucine, which appeared to have a stronger association with HOMA-IR among the obese. Of note, none of the P-values for interaction met the Bonferroni-corrected P-value threshold.
Discussion
Our findings are threefold: first, BMI is associated with broad alterations in metabolites representing multiple distinct biochemical pathways, including the metabolism of amino acids, the citric acid cycle, nucleotides, the urea cycle, and specific triacylglycerol species. Second, there was considerable overlap in metabolite profiles associated with BMI and other metabolic traits, including abdominal adiposity, insulin resistance, and dyslipidemia. By contrast, profiles of fasting glucose and blood pressure traits were relatively distinct. Lastly, obesity appeared to modify metabolic signatures of insulin resistance. Collectively, these data highlight the important interaction of obesity and metabolic health—and reveal with further granularity the intricateness of this relationship.
A number of recent studies have identified a range of metabolites as markers and possible effectors of metabolic disease [3,4,8,14–18]. We extend these analyses to a broad range of cardiometabolic phenotypes and present a more detailed analysis of metabolite profiles with respect to adiposity, dyslipidemia, metabolic disease, and blood pressure. Our group has previously examined metabolites in relation to these traits [5], however we now extend findings to more than double the previous sample size, with a much broader profiling platform that now includes negatively charged anions and lipid subspecies. A number of novel findings are worth noting. For example, mounting evidence supports the role of dysregulated branched-chain amino acid metabolism in human cardiometabolic disease [6,7,14,16,19–21], with experimental studies demonstrating activation of the mammalian target of rapamycin (mTOR) by branched chain amino acids leading to insulin resistance [14,22,23]. Our study demonstrates not only that branched-chain and aromatic amino acids are associated with BMI, but also other essential and non-essential amino acids (alanine, asparagine, glycine, proline) and amino acid intermediates. The latter included multiple metabolites of the tryptophan pathway (indole propionate, kynurenic and quinolinic acid), cystathionine in the cysteine pathway, and 2-aminoadipic acid, a lysine metabolite, which has previously been associated with incident diabetes [24]. Prior studies have demonstrated lower glycine concentrations to be associated with metabolic traits, worse glucose disposal rate, and diabetes [14,16,25,26]. We now confirm the association of lower glycine concentrations with unfavorable metabolic traits, and moreover show that its derivative metabolite, dimethylglycine, is positively associated with metabolic traits.
We observed that most amino acids associated with BMI were also associated with abdominal adiposity, insulin-resistance, HDL, and triglycerides. However, two essential amino acids were negatively associated with fasting glucose but not with BMI: histidine, which has been shown to protect against diabetic deterioration in mice [27], and threonine, a precursor to glycine synthesis. Among non-amino acid metabolites, considerable overlap in metabolite profiles existed between BMI, abdominal adiposity, insulin resistance, and lipid traits, but not with elevated fasting glucose, highlighting two known distinct metabolic phenotypes. Whereas impaired fasting glucose reflects predominantly a defect in insulin secretion, impaired glucose tolerance is more closely related to insulin resistance [28,29] and clinical outcomes [30]. Our study provides further granularity with respect to differences in metabolite intermediaries between these two phenotypes.
Multiple additional metabolic pathways appear to be altered in relation to metabolic traits, including carbohydrate metabolism, ketosis, lipid metabolism, and bile acid metabolism [3,4,8,14,17,31–33]. Our findings affirm the effect of BMI and related traits on this wide spectrum of metabolic pathways, and extend previous findings in two ways. First, our broad-based platform was able to shed further light regarding multiple metabolic intermediates involved in related pathways. For example, we demonstrate novel alterations in multiple intermediates of the citric acid cycle, including aconitate, isocitrate, and alpha-ketoglutarate, all of which are positively associated with BMI, abdominal obesity, HOMA-IR, and triglycerides. A previous study demonstrated that alpha-ketoglutarate was reduced after weight loss in obese insulin-resistant women [34]. Further, the citric acid cycle enzyme NADP+-dependent isocitrate dehydrogenase is expressed in liver and adipocytes, where it produces NADPH by converting isocitrate to alpha-ketoglutarate, and thereby plays a critical role in adipogenesis, particularly with respect to triglycerides and cholesterol [35].
Previous studies have demonstrated an inverse association of long-chain acylcarnitines [32], positive association with medium-chain acylcarnitines [31], and also greater fatty acid saturation levels [3,4,18] with metabolic traits. We extend these findings and demonstrate a distinct TAG pattern associated with BMI, central adiposity, and HOMA-IR, with lower carbon length and saturation associated with worse metabolic profiles. This finding is in concert with prior lipidomic studies, that have demonstrated the association of low carbon number and double-bond content TAG subspecies with incident cardiovascular disease [36]. Whether TAG subspecies are markers or effectors of disease remains unclear, however it is known that triglycerides may modulate lipoprotein retention and thrombogenicity as potential underlying mechanisms of atherosclerosis [37]. Together, these findings suggest that profiling of specific TAG species beyond conventional assessment of bulk triglyceride concentrations may lend more granular insights into underlying pathways contributing to cardiometabolic disease.
Interestingly, we found no metabolites that predicted change in BMI over time above and beyond baseline BMI. This is probably due to the fact that baseline and subsequent BMI measurements are highly correlated. By contrast, we found multiple metabolites that predicted future change in fasting glucose measures, including aminoadipic acid, ornithine, isoleucine, 5-HIAA, and cotinine, each of which has support from previous studies linking metabolite to glucose homeostasis and diabetes mellitus: Our group has previously demonstrated an association of 2-aminoadipic acid with incident diabetes [24]. In this study, administration of 2-aminoadipic acid to mice lowered their fasting glucose levels. Further, in cell-based experiments, 2-aminoadipic acid enhanced insulin secretion from pancreatic beta cell lines. This suggests that 2-aminoadipic acid may be involved in glucose homeostasis and underlying pathology. Ornithine is an amino acid that is part of the urea cycle. In animal models of diabetes, insulin deficiency was associated with increased levels of ornithine along with high branched chain amino acid levels such as isoleucine, [38] and recapitulates hyperaminoacidemia seen in human metabolic disease that precedes diabetes [6]. Amino acid supplementation in humans directly promotes insulin resistance [23], and may be related to effects on insulin signaling in skeletal muscle [39]. 5-HIAA is a metabolite of serotonin, a neurotransmitter that plays a central role in glucose homeostasis via central autonomic pathways and peripheral release from gastrointestinal enterochromaffin cells [40]. Genetic serotonin receptor deletion is associated with diabetes mellitus, further substantiating the role of serotonin in diabetes development [41]. Lastly, cotinine is the predominant metabolite of nicotine, and can be used as a marker of smoking status [42], which in turn has been linked to incident diabetes in cohort studies [43]. Smoking has been shown to lead to impaired glucose tolerance acutely [44], and abdominal adiposity [45], both of which may lead to eventual diabetes mellitus.
It is well recognized that within a given BMI category, a range of metabolic health phenotypes can be observed. Whether or not 'metabolically healthy obesity' exists is controversial [46], but it is known that heterogeneity exists in the relation of metabolic profiles with metabolic health across weight categories [19]. Our findings highlight the complex interaction of obesity with metabolic health, and demonstrate notable differences in metabolic signatures of insulin resistance in obese compared with non-obese individuals. For example, the amino acids isoleucine and proline were more strongly associated with insulin resistance in the non-obese. By contrast, the association of saturated TAGs and both conjugated and deconjugated bile acids with insulin resistance was markedly more pronounced in obesity. The association of bile acids and metabolic traits has previously been described [3,47–49]. Differences in obese compared with non-obese may also relate to the rate-limiting enzyme of bile acid synthesis, cholesterol 7α-hydroxylase (CYP7A), whose basal expression is increased in an obesity mouse model via both glucose-stimulated epigenetic modifications and the insulin/FoxO1 pathway [50]. Alterations in gut microbiota may also contribute, as intestinal flora can modify bile acids by deconjugation, oxidation, dehydroxylation, and sulfation [51], and have been implicated in the pathogenesis of diabetes [52].
Several limitations deserve mention. The 217 metabolites analyzed in this study were selected as part of a targeted platform geared toward small molecules based on previous studies linking several metabolites to insulin resistance [6]. Other metabolite platforms can acquire profiles in a less biased manner, however advantages for our approach included highly specific analyte identification, and capability for absolute quantification of metabolites. It remains unclear whether metabolites represent biomarkers or actual mediators of metabolic disease, and causal inferences cannot be drawn from our observational study. There were 154 individuals with diabetes mellitus in our sample, and we thus we had limited power to examine the association of metabolites with metabolic traits among this subset. We were unable to ascertain relations of plasma metabolites to the adipose tissue secretome, as these data were not available in our study sample. Lastly, we performed metabolite profiles at one single time point, precluding the assessment of cardiometabolic disease and medication effects on longitudinal changes in metabolites.
In summary, body mass index was associated with a broad range of metabolic alterations. Metabolite profiling illuminated considerable overlap with abdominal adiposity, insulin resistance, and dyslipidemia, but not with fasting glucose or blood pressure traits. The association of metabolic profiles with IR was modified by the presence or absence of obesity, highlighting heterogeneity in the association of obesity and metabolic disease. Future studies are needed to identify which metabolites could have importance as causal risk factors for metabolic disease, and, in turn, which metabolites could serve as potential targets for therapy.
Supporting Information
S1 Fig. Distribution of baseline body mass index in the sample.
https://doi.org/10.1371/journal.pone.0148361.s001
(PDF)
S1 File. Supplemental information to manuscript.
Supplemental methods are included in text. Metabolites assayed and descriptions are summarized (Table A), and correlations of metabolic traits across subsequent examinations displayed (Table B). Results for the association of non-lipid (Table C) and lipid (Table D) with HOMA-IR in obese and non-obese individuals are summarized.
https://doi.org/10.1371/journal.pone.0148361.s002
(DOCX)
S1 Table. Full results of metabolite associations with cardiometabolic traits.
https://doi.org/10.1371/journal.pone.0148361.s003
(XLSX)
Acknowledgments
This work was partially supported by the National Heart, Lung and Blood Institute's Framingham Heart Study (Contract No. N01-HC-25195 and HHSN268201500001I), National Institutes of Health grants K23-HL116780 (JEH), 2R01-DK081572 (MGL, REG, TJW), a Boston University School of Medicine Department of Medicine Career Investment Award (JEH), and the Ellison Foundation (SC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author Contributions
Conceived and designed the experiments: JEH SC REG TJW. Performed the experiments: JEH SC EPR CBC RSV REG TJW. Analyzed the data: MGL AG MHC MK. Contributed reagents/materials/analysis tools: MGL EPR CBC REG TJW. Wrote the paper: JEH MGL AG SC MHC MK EPR CBC RSV REG TJW.
References
- 1. Go AS, Mozaffarian D, Roger VL, Benjamin EJ, Berry JD, Borden WB, et al. (2013) Heart disease and stroke statistics—2013 update: a report from the American Heart Association. Circulation 127: e6–e245. pmid:23239837
- 2. Fox CS, Pencina MJ, Meigs JB, Vasan RS, Levitzky YS, D'Agostino RB Sr. (2006) Trends in the incidence of type 2 diabetes mellitus from the 1970s to the 1990s: the Framingham Heart Study. Circulation 113: 2914–2918. pmid:16785337
- 3. Suhre KK, Meisinger CC, Döring AA, Altmaier EE, Belcredi PP, Gieger CC, et al. (2010) Metabolic footprint of diabetes: a multiplatform metabolomics study in an epidemiological setting. PLoS ONE 5: e13953–e13953. pmid:21085649
- 4. Würtz P, Mäkinen V-P, Soininen P, Kangas AJ, Tukiainen T, Kettunen J, et al. (2012) Metabolic signatures of insulin resistance in 7,098 young adults. Diabetes 61: 1372–1380. pmid:22511205
- 5. Cheng S, Rhee EP, Larson MG, Lewis GD, McCabe EL, Shen D, et al. (2012) Metabolite profiling identifies pathways associated with metabolic risk in humans. Circulation 125: 2222–2231. pmid:22496159
- 6. Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, et al. (2011) Metabolite profiles and the risk of developing diabetes. Nat Med 17: 448–453. pmid:21423183
- 7. Stancakova A, Civelek M, Saleem NK, Soininen P, Kangas AJ, Cederberg H, et al. (2012) Hyperglycemia and a common variant of GCKR are associated with the levels of eight amino acids in 9,369 Finnish men. Diabetes 61: 1895–1902. pmid:22553379
- 8. Ferrannini E, Natali A, Camastra S, Nannipieri M, Mari A, Adam K-P, et al. (2013) Early metabolic markers of the development of dysglycemia and type 2 diabetes and their physiological significance. Diabetes 62: 1730–1737. pmid:23160532
- 9. Kannel WB, Feinleib M, McNamara PM, Garrison RJ, Castelli WP (1979) An investigation of coronary heart disease in families. The Framingham offspring study. Am J Epidemiol 110: 281–290. pmid:474565
- 10. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC (1985) Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 28: 412–419. pmid:3899825
- 11. Preis SR, Massaro JM, Robins SJ, Hoffmann U, Vasan RS, Irlbeck T, et al. (2010) Abdominal subcutaneous and visceral adipose tissue and insulin resistance in the Framingham heart study. Obesity (Silver Spring) 18: 2191–2198. pmid:20339361
- 12. Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, et al. (2011) Metabolite profiles and the risk of developing diabetes. Nat Med 17: 448–453. pmid:21423183
- 13. Rhee EP, Cheng S, Larson MG, Walford GA, Lewis GD, McCabe E, et al. Lipid profiling identifies a triacylglycerol signature of insulin resistance and improves diabetes prediction in humans. J Clin Invest 121: 1402–1411. pmid:21403394
- 14. Newgard CB, An J, Bain JR, Muehlbauer MJ, Stevens RD, Lien LF, et al. (2009) A branched-chain amino acid-related metabolic signature that differentiates obese and lean humans and contributes to insulin resistance. Cell Metab 9: 311–326. pmid:19356713
- 15. Fizelova M, Miilunpohja M, Kangas AJ, Soininen P, Kuusisto J, Ala-Korpela M, et al. (2015) Associations of multiple lipoprotein and apolipoprotein measures with worsening of glycemia and incident type 2 diabetes in 6607 non-diabetic Finnish men. Atherosclerosis 240: 272–277. pmid:25818853
- 16. Palmer ND, Stevens RD, Antinozzi PA, Anderson A, Bergman RN, Wagenknecht LE, et al. (2015) Metabolomic profile associated with insulin resistance and conversion to diabetes in the Insulin Resistance Atherosclerosis Study. J Clin Endocrinol Metab 100: E463–468. pmid:25423564
- 17. Wurtz P, Wang Q, Kangas AJ, Richmond RC, Skarp J, Tiainen M, et al. (2014) Metabolic signatures of adiposity in young adults: Mendelian randomization analysis and effects of weight change. PLoS Med 11: e1001765. pmid:25490400
- 18. Wurtz P, Makinen VP, Soininen P, Kangas AJ, Tukiainen T, Kettunen J, et al. (2012) Metabolic signatures of insulin resistance in 7,098 young adults. Diabetes 61: 1372–1380. pmid:22511205
- 19. Batch BC, Shah SH, Newgard CB, Turer CB, Haynes C, Bain JR, et al. (2013) Branched chain amino acids are novel biomarkers for discrimination of metabolic wellness. Metabolism 62: 961–969. pmid:23375209
- 20. Würtz P, Soininen P, Kangas AJ, Rönnemaa T, Lehtimäki T, Kähönen M, et al. (2013) Branched-chain and aromatic amino acids are predictors of insulin resistance in young adults. Diabetes Care 36: 648–655. pmid:23129134
- 21. Floegel A, Stefan N, Yu Z, Muhlenbruch K, Drogan D, Joost HG, et al. (2013) Identification of serum metabolites associated with risk of type 2 diabetes using a targeted metabolomic approach. Diabetes 62: 639–648. pmid:23043162
- 22. Tremblay F, Marette A (2001) Amino acid and insulin signaling via the mTOR/p70 S6 kinase pathway. A negative feedback mechanism leading to insulin resistance in skeletal muscle cells. J Biol Chem 276: 38052–38060. pmid:11498541
- 23. Krebs M, Krssak M, Bernroider E, Anderwald C, Brehm A, Meyerspeer M, et al. (2002) Mechanism of amino acid-induced skeletal muscle insulin resistance in humans. Diabetes 51: 599–605. pmid:11872656
- 24. Wang TJ, Ngo D, Psychogios N, Dejam A, Larson MG, Vasan RS, et al. (2013) 2-Aminoadipic acid is a biomarker for diabetes risk. J Clin Invest 123: 4309–4317. pmid:24091325
- 25. Thalacker-Mercer AE, Ingram KH, Guo F, Ilkayeva O, Newgard CB, Garvey WT (2014) BMI, RQ, Diabetes, and Sex Affect the Relationships Between Amino Acids and Clamp Measures of Insulin Action in Humans. Diabetes 63: 791–800. pmid:24130332
- 26. Wang-Sattler R, Yu Z, Herder C, Messias AC, Floegel A, He Y, et al. (2012) Novel biomarkers for pre-diabetes identified by metabolomics. Mol Syst Biol 8: 615. pmid:23010998
- 27. Lee YT, Hsu CC, Lin MH, Liu KS, Yin MC (2005) Histidine and carnosine delay diabetic deterioration in mice and protect human low density lipoprotein against oxidation and glycation. Eur J Pharmacol 513: 145–150. pmid:15878720
- 28. Festa A, D'Agostino R Jr., Hanley AJ, Karter AJ, Saad MF, Haffner SM (2004) Differences in insulin resistance in nondiabetic subjects with isolated impaired glucose tolerance or isolated impaired fasting glucose. Diabetes 53: 1549–1555. pmid:15161760
- 29. Carnevale Schianca GP, Rossi A, Sainaghi PP, Maduli E, Bartoli E (2003) The significance of impaired fasting glucose versus impaired glucose tolerance: importance of insulin secretion and resistance. Diabetes Care 26: 1333–1337. pmid:12716784
- 30. (1999) Glucose tolerance and mortality: comparison of WHO and American Diabetes Association diagnostic criteria. The DECODE study group. European Diabetes Epidemiology Group. Diabetes Epidemiology: Collaborative analysis Of Diagnostic criteria in Europe. Lancet 354: 617–621. pmid:10466661
- 31. Huffman KM, Shah SH, Stevens RD, Bain JR, Muehlbauer M, Slentz CA, et al. (2009) Relationships between circulating metabolic intermediates and insulin action in overweight to obese, inactive men and women. Diabetes Care 32: 1678–1683. pmid:19502541
- 32. Gall WE, Beebe K, Lawton KA, Adam K-P, Mitchell MW, Nakhle PJ, et al. (2010) alpha-hydroxybutyrate is an early biomarker of insulin resistance and glucose intolerance in a nondiabetic population. PLoS ONE 5: e10883.
- 33. Weir JM, Wong G, Barlow CK, Greeve MA, Kowalczyk A, Almasy L, et al. (2013) Plasma lipid profiling in a large population-based cohort. J Lipid Res 54: 2898–2908. pmid:23868910
- 34. Campbell C, Grapov D, Fiehn O, Chandler CJ, Burnett DJ, Souza EC, et al. (2014) Improved metabolic health alters host metabolism in parallel with changes in systemic xeno-metabolites of gut origin. PLoS ONE 9: e84260. pmid:24416208
- 35. Koh H-J, Lee S-M, Son B-G, Lee S-H, Ryoo ZY, Chang K-T, et al. (2004) Cytosolic NADP+-dependent isocitrate dehydrogenase plays a key role in lipid metabolism. J Biol Chem 279: 39968–39974. pmid:15254034
- 36. Stegemann C, Pechlaner R, Willeit P, Langley SR, Mangino M, Mayr U, et al. (2014) Lipidomics profiling and risk of cardiovascular disease in the prospective population-based Bruneck study. Circulation 129: 1821–1831. pmid:24622385
- 37. Sambola A, Osende J, Hathcock J, Degen M, Nemerson Y, Fuster V, et al. (2003) Role of risk factors in the modulation of tissue factor activity and blood thrombogenicity. Circulation 107: 973–977. pmid:12600909
- 38. Sailer M, Dahlhoff C, Giesbertz P, Eidens MK, de Wit N, Rubio-Aliaga I, et al. (2013) Increased plasma citrulline in mice marks diet-induced obesity and may predict the development of the metabolic syndrome. PLoS ONE 8: e63950. pmid:23691124
- 39. Patti ME, Brambilla E, Luzi L, Landaker EJ, Kahn CR (1998) Bidirectional modulation of insulin action by amino acids. J Clin Invest 101: 1519–1529. pmid:9525995
- 40. Lam DD, Heisler LK (2007) Serotonin and energy balance: molecular mechanisms and implications for type 2 diabetes. Expert Rev Mol Med 9: 1–24.
- 41. Bonasera SJ, Tecott LH (2000) Mouse models of serotonin receptor function: toward a genetic dissection of serotonin systems. Pharmacol Ther 88: 133–142. pmid:11150593
- 42. Benowitz NL, Bernert JT, Caraballo RS, Holiday DB, Wang J (2009) Optimal serum cotinine levels for distinguishing cigarette smokers and nonsmokers within different racial/ethnic groups in the United States between 1999 and 2004. Am J Epidemiol 169: 236–248. pmid:19019851
- 43. Yeh HC, Duncan BB, Schmidt MI, Wang NY, Brancati FL (2010) Smoking, smoking cessation, and risk for type 2 diabetes mellitus: a cohort study. Ann Intern Med 152: 10–17. pmid:20048267
- 44. Frati AC, Iniestra F, Ariza CR (1996) Acute effect of cigarette smoking on glucose tolerance and other cardiovascular risk factors. Diabetes Care 19: 112–118. pmid:8718429
- 45. Berlin I (2008) Smoking-induced metabolic disorders: a review. Diabetes Metab 34: 307–314. pmid:18468932
- 46. Kramer CK, Zinman B, Retnakaran R (2013) Are metabolically healthy overweight and obesity benign conditions?: A systematic review and meta-analysis. Annals of internal medicine 159: 758–769. pmid:24297192
- 47. Ahmad NN, Pfalzer A, Kaplan LM (2013) Roux-en-Y gastric bypass normalizes the blunted postprandial bile acid excursion associated with obesity. Int J Obes (2005) 37: 1553–1559.
- 48. Vincent RP, Omar S, Ghozlan S, Taylor DR, Cross G, Sherwood RA, et al. (2013) Higher circulating bile acid concentrations in obese patients with type 2 diabetes. Ann Clin Biochem 50: 360–364. pmid:23771134
- 49. Ho JE, Larson MG, Vasan RS, Ghorbani A, Cheng S, Rhee EP, et al. (2013) Metabolite profiles during oral glucose challenge. Diabetes 62: 2689–2698. pmid:23382451
- 50. Li T, Francl JM, Boehme S, Ochoa A, Zhang Y, Klaassen CD, et al. (2012) Glucose and insulin induction of bile acid synthesis: mechanisms and implication in diabetes and obesity. J Biol Chem 287: 1861–1873. pmid:22144677
- 51. Ridlon JM, Kang DJ, Hylemon PB (2006) Bile salt biotransformations by human intestinal bacteria. J Lipid Res 47: 241–259. pmid:16299351
- 52. Cani PD, Delzenne NM (2009) The role of the gut microbiota in energy metabolism and metabolic disease. Curr Pharm Des 15: 1546–1558. pmid:19442172