- Research article
- Open Access
- Open Peer Review
Cross-sectional examination of metabolites and metabolic phenotypes in uremia
BMC Nephrology volume 16, Article number: 98 (2015)
Although metabolomic approaches have begun to document numerous changes that arise in end stage renal disease (ESRD), how these alterations relate to established metabolic phenotypes in uremia is unknown.
In 200 incident hemodialysis patients we used partial least squares discriminant analysis to identify which among 166 metabolites could best discriminate individuals with or without diabetes, and across tertiles of body mass index, serum albumin, total cholesterol, and systolic blood pressure.
Our data do not recapitulate metabolomic signatures of diabetes and obesity identified among individuals with normal renal function (e.g. elevations in branched chain and aromatic amino acids) and highlight several potential markers of diabetes status specific to ESRD, including xanthosine-5-phosphate and vanillylmandelic acid. Further, our data identify significant associations between elevated tryptophan and long-chain acylcarnitine levels and both decreased total cholesterol and systolic blood pressure in ESRD. Higher tryptophan levels were also associated with higher serum albumin levels, but this may reflect tryptophan’s significant albumin binding. Finally, an examination of the uremic retention solutes captured by our platform in relation to 24 clinical phenotypes provides a framework for investigating mechanisms of uremic toxicity.
In sum, these studies leveraging metabolomic and metabolic phenotype data acquired in a well-characterized ESRD cohort demonstrate striking differences from metabolomics studies in the general population, and may provide clues to novel functional pathways in the ESRD population.
End stage renal disease (ESRD) is characterized by various metabolic disturbances linked to adverse outcomes, but the nature of these associations is poorly understood or even counterintuitive . For example, low serum albumin has consistently been linked to cardiovascular mortality in dialysis [2, 3], yet the association is not simply a result of inadequate nutrition, as intradialytic parenteral nutrition does not always improve survival in malnourished ESRD patients , and alternative pathways including inflammation are known to contribute . Dyslipidemia is also common in ESRD, including atherogenic increases in carbamylated and oxidized-LDL cholesterol. Clinical trials, however, have failed to demonstrate a survival benefit with statin therapy in ESRD, despite substantial reductions in LDL cholesterol [6, 7]. Indeed, several studies have identified an inverse association between cholesterol levels and uremic cardiovascular risk [8, 9]. Similar paradoxical associations, sometimes described as examples of ‘reverse epidemiology’ in ESRD, have been noted with body mass index (BMI) and blood pressure [8, 10–12] and survival on hemodialysis.
Metabolomic approaches enable high resolution interrogation of metabolic phenotypes. For example, studies have begun to highlight specific metabolites associated with obesity and diabetes risk in the general population, particularly branched-chain and aromatic amino acids, but also short chain acylcarnitines, the glutamate/glutamine ratio, and bile acids [13–23]. Ultimately, these findings may lead to a more refined understanding of underlying disease mechanisms. Whether these and other associations identified in the burgeoning metabolomics literature extend to the ESRD population is unknown. While initial applications of metabolomics to nephrology research have generated a broad view of uremia  and the immediate effects of the hemodialysis procedure [25, 26], no studies have explicitly examined the cross-sectional relationship between uremic metabolite alterations and metabolic or other clinical phenotypes captured in a given ESRD study population.
We recently performed liquid chromatography-mass spectrometry (LC-MS) based metabolite profiling on plasma obtained from incident hemodialysis subjects in the Accelerated Mortality on Renal Replacement (ArMORR) cohort, and highlighted baseline levels of long-chain acylcarnitines as markers of future cardiovascular death in a longitudinal study . Using this data set, we now leverage the rich clinical phenotyping available in ArMORR to examine cross-sectional associations between plasma metabolite levels and select metabolic phenotypes. Specifically, we sought to identify associations between metabolites and metabolic phenotype data in this well-characterized ESRD cohort, hypothesizing that these findings would differ from metabolomic associations derived from the general population. The aim of such work is to gain new insights into the functional metabolic pathways that underlie the poorly understood associations between clinical outcomes in ESRD and variables such as BMI, albumin, cholesterol, and blood pressure. Further, we examined the correlation between the subset of metabolites measured by our platform that are known uremic retention solutes and numerous clinical variables captured in ArMORR. Taken together, these results provide a new perspective on the metabolic disturbances that accompany and arise in ESRD.
The Accelerated Mortality on Renal Replacement (ArMORR) study is a prospective cohort study of 10,044 incident hemodialysis patients in any of 1,056 U.S. centers operated by Fresenius Medical Care North America between June 2004 and August 2005. Full details have been previously published [27–31]. All subjects underwent 1 year of follow-up except for those who died (15.2 %), voluntarily discontinued dialysis (5 %), underwent kidney transplantation (3 %), recovered renal function (4 %), or transferred to a dialysis unit outside the Fresenius system (12 %). Clinical data were prospectively collected by physicians at the point-of-care and included demographics, coexisting conditions, and routine laboratory studies. Plasma samples drawn at the beginning of a dialysis session within 14 days of initiation of outpatient hemodialysis, and that would otherwise have been discarded after routine clinical testing, were saved and stored in liquid nitrogen. We previously evaluated the associations between metabolite profile and mortality in the ArMORR cohort in a study of 100 randomly selected patients who died within the first year of dialysis from a cardiovascular cause, and 100 patients who were alive 1 year after starting dialysis, frequency matched for age, sex, and race . We included this entire group of individuals in the current study. This study was approved by the Massachusetts General Hospital IRB, which waived the need for informed consent because all personal identifiers were removed from the blood samples and clinical data before transfer to the investigators.
We applied two distinct LC-MS based methods to distinct plasma aliquots for each study subject. Amino acids, amino acid derivatives, urea cycle intermediates, nucleotides and other positively charged polar metabolites were profiled as previously described . Briefly, 10 μL of plasma were extracted with 90 μL of 74.9:24.9:0.2 vol/vol/vol acetonitrile/methanol/formic acid containing valine-d8 (Sigma-Aldrich; St Louis, MO). After centrifugation, supernatants underwent hydrophobic interaction chromatography using a 150 × 2.1 mm Atlantis hydrophobic interaction chromatography column (Waters, Milford, MA), and MS data were acquired on a 4000 QTRAP triple quadrupole mass spectrometer (AB SCIEX, Foster City, CA) using ESI and MRM in the positive ion mode. Organic acids, sugars, bile acids, and other negatively charged polar metabolites were profiled as previously described . Briefly, 30 μL of plasma were extracted with the addition of four volumes of 80:20 vol/vol methanol/water containing isotope-labeled inosine-15 N4, supernatants underwent chromatography on a 150 × 2.0 mm Luna NH2 column (Phenomenex; Torrance, CA), and MS data were acquired using a 5500 QTRAP triple quadrupole mass spectrometer (AB SCIEX; Foster City, CA) using ESI and MRM in the negative ion mode.
We first examined the association between plasma metabolites and pre-selected metabolic phenotypes: diabetes status, BMI, serum albumin, total cholesterol, and systolic blood pressure (SBP). We focused on these phenotypes because of their strong links to metabolism and their established, but incompletely understood, associations with ESRD mortality [33–39]. To investigate differences in metabolite profiles across these phenotypes, continuous outcomes measures (BMI, albumin, total cholesterol, and SBP) were divided into tertiles. Diabetes status was classified as “yes” or “no”.
Because of the number of metabolites measured by our LC-MS platform (166 total, Additional file 1: Table S1), many with significant inter-correlations, we performed partial least squares discriminant analysis (PLS-DA) to visualize the linear components that discriminate individuals across the categories (tertiles or class) of metabolic phenotypes. PLS-DA was performed on log transformed and auto-scaled (mean-centered and divided by the standard deviation of each metabolite) data using MetaboAnalyst 2.0 [40, 41]. This program performs PLS regression using the plsr function provided by the pls package in R, and classification and cross validation using the corresponding wrapper function offered by the caret package. Variable Importance in Projection (VIP) scores generated by this program estimate the importance of each variable in the projection used within the PLS model. A variable with a VIP score greater than 1 can be considered important in a given model. The metabolites with the highest VIP scores were further analyzed by comparing their levels across classes using Mann-U-Whitney and Kruskal-Wallis tests as appropriate.
Finally, heatmaps were created to represent Pearson correlation (r) and P-value matrices between available clinical variables in the cohort and the subset of metabolites measured by our platform that have previously been described as uremic retention solutes (http://www.uremic-toxins.org) [27, 42–45].
To account for multiple testing, we used a Bonferroni corrected significance threshold of P < 0.0003 (0.05/166) for comparisons of individual metabolites across phenotype class. For the correlation matrices, we used an adjusted significance threshold of P < 0.00007 (0.05/[29 metabolites* 24 clinical variables]). A sensitivity analysis was performed to gauge the influence of case–control status (i.e. whether or not a given individual died of a cardiovascular cause within one year of starting dialysis) on the associations described herein. Stratified analysis by case status did not alter the statistical significance of the models described. Therefore, we did not stratify models by mortality status and analyzed the entire group of 200 individuals together as a cohort. All analyses were performed using SAS software version 9.1.3 (SAS Institute) and MetaboAnalyst 2.0 software (www.metaboanalyst.ca).
Results and discussion
As shown in Table 1, the mean age of the study population was 69.5 years and 69 % of subjects were white. There was an equal representation of males and females and nearly half of the individuals had a history of diabetes or had diabetes listed as their cause of ESRD (49 %). The mean BMI was 26.5 kg/m2 (SD ±7.6 kg/m2) and the mean SBP was 144 mmHg (±27 mmHg). A minority of individuals reported a lipid disorder (12 %) and the median total cholesterol level was 162 mg/dL (quartile1-quartile3, 127–188 mg/dL). The median serum albumin level was 3.6 g/dl (3.2-3.8 g/dL).
Examination of metabolite profiles and select metabolic phenotypes
The PLS-DA approach allowed us to visualize and extract the metabolites that best separated individuals (Figs. 1,2,3,4 and 5; left panels) across phenotype tertile or class (Table 2). Because four of the clinical phenotypes we studied (BMI, serum albumin, total cholesterol, and SBP) are continuous measures, we recognize that tertile cut-offs do not demarcate distinct physiologic or pathophysiologic processes. Thus, plot overlap across tertiles was expected. Variable importance in projection (VIP) provided a score for each metabolite, ranking the metabolites according to their predictive power in the PLS model; the fifteen metabolites with the highest VIP scores for each plot are shown in Figs. 1, 2, 3, 4 and 5 (right panels), and the levels of these metabolites across tertiles (or class) of the phenotypes with corresponding test statistics are shown in Tables 3-7.
Values are median peak area for the metabolites (quartile 1, quartile 3)*P-value significant at the Bonferroni adjusted level of 3.0 × 10−4
For individuals with or without diabetes, glucose was one of the strongest discriminating metabolites, with higher levels among the individuals with diabetes (Fig. 1). Xanthosine-5-phosphate (XMP) and vanillylmandelic acid (VMA), however, had even higher VIP scores than glucose, with XMP levels higher and VMA levels lower among the individuals with diabetes (P = 5.0 × 10−7 and 1.5 × 10−6 for XMP and VMA, respectively). Other significant differences included higher levels of 3-phosphoglyceric acid (P = 4.9 × 10−5), methionine (P = 4.5 × 10−5), mevalonate-5-phosphate (P = 2.8 × 10−4), and ornithine (P = 3.2 × 10−4) among individuals with diabetes. Tyrosine, an aromatic amino acid, had a nominal association with diabetes status (P = 1.0 × 10−3). There were no significant associations between tyrosine levels, or any aromatic or branched-chain amino acid levels, with BMI (Table 4). In fact, none of these metabolites had even a nominal association with BMI in univariate comparisons (data not shown). Overall, the ability of metabolite profiles to discriminate individuals across tertiles of BMI was poor (Fig. 2). Among the top metabolites by VIP score, only 20:4 carnitine (arachidonyl carnitine), β-hydroxybutyrate, and benzoate had nominal P-values in univariate analysis (Table 4).
Values are median peak area for the metabolites (quartile 1, quartile 3)*P-value significant at the Bonferroni adjusted level of 3.0 × 10−4
Values are median peak area for the metabolites (quartile 1, quartile 3)*P-value significant at the Bonferroni adjusted level of 3.0 × 10−4
Thus, our data from an ESRD population do not recapitulate observations derived in the general population linking diabetes and obesity with elevations in branched chain and aromatic amino acids (except for the nominal association with tyrosine noted above), short chain acylcarnitines, the glutamate/glutamine ratio, or bile acids . One potential explanation is that kidney function (and dysfunction) has a direct effect on select metabolites. The interaction between kidney function and amino acid metabolism, in particular, has been closely examined. For example, the kidney is known to make a substantial contribution to whole body tyrosine appearance via intra-organ phenylalanine hydroxylation . Further, the metabolic acidosis that results from renal failure leads to increased leucine oxidation , and indeed, generalized muscle catabolism. The significant associations we did identify with diabetes status in ESRD, e.g. XMP (a purine breakdown product) and VMA (a catecholamine metabolite), have not been reported in the general population. To what extent these alterations reflect distinct pathophysiologic processes in uremia warrant further investigation.
For individuals in different tertiles of serum albumin, PLS-DA demonstrated moderate separation across groups (Fig. 3). Tryptophan and its downstream metabolite kynurenine were significantly higher in individuals with higher serum albumin levels (P = 1.8 × 10−4 and 2.5 × 10−4, respectively). Two additional tryptophan metabolites, indoxyl sulfate and indole lactate, had a trend for higher levels among individuals with higher serum albumin levels. Given recent observations that ascribe a functional role for tryptophan metabolism, specifically through indoleamine-2, 3-dioxygenase (IDO), in modulating the immune system and vascular tone [48, 49], it is an appealing link between metabolism, nutritional status, inflammation, and cardiovascular risk in ESRD. However, these data should be interpreted with caution. Because tryptophan and its catabolites are hydrophobic, their positive correlation with albumin may reflect their significant protein-binding  – our LC-MS method measures total, not free, plasma metabolite levels. A similar mechanism may underlie the trend for association between thyroxine and albumin levels. By contrast, creatine, the metabolite with the third highest VIP score, had a trend for higher levels among individuals with lower serum albumin levels. In a study that examined spent media from cultured muscle cells treated with mitochondrial respiratory chain inhibitors, as well as plasma obtained from individuals with respiratory chain diseases, creatine levels were reproducibly elevated . In cell culture, extracellular creatine was inversely correlated with the intracellular phosphocreatine/creatine ratio, suggesting that elevated plasma creatine may signal a low energetic state in tissues using the creatine phosphate shuttle. Interestingly, elevated plasma creatine and decreased plasma phosphocreatine levels both have a nominal association with increased 1-year cardiovascular mortality among incident dialysis patients in ArMORR . These observations raise the possibility that impaired mitochondrial respiration in muscle is linked to the pathogenesis of hypoalbuminemia in dialysis.
For individuals across tertiles of total cholesterol and SBP, PLS-DA again demonstrated moderate separation across tertiles (Figs. 4 and 5). For both analyses, long-chain acylcarnitines were among the metabolites with the highest VIP scores, with 18:2 carnitine (linoleylcarnitine, P = 9.0 × 10−7), 18:1 carnitine (oleoylcarnitine, P = 6.0 × 10−7), and 18:0 carnitine (stearoylcarnitine, P = 7.7 × 10−5) all significantly higher in individuals with lower total cholesterol levels (Table 6,) and linoleylcarnitine (P = 5.5 × 10−5) and oleoylcarnitine (P = 8.7 × 10−6) significantly higher in individuals with lower SBP (Table 7); select short-chain and medium-chain acylcarnitines were also significantly elevated among individuals with lower SBP. Like albumin, low cholesterol and blood pressure are established predictors of increased mortality in ESRD [8, 9, 11, 12]. These findings are consistent with our observation that long-chain acylcarnitines are associated with 1-year cardiovascular mortality in ESRD , and motivate further interest in long-chain acylcarnitines as markers, or even mediators, of altered lipid metabolism and cardiovascular function in uremia. Tryptophan, one of the strongest discriminators of albumin tertile, was also a top hit in regards to total cholesterol and SBP, with higher tryptophan significantly associated with lower levels for both. Finally, higher arginine levels were significantly associated with higher total cholesterol and higher deoxyadenosine levels were significantly associated with lower SBP.
Relation of uremic solutes with twenty-four clinical phenotypes
To examine the specificity of the long-chain acylcarnitine-phenotype associations described above, we looked more broadly at the correlation between all of the metabolites captured on our platform that are uremic retention solutes and 24 phenotypes captured in ArMORR. Correlations between previously reported uremic retention solutes measured by our LC-MS platform, including long-chain acylcarnitines, and a broad array of clinical phenotypes measured in ArMORR are shown in Fig. 6 (full results are shown in Additional file 1: Tables S2 and S3). The correlations for glucose, as measured by the mass spectrometer, provide an anchor for interpreting this data – LC-MS glucose had a strong positive correlation with glucose measured by the clinical laboratory (r = 0.80, P <10−15), and a modest negative correlation with serum sodium (r = −0.21, P = 7.0 × 10−3) that did not reach statistical significance, but is well recognized in the clinic (i.e. dilutional hyponatremia due to hyperglycemia). Aside from the correlation between the two glucose measures, the strongest correlations identified by this broader examination of uremic metabolites were indeed between long-chain acylcarnitines and lipid traits (total and LDL cholesterol) and blood pressure. Among lipid phenotypes, long-chain acylcarnitines such as oleoylcarnitine were more strongly correlated with total (r = −0.45, P = 2.6 × 10−9) and LDL cholesterol (r = −0.43., P = 1.7 × 10−8) than with HDL cholesterol (r = −0.15, P = 5.7 × 10−2) or triglycerides (r = −0.27, P = 5.0 × 10−4). Similarly, oleoylcarnitine was more strongly correlated with SBP (r = −0.38, P = 3.4 × 10−8) than diastolic blood pressure (r = −0.25, P = 3.0 × 10−3). Long-chain acylcarnitines were not significantly correlated with any non-lipid or non-blood pressure phenotype.
Among the other uremic retention solutes included in our analysis, only kynurenic and anthranilic acid (both tryptophan metabolites downstream of IDO) had more than one significant correlation with a clinical variable (Fig. 6). Kynurenic acid was significantly correlated with diastolic blood pressure (r = 0.32, P = 5.0 × 10−6), bicarbonate (r = −0.31, P = 1.1 × 10−5), and phosphorous (r = 0.38, P = 5.7 × 10−8), whereas anthranilic acid was significantly correlated with alkaline phosphatase (r = 0.29, P = 6.3 × 10−5) and ferritin (r = 0.35, P = 9.7 × 10−7), a frequently used marker of inflammation in hemodialysis. Notably, the other metabolites that had a nominal association with ferritin included two other tryptophan metabolites, kynurenic acid (r = 0.16, P = 2.6 × 10−2) and quinolinate (r = 0.20, P = 6.3 × 10−3), and malondialydehyde (r = 0.24, P = 9.0 × 10−4), an oxidative stress marker. These findings reinforce a potential link between the IDO pathway and uremic inflammation. Sorbitol (r = −0.34, P = 3.4 × 10−5), cytidine (r = −0.40, P = 1.1 × 10−6), and uracil (r = −0.42, P = 1.7 × 10−7) were significantly correlated with plasma glucose; urate (r = −0.35, P = 5.2 × 10−6) and malondialdehyde (r = 0.38, P = 5.7 × 10−7) were significantly correlated with triglycerides; hippurate (r = 0.31, P = 7.0 × 10−6) and oxalate (r = 0.36, P = 1.6 × 10−7) were significantly correlated with phosphorous; and indoxyl sulfate (r = −0.28, P = 6.8 × 10−5) was significantly correlated with bicarbonate. Indoxyl sulfate, which has been found to have several proinflammatory effects in model systems , had no correlation with ferritin (nor transferrin saturation, nor white blood cell count), but had a trend for a negative correlation with BNP (r = −0.44, P = 4.0 × 10−4).
Although significant work has been devoted to understanding the in vitro effects of uremic retention solutes, less is known about how many of these metabolites relate to clinical traits in humans. Thus, in addition to demonstrating the specificity of select long-chain acylcarnitine associations, our examination of uremic metabolites against clinical phenotypes provides a resource for future uremic toxin research.
Our study has several limitations. First, by examining cross-sectional relationships between metabolite levels and various clinical phenotypes, our study is unable to address causation. Instead, by applying parsimonious selection procedures and conservative adjusted significance thresholds, we seek to highlight notable associations for future investigations. Second, study subjects were selected in the context of a case–control study of 1-year cardiovascular mortality, raising the possibility of confounding by case status. However, as noted in the Methods, the statistical significance of the findings are not changed when analyzed stratified by case status. Nevertheless, we acknowledge that our sample is not a random selection of incident dialysis patients, potentially limiting the generalizability of our findings. Third, residual confounding from other sources is likely to influence the associations described, as multivariate adjustments were not pursued. However, because the relationship between various clinical phenotypes and metabolite levels in uremia is unknown, we did not want to obscure novel biological associations by statistical adjustment. Thus, we present these raw associations as a framework for interpreting and appropriately adjusting select findings in future metabolomics studies in ESRD. Finally, our results require replication in an independent sample, ideally selected randomly and including both incident and prevalent dialysis patients.
In this report, we describe the small molecule alterations that accompany various metabolic phenotypes in ESRD. For phenotypes like diabetes and obesity that have been examined using metabolomics in the general population, our data reinforce the notion that unique pathophysiologic processes arise at ESRD onset. Our finding that long-chain acylcarnitines are strongly and inversely correlated with cholesterol and blood pressure corroborate their potential value as markers of cardiovascular risk in ESRD, whereas the association of tryptophan levels across several phenotypes is of interest given recent studies that assign a functional role for tryptophan metabolism in inflammation and vascular biology. Ultimately, we hope that the breadth of the data presented herein serves as a springboard for investigating mechanisms of uremic toxicity and identifying therapeutic targets to improve the unacceptably high morbidity and mortality attributable to ESRD.
End stage renal disease
Low density lipoprotein
Body mass index
Liquid chromatography-mass spectrometry
- ArMORR cohort:
Accelerated Mortality on Renal Replacement cohort
Systolic blood pressure
Partial least squares discriminant analysis
Variable importance in projection
b-type natriuretic peptide
Kopple JD. How to reconcile conventional and altered risk factor patterns in dialysis patients. Semin Dial. 2007;20(6):602–5.
Kalantar-Zadeh K, Kilpatrick RD, Kuwae N, McAllister CJ, Alcorn Jr H, Kopple JD, et al. Revisiting mortality predictability of serum albumin in the dialysis population: time dependency, longitudinal changes and population-attributable fraction. Nephrol Dial Transplant. 2005;20(9):1880–8.
Kovesdy CP, Kalantar-Zadeh K. Why is protein-energy wasting associated with mortality in chronic kidney disease? Semin Nephrol. 2009;29(1):3–14.
Cano NJ, Fouque D, Roth H, Aparicio M, Azar R, Canaud B, et al. Intradialytic parenteral nutrition does not improve survival in malnourished hemodialysis patients: a 2-year multicenter, prospective, randomized study. J Am Soc Nephrol. 2007;18(9):2583–91.
Dukkipati R, Kopple JD. Causes and prevention of protein-energy wasting in chronic kidney failure. Semin Nephrol. 2009;29(1):39–49.
Fellstrom BC, Jardine AG, Schmieder RE, Holdaas H, Bannister K, Beutler J, et al. Rosuvastatin and cardiovascular events in patients undergoing hemodialysis. N Engl J Med. 2009;360(14):1395–407.
Wanner C, Krane V, Marz W, Olschewski M, Mann JF, Ruf G, et al. Atorvastatin in patients with type 2 diabetes mellitus undergoing hemodialysis. N Engl J Med. 2005;353(3):238–48.
Kalantar-Zadeh K, Block G, Humphreys MH, Kopple JD. Reverse epidemiology of cardiovascular risk factors in maintenance dialysis patients. Kidney Int. 2003;63(3):793–808.
Lowrie EG, Lew NL. Death risk in hemodialysis patients: the predictive value of commonly measured variables and an evaluation of death rate differences between facilities. Am J Kidney Dis. 1990;15(5):458–82.
Kopple JD, Zhu X, Lew NL, Lowrie EG. Body weight-for-height relationships predict mortality in maintenance hemodialysis patients. Kidney Int. 1999;56(3):1136–48.
Iseki K, Miyasato F, Tokuyama K, Nishime K, Uehara H, Shiohira Y, et al. Low diastolic blood pressure, hypoalbuminemia, and risk of death in a cohort of chronic hemodialysis patients. Kidney Int. 1997;51(4):1212–7.
Port FK, Hulbert-Shearon TE, Wolfe RA, Bloembergen WE, Golper TA, Agodoa LY, et al. Predialysis blood pressure and mortality risk in a national sample of maintenance hemodialysis patients. Am J Kidney Dis. 1999;33(3):507–17.
Newgard CB, An J, Bain JR, Muehlbauer MJ, Stevens RD, Lien LF, et al. A branched-chain amino acid-related metabolic signature that differentiates obese and lean humans and contributes to insulin resistance. Cell Metab. 2009;9(4):311–26.
Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, et al. Metabolite profiles and the risk of developing diabetes. Nat Med. 2011;17(4):448–53.
Menni C, Fauman E, Erte I, Perry JR, Kastenmuller G, Shin SY, et al. Biomarkers for type 2 diabetes and impaired fasting glucose using a nontargeted metabolomics approach. Diabetes. 2013;62(12):4270–6.
Floegel A, Stefan N, Yu Z, Muhlenbruch K, Drogan D, Joost HG, et al. Identification of serum metabolites associated with risk of type 2 diabetes using a targeted metabolomic approach. Diabetes. 2013;62(2):639–48.
Cheng S, Rhee EP, Larson MG, Lewis GD, McCabe EL, Shen D, et al. Metabolite profiling identifies pathways associated with metabolic risk in humans. Circulation. 2012;125(18):2222–31.
Wewalka M, Patti ME, Barbato C, Houten SM, Goldfine AB. Fasting Serum Taurine-conjugated Bile Acids are Elevated in Type 2 Diabetes: and Do Not Change with Intensification of Insulin. J Clin Endocrinol Metab. 2014;99(4):1442–51. doi:10.1210/jc.2013-3367. Epub 2014 Jan 16.
Gerhard GS, Styer AM, Wood GC, Roesch SL, Petrick AT, Gabrielsen J, et al. A role for fibroblast growth factor 19 and bile acids in diabetes remission after Roux-en-Y gastric bypass. Diabetes Care. 2013;36(7):1859–64.
Ridaura VK, Faith JJ, Rey FE, Cheng J, Duncan AE, Kau AL, et al. Gut microbiota from twins discordant for obesity modulate metabolism in mice. Science. 2013;341(6150):1241214.
Wang-Sattler R, Yu Z, Herder C, Messias AC, Floegel A, He Y, et al. Novel biomarkers for pre-diabetes identified by metabolomics. Mol Syst Biol. 2012;8:615.
McCormack SE, Shaham O, McCarthy MA, Deik AA, Wang TJ, Gerszten RE, et al. Circulating branched-chain amino acid concentrations are associated with obesity and future insulin resistance in children and adolescents. Pediatr Obes. 2013;8(1):52–61.
Moore SC, Matthews CE, Sampson JN, Stolzenberg-Solomon RZ, Zheng W, Cai Q, et al. Human metabolic correlates of body mass index. Metabolomics. 2014;10(2):259–69.
Atzler D, Schwedhelm E, Zeller T. Integrated genomics and metabolomics in nephrology. Nephrol Dial Transplant. 2014;29(8):1467–74.
Rhee EP, Souza A, Farrell L, Pollak MR, Lewis GD, Steele DJ, et al. Metabolite profiling identifies markers of uremia. J Am Soc Nephrol. 2010;21(6):1041–51.
Sato E, Kohno M, Yamamoto M, Fujisawa T, Fujiwara K, Tanaka N. Metabolomic analysis of human plasma from haemodialysis patients. Eur J Clin Invest. 2011;41(3):241–55.
Kalim S, Clish CB, Wenger J, Elmariah S, Yeh RW, Deferio JJ, et al. A plasma long-chain acylcarnitine predicts cardiovascular mortality in incident dialysis patients. J Am Heart Assoc. 2013;2(6), e000542.
Gutierrez OM, Mannstadt M, Isakova T, Rauh-Hain JA, Tamez H, Shah A, et al. Fibroblast growth factor 23 and mortality among patients undergoing hemodialysis. N Engl J Med. 2008;359(6):584–92.
Teng M, Wolf M, Lowrie E, Ofsthun N, Lazarus JM, Thadhani R. Survival of patients undergoing hemodialysis with paricalcitol or calcitriol therapy. N Engl J Med. 2003;349(5):446–56.
Teng M, Wolf M, Ofsthun MN, Lazarus JM, Hernan MA, Camargo Jr CA, et al. Activated injectable vitamin D and hemodialysis survival: a historical cohort study. J Am Soc Nephrol. 2005;16(4):1115–25.
Thadhani R, Tonelli M. Cohort studies: marching forward. Clin J Am Soc Nephrol. 2006;1(5):1117–23.
Rhee EP, Clish CB, Ghorbani A, Larson MG, Elmariah S, McCabe E, et al. Wang TJ. Gerszten RE. A combined epidemiologic and physiologic metabolomics approach improves chronic kidney disease prediction. J AM Soc Nephrol: Fox CS; 2013.
Kilpatrick RD, McAllister CJ, Kovesdy CP, Derose SF, Kopple JD, Kalantar-Zadeh K. Association between serum lipids and survival in hemodialysis patients and impact of race. J Am Soc Nephrol. 2007;18(1):293–303.
Kalantar-Zadeh K, Kilpatrick RD, McAllister CJ, Greenland S, Kopple JD. Reverse epidemiology of hypertension and cardiovascular death in the hemodialysis population: the 58th annual fall conference and scientific sessions. Hypertension. 2005;45(4):811–7.
Kalantar-Zadeh K, Kopple JD, Kilpatrick RD, McAllister CJ, Shinaberger CS, Gjertson DW, et al. Association of morbid obesity and weight change over time with cardiovascular survival in hemodialysis population. Am J Kidney Dis. 2005;46(3):489–500.
Bradbury BD, Fissell RB, Albert JM, Anthony MS, Critchlow CW, Pisoni RL, et al. Predictors of early mortality among incident US hemodialysis patients in the Dialysis Outcomes and Practice Patterns Study (DOPPS). Clin J Am Soc Nephrol. 2007;2(1):89–99.
Trivedi H, Xiang Q, Klein JP. Risk factors for non-fatal myocardial infarction and cardiac death in incident dialysis patients. Nephrol Dial Transplant. 2009;24(1):258–66.
Kovesdy CP, Kalantar-Zadeh K. Introduction: the reverse epidemiology controversy. Semin Dial. 2007;20(6):485.
Rhee CM, Leung AM, Kovesdy CP, Lynch KE, Brent GA, Kalantar-Zadeh K. Updates on the management of diabetes in dialysis patients. Semin Dial. 2014;27(2):135–45.
Xia J, Mandal R, Sinelnikov IV, Broadhurst D, Wishart DS. MetaboAnalyst 2.0--a comprehensive server for metabolomic data analysis. Nucleic Acids Res. 2012;40(Web Server issue):W127–33. doi:10.1093/nar/gks374. Epub 2012 May 2.
Xia J, Psychogios N, Young N, Wishart DS. MetaboAnalyst: a web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 2009;37(Web Server issue):W652–60.
Vanholder R, Abou-Deif O, Argiles A, Baurmeister U, Beige J, Brouckaert P, et al. The role of EUTox in uremic toxin research. Semin Dial. 2009;22(4):323–8.
Vanholder R, Baurmeister U, Brunet P, Cohen G, Glorieux G, Jankowski J. A bench to bedside view of uremic toxins. J Am Soc Nephrol. 2008;19(5):863–70.
Vanholder R, De Smet R, Glorieux G, Argiles A, Baurmeister U, Brunet P, et al. Review on uremic toxins: classification, concentration, and interindividual variability. Kidney Int. 2003;63(5):1934–43.
Duranton F, Cohen G, De Smet R, Rodriguez M, Jankowski J, Vanholder R, et al. Normal and pathologic concentrations of uremic toxins. J Am Soc Nephrol. 2012;23(7):1258–70.
Tessari P, Deferrari G, Robaudo C, Vettore M, Pastorino N, De Biasi L, et al. Phenylalanine hydroxylation across the kidney in humans rapid communication. Kidney Int. 1999;56(6):2168–72.
Reaich D, Channon SM, Scrimgeour CM, Daley SE, Wilkinson R, Goodship TH. Correction of acidosis in humans with CRF decreases protein degradation and amino acid oxidation. Am J Physiol. 1993;265(2 Pt 1):E230–5.
Wang Y, Liu H, McKenzie G, Witting PK, Stasch JP, Hahn M, et al. Kynurenine is an endothelium-derived relaxing factor produced during inflammation. Nat Med. 2010;16(3):279–85.
Munn DH, Mellor AL. Indoleamine 2,3 dioxygenase and metabolic control of immune responses. Trends Immunol. 2013;34(3):137–43.
Jourde-Chiche N, Dou L, Cerini C, Dignat-George F, Vanholder R, Brunet P. Protein-bound toxins--update 2009. Semin Dial. 2009;22(4):334–9.
Shaham O, Slate NG, Goldberger O, Xu Q, Ramanathan A, Souza AL, et al. A plasma signature of human mitochondrial disease revealed through metabolic profiling of spent media from cultured muscle cells. Proc Natl Acad Sci U S A. 2010;107(4):1571–5.
Pletinck A, Glorieux G, Schepers E, Cohen G, Gondouin B, Van Landschoot M, et al. Protein-bound uremic toxins stimulate crosstalk between leukocytes and vessel wall. J Am Soc Nephrol. 2013;24(12):1981–94.
SK received support from the National Kidney Foundation Young Investigator award and NIH award KL2TR001100; EPR was supported by the Extramural Grant Program by Satellite Healthcare, a not-for-profit renal care provider, and NIH grant K08-DK-090142. REG received support from NIH 2R01DK081572 and NIH 1R01 HL098280; RT received support from NIH award K24 DK094872.
The authors declare that they have no competing interests.
SK and EPR designed the study, performed the statistical analysis, and wrote the manuscript. EPR and CBC performed the MS experiments. JJD, GO, ASM helped with sample processing, data analysis, and generating figures. REG aided with study design, data analysis, and manuscript drafting. RT established the cohort from which plasma samples were derived. All authors read and approved the final manuscript.
Metabolites measured within the sample set. Table S2. Pearson correlation coefficients for select metabolites and available clinical variables in the ArMORR cohort. Table S3. P-values from Pearson correlations using select metabolites and available clinical variables in the ArMORR cohort.
About this article
Cite this article
Kalim, S., Clish, C.B., Deferio, J.J. et al. Cross-sectional examination of metabolites and metabolic phenotypes in uremia. BMC Nephrol 16, 98 (2015) doi:10.1186/s12882-015-0100-y
- Uremic toxins