1. Introduction
Chronic kidney disease (CKD), defined by abnormalities of renal function and/or structure over 3 months with clinical consequences, is a global healthcare issue that has been predicted to become one of the leading global causes of death by 2040 [
1]. Furthermore, its prevalence has significantly increased in the last decades, so much so that up to 10-15% of adults worldwide are affected. The main burden of CKD is not, as it is perceived by many, dialysis, and kidney transplant; rather, accelerated aging and premature death, often associated with increased cardiovascular risk, are the most prominent concerns [
2].
One specific aspect of CKD, particularly in early stages, is that often times patients do not show overt clinical symptoms; therefore, laboratory tests are crucial to monitor the disease. In this regard, the assessment of glomerular filtration rate (GFR) is key to evaluate renal function and is usually estimated (eGFR) based on serum creatinine concentrations using different equations, of which CKD-EPI is presently recommended. However, creatinine levels begin to stand out only when the disease is well under way and have been shown to be less accurate under special circumstances. In addition, eGFR is of little help in assessing the risk of progression to end-stage kidney disease (ESKD). In fact, even though CKD prevalence based on eGFR measurements is very similar in Western countries, the incidence of patients requiring renal replacement therapy is very different among them [
3]. Available clinical evidence strongly indicates that patients could greatly benefit from novel biomarkers that are able to early identify patients at risk.
Renal function affects the levels of metabolites in blood; hence, it is likely that some metabolites could be helpful in estimating kidney function. In this regard, the kidney plays a crucial role in the metabolism of proteins and amino acids (AA), which are constantly filtered and reabsorbed [
4]. Quantitative analysis of free AA in blood and urine has been shown useful for the diagnosis and management of metabolic disturbances as well as for the assessment of nutritional status, tissue damage and renal function [
5]. Indeed, the advent of metabolomics, which identifies and quantitatively compares small-molecule metabolites between different groups, has allowed to identify several circulating AA and related compounds that may be associated with CKD [
6,
7,
8]. However, results are mostly based on small cohorts of renal patients and oftentimes a comparison between different stages of the disease is lacking.
The goal of this work was to determine plasma concentrations of AA and biogenic amines in a large group of healthy subjects, CKD patients and patients with ESKD, to identify differences in their metabolic profiles and analytes with the potential to become biomarkers of renal disease.
2. Materials and Methods
2.1. Subjects
Between 2017 and 2022, 820 Caucasian Spanish subjects were recruited, namely (i) 231 patients diagnosed with different CKD stages at the Nephrology Service of the Badajoz University Hospital (Spain); (ii) 325 patients with end-stage kidney disease (ESKD, eGFR< 15 ml/min/1.73 m2), who were enrolled at the Advanced CKD Unit of the Badajoz University Hospital and at three Dialysis Units (Fresenius clinic and Llerena and Zafra hospitals in Badajoz, Spain); and finally, (iii) 264 healthy volunteers who were recruited from primary care clinics in Soria, Spain. Transplantation, pregnancy or breastfeeding, active infection, cancer, or acute kidney injury were all considered exclusion criteria. All patients were over 18 years of age and gave written consent for their participation in the study, which was approved by the Ethics Committees of the Badajoz University Hospital (21/04/2022) and was carried out in accordance with the Declaration of Helsinki and its subsequent revisions.
Diagnostic and prognostic stratification of renal patients was conducted with the KDIGO classification and table of progression risk and the CONSORTIUM-CKD equation. Kidney function was assessed by the CKD-EPI equation. Clinical records were examined to retrieve clinical, demographic and biochemical data.
2.2. Sample collection and amines analyzed
Blood samples (3 ml) were drawn in EDTA tubes from each participant and plasma was immediately obtained by centrifuging 10 minutes at 3000 rpm. Plasma samples were then aliquoted (500 μl) and stored at -80 ºC until analysis. Essential AA (EAA): histidine, isoleucine, leucine, lysine, methionine, phenylalanine, threonine, tryptophan, and valine; and non-essential AA (NEAA): alanine, arginine, asparagine, aspartic acid, citrulline, glutamic acid, glutamine, glycine, proline, L-serine and tyrosine were determined. In addition, seven other related metabolites and biogenic amines, namely kynurenic acid, kynurenine, asymmetric dimethyl arginine (ADMA), D-serine, serotonin, gamma aminobutyric acid (GABA), acetylcholine and creatine were also analyzed.
Ten different ratios involving AA were calculated as described by Duranton et al. [
9]. EAA/NEE, Val/Gly and Alanine/Branched-chain amino acids (BCAA) ratios are used to assess the overall nutritional state; the kidney’s ability for interconversions of AA is measured with the Tyr/Phe, Ser/Gly and Arg/Cit ratios. The relative importance of both the kynurenine pathway and ADMA was assessed with the Kynurenine/Trp ratio and the Arg/ADMA ratio. Finally, the Fisher’s ratio (defined as the sum of BCAA, Leu, Ile and Val, divided by the sum of aromatic AA, Tyr and Phe), and the GSG index [Glu/(Ser+Gly)] were also calculated, as they have been observed to be altered in several pathologies and metabolic dysfunctions [
10,
11].
2.3. Chemical and reagents
All chemical reagents were purchased from Sigma Aldrich (St. Louis, MO). The water, methanol and acetonitrile (J.T. Baker, Innovagen SL) used for the mobile phases were all LC-MS grade. The stock solutions of Asparagine, Aspartate, Glycine, Serine (300 µM); Alanine, Citrulline, Glutamate, Glutamine, Histidine, Kynurenic acid, Kynurenine, Methionine, Phenilalanine, Proline, Threonine, Tryptophan, Valine (50 µM); GABA, Arginine, ADMA, Leucine, Isoleucine, Lysine, Serotonin, Tyrosine, Acetylcholine (5 µM); and 20 mM isotopically labeled d4-acetylcholine were prepared in water and kept in artificial cerebrospinal fluid (CSF), which consisted of 145 mM NaCl; 2.68 mM KCl; 1.4 mM CaCl2; 1.0 MgSO; 1.55 mM Na2HPO4; 0.45 mM NaH2PO4 pH 7.4 at −80ºC (adjusted with NaOH).
2.4. Sample preparation and LC/MS/MS analysis
The procedure to prepare the samples and the calibration standards, listed in Supplemental
Table S1, was the same, except that the samples started with a volume of 20 µl and the calibration standard with 5 µl. Acetonitrile (1:4) was then added to collect 20 µl of the supernatant once the samples were centrifuged at 13,000 rpm for 10 minutes. Subsequently, 10 µl of 100 mM Na2CO3, 10 µl of BzCl (2% v/v in acetonitrile) and 10 µl of the internal standard were incorporated. To balance the concentration of organic content, 50 µl of water were added. The calibration standards and the internal standard were frozen at -20°C in aliquots to prevent multiple freeze/thaw cycles. An aliquot of internal standard was thawed on the day of use and a fresh benzoyl chloride (BzCl) solution was made daily.
Amines were analyzed by liquid chromatography coupled with tandem mass spectrometry LC/MS/MS. The technique was based on the method described by Wong et al. [
12]. BzCl derivatized samples were analyzed by LC-MS/MS using multiple reaction monitoring (MRM) dynamics with an Agilent 6410 TQ equipment. MRM conditions for all the metabolites are shown in Supplemental
Table S2). Five µl of sample (by triplicate) were injected into a HiP-ALS automatic injection module at room temperature. An ACE Excel 2 SuperC18 column (15 cm × 2.1 mm ID, 2 µm, 90 Å) kept at 27 °C was used for the separation of metabolites. Mobile phase A consisted of 10 mM 0.15% formic acid ammonium formate and mobile phase B, acetonitrile. The flow was set to 0.2 µl/min and the elution gradient was as follows: 5%; 0.0 min, 15% B; 0.01 min, 17% B; 0.5 min, 55% B; 14 min, 70% B; 14.5 min, 100% B; 18 min, 100% B; 19 min, 5% B; 19.1 min, 5% B; 24 min. An additional 10-minute period of column equilibration at 0% B was required to achieve reproducible chromatography. The required pressure over the gradient ranged from 117 to 254 bar. Electrospray ionization (ESI) was used in positive mode at 4000 V. The gas temperature and flow were 350°C and 11 L/min, respectively, with the nebulizer set to 15 psi. Automated peak integration was performed using Agilent MassHunter Quantitative Analysis for QQQ, version B.04.01. All peaks were visually inspected to ensure proper integration (
Figure 1).
2.5. Statistical analyses
Chi-square tests were utilized to compare categorical variables. Quantitative variables were compared between groups either by t-test / Mann-Whitney (2 groups) or by ANOVA / Kruskal-Wallis (>2 groups), depending on the normality of data distribution. Post-hoc tests were applied to examine differences between group pairs. Correlations between quantitative variables were assessed by the Spearman’s test.
In order to evaluate the impact of the different metabolic signatures, multivariate analyses were conducted with Umetrics SIMCA 18 software (Umeå/Malmö, Sweden) to create partial least squares-discriminant analysis (PLS-DA) and orthogonal PLS-DA (OPLS-DA) models, which were evaluated for fitness (R2) and predictive ability (Q2). This is a supervised dimension reduction analysis that allows visualization of how individuals within a group cluster with regard to their data-compressed principal components, i.e., it permits to stratify samples based on differences in profiles rather than individual variables. After removing outliers identified with PLS-DA, OPLS-DA models were constructed based on 80% of the population and validated in the remaining 20% (training and test sets). These models were then utilized to identify the most significantly altered metabolites between groups by calculating the Variable Importance in the Projection (VIP) score (values >1 are considered relevant).
To assess the potential use of these metabolites as biomarkers to discriminate either between controls and CKD patients or between CKD and ESKD, areas under receiving operating characteristic (ROC) curves (AUCs) were calculated with R software (proc package). AUCs for models containing traditional risk factors before and after the addition of the metabolic signatures were compared using the De Long’s test.
IBM SPSS v.22.0 (SPSS Inc., Chicago, IL, v.22.0) was utilized for the statistical analyses. Bonferroni correction for multiple testing was applied, resulting in a significance threshold of P <0.0016.
4. Discussion
CKD is a major public health concern whose burden has not declined at the same rate as other important non-communicable diseases [
13]. Early diagnosis of CKD would allow reducing cardiovascular complications as well as preventing fast progression to ESKD, whose patients are 100 times more likely to die of cardiovascular disease than the general population [
14]. It has been suggested that in many patients a low eGFR may reflect a variety of age-related coexisting comorbid conditions, and therefore be a better predictor of global health outcomes rather than of CKD progression outcomes [
15]. In this context, adding new biomarkers such as AA and biogenic amines profiling to the current evaluation of CKD could be useful to improving our existing prognostic tools.
In our study, the raw comparison of plasma levels showed great differences for all 28 compounds assayed between controls and CKD patients. Our results show that citrulline had the highest VIP score for the discrimination between them. Two large longitudinal studies, one American [
16] and one Korean [
17], including subjects without CKD at baseline, measured plasma metabolite concentrations during an 8-year follow-up. Both these studies also highlighted citrulline as being significantly associated with incident CKD. Citrulline metabolism occurs preferentially in proximal tubular cells [
18], therefore, renal injury could very well translate into the observed accumulation of plasma citrulline in CKD patients [
9,
19,
20]. In this regard, it has been shown that reduced kidney function is also associated with greater citrulline excretion, suggesting that elevated citrullinemia is due to metabolic causes and not just retention [
9]. This deficient citrulline biotransformation in the kidney may easily lead to decreased renal function as seen in animal models [
21], since it would imply a reduction of arginine bioavailability, an AA that is necessary for endothelial nitric oxide synthesis [
22].
Our findings revealed that the Kynurenin/Trp ratio showed the strongest correlation with the eGFR. This ratio, which has also been associated with CKD in the past [
17,
23], reflects the activity of indoleamine dioxygenase (IDO), which is the first step of the kynurenine pathway of tryptophan metabolism. The mechanism by which kynurenine accumulation and tryptophan declining plasma concentrations (higher ratios) may be associated with reduced kidney function is discussed below in the context of ESKD. The potential of this ratio to be a good measure of early CKD risk, a pressing need in this pathology, has been previously stressed by Lee et al., who reported an OR of 12.65 (6.55-24.44) for CKD prevalence and 3.20 (1.57-6.51) for CKD incidence [
17].
In general, differences in the measured metabolites between CKD and ESKD patients were less marked than those observed between controls and CKD patients. This could be partly explained by the large number of patients with CKD stage 4 in our population (n=120, 52% of the CKD group), who display eGFR values lower than 30 ml/min/1.73 m2, and therefore could be more similar to ESKD patients. In any case, the addition of our metabolic profile to a predictive model based on demographic, clinical and physiological characteristics improved up to 96.5% its ability to discriminate CKD from ESKD, an AUC value which is in the same range as that recently reported for a similar model with 100 metabolites (99.9%) [
19].
Despite the aforementioned low average eGFR in the CKD group, our findings showed 12 metabolites that were still significantly altered between ESKD and CKD. Interestingly, 10 of these, namely alanine, arginine, citrulline, glutamine, histidine, isoleucine, serine, tryptophan, valine and ADMA had also been shown to be able to discriminate between CKD and dialysis in a previous small study [
9]. Another recent metabolomic study reported that BCAA metabolism was downregulated in ESKD compared to CKD patients [
19]. We did not measure BCAA metabolites, but, in line with this report, levels of two of the three BCAA, valine and isoleucine, were also significantly reduced in our ESKD sample. In addition, our findings show that valine and leucine were the second and fourth most important metabolites to discriminate CKD from ESKD according to their VIP scores. There are several explanations for this finding. BCAA are anabolic, essential AA whose levels have been shown to decrease with metabolic acidosis, which worsens with CKD progression to ESKD because of an enhanced proteolysis and BCAA catabolism [
24]. Moreover, the levels of BCAA and other essential AA have also been shown to be lower in ESKD patients due to hemodialysis and reduced protein intake [
25].
The most relevant metabolite for CKD/ESKD discrimination in our study was tryptophan, which showed the highest VIP score. Tryptophan catabolism by IDO leading to kynurenin is known to increase with CKD progression and dialysis [
26,
27]. The upregulation of this pathway is supported by our aforementioned observation that the Kynurenin/Trp ratio displayed an inverse high correlation with eGFR. In the same line, two other metabolomic studies, although with only 32 and 77 patients [
9,
19], have also reported the increase of similar ratios (using different metabolites in the kynurenin pathway) in ESKD patients. Kynurenin is further biotransformed into several metabolites that play a relevant role in chronic inflammation, oxidative stress or apoptosis [
28,
29], which could be behind the increased inflammatory response and cell damage present in advanced CKD. Indeed, the modification of kynurenine pathway via inhibition of kynurenine hydroxylase has been shown to attenuate surgical complications in rodents [
30], which suggests that the downregulation of tryptophan catabolism to kynurenin might be a promising therapy to delay progression to ESKD.
This study has a number of strengths and limitations. A clear strong point was the sample size, which was far higher than that reported by the vast majority of studies analyzing biological concentrations of AAs in renal patients, and that allowed to capture the true variability of the population in terms of the concentrations of these compounds. Among the limitations, urine samples were not available for the control group and hence proteinuria values could not be obtained for these subjects. In the same manner, protein intake of participants was not assessed. However, even though a higher risk of malnutrition was observed for patients with advanced CKD (significantly lower Val/Gly and EAA/NEAA ratios), multivariate models showed how eGFR was a far better predictor of AA plasma concentrations than albumin values (data not shown), indicating that nutrition status by itself would not explain our findings.
In summary, we have shown how in a large population of healthy subjects and renal patients with a wide range of eGFR values, the use of a LC/MS-based targeted metabolomics approach may be useful to identify prognostic biomarkers of CKD risk and progression. Our results taken together indicate that plasma concentrations of AA and related amines have a high degree of predictive ability regarding renal function. In particular, citrulline and the tryptophan metabolic route to kynurenine showed the greatest potential as biomarkers, whose addition to the current evaluation of CDK risk and progression might improve the accuracy of the existing prognostic methods largely based on creatinine alone.