ChemHelper 4.4 Activator incl

Chronic kidney disease (CKD) results in the accumulation of metabolic waste products that are normally cleared by the kidney, known as uremia. Many of these waste products are from bacteria metabolites in the gut. Accumulation of uremic toxins in plasma and tissue, as well as the gut-plasma-tissue metabolic axis are important for understanding pathophysiological mechanisms of comorbidities in CKD. In this study, an untargeted metabolomics approach was used to determine uremic toxin accumulation in plasma, liver, heart and kidney tissue in rats with adenine-induced CKD. Rats with CKD were also given AST-120, a spherical carbon adsorbent, to assess metabolic changes in plasma and tissues with the removal of gut-derived uremic toxins. AST-120 decreased >55% of metabolites that were increased in plasma, liver and heart tissue of rats with CKD. CKD was primarily defined by 8 gut-derived uremic toxins, which were significantly increased in plasma and all tissues. These metabolites were derived from aromatic amino acids and soy protein including: indoxyl sulfate, p-cresyl sulfate, hippuric acid, phenyl sulfate, pyrocatechol sulfate, 4-ethylphenyl sulfate, p-cresol glucuronide and equol 7-glucuronide. Our results highlight the importance of diet and gut-derived metabolites in the accumulation of uremic toxins and define the gut-plasma-tissue metabolic axis in CKD.

The National Kidney Foundation defines chronic kidney disease (CKD) as the presence of kidney damage or a progressive decline in renal function lasting for three or more months. The prevalence of CKD is estimated to be upwards of 14% in the United States2,,. CKD is divided into 5 stages, which are classified based on glomerular filtration rate. In end stage renal disease, patients require dialysis to sustain life, and the resulting condition is known as the “residual syndrome”. Although biochemical mechanisms of the residual syndrome are not completely understood, uremia is considered a major cause of this condition.

Uremia is the accumulation of solutes that are normally cleared by the kidneys in CKD. At high concentrations, the accumulating solutes are considered uremic toxins as they begin to induce the residual syndrome. Many of these highly concentrated toxins are gut-derived metabolites generated via protein fermentation by bacteria in the colon. Gut-derived aromatic amino acid (AAA) metabolites, such as indoxyl sulfate and p-cresyl sulfate, have been extensively studied in CKD,,,. Indoxyl sulfate and p-cresyl sulfate play a significant role in cardiovascular disease, which is the leading cause of death in patients with CKD accounting for 45% of mortalities,,. Indoxyl sulfate concentration is correlated with aortic calcification and vascular stiffness in patients with CKD. High plasma levels of unbound p-cresyl sulfate increases the risk of cardiovascular and all-cause mortality in hemodialysis patients.

AST-120 is a spherical carbon adsorbent that has been used in Japanese CKD patients since 1991 to remove the precursors of uremic toxins produced in the gut. This compound is formulated to bind small molecules (100 to 10,000 Da) while larger molecules such as digestive enzymes and hormones are left unbound. Adsorption of gut-derived uremic precursors to AST-120 decreases circulating levels of these toxins. In clinical trials, AST-120 has been shown to mitigate symptoms of uremia in CKD,19.

A number of studies have quantified individual uremic toxin levels in patients with CKD. In 2003, the European Uremic Toxin Work Group compiled mean and median serum or blood levels of 90 uremic toxins using meta-analysis. Subsequent updates were published in 2007 and 2012, adding 14 and 56 uremic toxins, respectively,. This is an effective method for determining mean patient uremic toxin levels; however, the uremic toxins reported are a result of targeted analysis, which can limit the identification of novel toxins. Recent advances in mass spectrometry allow for identification of novel uremic toxins using an untargeted metabolomics approach.

Metabolomics is the study of all metabolites generated in a biological system. Thousands of metabolites can be detected using untargeted metabolomics without prior knowledge of their composition. In the setting of CKD, metabolomics can provide uremic signatures directly associated with decreased renal function. This information can be used to determine mechanisms of kidney disease initiation and progression. Metabolomics has been applied to diabetic nephropathy, polycystic kidney disease, acute kidney injury and renal cell carcinoma in patients. Recently, clinical and animal studies have begun to identify plasma uremic toxins using untargeted metabolomics,,. These metabolites have been identified using both NMR and LC-MS techniques. Advances in LC-MS using ultra performance liquid chromatography (UPLC) coupled with quadrupole time-of-flight mass spectrometry (UPLC-QToF/MS) allows for high throughput and highly sensitive metabolomics analysis. The advent of UPLC has resulted in increased chromatographic resolution and improved compound separation. Current QToF/MS techniques allow alternating scans of low and high collision energy for simultaneous acquisition of precursor and fragment ions in a single run, thereby improving workflow.

Plasma uremic toxin levels have been correlated with CKD complications, most notably are the cardiovascular effects including endothelial dysfunction and heart failure,. Metabolic changes in the kidney are known to occur as a result of kidney damage, and the accumulation of uremic toxins has been shown to further kidney injury,. However, effects of accumulating toxins on metabolic changes in other major organs such as the heart and liver have not been assessed. Therefore, the purpose of this study was to determine the metabolic changes in plasma, liver, heart and kidney tissue as well as the metabolic axis between these matrices in rats with CKD. We also assessed the plasma and tissue metabolic effects of removing gut-derived uremic toxins using AST-120 to further understand the gut contribution to the uremic condition in both plasma and tissues in CKD.

Validation of Experimental Model

The extent of CKD was assessed by plasma creatinine and urea levels (). Creatinine levels were 8.0-fold and 5.8-fold higher in CKD and CKD+AST-120 groups, respectively, compared to controls (p < 0.05). CKD and CKD+AST-120 rats had a 6.8 and 5.1-fold increase in plasma urea levels compared to controls. There were no significant differences in weights between each group.

Physical and biochemical characteristics of control, CKD and CKD+AST-120 rats.

Principle Component Analysis (PCA) and the Removal of Metabolites by AST-120 in Biological Matrices

To determine the metabolic changes in CKD with and without the administration of AST-120, samples were analyzed by ultra performance liquid chromatography with quadrupole time-of-flight mass spectrometry (UPLC-QToF/MS) in positive and negative electrospray ionization (ESI) modes. All samples, including pooled samples, were evaluated by multivariate analysis. Representative UPLC-QToF/MS chromatograms as well as an overlay of extracted ion chromatograms for all metabolites identified in negative and positive ESI modes are demonstrated by effective chromatographic resolution (Supplementary Figure S1). Pooled sample injections demonstrated reproducibility and were clustered in the center of principle component analysis plots (Supplementary Figure S2 and S3). The first component in plasma negative mode PCA did not show clustering of experimental groups. When using second and third components, CKD rat plasma samples appeared as a separate cluster compared to control and CKD+AST-120 groups (). The control and CKD+AST-120 samples were also grouped together in this analysis. The majority of metabolites were increased in CKD relative to control and CKD+AST-120 groups. Using XCMS Online multi-group experiment, 26 metabolites were significantly increased in CKD plasma samples compared to control. Of those 26 metabolites, 24 were significantly reduced by the administration of AST-120 to rats with CKD (). The PCA for liver metabolites in negative mode also demonstrated clustering of CKD samples; however, CKD+AST-120 showed some separation from control samples. Further analysis revealed that AST-120 removed 60.7% of metabolites that were significantly increased in CKD rat livers (). Positive ESI did not show marked differences between groups for plasma and liver samples. The metabolic profile in heart tissue was distinctly different between control and CKD rats with complete separation of samples in the PCA for both negative and positive ESI modes (C and A). CKD+AST-120 rat samples were clustered between control and CKD groups. AST-120 significantly decreased 58 of 93 metabolites in heart tissue that were increased in CKD (G and C). Metabolites between control, CKD and CKD+AST-120 kidney tissue were mutually exclusive in the PCA for negative ESI mode showing the greatest separation between all matrices analyzed (). In positive ESI mode, control and CKD components were well separated with slight overlap between CKD and CKD+AST-120 groups. Kidney tissue had the largest number of significantly increased metabolites in CKD compared to control. However, AST-120 reduced less than 40% (90 of 233) of metabolites that were significantly increased in CKD compared to control kidneys (H and D).

Principle component analysis and metabolites increased in CKD and decreased by AST-120 in negative ESI mode.

Principle component analysis and metabolites increased in CKD and decreased by AST-120 in positive ESI mode.

The metabolites responsible for the differences between control and CKD rats were assessed by OPLS-DA in plasma, liver, heart and kidney. Plasma metabolites in control and CKD samples were well described by OPLS-DA (R2(Y) = 0.92) with high predictive ability (Q2(Y) = 0.75). The majority of metabolites in plasma samples were increased in CKD relative to control. Many of these metabolites were AAA derivatives including phenyl sulfate, indoxyl sulfate, p-cresyl sulfate, 4-ethylphenyl sulfate, hippuric acid and p-cresol glucuronide ( and ). Tissue samples from rats with CKD also had high levels AAA derivatives. Few metabolites were significantly higher in control plasma samples compared to CKD. However, tryptophan and many fatty acids including palmitic acid, oleic acid, linoleic acid, arachidonic acid, eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) were present at higher levels in control livers. CKD rat livers were described by increases in equol sulfate and taurocholic acid compared to control livers along with a similar AAA derivative profile as CKD plasma in OPLS-DA (, R2(Y) = 0.98, Q2(Y) = 0.87). In heart tissue samples run in negative ESI mode, metabolites in the control group were defined by lysophosphatidylcholines and lysophosphoethanolamines (, R2(Y) = 0.97, Q2(Y) = 0.89). Heart metabolites that were increased in CKD included AAA derivatives as well as pantothenic acid and glutathione, which also appeared in positive ESI mode (G and C). Positive ESI OPLS-DA demonstrated increased levels of L-carnitine and L-carnitine derivatives in heart tissue of control rats compared to CKD. Rat kidneys had the greatest separation in negative and positive ESI mode OPLS-DA plots, and showed the best fit and predictive ability compared to all other biological matrices (, R2(Y) = 1.0, Q2(Y) = 0.95, , R2(Y) = 0.99, Q2(Y) = 0.92). CKD kidneys had higher levels of AAA derivatives and adrenic acid. Both positive and negative ESI mode demonstrated increased levels of lysophosphatidylcholines and lysophosphoethanolamines in CKD. Control kidneys had higher levels of similar fatty acids seen in control livers including linoleic acid and arachidonic acid. DHA glychocholic acid, L-carnitine and L-carnitine derivatives were also increased in control kidneys compared to CKD.

Summary of metabolites altered in plasma, liver, heart and kidney tissue in rats with CKD compared to control and CKD+AST-120.

Common Significant Metabolites

A second order analysis in metaXCMS was used to determine common metabolic differences in plasma, liver, heart and kidney tissue between control and CKD rats. All tissues shared 8 common metabolites, which were significantly increased in CKD compared to control rats (). All 8 metabolites were gut derived uremic toxins and included the AAA metabolites: phenyl sulfate, indoxyl sulfate, p-cresyl sulfate, hippuric acid, pyrocatechol sulfate and p-cresol glucuronide as well as 4-ethylphenyl sulfate and equol 7-glucuronide (, ). These 8 metabolites were significantly decreased in CKD+AST-120 rats compared to controls. Plasma levels of these metabolites also showed a significant correlation with liver, heart and kidney levels (, Supplementary Figure S4). In positive ESI mode, L-carnitine and many of the L-carnitine derivatives found to be significantly different between control and CKD rats in pairwise comparisons were shared between heart and kidney tissue (). Metabolite pathway analysis demonstrated significant metabolic disturbances in a number of pathways involved with lipid, bile acid and amino acid metabolism (see Supplementary Table S6 for details).

Meta-analysis of significantly different metabolites between control and CKD in biological matrices.

Common significantly different metabolites found in all biological matrices.

In this study, rat plasma, liver, heart and kidney metabolic signatures were evaluated in CKD with and without the administration of AST-120. Rat plasma metabolic changes in CKD have been assessed in previous studies with the use of AST-120,. Renal metabolic profiles in negative ESI mode have also been evaluated in rats with CKD; however, this study is the first to assess tissue metabolic signatures in heart, liver and kidney including the metabolic axis between plasma and tissues. We demonstrate that the uremic condition causes a number of metabolic changes in both plasma and tissues in CKD and that many of these can be mitigated with AST-120.

The accumulation of uremic toxins in plasma and their contribution to the residual syndrome can be deleterious to patients with CKD. Non-invasive plasma sampling is crucial for evaluating the uremic condition in patients; however, to further our understanding of the metabolic interaction between plasma and tissues in CKD, tissue metabolic signatures must be directly assessed. Interestingly, less than 14% of metabolites that were significantly different between control and CKD rats in liver, heart and kidney tissues were also significantly different in plasma. Therefore, tissue metabolic signatures of CKD were weakly represented in plasma samples.

Rats with CKD had reduced kidney function that was not significantly affected by AST-120. The majority of plasma metabolites increased in CKD were reduced by administration of AST-120, suggesting that a number of these metabolites are derived from the gut or produced from the effects of gut-derived uremic toxins on in vivo metabolic processes. Although kidney function was not altered by AST-120, more than 50% of metabolites were reduced in both liver and heart tissue metabolic profiles with AST-120 administration to rats with CKD. Therefore, a number of tissue metabolic processes are affected by gut-derived metabolites.

The 8 metabolites that were significantly increased in all biological matrices were also significantly reduced by AST-120. P-cresyl sulfate, p-cresyl glucuronide, indoxyl sulfate, hippuric acid and phenyl sulfate are known gut-derived uremic toxins. Gut precursors of p-cresyl sulfate, indoxyl sulfate, p-cresyl glucuronide and phenyl sulfate undergo sulfation or glucuronidation, which occurs through phase II drug metabolizing enzymes,. The major site of phase II drug metabolism is in the liver and therefore, hepatic metabolism of gut precursors plays a major role in the accumulation of these toxins. Gut-derived toxins can also affect hepatic drug metabolizing enzymes. Cytochrome P450 3A4 metabolizes approximately 40% of drugs on the market and indoxyl sulfate levels have been correlated with decreased levels of the CYP3A4 endogenous metabolite, 4β-hydroxycholesterol,. Once in the blood stream, these uremic toxins are also highly protein bound and both p-cresyl sulfate and indoxyl sulfate renal clearance is mediated by organic anion transporters (OAT) 1 and 3 in the kidney,. Only the free fraction of these uremic toxins can exhibit biological effects or undergo transport into cells. In OAT1 knockout mice, indoxyl sulfate, p-cresyl sulfate, phenyl sulfate and pantothenic acid accumulate in plasma. Therefore, increased levels of these toxins in both plasma and tissues are likely a result of decreased OAT1 function in the kidneys of rats with CKD.

In this study, we demonstrate for the first time that both p-cresyl sulfate and indoxyl sulfate accumulate in heart tissue and can be mitigated by AST-120 in vivo. P-cresyl sulfate is a gut-derived bacterial metabolite produced from tyrosine. In a prospective clinical study, free p-cresyl sulfate levels predicted cardiovascular mortality and all-cause mortality in dialysis patients. The effect of p-cresyl sulfate on cardiac toxicity and dysfunction were recently published in rats with CKD suggesting NADPH oxidase activation and ROS as possible mechanisms in p-cresyl sulfate mediated cardiac apoptosis. Indoxyl sulfate is formed from tryptophan and has been shown to cause pro-fibrotic, pro-inflammatory and pro-hypertrophic affects through p38 MAPK, p42/44 MAPK, and NF-κB pathways in both rat cardiac myocytes and fibroblasts in vitro. In patients, indoxyl sulfate has been associated with the first heart failure event, and removal of indoxyl sulfate by AST-120 has been shown to improve cardiac function,. Prevention of cardiovascular disease in CKD is essential as these patients are more likely to die of cardiovascular disease than progress to requiring dialysis. P-cresyl sulfate and indoxyl sulfate are increased 361-fold and 438-fold, respectively, in rat CKD heart tissue compared to control (). Therefore, accumulation of these uremic toxins in the heart may provide evidence for potential mechanisms of cardiovascular disease development in CKD.

Pyrocatechol sulfate is a dietary phenolic metabolite that has not been previously described as a uremic toxin. Equol 7-glucoronide and 4-ethyl phenyl sulfate are also produced from gut bacteria but are specifically derived from soy protein. Standard rat chow includes soy as the main source of protein and the diet used in this study contained daidzein and genistein aglycone equivalents between 350 and 650 mg/kg, which are metabolized by gut bacteria into precursors of equol 7-glucuronide and 4-ethylphenyl sulfate, respectively,. These uremic toxins have not been identified in patients with CKD. Intestinal bacteria produce equol in 20–30% of Western populations and 40–60% of Asian populations,. Consequently, equol may only accumulate in patients who both consume soy protein and have the gut bacteria to convert daidzein to equol. Although all patients with CKD will have high AAA derived uremic toxins levels, accumulation of equol 7-glucuronide and 4-ethyl phenyl sulfate demonstrate the importance of diet on the production of uremic toxins in CKD. The 8 common metabolites identified had the highest variable importance in projection (VIP) values defining CKD in plasma and tissue S-plots ( and ). Therefore, gut-derived uremic toxins provide the most definitive metabolic signatures of both plasma and tissue in CKD. Tissue accumulation of these gut-derived uremic toxins in patients is unknown; however, in this study, we demonstrate that plasma levels of the 8 gut-derived metabolites were highly correlated with tissue levels (). This suggests that the non-invasive measurement of these 8 gut-derived uremic toxins in plasma may be used as a surrogate for tissue levels in CKD.

In control animals, higher levels of tryptophan were observed, which has also been demonstrated previously. Indoxyl sulfate is formed from tryptophan and elevated levels of indoxyl sulfate have been shown to reduce tryptophan binding to serum albumin resulting in increased levels of free tryptophan. Increasing free tryptophan levels subsequently increases tryptophan availability for uptake into tissues for metabolism, thereby reducing plasma levels. This effect is mitigated in patients with the administration of AST-120. In liver tissue, free fatty acids were decreased in animals with CKD. These included palmitic acid, oleic acid, linoleic acid, arachidonic acid, EPA and DHA ( and Supplementary Table S5). All of these fatty acids are components in the rat diet. Metabolite pathway analysis of these fatty acids demonstrates their involvement in a number of pathways; however, for the majority of fatty acid pathways, only one metabolite was matched. Palmitic acid was matched in the fatty acid biosynthesis, fatty acid metabolism and fatty acid elongation in the mitochondria pathways. AST-120 recovered hepatic levels of free palmitic acid suggesting that gut-derived uremic toxins may be involved in perturbing these pathways. Indeed, a previous study identified increased palmitic acid levels in feces of rats with adenine-induced CKD. Therefore, future targeted studies are necessary to determine the contribution of reduced dietary absorption and metabolic pathway disturbances in the liver resulting in decreased hepatic free fatty acids levels in CKD. The depletion of both EPA and DHA in kidney tissue of rats with CKD is demonstrated in this study and others,. Although plasma levels of these fatty acids were not detected, patients on dialysis have lower plasma levels of EPA and DHA, which are unaltered in non-dialysis CKD patients. A recent clinical study supplementing CKD patients with omega-3 fatty acids showed increases in specialized pro-resolving lipid mediators, which can inhibit pro-inflammatory cytokine production and potentially decrease low grade inflammation in CKD. AST-120 recovered levels of a number of free fatty acids in animals with CKD (see Supplementary Table S5). However, the mechanisms of AST-120 ameliorating the reduction in free fatty acid levels for rats with CKD are unknown. Therefore, targeted metabolomics studies assessing compounds in these metabolic pathways are necessary to understand the mechanisms resulting in metabolic disturbances in CKD.

Fatty acids are an important energy source for cardiac tissue, which can be utilized via β-oxidation. It is estimated that 50–70% of ATP generated in the adult heart is a result of β-oxidation. Although saturated fatty acids such as palmitic acid undergo β-oxidation, the majority of β-oxidation occurs through unsaturated fatty acids including oleic acid and arachidonic acid. L-carnitine mediated transport of long chain fatty acids into the mitochondria is the rate-limiting step in fatty acid oxidation. In this study, changes in heart fatty acid levels were not consistent; however, we demonstrate decreased levels of L-carnitine and L-carnitine derivatives in heart tissue from rats with CKD. Cardiac tissue cannot synthesize L-carnitine, which must be obtained through plasma from dietary sources or production in the liver and kidney. L-carnitine is cleared through the dialysis membrane resulting in decreased levels for dialysis patients; however, levels in CKD patients are not correlated with kidney function. In humans, 75% of L-carnitine is obtained from the diet,57. The standard rat diet used in this study was not supplemented with L-carnitine. As a result, it is possible that synthesis of L-carnitine is decreased in rats with CKD; however, this effect is not seen in patients due to substantial dietary L-carnitine consumption. AST-120 recovered levels of L-carnitine in CKD rats and therefore gut-derived uremic toxins may affect the synthesis of L-carnitine (see Supplementary Table S5).

In summary, gut-derived metabolic signatures most prominently defined plasma, liver, heart and kidney tissues in rats with CKD. The majority of these compounds were AAA metabolites. Increased levels of soy protein metabolites 4-ethyl phenyl sulfate and equol-glucuronide, which have not been described in patients with CKD, emphasize the potential for dietary effects on individual uremic toxin profiles. This study also highlights the gut-plasma-tissue metabolic axis, directly demonstrating the accumulation of gut-derived uremic toxins in multiple major organs of rats with CKD.

Animal Models

A total of 30 male Wistar rats weighing 150 g were obtained from Charles River Laboratories, Inc. (Wilmington, MA) for this study. Animal care and experimental protocols and procedures were approved by the Western University Animal Care Committee and conducted in accordance with the Canadian Council on Animal Care. Rats initially weighed 150 g and after 5 days of acclimatization were divided into two groups: control (n = 10) and CKD (n = 20). The CKD group received 0.7% adenine supplemented into rat chow ad libitum for 5 weeks to induce kidney disease. Food was weighed daily for the CKD group and pair-fed as standard chow to control animals. After 5 weeks, the CKD group was equally divided and the diet was changed to standard chow (CKD, n = 10) or standard rat chow supplemented with 8% AST-120 (CKD+AST-120, n = 10) for an additional 3 weeks. All animals were weighed daily throughout the study. Three rats from the CKD and CKD+AST-120 groups experienced significant weight loss and were euthanized prior to the end of the study. At the conclusion of the study, control (n = 10), CKD (n = 7) and CKD+AST-120 (n = 7) rats were sacrificed and blood was obtained in heparinized tubes to acquire plasma. Liver, heart and left kidneys were removed, snap-frozen in liquid nitrogen and stored at −80 °C. Plasma creatinine and urea were determined by the London Laboratory Services Group (London, ON, Canada) using standard methods.

Sample Preparation and Batch Setup

Plasma samples were thawed and protein was precipitated with the addition of ice-cold acetonitrile (3:1, acetonitrile:plasma). Liver, heart and kidney tissue (100 mg) were homogenized in 250 μL of ice-cold acetonitrile. All samples were placed on ice for 20 minutes and centrifuged at 14,000 g for 5 minutes. The supernatant from tissue samples was diluted to 80% acetonitrile. Quality control samples were generated for each matrix by creating a pooled sample. Acetonitrile used for protein precipitation contained isatin (5 μg/mL) and flurazepam (25 ng/mL) as internal standards. Samples from individual rats were transferred to vials and 1 μL was injected in triplicate for each vial. Injections were randomized to reduce error associated with instrument drift and quality control pool samples were run every 6 injections. All samples were run in a single batch for each biological matrix.

Chromatographic Separation and Mass Spectrometry

Metabolites were separated using a Waters ACQUITY UPLC HSS T3 column (1.8 μm particle size, 100 mm × 2.1 mm). The column was maintained at 45 °C in a Waters ACQUITY UPLC I-Class system (Waters, Milford, MA). The mobile phase flow was set to 0.45 ml/min and consisted of water (A) and acetonitrile (B), both containing 0.1% formic acid. The UPLC conditions were as follows: 0–2 mins, 1–60% B; 2–6 mins 60–85% B; 6–8 mins 85–99% B; 8–10 mins 99–1% B. Mass spectrometry was performed using a Waters XevoTM G2S-QTof/MS. Metabolites were measured separately in positive and negative ESI modes. Capillary voltage and cone voltage were set at 2 kV and 40 V, respectively. The source temperature was 150 °C. The desolvation gas flow was set to 1200 L/h at 600 °C and the cone gas flow was 50 L/h. The data was acquired in centroid mode using an MSE method, which allows for both MS and MS/MS fragmentation in a single run. Functions 1 and 2 of the MSE method acquired data with a 0.1 s scan time in the range of 50–1200 m/z, which was calibrated using sodium formate clusters prior to running samples. Collision energy was set to 0 V and ramped from 15–50V for functions 1 and 2, respectively. Function 3 acquired lockspray to ensure mass accuracy. Leucine-enkephalin (500 ng/mL) was used as the lockmass set at a flow rate of 10 μL/min, measured every 10 s and averaged over 3 scans.

Data Analysis

Mutlivariate Data Analysis

Mutlivariate analysis of LC-MS data was achieved by Waters Markerlynx with EZinfo 2.0 (Umetrics, Umeå, Sweden) software packages. Peak intensities were normalized to total marker intensity in Markerlynx and subsequently transferred to EZinfo. Pareto scaling was used to dampen the selection of features with the highest variance. EZinfo was used for principle component analysis (PCA) of plasma, liver, heart and kidney, as well as orthogonal partial least squares discriminant analysis (OPLS-DA) between control and CKD rats for each matrix.

Univariate Data Analysis

XCMS Online (https://xcmsonline.scripps.edu) was used for univariate data analysis. Raw data files were initially converted to mzData files and the lockmass was removed using massWolf 4.3.1 and the chemhelper package in R version 3.2.0. Function 1 MS data was uploaded to XCMS Online and the data was processed as a multi-group experiment using the default UPLC with G2S MS parameters. Briefly, the centWave method was used for feature detection (m/z tolerance = 15 ppm, min peak width = 2 s, max peak width = 25 s). The Obiwarp method was used for retention time correction (profStep = 0.5) and chromatograms were aligned using the following parameters: mzwid = 0.1, minfrac = 0.5, bw = 2. Metabolite features were statistically analyzed using one-way ANOVA with Tukey’s post-hoc test. Isotopes and adducts were annotated using CAMERA (m/z absolute error = 0.015, ppm = 5) and arranged into feature groups. Features for each biological matrix were filtered for isotopes. The feature with the maximum intensity in each feature group and present in at least 80% of pooled sample injections was chosen for statistical analysis. Highly significant metabolites were selected to assess differences between biological matrices in control, CKD, and CKD+AST-120 rats (p < 0.01, q < 0.001, maximum intensity > 1500).

Pairwise metabolite differences between control and CKD rats that were common across all biological matrices were determined using metaXCMS. Highly significant metabolites for each tissue were aligned using a 20 second retention time tolerance and a 0.01 m/z tolerance.

Metabolite Identification

Metabolites that were considered highly significant between control and CKD rats from multivariate and univariate analysis were searched in METLIN and Human Metabolome Database (HMDB) metabolomics databases. Fragmentation patterns for each metabolite were compared to putative database compound fragmentation using MassFragment®. Compound standards were purchased or synthesized, and retention time and fragmentation pattern were compared to confirm the metabolite identity. Pyrocatechol sulfate was identified by incubation of plasma with β-glucuronidase from Helix pomatia, which has ≥10,000 units/g sulfatase activity. The sulfatase product fragmentation pattern was confirmed as pyrocatechol to a purchased standard.

Metabolite Pathway Analysis

Pathway analysis of identified metabolites was performed in MetaboAnalyst 3.0,. GlobalTest was used for pathway enrichment analysis and relative-betweeness centrality was selected for pathway topology analysis.

This work was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Canadian Foundation for Innovation (CFI). The authors are grateful to Masaki Fujioka and Fuyu Nishijima from Kureha Corporation for helpful discussions and supplying AST-120 for our studies.

Author Contributions T.J.V., A.H., D.A.F., A.T. and B.L.U. contributed to the experimental design. Animal experiments and procedures were conducted by T.J.V., A.H., D.A.F. and A.T. UPLC-QTof/MS experiments were conducted by T.J.V., A.H., D.A.F., A.T., A.S.K. and B.L.U. Compound standards were synthesized by P.K. and M.A.K. Data and statistical analysis were performed by T.J.V., A.H., D.A.F., A.T., A.S.K., L.E.M., S.E.N. and B.L.U. Reagents and analytical tools were contributed by B.L.U. The manuscript was drafted by T.J.V., with additions and revisions by all authors.