2. Materials and Methods
2.1. Study Design and Participants
This is a prospective, observational, single-center study of 50 adult patients with stage 5 CKD requiring renal replacement therapy who received a deceased donor KT, with a follow-up over six months post-transplant. Inclusion criteria included patients over 18 years old selected to receive a compatible deceased donor KT, while exclusion criteria included patients with advanced CKD not requiring renal replacement therapy, those on immunosuppressive drugs before the transplant, and those receiving combined transplants.
2.2. Data Collection
Clinical variables of the donor and recipient, along with baseline analytical data, were collected pre-transplant (T0). Subsequent clinical outcomes, analytical variables, and metabolomic samples were collected at two time points: T1 (one-week post-transplant) and T2 (six months post-transplant). For T1, samples were specifically collected on either day 7 or day 8 post-transplant to ensure that, in cases where patients required hemodialysis, the samples would be taken before the dialysis session to avoid any potential alterations caused by the hemodialysis process.
Blood samples for metabolomic analysis were collected and processed by centrifugation. The resulting plasma was then stored at -80°C until the completion of the study.
2.3. Metabolomic Analysis
For the metabolomic analysis, 190 μL of cold methanol was added to 10 μL of plasma to precipitate proteins and extract compounds. The mixture was vortexed for 30 seconds and incubated at -20ºC for 30 minutes. Subsequently, the samples underwent double centrifugation at 9,000 rpm for 10 minutes at 4ºC. From the resulting supernatant, 90 μL was transferred to a 96-well plate for liquid chromatography-mass spectrometry and 10 μL of an internal standard mix solution containing phenilanaline-d5, caffeine-d9, reserpine and leucine enkephalin at 20 μM, was added. Quality control (QC) samples were prepared by pooling 5 μL from each serum sample. A blank was prepared by replacing serum with ultrapure water to identify potential artifacts from tubes, reagents, and other materials.
Samples were analyzed on an Ultra-Performance Liquid Chromatography (UPLC) system coupled to a Q-Exactive Plus spectrometer. The chromatographic separation was performed using an Xbridge BEH Amide column (150 × 2 mm, 2.5 μm particle size; Waters). The autosampler was set to 4°C, and a flow rate of 95 μL/min was used with water, 20 mM ammonium acetate, and 20 mM ammonium hydroxide at pH 9.4, as Mobile Phase A and acetonitrile as Mobile Phase B. To minimize intra-batch variability and enhance reproducibility, the samples were injected in a random order, and a QC sample was analyzed every six plasma samples. Blank analyses were conducted at the end of the sequence.
Electrospray ionization was employed in both positive and negative modes, in full mass acquisition, with a resolving power of 140,000 and three events: mass ranges from 70 to 300 Da, 300 to 560 Da and from 350 to 1200. Data were acquired in the centroid mode. Raw data were converted to mzXML format using MSConvert and then processed in EI-MAVEN software (v0.11.0) by using an in-house polar compound library to generate a peak table containing m/z, retention times, and intensities of metabolites detected.
Data quality (reproducibility, stability) was evaluated by means of the internal standard’s stability and the QC’s coefficients of variation (CVs). Molecular features with CVs > 30% were removed from the data matrix. Additionally, a normalization method (LOESS) was used to eliminate intra-batch variability due to technical differences. Data from positive and negative modes were merged and filtered for the statistical analysis.
2.4. Statistical Analysis
Initially, descriptive statistical analyses were performed on the baseline variables of both donors and recipients, as well as on basic laboratory values and clinical outcomes. Subsequently, a statistical analysis was conducted on the metabolomic data obtained at time point T1, using the outcome variable of serum creatinine levels greater than 1.5 mg/dL at six months. The metabolomic dataset comprised 226 molecules from 48 patients (initially 50 patients; one patient was excluded due to early graft loss and another due to a sample error at T1). The categorical response variable (creatinine levels at six months post-transplant) had two levels: 26 patients with creatinine ≤ 1.5 mg/dL and 22 patients with creatinine > 1.5 mg/dL.
For the metabolomic analysis, a Partial Least Squares-Discriminant Analysis (PLS-DA) model was first constructed using all available metabolites. The model utilized seven principal components, and the optimal number of latent variables was selected through Leave-One-Out Cross-Validation (LOOCV) with the aim of minimizing classification error. A Variable Importance in Projection (VIP) score threshold of 1.5 was applied, identifying 22 polar metabolic variables as significant contributors to the initial model.
A second PLS-DA model was then constructed using only the 22 selected variables obtained in the first model. This final model utilized two principal components, with the number of latent variables again optimized using LOOCV to ensure the lowest classification error. Model validation was conducted using permutation testing with 1000 iterations, confirming that the R² coefficient obtained with the original response variable was greater than that obtained with the permuted response variable in 95% of cases. This provided an empirical p-value for model significance, adding robustness to our findings.
Performance metrics of the final model were evaluated based on LOOCV predictions, including the Receiver Operating Characteristic (ROC) curve, Area Under the Curve (AUC), accuracy, sensitivity, and specificity. These metrics provided a comprehensive assessment of model performance. Additionally, coefficients and 95% confidence intervals (CI) for each variable were obtained using the Jack-Knifing resampling technique. Only significant metabolomic variables were proposed as the most relevant for clinical interpretation, highlighting their potential significance in relation to the outcome variable.
3. Results
The study cohort consisted of 50 KT recipients with a mean age of 53.98 years, predominantly male (70%) and Caucasian (94%). Cardiovascular risk factors were prevalent: 82% had hypertension, 70% dyslipidemia, and 28% were classified as obese (Body Mass Index, (BMI)>30). The primary etiologies of CKD were hereditary/familial causes (26%) and unclassified causes (24%).
Regarding renal replacement therapy (RRT) prior to transplantation, 74% were on hemodialysis and 26% were on peritoneal dialysis. The mean duration of RRT was 3.44 years. Most recipients had residual diuresis of less than 500 ml (64%). Sensitization was observed in 18% of recipients, with 6% being highly sensitized (calculated Panel Reactive Antibody (cPRA) >98%). The mean Estimated Post-Transplant Survival (EPTS) score was 41.96%.
Donors had a mean age of 50.58 years, with 70% being female. Most donors were from donation after brain death (DBD) (58%), and 32% were expanded criteria donors (ECD). The mean Kidney Donor Profile Index (KDPI) was 52.66, and the mean cold ischemia time was 17.30 hours. Full demographic and clinical details are summarized in
Table 1.
A semi-targeted metabolomic approach identified 226 polar compounds. The results of the final PLS-DA model are visualized in
Figure 1. The scores plot (
Figure 1A) illustrates the projection of the observations (patients) onto the two components. These two components contributed to the 26.6% of the total X-variance (Component 1: 17.1%; Component 2: 9.5%). The two groups of subjects show robust clustering based on the two creatinine levels, although partially overlapped.
Figure 1B shows the representation of loadings explaining the separation of the two groups of patients. The variables used for the PLS-DA model, along with their chemical formula, spectrometric properties, as well as the standardized coefficients, confidence interval, and associated p-values (Jack-Knifing resampling) for the variable Cr >1.5 mg/dl are shown in
Table S1 (Supplementary Material). Those metabolites showing a significant p-value have been colored in red in
Figure 1B. We found that the following metabolites presented a significant standardized coefficient: 5-Methylcytosine, 3-Indolepropinoic acid, Homocarnosine, Phenylalanine, Choline, Xanthosine, L-Kynurenine, Kynurenic acid, Lipoic Acid, Histidine, Guaiacol and 3-Methylindole.
The model performance metrics after cross-validation are as follows. The model achieved an AUC of 0.958. The overall accuracy was 0.875, with a sensitivity of 0.8636 and a specificity of 0.8846. The model correctly identified 19 true positives and 23 true negatives, while misclassifying 3 false positives and 3 false negatives. The ROC curve, depicted in
Figure 1C, further illustrates the model’s performance, with the curve approaching the top left corner.
4. Discussion
The evolution of kidney graft function in the initial months post-transplant correlates with long-term outcomes. Delayed graft function (DGF) and acute rejection are among the most significant, though not isolated, determinants of renal disease progression [
6]. Monitoring involves serial assessments of renal function, urine sediment analysis, proteinuria, antibody measurements, and microbiological tests, with renal biopsy being the most accurate but invasive method [
7], leading to a growing interest in identifying non-invasive biomarkers that can predict KT outcomes earlier and more reliably.
Omics technologies, particularly metabolomics, have emerged as promising tools for biomarker discover. In KT, omics studies have primarily aimed to diagnose acute rejection [
8,
9,
10]. However, few studies have examined other events like DGF [
11] or overall progression [
3]. A notable study by Ihsan Yozgat et al. [
2] used untargeted metabolomics of urine samples to develop a PLS-DA model, correlating a panel of marker candidates with serum creatinine levels above 1.5 mg/dl one-year post-transplant, achieving an AUC of 0.9 for the one-week post-transplant model. This study highlighted the use of a metabolite panel for better discrimination but noted the challenge of collecting urine samples from DGF patients.
Our study used semi-targeted plasma metabolomics at one-week post-transplant, employing a precise in-house compound library. Employing UPLC system coupled to a Q-Exactive Plus spectrometer and a PLS-DA model for statistical analysis, we achieved excellent performance metrics with an AUC of 0.958 following cross-validation.
We identified five statistically significant metabolites with negative standardized coefficients, indicating their association with better graft function and suggesting a protective role. These metabolites are 3-methylindole (skatole), guaiacol, histidine, 3-indolepropionic acid, and α-lipoic acid.
3-methylindole, or skatole, is a metabolite derived from tryptophan metabolism by the indole pyruvate pathway facilitated by gut microbiota through the enzyme tryptophan dehydrogenase (TrpD) [
12] (
Figure 2). Sin-Hyoung Hong et al. demonstrated its protective effects by reducing endoplasmic reticulum stress, oxidative stress, lipogenesis, caspase activity, and apoptosis in hepatocytes [
13]. Although its direct effect on the kidney is not established, a review suggests potential benefits on CKD through intestinal homeostasis, due to the link between dysbiosis and CKD [
14]. However, high concentrations of skatole are associated with colorectal cancer and pneumotoxicity [
15], indicating the need for further investigation.
Guaiacol, or 2-methoxyphenol, is a phenolic compound known for its antioxidant properties, as it can donate electrons and hydrogen atoms to neutralize free radicals and prevent lipid peroxidation [
16]. Its antioxidant activity includes reversible inhibition of myeloperoxidase, an enzyme that induces oxidative stress and inflammation, leading to cardiovascular disease complications [
17,
18]. Although guaiacol’s relationship with renal disease is not well-described, its ability to reduce oxidative stress and inflammation may protect kidney graft function, as suggested by our study.
Histidine is an amino acid with an imidazole side chain that provides proton buffering, metal ion chelation, and antioxidant properties. It scavenges reactive oxygen species (ROS) and reactive nitrogen species (RNS) and sequesters advanced glycation end-products and advanced lipoxidation end-products [
19]. Lower levels of histidine are observed in CKD patients [
20], and it is used to treat CKD-related anemia [
21]. Its protective properties have led to its inclusion in organ preservation solutions like Custodiol
®, which is used in our patients.
3-Indolepropionic acid (IPA) is a tryptophan metabolite produced by gut microbiota through the enzyme tryptophanase (TNA) [
12] (
Figure 2). IPA exhibits potent antioxidant activities, scavenging hydroxyl radicals and improving gut barrier function [
22]. It downregulates fibrosis-related genes and inflammatory factors, offering significant renal protection. Additionally, IPA competes with toxic substances like indoxyl sulfate for organic anion transporters, reducing their accumulation in renal tubules and mitigating renal injury [
23]. Finally, IPA is recognized as a potential biomarker for protection against CKD [
24], and our study associates it with the preventing adverse renal graft outcomes.
α-Lipoic acid is a dithiol compound acting as a cofactor for several enzymes serving as an energy modulator. It exchanges thiol groups with other thiol-containing metabolites, making it a redox modulator and a potent antioxidant by reducing ROS and RNS, chelating metals, and regenerating antioxidants like vitamins C and E. It also has anti-inflammatory activities, targeting NF-κB and decreasing inflammatory cytokines [
25]. Its protective effects are documented in various renal diseases, including diabetic nephropathy, sepsis-associated acute kidney failure, nephrotoxicity-induced renal failure, and ischemia-reperfusion (I/R) injury [
26]. In I/R injury, α-lipoic acid combats oxidative stress, downregulates Na-K-ATPase and nitric oxide synthase, mitigates neutrophil infiltration, inhibits inflammation, and suppresses endothelin-1 upregulation [
27]. It is used as a therapeutic strategy in various solid organ transplants [
25] and could serve as a biomarker for protecting against adverse renal graft outcomes.
Conversely, our study revealed seven statistically significant metabolites with positive standardized coefficients, indicating their association with worse kidney graft outcomes. These metabolites are homocarnosine, 5-methylcytosine, xanthosine, choline, phenylalanine, kynurenic acid and L-kynurenine.
Homocarnosine, a dipeptide of gamma-aminobutyric acid and histidine, has been associated with various biological and pathological brain events. Its role in kidney diseases is less studied, though increased levels are observed in diabetic kidney disease compared to those with type 2 diabetes mellitus without nephropathy [
28], and in renal tissue after ischemia [
29], which may be related to the I/R injury observed in our KT recipients.
5-Methylcytosine (5-mC) is a methylated form of cytosine where a methyl group is attached to carbon at position 5. DNA methylation has been extensively studied in the context of KT and involves adding methyl groups to cytosines in CpG dinucleotides, which cluster in CpG islands, resulting in transcriptional silencing of affected genes [
30]. Although methylation is relatively stable, it can be reversed by TET enzymes, converting it to 5-hydroxymethylcytosine in an oxygen-dependent manner. Hypoxia during KT reduces TET activity, leading to 5-mC accumulation and downregulation of genes involved in suppressing kidney fibrosis and injury [
31]. This hypermethylation can predict chronic allograft injury [
31], as observed in our study, where elevated 5-mC levels predicted worse outcomes at six months.
Xanthosine is a purine nucleoside identified as a prognostic marker for CKD in several studies, often used alongside other metabolites, such as kynurenic acid, kynurenine and choline, which was also found in our study, to estimate CKD progression [
32,
33,
34]. The mechanism by which xanthosine increases CKD risk and its relationship with KT remains unclear. However, purine metabolism dysregulation in CKD, nucleotide breakdown [
35] in I/R injury, and the association of other purine metabolites like hypoxanthine as markers of ischemia [
36] could be related. An in vitro study showed that hypoxia significantly decreased purine nucleotide phosphorylase activity [
37], the enzyme required for converting xanthosine to xanthine, possibly explaining xanthosine accumulation in these patients.
Choline, an essential vitamin found in phospholipids, is involved in methylation reactions like converting homocysteine to methionine via its oxidation to betaine. It is linked to CKD progression through its conversion to trimethylamine-N-oxide, promoting renal fibrosis [
34,
38]. However, its association with KT appears to be more closely related to I/R injury. Choline is a marker of tissue ischemia linked to phospholipid degradation, exhibiting a biphasic release pattern with a peak immediately following global ischemia and a delayed phase during reperfusion [
39]. Lindeman J et al. observed choline uptake in kidney graft recipients from living or deceased donors without DGF, while those with DGF showed choline release, partially masked by its conversion to betaine, potentially associated with the role of betaine in methylation processes occurring in these patients [
40]. These findings align with our study, as DGF is related to posterior renal function and the deleterious effects of elevated choline levels on long-term renal function.
Phenylalanine is an essential amino acid that serves as a precursor for various compounds. Its primary metabolic pathway is the conversion to tyrosine, which can then follow multiple metabolic routes, including the synthesis of catecholamines, melanin, and energy metabolism through its deamination and conversion to acetoacetate and fumarate, intermediates of the Krebs cycle. The kidney is the major producer of tyrosine in the body, as it converts phenylalanine to tyrosine via phenylalanine hydroxylation. Consequently, impaired renal function leads to the accumulation of phenylalanine and a deficiency of tyrosine [
41], both linked with CKD progression and fibrosis [
42] due to their roles in energy metabolism. In KT, I/R injury decreases phenylalanine hydroxylase activity, leading to phenylalanine buildup [
43] and poorer graft outcomes.
The kynurenine pathway (KP) is the main route for tryptophan catabolism, essential for energy metabolism and immune modulation (
Figure 2). Excessive activation of this pathway is linked to various kidney diseases, including AKI, CKD, KT progression, I/R injury, infections in KT, and graft rejection. Elevation of these metabolites is linked to inflammation, with pro-inflammatory cytokines like IFN-γ upregulating indoleamine 2,3-dioxygenase (IDO), the rate-limiting enzyme of KP, and downregulate enzymes such as kynureninase and quinolinate phosphoribosyltransferase, leading to the accumulation of KP intermediates and disrupting de novo NAD+ biosynthesis [
44,
45]. Elevated KP metabolites, such as kynurenine and kynurenic acid (KYNA) correlate with poor renal outcomes and renal fibrosis [
34,
46]. Notably, kynurenine had the highest score as a biomarker for CKD in a large population analysis [
47].
In KT, increased tryptophan degradation and higher kynurenine levels correlate with worse kidney function [
48], can predict graft rejection [
49], and are associated with infections due to the inflammation present in both conditions, raising questions about their specificity as rejection markers. A study by Dharnidharka et al. [
50] found that the serum kynurenine/tryptophan (KYN/TRP) ratio was higher in acute graft rejection. However, peripheral blood CD4-ATP levels were useful in the differential diagnosis, being lower in the infection group. Our study supports the importance of monitoring KP metabolites like Kynurenine and KYNA as potential prognostic biomarkers and therapeutic targets in CKD and KT patients.
This study has several limitations. Firstly, the small sample size, due to the high cost and technical complexity of metabolomic studies, limits the ability to perform subgroup analyses, such as distinguishing between recipients of brain-dead versus cardiac-death donors, with the aim of not losing statistical power. Secondly, metabolomic studies are inherently exploratory and hypothesis-generating, serving as a foundation for future, more comprehensive research that requires validation in larger cohorts. Thirdly, although the long-term evolution of the renal graft has been correlated in other studies with renal function at six months, the short follow-up time in our study does not allow us to ensure this long-term evolution, so further investigation in this aspect would be interesting. Lastly, the use of a semi-targeted approach with an in-house library of compounds to limit irrelevant results and improve their characterization may have limited the discovery of other unknown metabolites.