Patients with End-Stage Renal Disease (ESRD) require dialysis or a kidney transplant for survival. Surprisingly, the incidence of ESRD due to diabetes has increased over the last 20 years despite improving hyperglycemia control and increased renoprotective drug use. To counteract this trend, it would be helpful to identify diabetic patients who are at risk for developing ESRD before they exhibit symptoms of renal disease.
Accordingly, there is an urgent need for biomarkers that are predictive of ESRD in diabetic patients.
The present invention is based, in part, on the discovery of biomarkers that are predictive of renal disease in patients who have diabetes. Accordingly, in certain embodiments, the invention described herein relates to a method of predicting risk of developing renal disease in a patient who has diabetes. The method of predicting risk comprises the steps of a) determining the levels of at least three metabolites selected from the group consisting of pseudouridine, C-glycosyltryptophan, myoinositol, threitol, p-cresol sulfate, 2-hydroxyisovalerate, 2-hydroxyisocaproate, glutaryl carnitine, N2, N2-dimethylguanosine, phenylacetylglutamine, arabitol, gulono-1,4-lactone, erythritol, erythronate, N4-acetylcytidine, urate, 2-hydroxyisocaproate, 2-oxoisoleucine, and 2-oxoisocaproate in a sample (e.g., a serum sample, a plasma sample) taken from the patient; b) comparing the levels of the metabolites in the sample from the patient to control levels of the metabolites; and c) predicting that the patient is at risk for developing renal disease when the levels of the metabolites in the sample from the patient are significantly higher than the control levels of the metabolites.
In other embodiments, the invention relates to a method of identifying a patient who has diabetes as being in need of a therapy to prevent or delay the onset of a renal disease. The method of identifying a patient comprises the steps of a) determining the levels of at least three metabolites selected from the group consisting of pseudouridine, C-glycosyltryptophan, myoinositol, threitol, p-cresol sulfate, 2-hydroxyisovalerate, 2-hydroxyisocaproate, glutaryl carnitine, N2, N2-dimethylguanosine, phenylacetylglutamine, arabitol, gulono-1,4-lactone, erythritol, erythronate, N4-acetylcytidine, urate, 2-hydroxyisocaproate, 2-oxoisoleucine, and 2-oxoisocaproate in a sample (e.g., serum sample, plasma sample) taken from the patient; b) comparing the levels of the metabolites in the sample from the patient to control levels of the metabolites; and c) implementing a therapy to prevent or delay the onset of a renal disease in the patient when the levels of the metabolites in the sample from the patient are significantly higher than the control levels of the metabolites.
The methods described herein are useful for identifying diabetic patients who are at an increased risk for developing renal disease as early as 5-10 years prior to the occurrence of clinical symptoms of the disease.
The incidence of End-Stage Renal Disease (ESRD) due to type 2 diabetes (T2D) increased over the last 20 years despite improving hyperglycemia control and increased renoprotective drugs use. (1) Clearly, a better understanding of the determinants responsible for progression to ESRD in T2D is urgently needed if this “epidemic” is to be contained.
Recently developed platforms for global metabolomic profiling are capable of examining hundreds of metabolites, so they are excellent tools to study complex metabolic alterations associated with progression of diabetic nephropathy.(2, 3) Reliable metabolomic data can be obtained with liquid or gas chromatography coupled with mass spectrometry (LC/GC-MS) or NMR spectroscopy. Among those, MS-based platforms are the most sensitive.(2, 4-6)
One of the hallmarks of progression to ESRD is plasma accumulation of certain metabolites, the so-called uremic solutes.(7-10) However, it is becoming apparent that increase in the levels of uremic solutes in blood may be more than a simple reflection of impaired kidney function.(11-13) The kidney is a key organ involved in the handling of major biochemical classes of metabolites. Kidney function includes filtration of metabolites via glomeruli, followed by their tubular secretion/reabsorption and synthesis/degradation in various components of the renal parenchyma. At present it is unclear whether elevated levels of uremic solutes precede or follow renal impairment. For example, elevated plasma concentration of uremic solutes may contribute to glomerular as well as tubular damage in diabetic nephropathy, and damage to those two components have been demonstrated in early nephropathy.(14, 15) Various alterations of certain biochemical classes of metabolites (amino acids, in particular) have been also reported in the associations with insulin resistance, type 2 diabetes or chronic kidney injury per se.(16-19)
To date, few metabolomic studies focusing on diabetic nephropathy have been performed in experimental models (20, 21) or in humans.(22-25) Nevertheless, the comparisons were either cross-sectional or focused on albuminuria progression rather than on the kidney failure, the ultimate outcome of the diabetic nephropathy.(22-25)
A description of example embodiments of the invention follows.
In certain embodiments, the present invention relates to a method of predicting risk of developing renal disease in a patient who has diabetes, comprising the steps of a) determining the levels of at least four metabolites selected from the group consisting of pseudouridine, C-glycosyltryptophan, myoinositol, threitol, p-cresol sulfate, 2-hydroxyisovalerate, 2-hydroxyisocaproate, glutaryl carnitine and N2, N2-dimethylguanosine in a sample taken from the patient; b) comparing the levels of the metabolites in the sample from the patient to control levels of the metabolites; and c) predicting that the patient is at risk for developing renal disease when the levels of the metabolites in the sample from the patient are significantly higher than the control levels of the metabolites.
As used herein, “patient” refers to a mammal (e.g., human, horse, cow, dog, cat). Preferably, the patient is a human.
In some embodiments, the patient has type 2 diabetes. In other embodiments, the patient has type 1 diabetes.
The methods disclosed herein are useful for predicting the risk of developing renal disease in diabetic patients before the onset of symptoms of a renal disease. Accordingly, in some embodiments, the patient does not have symptoms of a renal disease. In a particular embodiment, the patient has normal renal function. In another embodiment, the patient has mildly impaired renal function.
In one embodiment, the methods disclosed herein are useful for predicting the risk of diabetic nephropathy in a patient who has diabetes. In a further embodiment, the methods disclosed herein are useful for predicting the risk of ESRD in a diabetic patient.
The methods disclosed herein comprise the step of determining the levels of at least three metabolites selected from the group consisting of pseudouridine, C-glycosyltryptophan, myoinositol, threitol, p-cresol sulfate, 2-hydroxyisovalerate, 2-hydroxyisocaproate, glutaryl carnitine, N2, N2-dimethylguanosine, phenylacetylglutamine, arabitol, gulono-1,4-lactone, erythritol, erythronate, N4-acetylcytidine, urate, 2-hydroxyisocaproate, 2-oxoisoleucine, and 2-oxoisocaproate in a sample taken from the patient.
The sample that is taken from the patient can be any suitable bodily fluid sample (e.g., blood, plasma, serum, spinal fluid, lymph fluid, urine, amniotic fluid). In a particular embodiment, the sample is a plasma sample. In another embodiment, the sample is a serum sample.
The term “at least three metabolites” encompasses any combination of three or more (e.g., four, five, six, seven, eight, nine, ten, eleven, twelve, thirteen, fourteen, fifteen, sixteen, seventeen, eighteen or nineteen) metabolites selected from the group consisting of pseudouridine, C-glycosyltryptophan, myoinositol, threitol, p-cresol sulfate, 2-hydroxyisovalerate, 2-hydroxyisocaproate, glutaryl carnitine, N2, N2-dimethylguanosine, phenylacetylglutamine, arabitol, gulono-1,4-lactone, erythritol, erythronate, N4-acetylcytidine, urate, 2-hydroxyisocaproate, 2-oxoisoleucine, and 2-oxoisocaproate. Thus, any combination of three or more of these metabolites can be used to predict the risk of developing a renal disease. Preferably, the “at least three metabolites” include pseudouridine and/or C-glycosyltryptophan. More preferably, the “at least three metabolites” includes a nucleotide derivative (e.g., pseudouridine, N2, N2-dimethylguanosine), an amino acid derivative (e.g., C-glycosyltryptophan), a polyol (e.g., myoinositol, threitol), a phenyl compound (e.g., p-cresol sulfate), a branched amino acid derivative (e.g., 2-hydroxyisovalerate, 2-hydroxyisocaproate) and a branched chain acylcarnitine (e.g., 2-hydroxyisocaproate, glutaryl carnitin). In a particular embodiment, the levels of all metabolites in the group consisting of pseudouridine, C-glycosyltryptophan, myoinositol, threitol, p-cresol sulfate, 2-hydroxyisovalerate, 2-hydroxyisocaproate, glutaryl carnitine and N2, N2-dimethylguanosine are determined.
Suitable techniques and reagents for detecting levels of metabolites in a sample from a patient are known in the art. For example, levels of metabolites in a sample can be determined using mass spectrometry (MS) or NMR spectroscopy.
In a particular embodiment, levels of metabolites in a sample are determined using mass spectrometry (e.g., ESI MS, MALDI-TOF MS, tandem MS (MS/MS)). In general, mass spectrometry involves ionizing a sample containing one or more molecules of interest, and then m/z separating and detecting the resultant ions (or product ions derived therefrom) in a mass analyzer, such as (without limitation) a quadrupole mass filter, quadrupole ion trap, time-of-flight analyzer, FT/ICR analyzer or Orbitrap, to generate a mass spectrum representing the abundances of detected ions at different values of m/z. See, e.g., U.S. Pat. No. 6,204,500, entitled “Mass Spectrometry From Surfaces;” U.S. Pat. No. 6,107,623, entitled “Methods and Apparatus for Tandem Mass Spectrometry;” U.S. Pat. No. 6,268,144, entitled “DNA Diagnostics Based On Mass Spectrometry;” U.S. Pat. No. 6,124,137, entitled “Surface-Enhanced Photolabile Attachment And Release For Desorption And Detection Of Analytes;” Wright et al., “Protein chip surface enhanced laser desorption/ionization (SELDI) mass spectrometry: a novel protein biochip technology for detection of prostate cancer biomarkers in complex protein mixtures,” Prostate Cancer and Prostatic Diseases 2: 264-76 (1999); and Merchant and Weinberger, “Recent advancements in surface-enhanced laser desorption/ionization-time of flight-mass spectrometry,” Electrophoresis 21: 1164-67 (2000), each of which is hereby incorporated by reference in its entirety, including all tables, figures, and claims.
Levels of metabolites can also be determined using a Metabolon (Durham, N.C.) MS platform, as described herein.
In some embodiments, the sample from the patient is subjected to a liquid chromatography (LC) or gas chromatography (GC) purification step prior to mass spectrometry (LC/GC-MS). Methods of coupling liquid chromatography techniques to MS analysis are known in the art.
The methods disclosed herein further comprise the steps of comparing the levels of the metabolites in the sample from the patient to control levels of the metabolites and predicting that the patient is at risk for developing renal disease when the levels of the metabolites in the sample from the patient are significantly higher than the control levels of the metabolites.
A control level for a given metabolite can be obtained, for example, from a sample, or collection of samples, taken from diabetic patients who did not develop renal disease. Alternatively, a control level for a given metabolite can be based on a suitable reference standard. The reference standard can be a typical, normal or normalized range of levels, or a particular level, of a metabolite. The standards can comprise, for example, a zero metabolite level, the level of a metabolite in a standard cell line, or the average level of a metabolite previously obtained for a population of normal human controls. Thus, the methods disclosed herein do not require that the level of a metabolite be assessed in, or compared to, a control sample.
In accordance with the invention, a patient is predicted to be at risk for developing renal disease when the levels of the metabolites in the sample from the patient are significantly higher than the control levels of the metabolites. A statistically significant difference (e.g., an increase, a decrease) in the level of a metabolite between two samples, or between a sample and a reference standard, can be determined using an appropriate statistical test(s), several of which are known to those of skill in the art. In a particular embodiment, a t-test (e.g., a one-sample t-test, a two-sample t-test) is employed to determine whether a difference in the level of a metabolite is statistically significant. For example, a statistically significant difference in the level of a metabolite between two samples can be determined using a two-sample t-test (e.g., a two-sample Welch's t-test). A statistically significant difference in the level of expression of a gene between a sample and a reference standard can be determined using a one-sample t-test. Other useful statistical analyses for assessing differences in gene expression include a Chi-square test, Fisher's exact test, and log-rank and Wilcoxon tests.
As used herein, patient who is “at-risk” for developing renal disease has a level of risk for developing a renal disease that is higher than the level of risk for an individual represented by the relevant baseline population.
In other embodiments, the invention relates to a method of identifying a patient who has diabetes as being in need of a therapy to prevent or delay the onset of a renal disease, comprising the steps of a) determining the levels of at least four metabolites selected from the group consisting of pseudouridine, C-glycosyltryptophan, myoinositol, threitol, p-cresol sulfate, 2-hydroxyisovalerate, 2-hydroxyisocaproate, glutaryl carnitine and N2, N2-dimethylguanosine in a sample taken from the patient; b) comparing the levels of the metabolites in the sample from the patient to control levels of the metabolites; and c) implementing a therapy to prevent or delay the onset of a renal disease in the patient when the levels of the metabolites in the sample from the patient are significantly higher than the control levels of the metabolites.
As defined herein, “therapy” is the administration of a particular therapeutic or prophylactic agent to a subject (e.g., a non-human mammal, a human), which results in a desired therapeutic or prophylactic benefit to the subject (e.g., prevention or delay in the onset of a renal disease).
A suitable therapy for preventing or delaying the onset of a renal disease in a patient can be readily determined by a skilled medical professional (e.g., a physician, such as a nephrologist), taking into account various factors including, but not limited to, the patient's age, weight, medical history, and sensitivity to drugs. Exemplary therapies for preventing or delaying the onset of a renal disease include, for example, administration of drugs to treat hypertension, dietary changes, exercise, weight loss, glycemic control, proteinuria therapies, and albuminuria therapies, among others at least three.
Preferably, the therapy is implemented early enough to prevent or delay the onset of renal disease in the patient. In some embodiments, the therapy is implemented before the patient shows any clinical symptoms of renal disease. According to the invention, the methods disclosed herein can identify diabetic patients who are at risk for developing renal disease about 5-10 years, or more, before the occurrence of clinical symptoms of renal disease.
Study Groups and Methods:
Study Group
Between 1991 and 1995, a cohort with T2D was recruited into the Joslin Study of the Genetics of Type 2 Diabetes and Kidney Complications (half with normoalbuminuria and half with microalbuminuria or proteinuria). The cohort was followed-up until 2004 for the occurrence of ESRD or death unrelated to ESRD. Details of the recruitment, examination, and follow-up of this cohort were already published elsewhere.(26) Among 410 individuals, 59 developed ESRD and 84 died without ESRD as ascertained by United States Renal Data System(1), National Death Index(76) and medical records review.
For this nested case-control study, 40 out of 56 incident cases of ESRD who had sufficient stored plasma samples at baseline available were selected. A group of 40 subjects from among those who survived and were without ESRD as of the end of follow-up were selected as controls. They were grouped-matched with cases with regard to gender, age, and baseline eGFR. For comparison of cases and controls baseline clinical characteristics were used as reported previously.(26)
As a small reproducibility study, ten study subjects (balanced by caseness status) had plasma samples selected 2.2+0.8 years after baseline. All plasma samples were stored at −70° C. until analysis. The study protocol and informed consent procedures were approved by the Joslin Diabetes Center Institutional Review Board.
Global Metabolomics Profiling
All plasma samples (80 baseline and 10 from the early follow up timepoint) were subjected to global metabolomic profiling (Metabolon, Inc, Durham, N.C.). Sample extracts were prepared to recover a wide range of chemically diverse metabolites. The LC/MS portion of the platform incorporates a Waters Acquity UPLC system and a Thermo-Finnigan LTQ mass spectrometer, including an electrospray ionization (ESI) source and linear ion-trap (LIT) mass analyzer. The GC column is 5% phenyl dimethyl silicone and the temperature ramp is from 60° to 340° C. in a 17-minute period. All samples were then analyzed on a Thermo-Finnigan Trace DSQ fast-scanning single-quadrupole mass spectrometer using electron impact ionization. Biochemicals were subsequently identified by comparison to library entries of over 2400 purified standards for distribution to both the LC/GC-MS platforms for determination of their analytical characteristics.
After peak identification and quality control filtering (signal greater than 3 times background, retention index within a pre-specified platform-dependent window, ion quantification to library match within 0.4 m/z, and the MS/MS forward and reverse scores), the metabolites' relative concentrations were obtained from median-scaled day-block normalized data for each compound. Repeated samples from the same individual were run the same day. The metabolite was defined as common, if it was present in at least 80% of the individuals in the study group, and as stable over time, when Spearman correlation coefficient between two measurements taken from the same individual was ≧0.4. For detailed information regarding all metabolites detected by global profiling, please refer to Supplemental Table 1.
Uremic Solutes
The European Uremic Toxins (EUTox) Work Group, initiated in 1999, consists of 24 European Research Institutes and provides the most comprehensive encyclopedic list of systematically and critically reviewed uremic solutes/toxins.(7, 9) Metabolites measured with the global profiling were classified as uremic solutes/toxins based on the EUTox list prepared in 2003, revisited in 2012 as well as based on selected relevant other publications.(7-9, 11, 13, 32) Seventy eight uremic solutes are available in the Metabolon library. For detailed information of the detectable uremic solutes in this study, please see Supplemental Table 2.
Targeted Quantitative Measurements of Metabolites
Quantitative measurements of phenyl and indole compounds (p-cresol and indoxyl sulfate, phenylacetylglutamine and hippurate) were performed as a collaborative effort with Dr. T. W. Meyer's laboratory by stable isotope dilution LC/MS/MS using p-cresol-d8-sulfate (synthesized from p-cresol d831 (Cambridge Isotopes)), indoxyl-2,4,5,6,7-d5-sulfate (Isosciences), and N-benzoyl-d5-glycine and Nα-(phenyl-d5-acetyl)-L-glutamine (both C/D/N Isotopes) as internal standards. LC was performed on a Kinetex, C18 column maintained at 30° C. MS was performed on an Agilent 6430 Triple Quadrupole mass spectrometer with electrospray ionization in the negative mode (Agilent Technologies). Solute concentrations were calculated using manufacturer's software (MassHunter Quant). Ion transitions used for quantitation were m/z 187/107 for p-cresol sulfate, m/z 212/80 for indoxyl sulfate, m/z 178.1/134.2 for hippurate, and m/z 263.2/145.1 for phenylacetylglutamine with corresponding transitions for the deuterated internal standards. Spiked recoveries for all the metabolites were between 93 to 108%.
Quantitative measurements of modified nucleotide derivatives and amino acids were performed as a collaborative effort with Dr. Pennathur's laboratory. Isotopically labeled 13C6 tyrosine and 13C6 phenylalanine (Cambridge Isotope Laboratories, Andover, Mass., USA), 13C1 15N2 uracil (Movarek Biochemicals and Radiochemicals Brea, Calif., USA) were used to quantify the respective authentic compounds. Authentic pseudouridine (Santa Cruz Biotechnology, Dallas, Tex., USA) and N2N2-dimethylguanosine (Biolog Life sciences Institute, Bremen, Germany) were utilized as standards to construct calibration curves. Plasma samples were subjected to protein precipitation with 1:1 v/v 100% acetonitrile after addition of isotopically labeled internal standards. The supernatant containing the metabolites was subjected to liquid chromatography-electrospray ionization tandem mass spectrometry (LC-ESI-MS/MS) and metabolites were quantified using multiple reaction monitoring (MRM) in the MS/MS positive ion acquisition mode. For LC/ESI/tandem MS experiments, we utilized an Agilent 6410 triple Quadrupole MS system equipped with an Agilent 1200 LC system and an ESI source. Luna C-18 column was used for the LC separation. The mobile phase is 0.1% Formic acid (Solvent A) and acetonitrile with 0.1% formic acid (Solvent B). The following transitions were monitored for each of the analyte: a) m/z 113 to m/z 70 and m/z 116 to m/z 71 for uracil and 13C115N2 uracil; b). m/z 166/120 and m/z 172/126 for phenylalanine and 13C6 phenylalanine; c) m/z 182/136 and m/z 188/142 for tyrosine and 13C6 tyrosine, and d) m/z 245/209 and m/z 312/180 for pseudouridine and N2,N2-dimethylguanosine respectively. For analytes with available corresponding isotopic standard the quantification was computed by the ratio of the peak area of the analyte compared with the known amount of the isotopic standard. For pseudouridine and N2,N2-dimethylguanosine calibration curves were created with the authentic nucleoside spiked in the biological matrix. Extraction efficiency and quality control were monitored with the peak area of the labeled stable isotopes in control plasma samples spaced throughout the run to monitor intra and inter run efficiency. Amino acids. Quantification of amino acid levels were performed using the “EZ:faast” kit for free amino acid analysis from Phenomenex (Torrance, Calif., USA) with Norvaline added as an internal standard as described previously using an Agilent 6890 GC-MS system.(79) For the results of amino acids measured on two platforms, please see the Supplemental Table 3.
Myo inositol measurements were performed as a collaborative effort with Dr G. Berry's laboratory (Boston Childrens Hospital, Boston, Mass.). Myo-inositol was purchased from Sigma-Aldrich (USA). The internal standard (IS), [2H6]-myo-inositol, was purchased from Cambridge Isotope Laboratories, Inc. (MA, USA). The separation of samples and standards was performed with Shimadzu Prominence high performance liquid chromatography (HPLC) unit and a SUPELCOGEL Pb (300×7.8 mm; 5 μm; Supelco UK) column, with the column temperature set at 85° C. The HPLC was coupled to a 5500 QTRAP hybrid dual quadrupole ion trap mass spectrometer (AB Sciex, Ontario, Canada) operating in negative-ion mode. The MRM transitions monitored were 179/87 for myo-inositol and 185.1/89 for [2H6]-myo-inositol. Samples were run in duplicate and measurements with coefficient variation (CV)<10% were considered for analysis.
Urate was measured in the clinical laboratory of the Joslin Diabetes Center, Boston, Mass. using an enzymatic colorimetric assay kit and read on the Roche Modular P Chemistry analyzer (Roche Diagnostics).(80) The inter-assay CV range is 1.9% in this laboratory.
Data Analysis
Differences in clinical characteristics between the two study groups were tested by analysis of variance for continuous variables and chi-square test for categorical variables. Preliminary data cleaning included investigations of detectability, batch effects and outliers at the metabolite- and individual-levels (heatmap, principal component procedures; data not shown).(29, 77) Volcano plots were generated for common metabolites based on the fold difference between the outcome groups and the p-value obtained in a general linear model. Adjustment for multiple comparison was performed with a positive false discovery rate (pFDR) q value <0.05 for significance.(78) The effect of a metabolite was estimated using logistic regression. After transformation of the metabolite concentration to normally distributed ranks, its effect on risk of progression was expressed as the odds ratio for a one standard deviation difference.(26) Clinical covariates and metabolites measured quantitatively were transformed to their (base 10) logarithms for the logistic analysis. Non-common metabolites were analyzed as categorical variables by chi-square, but this analysis did not result in identifying additional significant metabolites (data not shown).
Correlations between the continuous variables were examined with Spearman rank correlation. Data reduction was carried out with hierarchical cluster analysis using the Ward method based on values transformed to normally distributed ranks. Top metabolites from the multivariate volcano plots analysis were included. More than a half of the detected metabolites lacked pathway identifiers, preventing us from a comprehensive canonical pathway analysis. Data analysis was performed with SAS 9.3 and JMP Pro 9.0.0 softwares (Cary, N.C.).
Results
Study Groups and their Characteristics:
A cohort with T2D patients attending the Joslin Clinic was recruited into the Joslin Study of the Genetics of Kidney Complications. Of the 509 patients examined between 1992 and 1996. 410 were followed until the end 2004. During 8-12 years of follow-up 59 (14.4%) patients developed ESRD, 84 (20%) died without progressing to ESRD and 267 (65.1%) remained alive without progressing to ESRD. Details of the follow-up study were already published.(26)
For the present nested case-control study, 40 patients who developed ESRD (cases of progressors to ESRD) were selected and matched them with 40 patients who were alive as of 2004 without ESRD (controls for non-progressors). Of the 80 patients, 75 identified themselves as Caucasians of European origin. Baseline characteristics of the two selected study groups are summarized in Table 1. The groups were very similar with regard to most clinical characteristics. Progressors, however, had higher urinary albumin excretion and slightly lower eGFR. Despite the differences noted in median AER and mean eGFR, there was substantial overlap of the distributions in the two study groups. At baseline the majority of progressors and non-progressors were in CKD stage 2. CKD stage 3 was present in 7% of controls and 22% of cases, respectively. Overall the distribution of CKD stages was not statistically different between the study groups. 87% of non-progressors had annual eGFR decrease less than 3.5 ml/min/1.73 m2. The median (25th, 75th percentile) decrease was −1.95 (−3.2, −0.8) ml/min/1.73 m2 and the slope was determined based on the serial creatinine measurements over 7.6 (6.5-12.4) years. The study groups did not also differ regarding baseline plasma levels of parathormone.
aBonferroni corrected.
bUntil ESRD development.
Results of Global Metabolomic Analysis:
Baseline plasma samples from study subjects were run on the Metabolon platform against a library containing mass spectra for 2400 chemically identified metabolites. A total of 445 named metabolites were detected. Of these, 183 belonged to drug-related metabolites or had low detectability and were excluded from further analysis. The remaining 262 were detected in at least 80% of the study subjects and were designated as “common.” The assay was repeated in the follow-up plasma sample taken one to three years later for a subgroup of study patients to estimate intra-individual variation for each metabolite. The results by biochemical classes are summarized in
Fold differences between the mean plasma concentration in progressors to ESRD and non-progressors were analyzed for the 262 common metabolites (Table 2). In the analysis of fold difference, 49 (19%) were significant after adjusting for multiple comparisons (q value <0.05). Among the subset of 119 metabolites stable over time, 28 (23%) differed significantly from 1.0. These results are summarized in Table 2 according to biochemical class and whether the metabolite has been identified as a uremic solute. The proportion of metabolites associated with ESRD progression was similar regardless of whether the metabolite was stable over time. Nevertheless, further analysis focused on the stable ones because their involvement in a prolonged process such as progression to ESRD is more plausible. Among the biochemical classes, the proportion associated with ESRD was high for amino acids and their derivatives, carbohydrates, and modified nucleotides (40%, 42%, and 57%, respectively), intermediate (16%) for other metabolites, and low (4%) for lipids (see
aValues of fold difference were significantly different from 1.0 at a q-value o0.05. See also FIGS. 2a and b.
bDetectable in plasma of X80% of patients.
cCommon and stable over time, that is, Spearman's correlation coefficient X0.4 between measurements taken 1-2 years apart from the same individual. Out of 29 stable metabolites associated with ESRD, five were not examined further.
Progression to ESRD According to Plasma Concentration of Uremic Solutes:
The fold difference between progressors and non-progressors and its q value for each of these metabolites are plotted in
Among the amino acid-derived uremic solutes associated with progression to ESRD were two, p-cresol sulfate and phenylacetylglutamine, produced by the gut microbiome. Their effects on the risk of progression to ESRD were strong. For example, the odds ratio for progression to ESRD for a one standard deviation increase in plasma p-cresol sulfate concentration was 2.3 (95% CI; 1.3, 3.9) in univariable analysis. The effect of phenylacetylglutamine was similar, but slightly less than that of p-cresol sulfate. Of the six polyol derived uremic solutes significantly associated with the risk of progression to ESRD, myo-inositol was the one most strongly associated, odds ratio: 3.2 (95% CI; 1.7, 5.9). Of the four nucleotide-derived uremic solutes significantly associated with the risk of progression to ESRD, three are derived from degradation of RNA and the fourth (urate or uric acid) is derived from degradation of DNA. The strongest association with progression to ESRD was for pseudouridine, odds ratio 7.8 (95% CI; 3.1, 19).
Progression to ESRD According to Plasma Concentration of Amino Acids and their Derivatives:
Thirty-nine metabolites representing amino acids or their derivatives were common and stable over time and 14 of them (29%) were associated with risk of progression to ESRD. The fold difference between progressors and non-progressors and its q-value for each of these metabolites are plotted in
The effects of the remaining 12 amino acids on risk of ESRD were estimated with logistic regression
A few amino acid derivatives were positively associated with progression to ESRD. C-glycosyltryptophan was elevated and the most significantly different between progressors and non-progressors among the metabolites shown in
Among other metabolites associated with ESRD risk, there were two lipids, dihomo-linolenate (20:3n3) and docosapentaenoate (n3 DPA; 22:5n3). Both are ω6 fatty acids and were negatively associated with risk of progression to ESRD. Gluconate and two metabolites involved in vitamin B metabolism: pantothenate and N1-methyl-2-pyridone-5-carboxamide were positively associated with progression to ESRD. Data for these five metabolites were not shown.
Reduction of Redundant Data
Among metabolites significantly associated with progression to ESRD, some may reflect on shared underlying biology. To evaluate the potential redundancy among these metabolites, a Spearman rank correlation matrix was created. Approximately two-thirds of the metabolites were strongly correlated (data not shown).
A cluster analysis revealed clusters (see
Targeted Quantitative Metabolites Measurements
To validate these findings, targeted quantitative measurements of nine common and stable metabolites over time identified on the Metabolon global biochemical profiling platform was performed. Five out of nine metabolites were the ones strongly associated with progression to ESRD. Among the remaining four, two metabolites, indoxyl sulfate and hippurate were uremic solutes reported in the literature, but not significantly associated with progression in the global profiling study. The remaining two measured were allantoin and uracil, linked in the metabolic pathways related to urate and pseudouridine, respectively. Allantoin was not sufficiently detectable by the quantitative methods employed in this study (data not shown). For each of the other metabolites, association with progression to ESRD was consistent with the result obtained via global profiling. Also the association remained after adjustment for AER, eGFR and HbA1c (Table 3). Odds ratio for ESRD progression was 1.7 (95% CI; 1.0, 2.8) per one standard deviation of the logarithmically transformed metabolite, p-cresol sulfate and phenylacetylglutamine and higher than 2.5 for every other significant metabolite (pseudouridine, p-cresol sulfate, myoinositol and urate). It was confirmed that hippurate and indoxyl sulfate were not significantly associated with ESRD in this study. Logistic regression of the pseudouridine/uracil ratio did not improve the discrimination ability of the model (data not shown).
aTo convert urate (uric acid) mmol/l concentrations to mg/dl, divide the values by 59.5. Concentration of the uric acid 275 mmol/l is equal to 4.6 mg/dl.
This study was performed in subjects with T2D, the majority of whom had normal renal function at baseline. Half progressed to ESRD and half did not during a decade of follow-up. In baseline plasma, progressors could be distinguished from non-progressors by high concentrations of metabolites referred to as uremic solutes and low concentrations of certain amino acids and their derivatives. This is the first demonstration that abnormal plasma concentrations of certain metabolites are associated with risk of progression to ESRD at a very early stage of diabetic nephropathy.
Uremic solutes, as catalogued by EUTox group,(7, 9) comprise compounds of different biochemical classes: amino acid derivatives, certain alcohol/polyols and modified nucleosides among them. Of the 18 known uremic solutes that were detected as common and stable metabolites with the Metabolon platform, 12 were elevated in subjects who progressed to ESRD. On the basis of results of multiple cross-sectional studies (7-10) the immediate interpretation of these findings might be that the increased concentration of uremic solutes was due to significant impairment of renal function in subjects who progressed to ESRD during 8-12 years of follow-up. However, the majority of progressors at baseline were in CKD stages 1 and 2 and had only slightly lower baseline eGFR than non-progressors (75±19 ml/min vs. 87±23 ml/min). In multivariable analyses, the odds ratios of progression to ESRD for specific uremic solutes were only minimally affected by adjustment for eGFR, AER and HbA1c.
Phenyl compounds, such as p-cresol sulfate and phenylacetylglutamine are the most extensively studied solutes known to increase in the uremic state.(8, 11, 13) These solutes can be toxic to endothelial cells and can contribute to increased risk of cardiovascular complications in patients with renal impairment.(30, 31) In humans these metabolites are exogenous and are produced by intestinal bacterial flora before they are absorbed into plasma and excreted through the kidney.(8, 11, 13, 32) Evidence confirming the microbiome as a source for these solutes was recently provided in a study of ESRD subjects with and without a colon.(11) In this study, high plasma concentrations of these solutes were associated with progression to ESRD.
Three additional solute derivatives of amino acids that are synthesized in the gut, phenol sulfate, indoleacetate and 3-indoxyl sulfate, were elevated in plasma of progressors compared to non-progressors, however, the differences did not reach statistical significance (see supplemental Table 2). All these findings support the notion that the gut microbiome (11, 13, 32) might control plasma levels of amino acid-derived uremic solutes, and their high levels increase the risk of progression to ESRD in subjects with T2D. The increase in circulating levels of these solutes also may be attributed to dietary factors.(34)
Elevated plasma concentrations of myo-inositol and other polyols were strongly associated with progression to ESRD in this study. These metabolites accumulate in plasma in the uremia state and during acute kidney injury.(9, 35, 36) The kidney is the most important organ for whole body metabolism of myo-inositol, as both the synthetic and degradative enzymes, L-myo-inositol-1-phosphate synthetase and myo-inositol oxygenase (MIOX), are highly expressed in the renal parenchyma.(37, 38) Myo-inositol can be obtained also from dietary sources. Regulation of its plasma levels is complex, involving glomerular filtration, reabsorption in proximal tubules in competition with glucose transport, apparent conversion to chiro-inositol in tissues and, finally, catabolic degradation of myo-inositol by MIOX, a protein that is upregulated by hyperglycemia.(36, 38-44) The high concentration gradients of myo-inositol in certain tissues such as kidney and the vascular endothelium are maintained by the osmoregulatory sodium/myo-inositol cotransporter 1 (SMIT1)(39, 40) Glucose/myo-inositol imbalance was demonstrated to induce proliferative and pro-fibrotic response in the proximal tubules in vitro and to alter the immune responses in leukocytes.(45, 46) Increased urinary excretion of myo-inositol and inositol imbalance in muscle tissue (high myo/chiro-inositol ratios) were reported in subjects with T2D.(42-44) In this study chiroinositol levels were hardly detectable (data not shown).
Plasma concentrations of several nucleotide derivatives that are considered to be uremic solutes were also strongly associated with progression to ESRD in this study. Among these derivatives, elevated concentration of pseudouridine in plasma was the strongest and most statistically significant predictor of progression. Pseudouridine belongs to the group of modified nucleosides that are regarded as indicators of whole-body RNA turnover.(47) These metabolites are increased in patients with malignancies,(48) and with uremia.(8, 49-51) Pseudouridine is synthesized from uracil(52, 53) and constitutes an end product, as it is not catabolized in humans.(54) Trace study in humans with radiolabeled pseudouridine showed that the kidney handling includes glomerular filtration and also tubular reabsorption.(55) High plasma levels of two other pyrimidine derivatives that are closely correlated with plasma levels of pseudouridine also increased the risk of ESRD progression (N2, N2-dimethylguanosine, and N4-acetylcytidine). These three nucleotides increase in plasma as urinary excretion decreases and accumulate significantly in the uremic state.(51)
Urate (uric acid) is a metabolite of purine metabolism. In this study, increased plasma concentrations were associated with progression to ESRD. Urate is another compound known to accumulate in the uremic state.(8, 9) Its increase, however, is disproportionally small due to compensatory mechanisms including increased enteric excretion, decreased production(56) and possibly altered tubular handling.(57) In a recent study, elevated plasma level of urate was a strong predictor of early renal function decline during follow-up of a large cohort of subjects with T1D.(58)
Kidney protein turnover, as compared with muscle and splanchnic turnover, is characterized by the highest rates of protein synthesis and amino acid oxidation, mainly in the tubulointerstitium.(17, 66) Depletion of the circulating pools of branched chain amino acids and tryptophan are known phenomena accompanying advanced chronic kidney disease.(10, 16, 17) On the other hand, increase in those amino acids was shown to predict development of T2D. (19) Interestingly, this study revealed that not only branched and aromatic, but also all other essential amino acids and their derivatives were lower in subjects who progressed to ESRD than in those who were non-progressors. In an experimental model of acute kidney injury, one of the strongest metabolic responses to nephrotoxins was massive excretion of all essential amino acids.(35, 67) It needs to be determined whether impaired tubular reabsorption contributes to the decreased levels of amino acids in early diabetic nephropathy. Amino acids that were decreased in subjects at risk of ESRD in this study all have a neutral charge. They are handled by the B0AT1 cotransporter responsible for the luminal influx, the heterodimeric exchanger and the facilitated diffusion transporter TAT1 in charge of the basolateral efflux and possibly by other transporters.(68, 69)
Among amino-acid derivatives, C-glycosyltryptophan showed a different pattern of association than the others. Its plasma concentration was the highest in progressors when compared with non-progressors. After pseudouridine it had the second highest fold difference between the study groups. Plasma concentrations of both were very highly correlated. Interestingly, they both carry a C-glycosylation linkage, a rare type of posttranslational protein modification.(70) Increased expression of the proteins containing certain forms of C-glycosylated tryptophan in the aortic vessels have been reported in the diabetic rats.(71) C-glycosylated tryptophan also correlates with eGFR.(72)
Analysis of acylcarnitines revealed that the increased concentrations were independently associated with risk of progression to ESRD. Acylcarnitines are filtered through the kidney and about 75% are excreted into urine.(73) Serum acylcarnitines deriving from lipid and amino acids are inversely correlated with GFR in individuals with normal as well as with impaired renal function.(18, 74, 75) Acylcarnitines transport is regulated by organic carnitine transporters in the kidney.(61) In this study, amino acid-deriving (but not lipid-deriving) acylcarnitines were increased in the subjects at risk. They are generated via beta-oxidation of the branched chain amino acids. Those amino acids and their intermediate keto acid derivatives were also depleted in this study (2-oxoisoleucine, 2-oxoisocaproate).
Supplemental Tables 1A-1D.
Analytical and intraindividual performance of the 445 metabolites detected by global metabolomic profiling in the study subjects stratified by the biochemical classes (1A—amino acids, 1B—carbohydrates, 1C—lipids, 1D—metabolites that belong to other than the major three classes. Metabolites are stratified by their detectability and subsequently presented in the alphabetic order. Drug related metabolites are not displayed.
The relevant teachings of all patents, published applications and references cited herein are incorporated by reference in their entirety.
While this invention has been particularly shown and described with references to example embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the scope of the invention encompassed by the appended claims.
This application claims the benefit of U.S. Provisional Application No. 62/103,709, filed on Jan. 15, 2015. The entire teachings of the above application are incorporated herein by reference.
This invention was made with government support under Grant Nos. P30DK036836, DK41526, and DK67638 from National Institutes of Health (NIH). The government has certain rights in the invention.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US16/13661 | 1/15/2016 | WO | 00 |
Number | Date | Country | |
---|---|---|---|
62103709 | Jan 2015 | US |