- Research article
- Open Access
- Open Peer Review
Plasma metabolites and lipids associate with kidney function and kidney volume in hypertensive ADPKD patients early in the disease course
BMC Nephrologyvolume 20, Article number: 66 (2019)
Autosomal dominant polycystic kidney disease (ADPKD) is the most common hereditary kidney disease and is characterized by gradual cyst growth and expansion, increase in kidney volume with an ultimate decline in kidney function leading to end stage renal disease (ESRD). Given the decades long period of stable kidney function while cyst growth occurs, it is important to identify those patients who will progress to ESRD. Recent data from our and other laboratories have demonstrated that metabolic reprogramming may play a key role in cystic epithelial proliferation resulting in cyst growth in ADPKD. Height corrected total kidney volume (ht-TKV) accurately reflects cyst burden and predicts future loss of kidney function. We hypothesize that specific plasma metabolites will correlate with eGFR and ht-TKV early in ADPKD, both predictors of disease progression, potentially indicative of early physiologic derangements of renal disease severity.
To investigate the predictive role of plasma metabolites on eGFR and/or ht-TKV, we used a non-targeted GC-TOF/MS-based metabolomics approach on hypertensive ADPKD patients in the early course of their disease. Patient data was obtained from the HALT-A randomized clinical trial at baseline including estimated glomerular filtration rate (eGFR) and measured ht-TKV. To identify individual metabolites whose intensities are significantly correlated with eGFR and ht-TKV, association analyses were performed using linear regression with each metabolite signal level as the primary predictor variable and baseline eGFR and ht-TKV as the continuous outcomes of interest, while adjusting for covariates. Significance was determined by Storey’s false discovery rate (FDR) q-values to correct for multiple testing.
Twelve metabolites significantly correlated with eGFR and two triglycerides significantly correlated with baseline ht-TKV at FDR q-value < 0.05. Specific significant metabolites, including pseudo-uridine, indole-3-lactate, uric acid, isothreonic acid, and creatinine, have been previously shown to accumulate in plasma and/or urine in both diabetic and cystic renal diseases with advanced renal insufficiency.
This study identifies metabolic derangements in early ADPKD which may be prognostic for ADPKD disease progression.
HALT Progression of Polycystic Kidney Disease (HALT PKD) Study A; Clinical www.clinicaltrials.gov identifier: NCT00283686; first posted January 30, 2006, last update posted March 19, 2015.
Autosomal dominant polycystic kidney disease (ADPKD) is the most common hereditary kidney disorder and affects 1/1000 individuals. It is characterized by gradual enlargement of numerous cysts in the kidneys over decades, and the disease process begins long before loss of estimated glomerular filtration rate (eGFR) occurs. There are at least three definite genetic causes of ADPKD. The majority of ADPKD cases (~ 75%) are caused by mutations in polycystin 1 , and second most common (~ 15%) are mutations in polycystin 2 (PKD2) . Recently a third causative gene in ADPKD and autosomal dominant polycystic liver disease (ADPLD) was identified to be GANAB, responsible for 0.3% of all ADPKD . Mutations in GANAB result in a defect in the maturation of PKD1 such that it fails to localize in the plasma membrane . PKD1 binds to PKD2  and this protein complex signals tubular morphogenesis through the formation of an ion-channel . When GANAB is mutated and PKD1 maturation is mostly blocked, then PKD1 doesn’t interact correctly with PKD2 and PKD2 fails to localize in the cilia . This leaves 5–10% of ADPKD patients with no detectable mutation after DNA sequencing of their PKD1 and PKD2 genes .
The course of ADPKD is variable depending not only on which gene is mutated  and the strength of the mutation , but also on environmental factors such as dietary sodium intake  and smoking exposure , both of which contribute to disease severity and progression. A high sodium intake could contribute to an increased rate of disease progression in ADPKD via various possible mechanisms  and not just by its effects on hypertension. Smoking has been found to be associated with established proteinuria and this is highly associated with hypertension and a more severe renal cystic phenotype . Hypertension is common and occurs early in ADPKD, prior to loss of kidney function, and is associated with more rapid disease progression and increased kidney size [11, 12]. With the recognition that ADPKD is at least in part a metabolic disease as has become evident with the discovery that glucose, histidine, glutamine, and arginine metabolic pathways are reprogrammed [13,14,15,16], there is the opportunity to investigate such pathways to develop new therapies for this disease.
With the continuing advent of new therapeutic approaches, identifying patients at risk for progression to ESRD has become paramount. The Consortium for Radiologic Imaging Studies of Polycystic Kidney Disease (CRISP) study established ht-TKV as assessed by magnetic resonance imaging as an accurate and reliable measure of renal cyst burden, which also strongly predicts those patients who will progress to advanced chronic kidney disease (CKD) . Ht-TKV is now approved by the FDA as both a prognostic and a clinical trial enrichment imaging biomarker for ADPKD . We asked whether plasma metabolomic analyses can further identify circulatory signatures which are indicative of increased ht-TKV or kidney function (eGFR) in individuals at risk for progression to ESRD early in the course of disease. Such studies will shed light on mechanisms of disease, including metabolic reprogramming. The HALT A clinical trial was designed and developed to test the hypothesis that rigorous blood pressure control (< 110/75 mmHg) or combined angiotensin converting enzyme inhibitor and angiotensin receptor blocking therapy would slow disease progression defined as the rate of increase in ht-TKV and the chronic slope of decline in eGFR. Importantly, low blood pressure control was found to have a significant protective impact on both the rate of ht-TKV increase as well as the chronic eGFR slope. For this study, we now have utilized plasma from the baseline visit from participants in the HALT-A randomized clinical trial [12, 19], after stopping all antihypertensive medications for at least two weeks, in an attempt to correlate expression of specific metabolites with simultaneously measured eGFR and ht-TKV.
Ethics approval and consent to participate
IRBs at each participating center approved the HALT-A clinical trial. These included University of Colorado Health Sciences Center, Denver (Aurora), Colorado; Emory University School of Medicine, Atlanta, Georgia; University of Kansas Medical Center, Kansas City, Kansas; Tufts Medical Center, Boston, Massachusetts; Beth Israel Deaconess, Boston, Massachusetts; Cleveland Clinic, Cleveland, Ohio; Mayo Clinic, Rochester, Minnesota. All subjects provided written consent, and parental consent was required and obtained for all participating subjects between the ages of 15 and 18 years of age at each participating site. Two Data Coordinating Centers were involved in the HALT-A clinical trial: Washington University, St. Louis Missouri and the University of Pittsburgh, Pittsburgh, Pennsylvania. IRB study approval was obtained at each Data Coordinating Center.
This study was approved by the Institutional Review Boards both at the University of California, Davis, and at the University of Chicago. The HALT Progression of Polycystic Kidney Disease (HALT PKD, NCT00283686, http:// clinicaltrials.gov) Study A is a double-blind, placebo-controlled trial consisting of 588 hypertensive participants (15 to 49 years of age, with eGFR > 60 ml/1.73 m2) . Briefly, subjects provided written consent and were enrolled, after which antihypertensive medications were withdrawn for two weeks. Patients who demonstrated elevations in blood pressure during washout were placed on labetolol therapy. At baseline, blood was obtained for centralized eGFR measures, MR imaging performed for ht-TKV measurement and plasma and urine collected and appropriately stored for future analyses. Plasma samples were handled immediately at each participating clinical center and stored within 30 min at − 80 degrees. Samples were shipped in bulk to the NIDDK repository for storage prior to this study. Subjects were randomized in a 2 × 2 factorial design to determine whether rigorous blood pressure control (< 110/75 mm HG vs. < 130/80 mmHg) or additive inhibition of the renin-angiotensin-aldosterone-system (RAAS) impacted disease progression. The primary outcomes of our study were baseline ht-TKV and eGFR measures. We analyzed baseline fasting plasma samples obtained from a subset of 277 Caucasian ADPKD patients who were known to have mutations in either the PKD1 and PKD2 genes.
Non-targeted metabolomics analysis
Samples and standards were analyzed as described . Peaks were deconvoluted and detected using Leco ChromaTOF software and were matched to FiehnLib mass spectra and retention time libraries. The Binbase software was used for post-curation and peak replacements, and the sum of intensities for all known compounds was used to normalize data.
Non-targeted lipidomics analysis
Samples and standards were analyzed as described . The sample temperature was maintained at 4 °C.
Metabolomic and lipidomic analyses were performed on all 277 plasma samples by gas chromatography-time of flight mass spectrometry (GC-TOF/MS). Patient samples were analyzed along with 35 (metabolomics) and 36 (lipidomics) splits of the “reference” (pooled) quality control sample placed evenly in between patient samples for quality control assessment. Run order of patient samples was randomly determined. Note that for the lipidomics data, one patient sample showed anomalous behavior and was markedly dispersed from the rest of the patient samples (i.e., contained values above and below 3 standard deviations from the mean) and thus removed from further analysis.
A total of 116 plasma metabolites and 101 lipids were identified, quantitated, and included in downstream statistical analysis. Also 174 unknown plasma metabolites were analyzed separately in the same manner described below. For lipidomic data, lipids quantified in fewer than 50% of samples were discarded from downstream analysis to reduce the bias that could be induced by imputation due to the detection limit. Unobserved values for any remaining undetected lipid below the pre-defined detection limit were imputed as one-half of the lipid-specific minimum of the observed values prior to normalization and statistical inference tests.
For between-sample normalization, the intensity values for each sample were summed, and the median value of the sums across all samples were determined. The intensity values of each sample were then scaled such that the sum of the scaled intensities equaled the median value of all samples. Thus, the sum of the scaled intensities was the same for all samples. Normalized intensity values were then log2 transformed to reduce the influence of extreme values and to meet the homogeneity of variance assumption.
Prior to assessing the association between metabolite intensity levels and disease severity (baseline eGFR and ht-TKV), we assessed the distributions of eGFR and ht-TKV and then log transformed them to meet normality and homoscedasticity of linear regression models. To identify individual metabolites whose intensities are significantly associated with eGFR and ht-TKV, we performed association analyses using linear regression with each metabolite signal level as the primary predictor variable and baseline eGFR and ht-TKV as the two outcomes of interest, analyzed separately. All analyses accounted for the covariates attributed to the HALT-A study design including age, sex, PKD genotype (PKD1 or PKD2), and participating clinical site of sample collection. Other HALT covariates, including baseline systolic and diastolic blood pressure, weight, body mass index (BMI), height-adjusted liver volume, and urine levels for sodium, potassium, and creatinine, were evaluated for inclusion in multiple linear regression models to determine whether any of those covariates are influential in estimating the magnitude of association of metabolite levels with the outcomes. Variable inclusion and exclusion were determined based on Akaike weight-based relative importance (RI) scores , with a suggestive threshold of 0.80 in conjunction with likelihood ratio tests (LRT) at p-value < 0.05. Multicollinearity was inspected for possible collinearity issues using variance inflation factor (VIF) score with a suggestive threshold of 2.5. We investigated potential effect modifiers such as sex and PKD genotype by adding interaction terms in regression models. When an interaction term reached statistical significance, post-hoc analysis was performed stratified by the effect modifier. Five missing values for BMI were imputed with the mean of the observed values. We used the Storey’s false discovery rate (FDR) q-values to correct for multiple testing . All statistical analyses were conducted in R (version 3.4.) language and environment.
Reproducibility of metabolomics and lipidomics data
Due to the large number of samples (277; Table 1) and the necessity to process the samples over a prolonged time period of several days, we first investigated reproducibility and time-dependent variation using the repeated measurements from the 35 and 36 technical replicates, for metabolomics and lipidomics respectively of the reference quality control sample which were inserted throughout the GC-TOF/MS analysis to assess for systematic shifts among samples. Boxplots of the intensity values of the reference sample replicates revealed no systematic time-dependent pattern in the magnitude or distribution of intensities (Additional file 1: Figure S1 and S2). In order to provide a direct quantitative measure of technical between-sample variability, we calculated the coefficient of variation (CV) across the 35 and 36 measurements for each metabolite and lipid in the reference sample. The CV values of the 35–36 replicates ranged from 0.006 to 0.445 with mean = 0.051 and SD = 0.045 for metabolites and from 0.232 to 0.250 with mean = 0.241 and SD =0.004 for lipids. These low CV values demonstrate adequate technical reproducibility of sample analyses during the course of the GC-TOF/MS experiment.
Selection of accountable covariates in association with metabolite signals of baseline eGFR and ht-TKV
To determine which covariates should be included in multiple regression models, Akaike weight-based RI scores were calculated as a merit for importance of inclusion of each covariate in a model for all metabolites (ranged from 0.29 to 0.905 for eGFR and from 0.28–0.812 for ht-TKV) and final variable selection was carried out via likelihood ratio test (LRT). Age, genotype, sex, site of sample collection, and BMI were significantly (p < 0.05) influential in association between metabolites and outcomes, eGFR and ht-TKV (Fig. 1) and thus we retained these covariates in the final models. Although the site of sample collection was not significant for ht-TKV at p-value < 0.05, we decided to include it in the final model to account for possible study site differences. The variable selection procedures led to the following final models for association analysis explicitly:
eGFR = α + β1*M_Signal + β2*Age + β3*Genotype + β4*Sex + β5*Site + β6*BMI.
ht-TKV = α + β1*M_Signal+ β2*Age + β3*Genotype + β4*Sex + β5*Site + β6*BMI.
where M_ Signal is the signal level of a given metabolite/lipid.
We further investigated potential multicollinearity issues among the selected variables in the above models and found no presence of collinearity.
Identification of plasma metabolites and lipids associated with baseline eGFR and ht-TKV
After accounting for the covariates discussed above, we discovered that 20 metabolites significantly associated with eGFR at p < 0.05. Of those, 12 metabolites remained significant at FDR q-value < 0.05 (Table 2). All of the significant metabolites exhibited negative associations with eGFR, meaning that the metabolite levels were increased with decreasing eGFR levels (Fig. 2). Of note, indole-3-lactate is an exogenously produced microbial metabolite, while pseudo-uridine is a known uremic toxin (see Discussion).
We next evaluated metabolites that associated with ht-TKV. Ten metabolites were found to be significantly associated with ht-TKV at p < 0.05; however, none of them maintained their significance at FDR q-value < 0.05 (Table 3). Four metabolites (taurine, tyrosine, glutamic acid, and 2-hydoxybutanoic acid) were negatively correlated with ht-TKV, while the others exhibited a positive correlation (Fig. 2).
In analyzing the lipidome, we discovered 11 lipids which were significantly associated with eGFR at p < 0.05. However, no lipids remained significant at FDR q-value < 0.05. Acylcarnitines and triglycerides exhibited negative associations with eGFR, and phosphatidylcholines and diglyceride exhibited positive associations (Table 4). Upon analysis of the relationship between the lipidome and ht-TKV, we found that five triglycerides and one sphingomyelin were significantly associated with ht-TKV at p < 0.05, however, of these only two triglycerides maintained their significance at FDR q-value < 0.05 (Table 5). All six of these lipids exhibited positive associations with ht-TKV.
Assessment of effect modification by sex or PKD genotype
We further investigated whether the relationship between metabolite level and eGFR or ht-TKV differ depending on patient’s sex or their PKD genotype, to evaluate the possibility that either could have a significant role in metabolomics signature expression in association with eGFR or ht-TKV. To accomplish this, we added an interaction term to the model depicted by equation 1) above for eGFR and to the model depicted by equation 2) above for ht-TKV. For eGFR, none of the metabolites showed significant effect modification by sex or genotype at FDR q-value< 0.05 (Additional file 2: Table S1). For significant metabolites at p < 0.05, the direction of the metabolite-eGFR association was generally reversed by sex. Positive metabolite-eGFR associations for females were altered to negative associations for males, and vice-versa; similar data were demonstrated for metabolite-PKD genotype associations (Fig. 3). By contrast, creatinine exhibited negative associations for both males and females, but at different degrees: the magnitude of association was 3-fold greater in males compared to females (Fig. 3 (a), Additional file 2: Table S1).
For ht-TKV, none of the metabolites showed significant effect modification by sex or genotype at FDR q-value < 0.05 (Additional file 2: Table S2). The majority of the metabolites that displayed significant effect modification by sex or genotype at p < 0.05 modified the direction of the metabolite-ht-TKV association (Fig. 3). Pseudo-uridine and myo-inositol exhibited positive association regardless of sex, but their magnitude of association differed depending on sex. The degree of association for pseudo-uridine was 7-fold greater in males compared to females, and for myo-inositol was 15-fold greater in males than females (Fig. 4(a), Additional file 2: Table S2).
For both eGFR and ht-TKV, no lipids showed significant effect modification by sex or genotype at FDR q-value < 0.05 (Additional file 1: Tables S1 & S2). The majority of those lipids significant at p-value < 0.05 had the directionality of their associations with eGFR (Fig. 3) and ht-TKV (Fig. 4) altered by sex and genotype, respectively. It should be noted that none of the metabolites and lipids reached FDR-adjusted significance after taking multiple testing into account and thus further investigation is needed to validate whether sex or genotype are necessary to treat them as effect modifiers in the associations with eGFR or ht-TKV for metabolites or lipids.
Association of metabolites and lipids with eGFR and ht-TKV are exclusively independent
In the above results, we have identified metabolites and lipids associated with eGFR and ht-TKV. As eGFR and ht-TKV are clinically correlated traits for disease severity, we decided to perform further analysis using those outcomes as covariates in an attempt to identify metabolites associated only with eGFR or ht-TKV exclusively independent of each other. To accomplish this, ht-TKV was added as an additional covariate to the model depicted by equation 1) for eGFR as previously described, while eGFR was added as an additional covariate to the model depicted by equation 2) for ht-TKV. When controlled for ht-TKV as a covariate, four metabolites were determined to be significantly associated with eGFR at FDR q-value < 0.05, including pseudo-uridine, indole-3-lactate, creatinine, and uric acid (Additional file 2: Table S3); all of which were among the twelve metabolites previously identified as significant at FDR q-value< 0.05 for eGFR without controlling for ht-TKV (Table 2). Directionality of the eGFR-metabolite association remained unchanged compared to the ht-TKV unadjusted analysis results. All significant metabolites exhibited negative associations with eGFR independent of ht-TKV, with the exception of pipecolinic acid which demonstrated a positive relationship with eGFR. No lipids were determined to be significant at the FDR q-value< 0.05 level (Additional file 2: Table S3).
When evaluating ht-TKV while controlling for eGFR as a covariate, only taurine was found to be significant at p < 0.05, although it was not significant at FDR q-value < 0.05 (Additional file 2: Table S3). Taurine was also found to be associated with ht-TKV in the previous analysis that did not control for eGFR as a covariate (Table 3). Adjustment for eGFR did not alter the directionality of association between the metabolite and ht-TKV. In addition, triglyceride (51:3) was significant at FDR q-value< 0.05, as previously found without controlling for eGFR (Table 5), while the directionality of the association remains unchanged regardless of controlling for eGFR as a covariate in the analysis. The eGFR- and ht-TKV unique metabolites therefore, suggest that even those two highly correlated traits may be mediated, at least in part, by different metabolic pathways in kidney function.
Many metabolic disorders have been associated with ADPKD , including hypocitraturia [25, 26], hypomagnesiuria [25, 27, 28], hyperuricemia and gout , hypophosphataemia , altered arginine metabolism , abnormal lipid metabolism [31, 32] and nephrolithiasis [25, 27]. Recently the Warburg effect metabolic change that is characteristic of cancer has been identified in ADPKD .C-Glucose tracking in a kidney specific PKD1-null mouse model found that cystic kidneys had increased aerobic glycolysis, and transcriptomics of human ADPKD kidneys revealed gluconeogenesis enzymes were downregulated while some glycolytic enzymes were upregulated . For this report, we asked which plasma metabolites can serve as indicative of baseline eGFR and ht-TKV in the HALT A cohort of hypertensive ADPKD patients with relatively preserved (> 60 ml/min) kidney function. We found several metabolites that significantly correlated with eGFR and two lipids that correlated with ht-TKV at FDR q-value< 0.05. All of the metabolites which correlated with eGFR had plasma concentrations inversely correlated with the patient’s eGFR. Since all of the significant (FDR q < 0.05) metabolites showed negative associations with eGFR, it is possible that their increased presence is due merely to decreased filtration or tubular reabsorption. However, in light of previous findings of these metabolites in uremic serum, it is possible that disease specific changes are related to the utilization or transport of these metabolites. It is of great interest that these metabolites were seen to be increased in individuals with relatively preserved kidney function (Table 1), which suggests that there may be an early disruption of their transport or utilization in patients with kidney disease considered to be early with normal kidney function.
Non-targeted metabolomics is a powerful tool to generate hypotheses and discover novel small molecular biomarkers that may have therapeutic impact. However, there are limitations to this technique , some of which are apparent in this study. In most GC/MS-based metabolomics studies, it is common to come across unidentified metabolites for which a mass/charge ratio is present but structural identification is lacking. Many of these compounds can be identified by nuclear magnetic resonance or other advanced target-specific techniques, but such analysis requires prior knowledge of specific targets as well as a substantial investment and considerably larger quantities of starting material . In this study, one such metabolite, 191801, was found to be highly correlated with eGFR (FDR q < 0.0001) even after controlling for ht-TKV as a covariate (Additional file 2: Table S3). While we can infer properties and the possible origin of this metabolite based on other studies in the analyzing laboratory (i.e. that it is likely of microbial origin), structural information on this metabolite at this point is lacking.
We are the first to examine the metabolome of plasma from ADPKD patients with intact kidney function. Other ADPKD metabolomics studies have been done on mouse urine [36, 37], mouse kidneys  and mouse cell culture media . Of the twelve identified plasma metabolites which we found to be associated with eGFR, pseudo uridine and uric acid have long been known to be present in uremic serum  and these compounds have more recently been shown to be markers of high risk or progression to ESRD in type 2 diabetes . Higher baseline serum uric acid in ADPKD patients has also been retrospectively associated with a high risk of ESRD  and uric acid is increased in the urine of a juvenile cystic mouse model . One study found an association between increased serum uric acid and decreased creatinine clearance rate in ADPKD, without any evidence of increased uric acid production , a second study found ADPKD patients had a higher incidence of hyperuricemia than patients with CKD or normal renal function despite all groups having normal rates of uric acid excretion  and a third study found ADPKD patients with normal renal function had elevated serum uric acid . Neither study found evidence of increased uric acid synthesis. The CRISP study found that higher baseline serum uric acid was not associated with declines in eGFR or increases in ht-TKV in ADPKD patients over a 6-year period  so the mechanism behind a higher incidence of hyperuricemia in ADPKD remains unresolved. We are currently undergoing analyses to determine if any of these 12 metabolites at baseline in ADPKD are predictive of faster disease progression. In addition, urine metabolites are being evaluated for a future study, as it is clear that there is high biomarker potential in this biofluid [8, 37, 43, 44].
Several metabolites, including indole-3-lactate and likely the unknown compound 191801, are microbial metabolites, the former being the end-product of tryptophan metabolism produced by gut microbes , and the latter is found in wine, microalgae, fungi and yeast, as well as in human urine and plasma and is therefore likely a microbial metabolite (please see https://binvestigate.fiehnlab.ucdavis.edu/#/bin/191801 for further analytical information on the metabolite 191801 as well as in what species it has been observed). Indole-3-lactate and indole-3-acetate are both present in the urine of a juvenile cystic mouse model . Meanwhile creatinine, which is a component of the CKD-EPI eGFR equation, was found to be significantly associated with eGFR levels, which serves as a supportive positive internal control for the analysis presented here.
We found significant associations between two large chain length triglycerides and baseline values of ht-TKV in ADPKD patients. This kind of association between specific triglycerides and a marker of progression risk has not been previously described in early ADPKD patients with relatively intact kidney function. Lipid dysregulation is associated with ADPKD in mouse models, including impaired fatty acid oxidation in the kidney  and upregulation of apolipoprotein-related genes in PKD-1 null kidneys . Fatty acid metabolism was found to be altered in ADPKD patients both with normal renal function (HALT A) and decreased renal function (HALT B) . A comparison between ADPKD with normal kidney function and healthy patients found only 3 components of metabolic syndrome were increased in ADPKD patients, namely hypertension, fasting glycemia and abdominal obesity . A different study of young ADPKD patients with normal renal function found increased insulin resistance, inflammation and low density lipoprotein (LDL) in ADPKD . A study in Southern India found ADPKD patients of a wide variety of ages had moderately increased plasma triglycerides, LDL and very low density lipoprotein (VLDL), and a ~ 50% reduction in high density lipoprotein (HDL) compared to healthy controls . The Consortium for Radiologic Imaging Studies of Polycystic Kidney Disease (CRISP) cohort 6-year study found baseline serum HDL-cholesterol levels had a protective effect against increased kidney growth and faster declines in GFR . Hypertriglyceridemia is not a feature of ADPKD in contrast to CKD (reviewed in ), which is also characterized by changes in acylcarnitine levels (reviewed in ). We are currently undergoing analyses to determine if the plasma triglyceride concentrations at baseline in ADPKD are predictive of faster disease progression.
We also evaluated the metabolites when accounting for sex and genotype as possible effect modifiers in their association with outcome. We found evidence that metabolite-outcome associations might be altered by sex or by genotype. However, none of the metabolites were found to be significant at FDR q < 0.05 when controlling for multiple testing so neither sex nor PKD genotype are likely a differentiating factor in the association of all or most metabolites with those outcomes, but nonetheless further investigation is warranted.
In summary, we have shown the presence of metabolites and lipids which highly correlate with eGFR and ht-TKV in ADPKD patients with normal kidney function. These metabolites and lipids can be useful for uncovering metabolic derangements, as well as reprogramming, which appear to occur early in disease. With further evaluation we may find additional predictors of early kidney dysfunction that are potentially less costly than measuring ht-TKV, which is currently the only measurement known to predict the rate of disease progression. Perhaps these metabolites will allow for prognostication and stratification of patients into slow and rapid progressors to ESRD in ADPKD.
This study identifies metabolic derangements in early ADPKD which may be prognostic for ADPKD disease progression. Such data will be useful for future studies designed to predict outcome of disease based on such early metabolic changes and will likely lead to new therapeutic paradigms for a disease with quite limited therapeutic options.
Autosomal dominant polycystic kidney disease
Body mass index
Chronic kidney disease
Consortium for Radiologic Imaging Studies of Polycystic Kidney Disease
N-Cyclohexyl-N′-dodecanoic acid urea
Coefficient of variation
Estimated glomerular filtration rate
End stage renal disease
Storey’s false discovery rate
Gas chromatography-time of flight/mass spectrometry
High density lipoprotein
Height corrected total kidney volume
Low density lipoprotein
Likelihood ratio tests
Methyl tert-butyl ether
Variance inflation factor
Very low density lipoprotein
Consortium TEPKD. The polycystic kidney disease 1 gene encodes a 14 kb transcript and lies within a duplicated region on chromosome 16. Cell. 1994;77:881–94.
Mochizuki T, Wu G, Hayashi T, et al. PKD2, a gene for polycystic kidney disease that encodes an integral membrane protein. Science (New York, NY). 1996;272:1339–42.
Porath B, Gainullin VG, Cornec-Le Gall E, et al. Mutations in GANAB, encoding the glucosidase IIalpha subunit, cause autosomal-dominant polycystic kidney and liver disease. Am J Hum Genet. 2016;98:1193–207.
Qian F, Germino FJ, Cai Y, Zhang X, Somlo S, Germino GG. PKD1 interacts with PKD2 through a probable coiled-coil domain. Nat Genet. 1997;16:179–83.
Boletta A, Qian F, Onuchic LF, et al. Polycystin-1, the gene product of PKD1, induces resistance to apoptosis and spontaneous tubulogenesis in MDCK cells. MolCell. 2000;6:1267–73.
Braun WE, Abebe KZ, Brosnahan G, et al. ADPKD progression in patients with no apparent family history and no mutation detected by sanger sequencing. Am J Kidney Dis. 2018;71:294–6.
Chebib FT, Torres VE. Autosomal dominant polycystic kidney disease: Core curriculum 2016. Am J Kidney Dis. 2016;67:792–810.
Ganti S, Taylor S, Kim K, et al. Urinary acylcarnitines are altered in kidney cancer. IntJCancer. 2012;130:2791–800.
Chapman AB, Johnson AM, Gabow PA, Schrier RW. Overt proteinuria and microalbuminuria in autosomal dominant polycystic kidney disease. Journal of the American Society of Nephrology : JASN. 1994;5:1349–54.
Torres VE, Grantham JJ, Chapman AB, et al. Potentially modifiable factors affecting the progression of autosomal dominant polycystic kidney disease. Clinical journal of the American Society of Nephrology : CJASN. 2011;6:640–7.
Gabow PA, Johnson AM, Kaehny WD, et al. Factors affecting the progression of renal disease in autosomal-dominant polycystic kidney disease. Kidney Int. 1992;41:1311–9.
Schrier RW, Abebe KZ, Perrone RD, et al. Blood pressure in early autosomal dominant polycystic kidney disease. NEnglJMed. 2014;371:2255–66.
Flowers EM, Sudderth J, Zacharias L, et al. Lkb1 deficiency confers glutamine dependency in polycystic kidney disease. Nat Commun. 2018;9:814.
Rowe I, Chiaravalli M, Mannella V, et al. Defective glucose metabolism in polycystic kidney disease identifies a new therapeutic strategy. Nat Med. 2013;19:488–93.
Trott JF, Hwang VJ, Ishimaru T, et al. Arginine reprogramming in ADPKD results in arginine-dependent cystogenesis. Am J Physiol Renal Physiol. 2018.
Hwang VJ, Kim J, Rand A, et al. The cpk model of recessive PKD shows glutamine dependence associated with the production of the oncometabolite 2-hydroxyglutarate. American Journal of Physiology: Renal Physiology. 2015;309:F492–F8.
Yu ASL, Shen C, Landsittel DP, et al. Baseline total kidney volume and the rate of kidney growth are associated with chronic kidney disease progression in autosomal dominant polycystic kidney disease. Kidney Int. 2018;93:691–9.
Perrone RD, Mouksassi MS, Romero K, et al. Total kidney volume is a prognostic biomarker of renal function decline and progression to end-stage renal disease in patients with autosomal dominant polycystic kidney disease. Kidney international reports. 2017;2:442–50.
Torres VE, Abebe KZ, Chapman AB, et al. Angiotensin blockade in late autosomal dominant polycystic kidney disease. NEnglJMed. 2014;371:2267–76.
Fiehn O. Metabolomics by gas chromatography-mass spectrometry: combined targeted and untargeted profiling. Current protocols in molecular biology 2016;114:30.4.1–.4.2.
Cajka T, Smilowitz JT, Fiehn O. Validating quantitative untargeted Lipidomics across nine liquid chromatography-high-resolution mass spectrometry platforms. Anal Chem. 2017;89:12360–8.
Johnson JB, Omland KS. Model selection in ecology and evolution. Trends Ecol Evol. 2004;19:101–8.
Storey J. A direct approach to false discovery rates. J Royal Stat Soc B. 2002;64:479–98.
Mao Z, Xie G, Ong AC. Metabolic abnormalities in autosomal dominant polycystic kidney disease. Nephrol Dial Transplant. 2015;30:197–203.
Torres VE, Wilson DM, Hattery RR, Segura JW. Renal stone disease in autosomal dominant polycystic kidney disease. Am J Kidney Dis. 1993;22:513–9.
Masoumi A, Reed-Gitomer B, Kelleher C, Bekheirnia MR, Schrier RW. Developments in the management of autosomal dominant polycystic kidney disease. Ther Clin Risk Manag. 2008;4:393–407.
Nishiura JL, Neves RF, Eloi SR, Cintra SM, Ajzen SA, Heilberg IP. Evaluation of nephrolithiasis in autosomal dominant polycystic kidney disease patients. Clin J Am Soc Nephrol. 2009;4:838–44.
Idrizi A, Barbullushi M, Koroshi A, et al. Urinary tract infections in polycystic kidney disease. Med Arh. 2011;65:213–5.
Mejias E, Navas J, Lluberes R, Martinez-Maldonado M. Hyperuricemia, gout, and autosomal dominant polycystic kidney disease. Am J Med Sci. 1989;297:145–8.
Pavik I, Jaeger P, Kistler AD, et al. Patients with autosomal dominant polycystic kidney disease have elevated fibroblast growth factor 23 levels and a renal leak of phosphate. Kidney Int. 2011;79:234–40.
Allen E, Piontek KB, Garrett-Mayer E, Garcia-Gonzalez M, Gorelick KL, Germino GG. Loss of polycystin-1 or polycystin-2 results in dysregulated apolipoprotein expression in murine tissues via alterations in nuclear hormone receptors. Hum Mol Genet. 2006;15:11–21.
Menezes LF, Lin CC, Zhou F, Germino GG. Fatty acid oxidation is impaired in an orthologous mouse model of autosomal dominant polycystic kidney disease. EBioMedicine. 2016;5:183–92.
Chiaravalli M, Rowe I, Mannella V, et al. 2-Deoxy-d-glucose ameliorates PKD progression. J Am Soc Nephrol. 2016;27:1958–69.
Weiss RH, Kim K. Metabolomics in the study of kidney diseases. NatRevNephrology. 2011;8:22–33.
Wettersten HI, Weiss RH. Applications of metabolomics for kidney disease research: from biomarkers to therapeutic targets. Organogenesis. 2013;9.
Menezes LF, Zhou F, Patterson AD, et al. Network analysis of a Pkd1-mouse model of autosomal dominant polycystic kidney disease identifies HNF4alpha as a disease modifier. PLoS Genet. 2012;8:e1003053.
Taylor SL, Ganti S, Bukanov NO, et al. A metabolomics approach using juvenile cystic mice to identify urinary biomarkers and altered pathways in polycystic kidney disease. American Journal of Physiology: Renal Physiology. 2010;298:F909–F22.
Niwa T, Takeda N, Yoshizumi H. RNA metabolism in uremic patients: accumulation of modified ribonucleosides in uremic serum. Technical note Kidney Int. 1998;53:1801–6.
Niewczas MA, Sirich TL, Mathew AV, et al. Uremic solutes and risk of end-stage renal disease in type 2 diabetes: metabolomic study. Kidney Int. 2014;85:1214–24.
Helal I, McFann K, Reed B, Yan XD, Schrier RW, Fick-Brosnahan GM. Serum uric acid, kidney volume and progression in autosomal-dominant polycystic kidney disease. Nephrol Dial Transplant. 2013;28:380–5.
Hosoya T, Ichida K, Tabe A. Sakai O. A study of uric acid metabolism and gouty arthritis in patients with polycystic kidney. Nihon Jinzo Gakkai Shi. 1993;35:43–8.
Lai S, Mastroluca D, Matino S, et al. Early markers of cardiovascular risk in autosomal dominant polycystic kidney disease. Kidney Blood Press Res. 2017;42:1290–302.
Kind T, Tolstikov V, Fiehn O, Weiss RH. A comprehensive urinary metabolomic approach for identifying kidney cancer. AnalBiochem. 2007;363:185–95.
Kim K, Taylor SL, Ganti S, Guo L, Osier MV, Weiss RH. Urine metabolomic analysis identifies potential biomarkers and pathogenic pathways in kidney cancer. OMICS. 2011;15:293–303.
Wikoff WR, Nagle MA, Kouznetsova VL, Tsigelny IF, Nigam SK. Untargeted metabolomics identifies enterobiome metabolites and putative uremic toxins as substrates of organic anion transporter 1 (Oat1). J ProteomeRes. 2011;10:2842–51.
Klawitter J, Klawitter J, McFann K, et al. Bioactive lipid mediators in polycystic kidney disease. J Lipid Res. 2014;55:1139–49.
Pietrzak-Nowacka M, Safranow K, Byra E, Binczak-Kuleta A, Ciechanowicz A, Ciechanowski K. Metabolic syndrome components in patients with autosomal-dominant polycystic kidney disease. Kidney Blood Press Res. 2009;32:405–10.
Veeramuthumari P, Isabel W. Clinical study on autosomal dominant polycystic kidney disease among south Indians. International Journal of Clinical Medicine. 2013;2013:200–4.
Barter P. Lipoprotein metabolism and CKD: overview. Clin Exp Nephrol. 2014;18:243–6.
Hocher B, Adamski J. Metabolomics for clinical use and research in chronic kidney disease. Nat Rev Nephrol. 2017;13:269–84.
This work was supported by NIH grant 1R01DK107416-01A1 and by a grant from Dialysis Clinics Incorporated. Neither of the funding bodies played any role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Eligible participants in the HALT PKD trial were enrolled at seven clinical sites from February 2006 through June 2009. All the participants provided written informed consent.
Consent for publication
Partial funding was obtained from a grant from Dialysis Clinics Incorporated but DCI had no impact on the content of this article.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Boxplots of the distribution of Metabolite intensity levels (on log 2 scale) for all 35 reference samples of metabolomics data. Figure S2. Boxplots of the distribution of Lipid intensity levels (on log 2 scale) for all 36 reference samples of lipidomics data. (PPTX 103 kb)
Table S1. Significant (p < 0.05) effect modification of patient’s sex/genotype on metabolite/lipid-eGFR association. Table S2. Significant (p < 0.05) effect modification of patient’s sex/genotype on metabolite/lipid-ht-TKV association. Table S3. Plasma metabolites and lipids significantly (p < 0.05) associated with (a) eGFR, after controlling for ht-TKV as a covariate and with (b) ht-TKV, after controlling for eGFR as a covariate. (XLSX 19 kb)