The metabolomic differential plasma profile between dialysates. Pursuing to understand the mechanisms of citrate dialysate clinical benefits

1 Introduction

Metabolomics is the comprehensive study of endogenous small mass molecules (i.e., less than 1.5 kDa) referred to as metabolites, such as lipids, amino acids, vitamins, steroids, sugars, nucleotides, fatty acids, and organic acids. Unlike genomics, metabolites are dynamic, as extrinsic factors constantly modify them. Their study potentially provides a deeper insight into the human body’s response to various situations and delivers the information necessary to analyze their pathophysiological relevance (Roca et al., 2021). In the field of chronic kidney disease (CKD), it has been used to discover new uremic toxins, biomarkers that improve the stratification of the disease, and its interaction with other diseases such as diabetes, inflammatory conditions, stress, energy metabolism, kidney disease progression, or renal cancer, among others (Sharma et al., 2013; Duranton et al., 2014; Hallan et al., 2017; Hocher and Adamski, 2017; Oto et al., 2020; Zhu et al., 2021).

Moreover, with these methods and specifically in dialysis-dependent CKD (DD-CKD) patients, several metabolites have been correlated with cardiovascular disease and death, inflammatory parameters, nutritional status, cognitive function, hypoxia and oxidative stress, body mass wasting, uremic pruritus, and sleep cycle disturbance (Kurella Tamura et al., 2016; Hu et al., 2019; Kalantari and Nafar, 2019; Bolanos et al., 2021). It has also been used to compare dialysis techniques, but no study on dialysate effects has been performed (Sato et al., 2011). In fact, the dialysate composition is a relevant element of the treatment that has been understudied (Basile and Lomonte, 2015).

Currently, bicarbonate-based dialysate needs a buffer to prevent precipitation of bicarbonate salts with the bivalent cations, and acetate, at doses of 3–4 mmol/L, is the most used (Petitclerc et al., 2011). However, its use is associated with hemodynamic instability, increased oxidative stress, and an inflammatory cardiac profile (Gabutti et al., 2009; Grundström et al., 2013; Pérez-García et al., 2017; Dellepiane et al., 2019; Broseta et al., 2021). Therefore, several compounds have been tested to serve as alternatives to acetate. From them, citrate is being postulated as the preferred option because of its association with an improved hemodynamic tolerance, a reduction in cardiac inflammatory biomarkers, a better acid-base status, a reduction in vascular calcification, and the dose of anticoagulants needed, in most cases by unclear mechanisms (Gabutti et al., 2009; Kossmann et al., 2009; Daimon et al., 2011; Bryland et al., 2012; Kuragano et al., 2012; Rocha et al., 2014; de Sequera et al., 2015; Pérez-García et al., 2017; Lorenz et al., 2018; Dellepiane et al., 2019; Sequera et al., 2019; Trakarnvanich et al., 2019; Villa-Bellosta et al., 2019; Broseta et al., 2021). Unfortunately, all these mentioned metabolomic studies were performed in patients using acetate as dialysate buffer, while there is no data on the metabolic profile in patients using citrate or comparing both.

In this paper, we analyze the metabolic profile in plasma from DD-CKD patients treated with either citrate or acetate dialysates to investigate the metabolic variations induced by the hemodialysis treatment itself and the differences between both dialysates aiming not to demonstrate causal relations but to contribute to this lack of information by raising hypotheses on these unknown biochemical mechanisms.

2 Materials and methods2.1 Study design and participants

DD-CKD patients in chronic hemodialysis at Hospital Universitari i Politècnic La Fe in Valencia, Spain, were considered for inclusion in this unicentric, cross-over, prospective study. Inclusion criteria were being over 18 years old, prevalent (i.e., for at least 3 months), with a treatment scheme of 4-h sessions three times per week; while exclusion ones were having been hospital admitted or discharged within the previous month, being on treatment with a low calcium dialysate (1.25 mmol/L) or declining to give written informed consent.

Patients were followed-up for twenty-four dialysis sessions, twelve with each dialysate (i.e., Fresenius® ACF 3A5 acting as the acetate dialysate (AD), that contains 4 mmol/L of acetate, and Fresenius® SmartBag CA 211.5 as the citrate one (CD), with 1 mmol/L of citrate (Table 1)). Dialysis parameters (blood and dialysate flows and treatment time), medical treatment received, and dialysates’ components, other than the acidifier, remained unchanged during the study to avoid potential confounders. More details on the characteristics of the parameters of the dialysis treatments can be found in a previously published work (Broseta et al., 2021).

TABLE 1. Dialysate characteristics and compounds.

2.2 Blood samples collection and preparation

In the twelfth session, which was always a midweek one, patients’ blood samples were taken at four different moments: pre-dialysis; 60 min after the start of the session (mid-dialysis), the time at which metabolism would be saturated, as based on previous studies Hernández-Jaras et al., 1994; post-dialysis; and 30-min post-dialysis (rebound) by analogy to the urea rebound. They were drawn in ethylenediaminetetraacetic acid (EDTA) tubes and processed within 30 min of collection to avoid platelet activation, protein production, and the decomposition of thermolabile compounds.

Subsequently, a protein precipitation step is made by adding to the 50 µL of plasma samples a cold solvent of 150 µL of acetonitrile (ACN) with 0.1% formic acid (FA) and vortexed for 30 min at -20°C. Afterward, it was centrifuged at 13,000 g for 10 min at a temperature of 4°C to separate the cellular fraction from the plasma. The upper phase, corresponding to the latter, was recovered and aliquoted into Eppendorf’s tubes for subsequent freezing and storing at -80°C.

20 μL of this extract were transferred to a 96 well-plate for LC-MS analysis. Plasma was then diluted by adding 70 μL of H2O with 0.1% v/v FA and 10 μL of an internal standard solution containing phenylalanine-d5, caffeine-d9, leukine enkephaline, and reserpine in H2O:CH3OH (1:1, 0.1% v/v HCOOH) at 20 μM. Once the plate was ready, 10 μL of each sample were collected and prepared for quality control (QC). On its part, blank samples were prepared by replacing plasma with ultrapure H2O in the same extraction tube and following the same sample preparation process as the real plasma ones.

To avoid intra-batch variability and enhance quality and reproducibility, the scheme analysis of samples was performed by random injection order, and at least 5 QC were analyzed at the beginning of the sequence to condition the column and equipment, and at every 5–7 plasma samples to monitor and correct changes in the instrument response as well for filtration and identification purposes. Blank analysis was performed at the end of the sequence and used to identify artifacts from sampling.

2.3 Untargeted Metabolomics Based on UPLC-Q-ToF Mass Spectrometry.

The metabolomic analysis of the processed samples was carried out in the Analytical Unit of the Medical Research Institute Hospital La Fe by an Ultra-Performance Liquid Chromatography (UPLC) system coupled to an iFunnel Q-ToF Agilent 6,550 mass spectrometer (Agilent Technologies, CA, United States). Reversed-phase chromatographic separation was performed using a UPLC BEH C18 column (100 × 2.1 mm, 1.7 μm, Waters, Wexford, Ireland). Autosampler and column temperatures were set to 4 and 40°C, respectively, and the injection volume was 5 μL. Mobile phase A and Mobile phase B consisted of H2O and ACN, both containing 0.1% of FA. The gradient elution was 14 min at a flow rate of 400 μL/min. The mobile phase A (H2O 0.1% v/v HCOOH) was maintained at 98% for 1 min then decreased to 75% in 2 min, 50% in 3 min, and 5% in 3 more min. 95% of mobile phase B (CH3CN, 0.1% v/v HCOOH) was held for 3 min, and then a 0.55 min gradient was used to return to the initial conditions, which were held for 2.5 min for a total column recovery. Full scan MS data from 100 to 1,700 Da were collected in positive (ESI +) electrospray ionization mode. All reagents and chemicals were purchased from Sigma Aldrich (St. Louis, MO, United States). An in-depth review of the sample preparation methodology used by our group has already been published (Roca et al., 2021).

2.4 Data pre-processing

Data processing was done using an in-house R (v.3.6.1) processing script with XCMS and CAMERA packages for peak detection, noise filtering, alignment and normalization. Parameters selected were: CentWave method (noise = 1,000; ppm = 50; k = 5; I = 300; snthresh = 10; min = 6; max = 25; integrate = 2; nS = 4); matching peaks across samples (grouping): Density method (mzwid = 0.05; bw = 1; mF = 0.5; mS = 1); retention time correction (m = 0; e = 0); integration of samples: FillPeaks method; median fold change normalisation: medFC technique; and CAMERA isotopes: FindIsotopes.

The data matrix finally obtained was composed of molecular features consisting of two values: accurate mass (m/z) and retention time (min). To avoid bias, all the samples were processed simultaneously and analyzed in the same batch; therefore, no inter-batch correction was needed. Samples were randomly analyzed to avoid possible intra-batch variance, and internal standards were used to calculate possible variations or drifts during the analysis. Also, a locally estimated scatterplot smoothing (LOESS) normalization with QC samples was performed to normalize our data. The data were filtered according to the quality assurance criteria of coefficient of variation <30% in QC samples and if the percentage of 0 was greater than 60%.

2.5 Metabolomics statistical analysis

Firstly, a pre-selection of significant molecular features between groups was selected by a Volcano Plot carried out using an in-house script in the R platform, combining a Fold Change (FC) method with the significance of a paired Student’s t-test for normally distributed variables or Wilcoxon signed-rank test for skewed data after performing Shapiro-Wilks test for normality. From this, molecular features with a stronger combination of FC (|log2⁡FC|>1) and statistical significance (p-value < 0.05) in each comparison were selected after a false discovery rate adjustment by a Benjamini–Hochberg procedure.

A multivariate analysis was then carried out with the significant features selected in this previous analysis, by using the OMICS skin of the SIMCA software (Sartorious Stedim Biotech, Aubagne, France). Firstly, an exploratory unsupervised principal component analysis (PCA) was performed to extract as much information as possible and to try to recognize patterns of behavior, simplifying the variability of the data by looking for how they were distributed according to their similarities or differences and grouping them into principal components. In a second step, a supervised orthogonal projection to latent structures discriminant analysis (OPLS-DA) was used to determine the main discriminant variables responsible for the differences between groups. The validity and robustness of the models were evaluated by R2(Y) (goodness-of-fit) together with Q2(Y) (goodness-of-prediction), considering a Q2(Y) prediction ability higher than 0,5 as the acceptability threshold, a p-value of the analysis of variance testing of cross-validated predictive residuals (CV-ANOVA) analysis inferior to 0,05 as significant, and a 1,000-iterations permutation test with the new diagnoses R2 and Q2 significantly lower than the real ones.

The selected variables were obtained from the first thirty variables ordered by the variance in importance in projection values (VIP) with a score greater than 1, a jack-knife confidence interval that did not include zero, and a FC greater than 1.2. Each feature was also verified by extracting each ion chromatogram (m/z) in some QC raw data and checking for peak shape and retention time. For each comparison, the PCA score plot, OPLS-DA score plot, and the list of molecular features finally selected are presented.

2.6 Potential metabolites annotation

Once the discriminant variables were selected, tentative annotation of the metabolites was made using the METLIN database (https://metlin.scripps.edu/) and Human Metabolome Database (HMDB) (http://www.hmdb.ca/) by querying their m/z within a range of ±10 ppm. The following adducts were included: M + H and M + Na and M + H-H2O. A verification of the fragmentation spectra in MS/MS mode of the metabolites annotated in the previous step was performed by comparing each experimentally MS/MS spectra obtained from both data dependent (DDA) and data independent analysis (DIA) carried out in some QC, by comparing their MS/MS spectra with those presented in databases. Annotation of metabolites has been made according to the Chemical Analysis Working Group of the Metabolomics Standards Initiative (Sumner et al., 2007). Where identified metabolites (level 1) are those confirmed based on the agreement of their accurate mass (m/z), retention time, and MS/MS spectral with commercially available chemical standards. If unavailable, metabolites are putatively annotated (level 2) when their m/z and MS/MS spectra match with HMDB and Metlin databases or putatively characterized (level 3) if only their m/z coincide with those databases. Molecular features not annotated (level 4) represent less reliable annotation classifications.

3 Results3.1 Participants

Twenty-one included patients with a mean age of 62.25 ± 13.77 years (range 33.05–82.26) years, a body mass index of 25.19 ± 5.52 kg/m2, of which eleven (52%) were females, were included. Two (9.5%) were on high-flux hemodialysis, and nineteen (90.5%) were on post-dilution online hemodiafiltration, with a mean dialysis vintage of 151.48 ± 308.43 months. Twelve (57.1%) patients were dialyzed through an arteriovenous fistula and nine (42.9%) through a tunneled catheter. There were no statistically significant differences between dialysate groups in terms of Kt/V (2.06 ± 0.42 with AD vs. 2.15 ± 0.35 with CD), dialyzed blood volume (84.25 ± 13.79 L with AD vs. 86.41 ± 9.89 L with CD), or substitution volume (21.12 ± 7.95 L with AD vs. 20.9 ± 7.71 with CD). Individual measured values and clinical characteristics of each patient can be seen in Table 2.

TABLE 2. Clinical characteristics and dialysis parameters of the included patients.

3.2 Variations between measurement times

In a preliminary analysis, the pre-dialysis samples were found to be statistically significantly different compared to the mid-dialysis, post-dialysis, and rebound ones; however, there were no differences in the pairwise comparison of the mid-dialysis, post-dialysis, and rebound samples. This finding happened both when using AD and CD and reflects that 60 min after the beginning of the session, the metabolism has achieved its saturation and remains stable during the dialysis session, as well as that these changes in the metabolism are kept for 30 min after its termination, not arriving to the equilibrium. Thus, we only further analyzed the pre- and post-dialysis measurements considering that these would be more reliable than those of mid-dialysis or the rebound.

The Volcano Plot showed 43 significantly different variables between patients’ pre- and post-dialysis blood samples after dialysis performed with the AD (Figure 1A). With these 43 variables, PCA (Figure 1B) and OPLS-DA score plots (Figure 1C) were constructed, showing that the model discriminates between pre- and post-dialysis times. The model diagnostics were adequate, and the model was successfully validated. As we can observe in Figure 1C, the pre-dialysis samples showed quite a lot of variability between them, while at post-dialysis time, the variability decreased and behaved similarly. From this model, the most important discriminant variables were selected according to their VIP score (>1), a jack-knife confidence interval that did not include zero, and a FC greater than 1.2. A search in the databases was performed for annotation of those features as potential metabolites. The annotated variables are shown in Table 3 if they were found in the comparisons between times with either dialysate or in Table 4 if they were only found with AD or CD.

FIGURE 1. Metabolomic variations between pre- and post-dialysis measurements when acetate dialysate is used. (A) Volcano plot, (B) PCA score plot, (C) OPLS-DA score plot.

TABLE 3. Identification of the most discriminant metabolomic variables between pre- and post-dialysis with both dialysates.

TABLE 4. Identification of the most discriminant metabolomic variables between pre- and post-dialysis with differences between dialysates.

Similarly, the Volcano Plot (Figure 2A) comparing the pre- and post-dialysis blood samples using CD showed 42 significantly different variables. PCA (Figure 2B) and OPLS-DA (Figure 2C) score plots were also elaborated, and the model was validated. The pre-dialysis samples also had increased variability compared to those of the post-dialysis time. The most important discriminant variables were selected according to the previous procedure, and those annotated metabolites are described in Tables 3, 4, as previously mentioned.

FIGURE 2. Metabolomic variations between pre- and post-dialysis measurements when citrate dialysate is used. (A) Volcano plot, (B) PCA score plot, (C) OPLS-DA score plot.

3.3 Variations between dialysates

No metabolomic differences were found between pre-dialysis samples with AD vs. CD. On its part, when comparing post-dialysis and rebound measurements, only one variable was selected in the Volcano Plot (Figure 3), being the same in both timings; thus, a multivariate analysis was not performed. This variable and its identification are shown in Table 5. According to their potential clinical relevance discussed below, the relative intensities of the selected metabolic variables have been represented as violin plots (Figure 4).

FIGURE 3. Volcano plot of the metabolomic post-dialysis variations between acetate and citrate dialysates.

TABLE 5. Identification of the most discriminant metabolomic variables between AD and CD in post-dialysis and rebound measurements.

FIGURE 4. Violin plots of selected relative intensities according to clinical interest. *Statistically significant difference.

4 Discussion

CKD has an important impact on metabolism. Thus, metabolomics has been instrumental in identifying new biomarkers that can enhance our understanding of the mechanisms and pathways that underlie renal diseases and improve diagnosis and prognosis estimations (Hocher and Adamski, 2017). In that sense, hemodialysis, the most used kidney replacement therapy though it does not replace all healthy kidney functions, is associated with perturbations in plasma metabolic profiles and would add to the alterations already found in non-dialysis-dependent chronic kidney disease patients in comparison to healthy controls.

This hypothesis-generating untargeted metabolomics study contributed to the identification of different metabolites whose measured plasma concentrations were affected by the diffusive clearance observed from the hemodialysis sessions per se but also successfully identified some others whose concentrations were affected by the dialysate used. Concretely, 31 metabolic variables were identified, grouped by class, and interpreted within a potentially clinical real-life context with the aim to deepen the beneficial mechanistic effects shown by citrate versus acetate in clinical or phenotypical terms. But its results can also be further studied to identify novel uremic toxins, alternatives to improve the urea kinetics model to approach dialysis dosing, and the role of CKD in the high cardiovascular mortality, infectious diseases, and impaired cognitive function, among others.

Notably, only one metabolic variable was found to be significantly different between dialysates at post-dialysis and rebound moments, being 8 to 10 times higher when CD was used as dialysate instead of AD. These metabolic variables were putatively characterized as branched-chain amino acids (BCAAs) (i.e., D-leucine, L-leucine, L-alloisoleucine, L-isoleucine); compounds that are proteinogenic amino acids and have different metabolic routes, with leucine being ketogenic and isoleucine being both a glucogenic and ketogenic amino acid. BCAAs also play a key role in stress states, energy, and muscle metabolism; in fact, isoleucine deficiency is characterized by muscle tremors (Wilkinson et al., 2013; Duan et al., 2016; Plotkin et al., 2021). But its accumulation in blood and other body fluids can exert neurotoxic effects (Yudkoff et al., 2005; Murín and Hamprecht, 2008).

Isoleucine is synthesized from pyruvate and α-ketobutyrate, compounds concentrations that may be affected by the gain of acetate and citrate during dialysis (Hernández-Jaras et al., 1994), and catabolized to α-ketoglutarate which is oxidized and split into propionyl-CoA, which is converted into succinyl-CoA, a tricarboxylic acid (TCA) cycle intermediate which can be converted into oxaloacetate for gluconeogenesis (hence glucogenic); and acetyl-CoA, which can be fed into the TCA cycle or used in the synthesis of ketone bodies or fatty acids (Adeva-Andany et al., 2017; Rajendram et al., 2018).

On its part, leucine particularly stimulates protein synthesis, increases the reutilization of amino acids in many organs, and reduces protein breakdown; this is promoted because leucine also induces insulin release (Stipanuk, 2007; Yang et al., 2010; Yang et al., 2012; Duan et al., 2016). However, like other BCAAs, this effect is associated with insulin resistance in the long term (Yang et al., 2010; Lynch and Adams, 2014; Bloomgarden, 2018). Furthermore, given its ketogenic properties, leucine is an important source of calories and could be considered an even more important energy source than glucose (Newsholme et al., 2011). Moreover, leucine is also a major component of the subunits in ferritin and other ‘buffer’ proteins (Khan, 2018), and it is required in stress states such as surgery, trauma, cirrhosis, infections, or fever (Holeček, 2018), not only as a great energy source but also given its potential capacity to attenuate inflammatory responses (Kato et al., 2016; Xia et al., 2017; An et al., 2020). CKD and hemodialysis could also be considered stress states, although leucine’s role in these specific scenarios is yet to be studied.

The dialysis process affected the relative intensity of 17 metabolic variables, as significant reductions were identified when both dialysates were used. Among them, some uremic toxins were identified according to the European Uremic Toxin Working Group (Duranton et al., 2012). Phenols were putatively characterized, from which p-cresol, a metabolite of aromatic amino acid metabolism produced by intestinal microbiota (mainly enterobacteria), seems to be the more plausible isomer form. At concentrations encountered during uremia, p-cresol inhibits phagocyte function and decreases leukocyte adhesion to cytokine-stimulated endothelial cells (Brunet et al., 2003). It also has been reported to have several toxic effects (e.g., diminishing the oxygen uptake of rat cerebral cortex slices; increasing the free active drug concentration of warfarin and diazepam; being related to growth retardation in the weanling pig; altering cell membrane permeability, at least in bacteria; to induce LDH leakage from rat liver slices; inducing susceptibility to auditive epileptic crises; blocking cell K+ channels; inhibiting the release of platelet-activating factor in rat peritoneal macrophages; and altering the hepatocyte growth and increase aspartate aminotransferase release (Vanholder et al., 1999)). Although measured as relative intensity reduction, the FC of 0.16 in our results would represent 84% removal which differs from the 30% previously described in high-flux hemodialysis (Vanholder et al., 1999).

Other detected metabolic variables that are also uremic toxins were (Duranton et al., 2012): N1-Methyl-2-pyridone-5-carboxamide, which is a product of nicotinamide-adenine dinucleotide (NAD) degradation that produces inhibition of PARP-1, which in turn leads to failure of DNA repair (Rutkowski et al., 2003; de Lorenzo et al., 2013), and has also been related to colorectal cancers and pellagra (Creeke et al., 2007; Brown et al., 2016); and N-(3-acetamidopropyl)pyrrolidin-2-one, whose levels are increased in non-Hodgkin’s lymphoma (Hessels et al., 1991; van den Berg et al., 1986) and is a catabolic product of spermidine which has been identified as a biomarker of glomerular filtration rate decline. Moreover, creatinine and hydrogen phosphate were also significantly removed during dialysis.

Apart from these molecules, differences between measurements were found with one but not the other dialysate. AD affected plasma levels of three of the measured metabolites and CD, eleven. Two acylcarnitines, 3-Hydroxydecanoyl and 9-Decenoylcarnitine, were putatively characterized as a metabolic variable whose intensities significantly reduced after an AD dialysis. The general role of acylcarnitines is to transport acyl groups (organic acids and fatty acids) from the cytoplasm into the mitochondria so they can be broken down to produce energy. This process is known as beta-oxidation (Dambrova et al., 2022). Regarding their relation to CKD, serum acylcarnitines increase in CKD and DD-CKD patients due to their impaired renal excretion (Fouque et al., 2006) and have been associated with IgA nephropathy and diabetic nephropathy as potential biomarkers (Xia et al., 2019; Mu et al., 2022). In fact, its elevation could indicate mitochondrial dysfunction and seems associated with cardiovascular mortality in incident dialysis (Kalim et al., 2013). Moreover, some of them are also pro-apoptotic (Ferrara et al., 2005). According to a recent review, both 3-Hydroxydecanoyl and 9-Decenoylcarnitine would be classified as medium-chain acylcarnitines, which are somewhat less abundant than short-chain acylcarnitines (Dambrova et al., 2022). They are formed either through esterification with L-carnitine or through the peroxisomal metabolism of longer chain acylcarnitines (Ferdinandusse et al., 1999; Violante et al., 2019), and have been related to inherited disorders of fatty acid metabolism (Dambrova et al., 2022). Particularly, 9-decenoylcarnitine is elevated in the plasma of overweighted subjects (Kang et al., 2018) and decreased in patients with schizophrenia or familial Mediterranean fever (Kiykim et al., 2016; Cao et al., 2020). The study of acylcarnitines is an active area of research, and many novel roles in health and disease will likely be uncovered (Dambrova et al., 2022).

Among the eleven metabolic variables’ intensities that significantly differ when using CD, three uremic toxins were identified (i.e., Trimethylamine N-oxide, Phenol, and Hippuric acid) (Duranton et al., 2012). Trimethylamine N-oxide (TMAO) is an oxidation product of trimethylamine derived from the conversion of dietary intake of lecithin or carnitine by the intestinal microbiota (Tang et al., 2013). It is used by the body as an osmolyte to counteract the effects of increased concentrations of urea (Lin and Timasheff, 1994). Increased TMAO levels are associated with an increase in cholesterol deposition (Koeth et al., 2013) and risk of incident major adverse cardiovascular events (Tang et al., 2013); thus, the higher the clearance, the most beneficial the cardiovascular effect achieved (Mair et al., 2018).

Apart from these uremic toxins, when CD was used, a relative intensity significant decrease of two more metabolites (i.e., 2.3,4,5-Tetrahydro-2-pyridinecarboxylic acid and D-1-Piperideine-2-carboxylic acid) playing a key role in the lysine degradation (Chang and Charles, 1995; Matthews, 2020), which is an essential amino acid that is a necessary building block for proteins, plays a major role in calcium absorption, building muscle protein, recovering from stress conditions, and the production of hormones, enzymes, and antibodies (Singh et al., 2011).

Among the 31 metabolic variables, seven of them could not be annotated and remain unknown. One of them, specifically m/z 539.9853, results of great interest as it is the unique one that has a FC > 1 (2.87), reflecting that its production during dialysis with CD surpasses its clearance in the dialysate.

Moreover, some metabolic variables annotated have not been further discussed since we wanted to focus on those with a potential clinical implication or prognosis in CKD, the dialysis treatment, or the effects of acetate or citrate. Those are compounds that act as dietary components (i.e., piperidines, 2,5-Dihydro-2,4,5-trimethyloxazole, 2-Acetylpyrrolidine, Branched alkanes, 1-Amino-propan-2-ol, 1-Methyl-2-pyrrolecarboxaldehyde, 2-Acetylpyrrole, 4-Aminophenol, Hydroxybenzaldehydes, 3-(2-Furanyl)-2-propenal, Propionylcholine, Sesquiterpenoids, xi-1-(Propylthio)-1-propanethiol, Dialkyldisulfides, XanthotoxolHypoxanthines, Aminochrome o-semiquinone, 2-benzoxazolol, L-3-Aminodihydro-2(3H)-furanone, 1,2-Dihydroxy-3-keto-5-methylthiopentene, 1-(Methylsulfanyl)-1-oxopropan-2-yl acetate, Bergaptol, N-(Carbethoxyacetyl)-4-chloro-L-tryptophan, 7-Chloro-6-demethylcepharadione B, N(2)-phenylacetyl-L-glutaminate); found in food of plants origin (i.e. 4-Methyleneproline), or food aditives (i.e., Benzoic acid, Epsilon-caprolactam, 1-(4-Methoxyphenyl)-2-nitroethylene, Methyl n-formylanthranilate, 4-(1.1,3,3-Tetramethylbutyl)-phenol, Etaspirene, 1,6-Hexanedithiol); some drug compounds or their metabolites (i.e., morpholines, aminocaproic acid, Pregabalin, methacholine, arabinofuranosylguanine, 3′,4′-Dihydrodiol); and even other chemicals that represent environmental exposures or exposome (i.e., Benzyl alcohol, S-1-1′Bi-2-naphtol). Of particular interest among these groups are DL-2-Aminooctanoic acid, an amino acid that has been related to colorectal cancer (Brown et al., 2016); Neocaspirene, which is a monoterpenoid, lipidic constituent of the cellular membrane, but also a food additive; and L-isoleucyl-L-proline, which is an incomplete breakdown product of protein digestion or protein catabolism and has been related to asthma (Mattarucchi et al., 2012). These three metabolites were reduced in the dialysis session with either dialysate. On its part, adrenochrome, an oxidation product of adrenaline acquired by ingestion, only seemed statistically significantly different when CD was used.

This study has some limitations. Firstly, the analysis was only made in the positive electrospray ionization mode, which may have led to the non-identification of some metabolites that could have been measured in the negative mode. Secondly, the metabolomic analysis of the effluent was not performed; hence, this study failed to discriminate whether reduced metabolite blood relative intensities were due to their dialysis clearance or consumption in a particular metabolic pathway. In that sense, comparing relative intensity means is less meaningful than differential tendencies between times for each dialysate. Thirdly, our data only determines trends or changes in response intensities but not exact concentrations. Fourthly, if we had identified every annotated metabolite, we would have been able to continue to a validation phase that is missing. And finally, the chosen times for sample extraction were helpful to verify that 60 min after the start of the session, the metabolism had reached saturation and that 30 min after the dialysis session, it still did not attain the baseline status, and therefore, the dialysis effects last longer on the metabolomic profile than the usual urea rebound; but fail to provide data of greater interest, so more measurements would have given more information on the kinetics of these metabolites.

In conclusion, this hypothesis-generating untargeted metabolomic study provides essential data on the metabolic profile of DD-CKD patients, the effect of the dialysis technique on the metabolome, and some potential pathways that differ when CD or AD are used. There are five pillars on which further research may bring light. Firstly, we identified different already known uremic toxins; some of them cleared better when CD was used. Secondly, the CD also has a role in the lysine degradation pathway that was not found with AD. Thirdly, and on the contrary, AD showed an effect on acylcarnitines clearance not shown with CD. Fourthly, BCAAs increased intensities in CD post-dialysis measurement compared to AD ones. And lastly, the identification of a metabolic variable that we could not annotate was the only one whose intensity increased during the dialysis session and only when CD was used.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: and https://github.com/Broseta/Citrate-dialysate,Broseta/Citrate-dialysate.

Ethics statement

The studies involving human participants were reviewed and approved by the Hospital ClÃ-nic of Barcelona Clinical Research Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

Author contributions

JB, MR, and JH-J. conceived the study. RD, SB-G, AS, and PS-P, also contributed substantially to its conception. JB, LL, and AG-B acquired the data from the electronic health records and collected the blood samples. MR. performed the metabolomic analysis. JB. made the statistical analysis. JB, EC, and DR-E. made the figures and drafted and revised the paper. All authors have revised the drafts and approved the final one.

Acknowledgments

We would like to publicly thank all participating patients and the nursing staff of the Dialysis and Investigation Units of the Department of Nephrology of Hospital Universitari i Politècnic La Fe for their selfless and enthusiastic collaboration in this study. Part of the equipment used in this work was co-funded by the Generalitat Valenciana and European Regional Development Fund (FEDER) funds (PO FEDER of Comunitat Valenciana 2014 - 2020)

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Adeva-Andany M. M., López-Maside L., Donapetry-García C., Fernández-Fernández C., Sixto-Leal C. (2017). Enzymes involved in branched-chain amino acid metabolism in humans. Amino Acids 49 (6), 1005–1028. doi:10.1007/S00726-017-2412-7

PubMed Abstract | CrossRef Full Text | Google Scholar

An Y. H., Kim J., Kim H. J., Lim K. (2020). Effects of leucine-enriched essential amino acid supplementation on muscular fatigue and inflammatory cytokines in wheelchair basketball players. Phys. Act. Nutr. 24 (2), 38–46. doi:10.20463/PAN.2020.0013

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolanos C. G., Pham N. M., Mair R. D., Meyer T. W., Sirich T. L. (2021). Metabolomic analysis of uremic pruritus in patients on hemodialysis. PLoS One 16 (2), e0246765. doi:10.1371/JOURNAL.PONE.0246765

PubMed Abstract | CrossRef Full Text | Google Scholar

Broseta J. J., López-Romero L. C., Cerveró A., Devesa-Such R., Soldevila A., Bea-Granell S., et al. (2021). Improvements in inflammation and calcium balance of citrate versus acetate as dialysate buffer in maintenance hemodialysis: A unicentric, cross-over, prospective study. Blood Purif. 50 (6), 914–920. doi:10.1159/000513419

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown D. G., Rao S., Weir T. L., O'Malia J., Bazan M., Brown R. J., et al. (2016). Metabolomics and metabolic pathway networks from human colorectal cancers, adjacent mucosa, and stool. Cancer Metab. 4 (1), 11. doi:10.1186/S40170-016-0151-Y

PubMed Abstract | CrossRef Full Text | Google Scholar

Bryland A., Wieslander A., Carlsson O., Hellmark T., Godaly G. (2012). Citrate treatment reduces endothelial death and inflammation under hyperglycaemic conditions. Diab. Vasc. Dis. Res. 9 (1), 42–51. doi:10.1177/1479164111424297