Abstract
Chronic kidney disease (CKD), a major public health problem, is associated with decreased anabolic response to insulin contributing to protein-energy wasting. Targeted metabolic profiling of the response to oral glucose tolerance testing (OGTT) may help identify metabolic pathways contributing to disruptions to insulin response in CKD. Using targeted metabolic profiling, we examined plasma metabolome in 41 moderate-to-severe non-diabetic CKD patients with estimated glomerular filtration rate (eGFR)<60ml/min per 1.73m2 (38.9±12.7) and 20 healthy controls with normal eGFR (87.2±17.7) before and after 2h of 75g oral glucose load. Compared to controls, CKD participants had higher lactate: pyruvate (L:P) ratio both at fasting and after oral glucose challenge. Total energy production estimated through GTP:GDP ratio was impaired during OGTT despite similar fasting GTP:GDP ratio. CKD group had sustained elevation of vitamin B family members, TCA cycle metabolites, and purine nucleotides in response to glucose challenge. Metabolic profiling in response to OGTT suggests a broad disruption of mitochondrial energy metabolism in CKD patients. These findings motivate further investigation into insulin sensitizers in patients with non-diabetic CKD and their impact on energy metabolism.
Introduction
Chronic kidney disease (CKD) has been recognized as a public health burden worldwide with the estimated global prevalence of 14% affecting over 650 million people globally [1]. Insulin resistance (IR) is one of the very early metabolic alterations in CKD increasing with severity of kidney disease [2]. IR is considered a cardio-metabolic risk factor correlated with increased systemic inflammation [3] and oxidative stress [4]. Furthermore, impaired glucose tolerance has been closely associated with increased risk of all-cause and cardiovascular mortality [5, 6].
Multiple pathophysiologic features of CKD contribute to IR in CKD [7] including chronic inflammation [8–10], altered gut microbiome [11], oxidative stress [12–14], metabolic acidosis [15], and accumulation of uremic toxins [11]. However, evidence for the biologic basis for IR in humans with CKD is lacking.
Studies in animal models indicate the IR related to CKD is the result of changes in insulin signaling both at the receptor and post-receptor levels [16] leading to blunted glucose uptake, metabolism, and storage in target tissues [17]. Peripheral tissues such as skeletal muscle are the main sites impacted by IR in CKD and link these early metabolic alterations with functional decline and frailty [18]. Compromised mitochondrial bioenergetics and physical performance has been reported in patients with CKD [19]. Post-receptor defects observed in CKD may be the consequence of excess insulin stimulation and heightened inflammation leading to augmented degradation of the insulin receptor by the ubiquitin proteosome pathway [20].
Although tissue insensitivity to insulin is the main contributor to IR, changes in insulin secretion and degradation have also been established in CKD [21]. In particular, a recent study using the gold-standard high-dose hyperinsulinemic-euglycemic clamp to investigate peripheral insulin sensitivity showed disturbances in insulin clearance as a principle characteristic distinguishing patients with non-diabetic CKD to controls [22]. This study demonstrated profound alterations in the plasma metabolome response to insulin in patients with CKD compared to controls.
However, understanding of the biologic basis for IR in humans with CKD at more physiologic, endogenous levels of insulin, such as that elicited by oral glucose tolerance testing (OGTT) remain unknown.
Oral glucose tolerance tests provide estimates of IR [23] but also captures additional physiologic processes by stimulating gut-derived incretin hormones known to augment insulin secretion [24, 25]. Metabolomic profiling of the response to OGTT may reveal specific metabolites and metabolic pathways underlying impaired glucose homeostasis in CKD thus identifying therapeutic targets for improving insulin sensitivity in this vulnerable population. We performed targeted plasma metabolic profiling comparing non-diabetic patients with moderate-severe CKD to healthy controls during an oral glucose challenge. In this study, we evaluated if CKD is associated with impaired plasma metabolic response to oral glucose challenge compared to controls. Our goal was to identify specific metabolic alterations in non-diabetic CKD to better understand the metabolic and biological basis for glucose intolerance and IR in CKD.
Results
Characteristics of the study participants
The study included 41 CKD participants with a mean eGFR of 38.9±12.7 ml/min per 1.73 m2, and 20 control subjects with a mean eGFR of 87.2±17.7 ml/min per 1.73 m2. The mean age was 61.2±12.9 years, 44% were women, and 22% were self- reported Black (Table 1). Compared with controls, participants with CKD had higher body weight, lower fat-free mass, higher BMI, and daily calorie intake, and greater plasma inflammatory markers.
Chronic kidney disease was defined as estimated glomerular filtration rate <60 ml/min per 1.73 m2; controls as ≥60 ml/min per 1.73 m2. SD, standard deviation; IQR, interquartile range; CCB, calcium channel blocker; RAASi, renin-angiotensin-aldosterone system inhibitor; eGFR, estimated glomerular filtration rate; CRP, C-reactive protein.
Oral glucose challenge (OGTT) was associated with a significant reduction in a wide range of metabolites primarily amino acids, purine nucleotides, and dicarboxylic acids. After adjusting for age, sex, race/ethnicity, body weight, and LC-MS batch, 65% (58/88) of the detected plasma metabolites were significantly altered post oral glucose challenge in the overall cohort (Table 2 and Supplemental Table 1). The largest reductions in individual metabolite concentrations post OGTT were observed for linoleic acid (fatty acid metabolism), ADP (nucleotide/purine nucleotide), and nicotinamide (vitamin) with a percent change of -82%, -79%, and -63%, respectively (Table 2). Only 6% (5/88) of the detected metabolites significantly increased post OGTT compared to fasting state. These metabolites included erythrose (sugar), glucose (glycolysis/sugar), and kynurenate (tryptophan cycle) with a percent change of +53%, +28%, and +17% respectively (Supplemental Table 2).
Differences in plasma metabolites during OGTT, compared with the fasting state
CKD attenuates the plasma metabolome response to glucose challenge in metabolites from predominantly the vitamin B family, TCA cycle intermediates, and purine nucleotides. After adjustment, changes in 15% (13/88) of detected plasma metabolites were significantly altered between CKD and controls post oral glucose challenge (Table 3 and Supplemental Table 2). Overall, CKD was associated with higher plasma levels of these metabolites such as succinate, inositol, nicotinamide, and glucose in response to OGTT compared to control. A notable exception to this was kynurenate. In general, controls had on average a greater decline in plasma metabolites in response to the oral glucose challenge compared to participants with CKD. The largest difference in response to the glucose challenge comparing CKD and controls were observed in ADP (nucleotide/purine nucleotide), glycochenodeoxycholate (bile acid metabolism), and IMP (nucleotide/purine metabolism). The ratio in fold change comparing CKD to control for these metabolites were 3.7, 3.6, and 3.1 respectively (Table 3).
Differences in plasma metabolites response post oral glucose tolerance test by chronic kidney disease status.
The response to the oral glucose challenge demonstrates a broad disruption of amino acid and mitochondrial energy metabolism in CKD patients. In the entire cohort, oral glucose challenge impacted alpha-Linoleic acid metabolism (p-value=3.84 x 10-31), aminoacyl-tRNA biosynthesis (p-value=1.15 x 10-26) and arginine biosynthesis (p-value=2.66 x 10-25) (Figure 1A). In addition, metabolic pathways involving amino acid metabolism such as branch chain amino acids (valine, leucine, isoleucine), phenylalanine, and histidine metabolism were also impacted. Compared to controls, participants with CKD demonstrated significant aberrations in alpha-Linoleic acid metabolism (p-value=2.91 x 10-5), nicotinamide metabolism (p-value=5.18 x 10-5), and arginine biosynthesis (p-value=1.96 x 10-4) (Figure 1B).
(A) in the overall cohort (n=61) and (B) comparing CKD (n=41) with controls (n=20). The size and color of the nodes represent pathway impact value and p-value, respectively.
CKD is associated with greater lactate to pyruvate (L:P) ratio, an indicator of mitochondrial respiratory chain impairment
A lactate to pyruvate ratio greater than 30 in patients suggests a primary respiratory chain dysfunction [16]. L:P ratio at Fasting and post-OGTT was significantly higher in participants with CKD compared to controls (Figure 2A). The mean L:P ratio at fasting was 28.9 (IQR of 26.2, 31.2) for controls compared to 36.7 (IQR of 29.3, 42.3) among participants with CKD (p-value=0.0012). The high L:P ratio was sustained post OGTT with a mean of 32.8 (IQR of 30.7, 37.2) in CKD in contrast to 29.5 (IQR of 27.8, 32.7) in controls (p- value<0.0023) (Figure 2A). The difference in L:P ratio is mostly driven by an increase in pyruvate levels however, lactate (fold change (FC), of 1.15 vs 1.20) and pyruvate (FC of 1.08 vs 1.34) elevations were not meaningfully different in CKD compared to controls during OGTT (Supplemental Table 2).
Panel A: L:P ratio at fasting and OGTT comparing CKD to controls. Panel B: GTP:GDP ratio at fasting and during OGTT comparing CKD to controls. Data points represent mean and error bars represent SD. Statistical analysis using Mann Whitney 2 tailed test *p-value<0.05, **p-value<0.01.
Persons with CKD have similar global energy production as the controls indicated by equivalent GTP:GDP ratio but an impaired response to OGTT. To assess the level of global energy production and the changes after the glucose challenge, we looked at the GTP:GDP ratio during fasting and OGTT separately. The GTP:GDP ratio is known to reflect the ATP:ADP ratio [26]. After the oral glucose challenge, it is expected that the GTP and ATP levels will increase and ADP and GDP levels will decrease, thus resulting in a greater GTP:GDP and ATP:ADP ratios.
The median GTP:GDP ratio at fasting was similar in CKD and controls (p-value=0.056) (Figure 2B). As expected, the GTP:GDP ratio increased during OGTT compared to fasting in both groups with the median of 3.18 to 3.47 in controls (p-value=0.042) and 2.97 to 3.24 in CKD (p- value=0.007); however, the ratio did not increase as sharply with the median of 3.47(IQR of 3.05, 3.95) and 3.24 (IQR of 2.824, 3.48) in controls and CKD respectively (p-value=0.044) (Figure 2B). This difference is more evident when considering GTP:GDP ratio of each group separately. Similar fasting GTP:GDP ratios and lower GTP:GDP ratio at OGTT suggest an impaired insulin response in persons with CKD.
Insulin secretion and sensitivity measured during OGTT did not differ among participants with and without CKD
We evaluated plasma insulin concentrations measured at 0, 10, 20, 30, 60, 90, and 120 minutes after glucose ingestion and found no significant differences between CKD and controls (Figure 3A). Similarly, the Matsuda index estimation of insulin sensitivity was not significantly different between CKD compared to controls (p-value=0.08) (Figure 3B). The median Matsuda index in CKD and controls was 3.93 (IQR of 2.44, 5.63) and 3.70 (IQR of 2.30, 9.13) respectively.
A) Plasma insulin concentration. Data points represent means and error bars represent 95% CI. B) Matsuda index.
Changes in plasma metabolites in response to the oral glucose challenge is associated with CKD signature and inflammation markers. We performed weighted gene co-expression network analysis (WGCNA) to better understand how closely related metabolite groups associate with CKD. We identified 8 metabolite modules in response to glucose load. Of the total 8 metabolite modules, those indicated in blue and green modules were significantly correlated with CKD signature and known CKD-associated inflammation markers (Figure 4A, 4B, and 4C). Green module positively correlated with CKD status, plasma cystatin C, and plasma creatinine and negatively correlated with eGFR (Figure 4B). Some notable metabolites from the green module include uric acid, kynurenate, and malonic acid (Supplemental Table 3). The blue module also positively correlated with plasma cystatin C and plasma creatinine (Figure 4C). Both modules also positively correlated with plasma TNF-α (Figure 4C). Blue module predominantly consisted of amino acids such as serine, phenylalanine, asparagine, and arginine (Supplemental Table 3).
a) Identified metabolite modules. Correlation between plasma metabolite modules and b) plasma biomarkers of kidney disease and c) inflammation markers. Each row in the table corresponds to a module, and each column represents a kidney biomarker, as indicated. The correlation coefficient is color-coded as indicated in the color key legend (red = positive and green = negative correlation) “***” = p- value<0.001, “**” = p-value<0.01, “*” = p-value<0.05, “.” = p-value<0.10. (d) The number of metabolites in each of the metabolite-module.
Discussion
Using targeted metabolic profiling, we identified several attenuated biological and metabolic pathways in response to an oral glucose challenge in the overall cohort with evidence of marked heterogeneity by CKD status. First, we observed that compared to the fasting state OGTT resulted in a significant decrease in a wide range of metabolites including purine nucleotides and amino acids in the overall cohort. Second, patients with CKD demonstrate an attenuated plasma metabolome response to oral glucose challenge. Our prior findings in this cohort were confirmed that these changes occur in the absence of significant differences in insulin secretion and insulin sensitivity (Figure 3A and 3B) [22]. In response to the glucose challenge, we observed changes in 13 metabolites significantly differ in participants with CKD compared to controls, a significant portion of which are known uremic toxins associated in persons with CKD. Third, we observed an attenuated plasma response to OGTT particularly in vitamin B family members, mitochondrial energy metabolism, and purine metabolism. Lastly, we identified a set of metabolites significantly associated with CKD status and inflammatory markers. Together, our findings suggest CKD is associated with suppressed anabolic response to glucose challenge consistent with prior findings of impaired anabolic response to insulin in this population [52].
The application of metabolomics to the oral glucose tolerance testing demonstrated several insights into the metabolic response to a glucose challenge. First, OGTT resulted in a reduction in the concentrations of a broad range of plasma metabolites 2 hours post-glucose challenge compared to fasting suggesting an anabolic response to insulin. Second, we noted marked reduction in plasma levels of TCA cycle intermediates suggesting TCA cycle activation consistence with the known high energetic demands of protein translation and fatty acid synthesis [27] (Table 2). Third, we observed a decrease in all detected amino acids and their metabolites suggesting suppression of proteolysis and cataplerosis in the overall cohort (Table 2). In particular, we detected significant decreases in succinate, oxaloacetate, aconitate, GDP, malonate, and ADP consistent with the known action of insulin in stimulating glucose metabolism and inhibition of catabolic pathways such as lipolysis, proteolysis and ketogenesis (Table 2).
We observed significant heterogeneity by CKD status in the plasma metabolome response to glucose load revealing underlying broad disruption in energy metabolism. In animal models, CKD leads to a post-insulin receptor defect contributing to IR [28, 29]. Studies have shown that IR is a major contributor to altered metabolism in CKD with adverse consequences for lipid [30], carbohydrate [31], and protein metabolism [32] impacting global energy metabolism [20]. Our investigation identified CKD-associated disruptions in GDP, ADP, and succinate; all metabolites involved in the TCA cycle. We showed these metabolites were significantly higher in CKD compared to controls post-OGTT suggesting an impaired insulin response to activate the TCA cycle (Table 3). These findings are also consistent with previous study demonstrating diminished TCA cycle activity in patients with non-diabetic CKD [33]. Our findings underscore the importance of IR associated with metabolic defects in CKD and its potential contribution to mitochondrial dysfunction. This association between IR and mitochondrial energy metabolism was demonstrated in a recent study suggesting that muscle specific insulin receptor knockout in a murine model impaired mitochondrial respiration, decreased ATP production, and increased reactive oxygen species [34]. This may be driven by decreased complex 1 (NADH oxidase) driven mitochondrial respiration and supercomplex assembly which may subsequently manifest as an increase in NADH/NAD ratio.
One prominent feature of mitochondrial dysfunction we noted in CKD versus control was elevated lactate-to-pyruvate (L:P) ratio. The L:P ratio is known to be in near equilibrium with the NADH/NAD ratio with an elevated lactate-to-pyruvate (L:P) ratio serving as an indicator for impaired mitochondrial function, TCA cycle, and pyruvate decarboxylase activity; a highly regulated subunit of pyruvate dehydrogenase (PDH) [35–37]. A L:P ratio of >25 has been suggestive of respiratory chain dysfunction [38]. Increased levels of pyruvate, lactate, and L:P ratio has also been shown in patients with acute kidney failure [39]. Glycolysis products, lactate and pyruvate, were increased in both groups in response to the glucose challenge signaling glycolysis activation. Lactate and pyruvate however were also both more elevated in CKD compared to the controls suggesting reduced activity of pyruvate dehydrogenase complex (PDH) (Supplemental Table 2). PDH is a key insulin-mediated enzyme complex with activity in the muscle which is acutely stimulated by insulin [40]. Indeed, reduced activity of insulin- stimulated PDH has been shown in skeletal muscle of patients with type 2 diabetes [40, 41].
Reduced muscle PDH activity in patients with stage 4 and 5 CKD has also been reported as a contributing factor to muscle weakness [42]. We also observed an elevated (L:P>30) lactate: pyruvate ratio in CKD compared to controls both at fasting and post-OGTT suggesting a disorder in mitochondrial respiratory chain complex (Figure 2A). Taken together our findings indicate persons with CKD have impaired mitochondrial metabolism evidenced by metabolic profiles suggesting blunted activation of TCA cycle, redox imbalance, as well as disruption in the mitochondrial respiratory chain complex.
CKD was associated with disruption in purine and pyrimidine metabolism, a crucial source of necessary energy and cofactors needed for bioenergetics and biomolecular demands of metabolism in addition to serving as building blocks of DNA and RNA [43]. Purines are found in biomolecules such as ATP, GTP, cyclic AMP, NADH, and coenzyme A [44]. We noted several indicators of impairment in the purine nucleotide cycle. First, four metabolites involved in purine nucleotide cycle GDP, ADP, adenylosuccinate, and IMP all had attenuated reduction in response to the glucose challenge in CKD compared to controls suggesting an impairment in purine nucleotide cycle (Table 3). Fluctuations in ATP, ADP, and AMP and corresponding changes in GTP and GDP are regulated by purine nucleotide cycle, particularly in skeletal muscle [45, 46]. During muscle contraction in the skeletal muscle a high ATP:ADP and GTP:GDP ratio is maintained during high ATP usage [34]. The purine cycle also enhances glycolysis and the TCA cycle via enhancing the rate of glycolysis by activating phosphofructokinase (PFK) and production of fumarate serving as an anaplerotic substrate supporting the TCA cycle [47, 48]. Second, to assess disruptions in its regulatory role of energy molecules during glycolytic changes, we investigated changes in GTP:GDP ratios at fasting and during OGTT (where GTP:GDP ratio is expected to increase [26, 49]) in both groups. We found that patients with CKD had similar fasting GTP:GDP ratios compared to controls, but significantly lower GTP:GDP ratio post-OGTT suggesting an impairment in the purine nucleotide cycle interfering with energy generation. Purine nucleotide cycle has not been previously explored in kidney disease and warrants further investigation in patients living with CKD given its role in various metabolic pathways fundamental in skeletal muscle energy metabolism.
CKD was also associated with disruption of metabolites in the vitamin B family with glucose challenge. Vitamin B family members are essential for metabolism involved in major metabolic pathways acting as cofactors in catabolic and anabolic pathways[50]. We found a decrease in vitamin B family derivatives in repones to OGTT in both groups. Among vitamin B family members, the most significant decrease post OGTT was observed in nicotinamide (a form of vitamin B3), important in various oxidation/reduction reactions. Compared to controls where nicotinamide levels profoundly decreased in response to glucose challenge, the levels of nicotinamide in participants with CKD remained elevated (Table 3) consistent with our prior findings using insulin clamp testing [51]. Our data confirms findings from these other studies showing nicotinamide and other NAD+ catabolites accumulating in uremic patients suggesting an insulin mediated defect to use nicotinamide during anabolic reactions. Dynamic changes in other B vitamins were also detected post glucose challenge (Supplemental Table 2). Biotin (vitamin B7) needed for energy metabolism, fat synthesis, amino acid metabolism and glycogen synthesis [38] significantly dropped in the controls whereas a slight increase was observed in CKD post glucose challenge (Table 3). Similarly, pyroxidal phosphate, the active form of vitamin B6 needed for glycine and serine metabolism in humans [38]; pantothenic acid (vitamin B5) needed for carbohydrate, protein, and fat metabolism also remained elevated post glucose challenge in CKD compared to controls [52] (Supplemental Table 2). The relevance of these changes is underscored by the essential role these metabolites play in metabolic pathways influenced by insulin action, such as amino acid metabolism (Vitamin B6), TCA cycle (Vitamin B5), and various redox reactions (Vitamin B3). Together, elevated levels of vitamin B family members points toward an impairment in activating insulin mediated anabolic pathways involving accumulation of these cofactors in patients with CKD.
There were several metabolic profile changes in response to the glucose challenge characteristic of CKD [53] and levels of CKD associated inflammation markers identified using WGCNA analysis [54, 55]. First, a module consisting of malonate, leucic acid, IMP, kynurenate, uric acid and inositol had a significant positive correlation with plasma creatinine, plasma cystatin C, CKD status, and a significant negative correlation with eGFR (Figure 3B). This module also had a significant positive correlation with TNF-α levels (Figure 3C). Interestingly, uremic retention solutes and uremic toxins including inositol, malonate, and uric acid make up a considerable number of correlating metabolites in this module.
Inositol is a biologically relevant metabolite involved in energy metabolism and considered a uremic retention solute. It was attenuated post-OGTT in the CKD group compared to controls, and it’s also part of the green module correlating with CKD status and TNF-α levels. Inositol is a polyol produced from glucose-6-phoshate and a water-soluble low molecular weight solute playing an important role in a cellular signaling [56]. In humans, the kidney is the main organ responsible for biosynthesis and catabolism of myo-inositol [57, 58]. Serum inositol levels in as been negatively correlation with GFR [59–61] and implicated to have an adverse impact on renal progression. Among patients with type 2 diabetes, elevated plasma concentration of myo-inositol has been associated with greater likelihood of progression to end stage renal disease (ESRD) [62]. An in vitro study has also shown inositol as an inhibitor of mitochondrial fission by direct inactivation of 5’-AMP-activated protein kinase (AMPK) [63], a primary energy sensor and regulator of cellular metabolism [64]. Our study, confirms and builds on the previous studies confirming elevated myo-inositol levels in patients with non-diabetic CKD. Further studies are needed to investigate the role of altered inositol metabolism in defective mitochondrial dynamics contributing to mitochondrial dysfunction and cell damage in CKD.
Two other known uremic toxins characterizing differences in the response to glucose challenge in CKD compared to controls were malonate and uric acid. Elevated malonate, a potent inhibitor of cellular respiration, has been reported in ESRD patients[65, 66]. Malonate binds to the active site of succinate dehydrogenase (Complex II) of the electron transport chain without reacting thus competing with the usual substrate succinate and disrupting oxidative phosphorylation contributing to mitochondrial dysfunction in CKD [67]. High serum uric acid levels are very common in CKD patients [68] and has been shown to be correlated with greater IR, greater C- reactive protein levels, and lower eGFR in patients with CKD [69, 70]. Uric acid activates NADPH oxidase enzyme increasing intracellular oxidative stress, mitochondrial injury, and ATP depletion [71, 72]. The protein bound solute kynurenate was also significantly altered post glucose challenge in the green module (Table 3, Figure 3B, Figure 3C). Elevated serum kynurenate levels have been associated with reduced kidney function[73], greater TNF-α levels [74] and greater IR [75] and is a well-documented marker for kidney disease [73, 74, 76, 77]. In summary, WGCNA analysis suggests uremic toxins are strongly associated with disruption in response to glucose challenge in CKD compared to controls and linked to inflammatory markers in our cohort.
Our study has several notable strengths but also limitations. First, we assessed dynamic changes in response to glucose challenge and accounted for several potential confounding factors by adjusting for age, sex, race, body weight, and batch in our analysis. Second, we applied targeted metabolic profiling to assess a broad range of attenuated metabolic pathways impacted by physiological effects of insulin and identified its potential disturbed mechanisms in CKD. Our study was not without limitations. We relied on plasma metabolite levels, so we had no direct evidence of metabolic alterations inside the cells or any specific tissue. Differences in gut microbiome can influence circulating metabolites and the changes we observed upon glucose loading could be in part due to differences in intestinal microbiome. We can’t rule out the potential impact of residual confounding by differences in unmeasured characteristics between CKD and controls. Finally, we can’t precisely identify the impact of potential differences in the incretin induced secretion of insulin between our two groups. However, our data confirms prior published study from our cohort where we demonstrated no significant difference in the plasma insulin levels and insulin sensitivity during OGTT comparing CKD to controls (Figure 3A and 3B) [22]. Further studies are necessary to interrogate the association of reduced activation of TCA cycle, redox imbalance, and electron transport chain efficiency in human tissues especially in the skeletal muscle.
In summary, our results show that numerous plasma metabolites are altered in response to oral glucose challenge in CKD. We showed that response to glucose challenge in CKD is associated with disruption in plasma levels of TCA cycle intermediates, purine nucleotide cycle metabolites, and vitamin B family members relative to controls. This abnormal metabolic profile in response to glucose challenge adds to prior studies indicating the IR in CKD is predominantly characterized by a depressed anabolic response and disruption in energy metabolism. We also showed that disruptions in plasma metabolic profile in response to glucose load, in particular uremic toxins and retention solutes, are the main metabolites correlating with CKD signature and inflammatory biomarkers. Given recent evidence for clinical benefits of insulin sensitizers, there is an urgent need for future studies interrogating the pleotropic effects of insulin sensitizers on mitochondrial energy metabolism and anabolism in CKD.
Methods
Population
The Study of Glucose and Insulin in Renal Disease (SUGAR) is a cross-sectional study of glucose and insulin metabolism in moderate-to-severe nondiabetic CKD. From 2011– 2014, we recruited 98 participants from Nephrology and Primary Care clinics associated with the University of Washington and neighboring institutions in Seattle (Washington, USA). Among the 98 recruited participants, 95 had adequate plasma samples collected for metabolomics, and of these, 61 had plasma samples before and after oral glucose tolerance testing(41 CKD and 20 controls) [22]. Sixty-one participants were included from the SUGAR study where 41 of them were moderate to severe non-diabetic CKD patients (eGFR < 60 ml/min per 1.73 m2) and 20 healthy controls (eGFR > 60 ml/min per 1.73 m2). Plasma lipid and non-lipid metabolites were measured after overnight fasting, 2h after a 75-gram oral glucose load for OGTT and during the clamp as reported previously [52]. Plasma biomarkers of kidney function and inflammation were measured in the fasting blood. Albumin excretion rate was calculated from 24h urine sample. A detailed description of the study design and distribution of the clinical phenotypes between CKD patients and healthy controls has been reported earlier [22].
Metabolomics
Targeted metabolomics based in a liquid chromatography-tandem mass spectroscopy (LC-MS/MS) platform was performed at the Northwest Metabolomics Research Center as described previously [78]. The LC-MS/MS analyses were performed using an Agilent 1260 LC (Agilent Technologies) with a SeQuant ZIC-cHILIC column (150 × 2.1 mm × particle size 3.0 μm, Merck KGaA), which was coupled to a AB Sciex QTrap 5500 MS (AB Sciex) system equipped with a standard electro-spray ionization (ESI) source. The details for sample preparation, sample collection times and reagents has been described earlier [51].
Statistical analysis
Plasma non-lipid metabolite data were available from 61 subjects at fasting and 2h post-oral glucose load. Raw metabolite data were converted to ratio of 2h post-oral glucose load and corresponding fasting values and checked for batch effect of the ratio values. All clinical and metabolite data were checked for normality. To account for drift of sample preparation batches, raw metabolite data was normalized using Systematic Error Removal Using Random Forest (SERRF)[79]. To account for the correlation of measurements within participants, we examined the fold changes associated with the OGTT procedure via a linear mixed model with random intercepts, regressing the log transformed SERRF normalized metabolite on the sample type (during OGTT versus fasting sample), adjusting for age, sex, race/ethnicity, body weight and batch. To evaluate whether the effect of the oral glucose challenge procedure differed between CKD and non-CKD participants, we used linear mixed effect modeling with random intercepts where SERRF normalized metabolites were regressed on a sample type (OGTT vs fasting), CKD status, and their interaction, additionally adjusting for covariates listed above. Data are presented as means ± SD unless otherwise indicated. Statistical analysis was performed using R 3.6.1 for windows release [80]. P-values were adjusted for multiple comparisons using the “Benjamini Hochberg” appraoch, and an adjusted p-value <0.05 was considered significant for all analyses unless stated otherwise.
Pathway-associated metabolite sets enrichment and metabolite pathway analysis were performed using MetaboAnalystR v4.0 [81] with the KEGG human metabolite database [82]. Pathway analysis using SERRF normalized pre OGTT and post OGTT metabolite level was used to evaluate metabolic changes in response to OGTT in the entire cohort. Significantly altered metabolic pathway by CKD status was determined by comparing the changes in metabolite levels calculated by subtracting log-transformed post OGTT SERRF normalized values to pre OGTT time points.
Differences in insulin sensitivity (Matsuda index) and plasma metabolites including GTP:GDP ratio and lactate: pyruvate ratio between fasting and 2h post-oral glucose load were evaluated separately by Mann-Whitney U test using GraphPad Prism 9.0 (GraphPad Software, Inc., San Diego, CA).
Metabolite data were checked for excessive missing values by using the R package WGCNA’s [83] “goodSampleGenes” test. The correlation between kidney biomarkers was determined by Spearman correlation. The ratio of plasma metabolites post oral glucose to pre glucose challenge was used for correlation network analysis using WGCNA R package. A soft threshold approach was used with a power of 6 (based on scales free topology) in a WGCNA default unsigned network with dynamic tree cutting (deep split = 2) and a min Module Size = 3 as parameters for the dynamic tree cut function [84]. The module eigengene, defined as the first principal component (PC) of a module’s metabolite concentration matrix, was used to calculate the Spearman correlation between a metabolite module and kidney biomarkers.
Author Contributions
Conceptualization was contributed by AA, BR, MNH, LRZ, DR, and IHDB. Methodology was contributed by AA, MNH, LRZ, DR, and BR. Formal analysis was contributed by AA and MNH. Investigation was contributed by BR, MNH, and DR. Resources were contributed by IHDB, BR, and DR. Data curation was contributed by AA and MNH. Writing the original draft was contributed by AA, BR, and MNH. Writing the review and editing were contributed by AA, BR, MNH, LRZ, IHDB, BJB, JG, LRD, MC, DR. Visualization was contributed by AA and MNH. Supervision was contributed by IHDB, BR, BJB, and DR. Project administration was contributed by IHDB and DR. Funding acquisition was contributed by BR, IHDB, DR, and BJB.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Supplemental Figures
Differences in plasma metabolites in response to glucose load compared to fasting. Results are from a regression SERRF normalized metabolites on sample type (OGTT vs fasting) adjusted for age, sex, race, weight, and batch. A p-value of <0.05 was used to determine significance (in bold). The fold change is in response to oral glucose is the adjusted fold change associated with OGTT compared to fasting state (e.g., a fold change of 1.28 indicates a 28% increase in metabolite levels).
Differences in plasma metabolic response post glucose challenge by CKD status. Result from regression analysis using SERFF normalized on sample type (fasting vs OGTT) by CKD status adjusted for age, sex, race, weight, and batch. Fold changes represent changes in metabolite levels with fasting levels after glucose load. (e.g. a fold change of 0.92 indicates a reduction of %8 in metabolite level. P-value for interaction represents the heterogeneity in fold change by disease status.
The list of metabolites in each module from the WGCNA analysis.
Acknowledgements
This project was funded in part by an unrestricted grant from the Northwest Kidney Centers and R01DK087726, R01DK087726-S1, R01DK129793-01 (BR), R01DK087726, R01DK087726- S1, K01 DK102851 (JAA, IHDB), K23DK100533 (JLG), K24 DK096574 (TRZ), and P30 DK017047 (University of Washington Diabetes Research Center) and Dialysis Clinics Incorporated, C-4122 (BR).
Footnotes
Armin Ahmadi: arahmadi{at}ucdavis.edu; Nazmul Huda: mnhuda{at}ucdavis.edu; Brian Bennett:, brian.bennett{at}ucda.gov; Jorge Gamboa: Jorge.gamboa{at}vumc.org; Leila Zelnick: LZelnick{at}nephrology.washington.edu; Daniel Raftary: draftary{at}uw.edu; Ian de Boer: IDeboer{at}nephrology.washington.edu; Lucas R. Smith: lucsmith{at}ucdavis.edu; Maria Chondronikola: mchondr{at}ucdavis.edu; Baback Roshanravan: broshanr{at}ucdavis.edu
Conflict of interest statement: The authors have declared that no conflict of interest exists.
References:
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.
- 10.↵
- 11.↵
- 12.↵
- 13.
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.↵
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵