Abstract
Introduction Limited tools exist for predicting kidney function in long-term kidney transplant recipients (KTRs). Elabela and apelin are APJ receptor agonists that constitute the apelinergic axis, which is a recently discovered system regulating vascular and cardiac tissue, in opposition to renin-angiotensin-aldosterone.
Methods Longitudinal, observational cohort of 102 KTRs who maintained graft function ≥24 months, with no acute rejection history or current active or chronic infection. Serum apelin, elabela, fibroblast growth factor 23 (FGF-23) and α-Klotho were tested using enzyme-linked immunoassay and compared with a control group of 32 healthy volunteers.
Results Median (IQR) follow-up time was 83 (42, 85) months. Higher serum FGF-23 and elabela, but lower α Klotho concentrations were observed in KTRs. Most KTRs had stable trajectories of renal function. All candidate markers were significantly associated with mean two-year eGFR over follow-up, which itself was validated respective to ‘death with functioning graft’ censored dialysis requirement. Using a cross-validation approach, we demonstrated eGFR at initial visit as the most salient predictor of future renal function. Machine learning models incorporating both clinical and biochemical (candidate markers) assessments were estimated to explain 15% of variance in future eGFR when considering eGFR-independent predictions.
Conclusions Utilization of machine learning tools that incorporate clinical information and biochemical assessments, including serum amrkers of the apelinergic axis, may help stratify risk and aid decision making in the care of long term KTRs.
1. Introduction
Chronic kidney disease (CKD) is an irreversible disorder with a slowly progressive course that is tied to excess morbidity. Studies estimate its global prevalence at over 10% within the general population, with non-linear growth and even higher rates among the elderly1. Cardiovascular (CV) disease remains one of the leading causes of death across the whole spectrum of CKD. CV-related mortality is starkly elevated in end-stage CKD, especially in subjects requiring renal replacement therapy (RRT), but also remains significantly elevated after kidney transplant (KTx)2.
Apelin (APLN) and elabela (ELA; toddler or apela) are two endogenous ligands of the APJ, a G protein coupled receptor tied to several intracellular transduction pathways affecting adenylyl cyclase activity, ion gradients (including calcium shift) and nitric oxide synthesis3,4. APJ shows widespread expression in multiple tissues, including cardiac, vascular and renal organs1,5,6. The APJ receptor shows modest homology and comparable tissue distribution to angiotensin receptor type 1, which has led to speculation of physiological role as a counterbalance to the renin-angiotensin-aldosteron (RAAS) axis, alterations of which result in accelerated atherosclerosis and various organ disorders7–10.
APLN and ELA are 77- and 54-amino acid preproproteins cleaved by tissue proteases into shorter peptides, which are secretable11. While the C-terminal APLN sequence determines receptor binding, the N-terminal interacts with the APJ receptor, which may explain variability in tissue affinity of specific isoforms1,12. Most immunological assays target the C-terminal fragment of ELA, which precludes differentiation of specific isoforms. Immunoassay remains the only robust method of testing ELA concentrations13,14. Contrary to prior hypotheses suspecting an ELA-specific receptor, ELA has been demonstrated to act as a ligand of the APJ receptor, with near selective expression within the vascular endothelium of the kidney15–17.
In contrast to the early post-Ktx period, long-term kidney transplant recipients (KTRs) that maintain stable graft function represent a distinct population that is susceptible to CV and metabolic diseases, which represents the major risk factor for graft loss and/or death18. Developing prognostic tools for these patients is of particular interest due to the paucity of measures that are validated specifically in such populations. Due to its strong link with CV disease, the recently discovered apelinergic axis represents an important regulator of tissue homeostasis and potential source of biomarkers.
This study was undertaken to develop a predictive model for future renal function in stable, long-term KTRs. This analysis includes an exploratory investigation of the potential relevance of apelinergic markers in serum, including their association with clinical co-variates and definitive post-KTx outcomes. While models derived based on large population of KTRs have been previously published, the population of stable, long-term KTRs represent a unique group, for which the validity and relevance of published risk scores is uncertain.
2. Methods
2.1 Patients
This was a longitudinal, observational cohort study, which enrolled 102 consecutive patients under ambulatory care at the University Hospital in Kraków, Poland. The recruitment process had an annual timeframe between September 1, 2016 and October 31, 2017, which was pre-defined to last at least one year due to slow rate of recruitment. Subjects were considered for study participation only if they were long-term KTRs, which was defined as having maintained allograft function for at least 24 months. Episodes of acute rejection remain one of the major risk factors for chronic rejection19, therefore patients in whom history of acute rejection was recorded (whether cellular, humoral or vascular), were not included. We further excluded patients with infection, both active (ie signs or symptoms of any acute infection at presentation) or chronic (based on serological evidence in routine work-up, eg viral hepatitis, HIV). Due to altered mineral-bone imbalance, individuals with prior parathyroidectomy status or a history of malignancy (excluding non-melanoma skin cancer), were also excluded due to the presence of (potential) salient interfering factors. The recruitment process is illustrated by a flowchart (Figure S1). The inclusion and exclusion criteria are consistent with our earlier studies.
2.2 Follow-up
Over time, electronic medical care records were screened manually by four physicians (K.B., M.B., A.S., K.N.) to identify visits at relatively equally spaced intervals, with a pre-defined range of 3 to 6 months apart. The total follow-up time for renal function measures was comparable for most patients, with a median (IQR) timeframe of 23.6 (22.8, 24.4) months and range of 21.7-25.3 months. The occurrence of events of poor prognosis, such as allograft loss (defined as permanent dialysis transfer or re-transplantation listing) or death was determined through in person or telephone contact with the patient, family and/or dialysis center. This process was carried out successively until December 30, 2023 (censoring date).
2.3 Covariates
Clinical data, which included age, sex, body mass index (BMI), KTx characteristics (primary etiology of kidney disease, immunosuppresive treatment, history of delayed graft function) and co-morbidity (CV disease, hypertension, diabetes and dyslipidemia) were gathered. CV disease was defined as records of either atherosclerotic CV disease (ie, a composite of coronary artery disease (CAD), prior myocardial infarction (MI) or prior stroke) or heart failure (HF).
Nearly all patients were treated with triple immunosuppressive therapy (glucocorticoids, mycophenolate mofetil and a calcineurin inhibitor) and none had history of induction pre-treatment. Glomerulonephritis, autosomal dominant polycystic kidney disease and reflux nephritis were the most common primary causes of KTx. Information regarding death or dialysis/re-transplantation was determined based on in-person and/or telephone contact with the patient, family and/or dialysis center.
2.4 Biochemical testing
Peripheral venous cannulation was performed at presentation on the morning following an overnight fast into EDTA tubes, which were subsequently frozen at -70 C and stored until analysis. Standard biochemistry assays were performed on either of three automatic analyzers: Hitachi 917 (Hitachi, Japan); Modular P (Roche Diagnostics, Germany); Sysmex XE 2100 (Sysmex, Japan). For non-routine biochemical assays, testing was performed in batch, using commercially available immunoenzymatic kits according to manufacturer instructions. A healthy control group of 32 volunteers (16 female; 16 male), aged between 29 and 74 years (mean 50 years) was recruited using convenience sampling among willing physicians and medical personnel. A more detailed description of this control group is available elsewhere20.
Serum apelin was measured using Human Apelin ELISA Kit (EIAab Science, Wuhan, China), with a detectable level of 62.5 - 4000 pg/mL, sensitivity of 30 pg/mL, and intra- and inter-assay precision of <7% and <9%, respectively.
Serum elabela was measured using Cat. No. S-1508 ELABELA (Peninsula Laboratories, San Carlos, CA, USA), with a detectable level of 0-100 ng/mL, and intra- and inter-assay precision ∼10% and ∼15% (as reported21), respectively. No cross-reactivity with apelin has been reported.
Serum klotho was measured using Human Soluble α Klotho kit (IBL, Gunma, Japan), with a detectable level of 93.75 – 6 000 pg/mL, sensitivity of 6.15 pg/mL, and intra- and inter-assay precision of ∼3% and ∼7%, respectively.
Serum FGF-23 was measured using Human FGF-23 Intact ELISA kit (Immunotopics, San Clemente, USA), with a detectable level of 31.25-2000 pg/mL, sensitivity of 15 pg/ml, and intra- and inter-assay precision of ∼7% and ∼11%, respectively.
2.5 Bioethics statement
The conduct of this study was in adherence with the Declaration of Helsinki principles, as well as International Conference on Harmonization/Good Clinical Practice regulations. Study approval was obtained from the Bioethics Committee of the Jagiellonian University and all patients provided written informed consent prior to study inclusion.
2.6 Statistical analysis
Analyses were performed in R 4.3.2 (R Core Team, 2024, R Foundation for Statistical Computing, Vienna, Austria).
Nominal variables were summarized as counts and proportions (N, %). Variable distribution was assessed using density plots, the Shapiro Wilk test and skewness measures. Continuous variables were summarized as mean (standard deviation) or median (interquartile range), as deemed appropriate. Comparison across groups is performed using Welch’s t test (supplemented with permutation based testing) for continuous variables. The measure of association for two dichotomous variables was assessed with x2 test or Fisher’s exact test.
The nephro package was used to calculate estimated glomerular filtration (eGFR) according to the Chronic Kidney Disease-Epidemiology Collaboration (CKD-EPI) formula based on serum creatinine22. While the median may be conceptually more robust than the mean, due to very high correlation between both measures, the mean was selected.
Prognostic models were constructed using random forest, as a high performance machine learning model, using the ranger package within the tidymodels package. Due to generally low rates of missing data per variable (<10%), we utilized imputation using bagging techniques. Tuning was performed with grid search for hyperparameter selection using five times repeated tenfold cross-validation. Due to sample size considerations, feature selection was determined using an initial model with all potential covariates, further trimmed according to variable importance and expert knowledge. Permutation-based model breakdown techniques using the vip and DALEX package were utilized to analyze feature contribution to model prediction from a global and local perspective, respectively.
Statistical tests were two-tailed and P value < 0.05 was deemed statistically significant.
3. Results
3.1 Baseline clinical characteristics of kidney transplant recipients
This was a cohort of middle-aged adults, with mean (SD) age of 50.8 (14.5) years, and male predominance (N=72, 70.6%). Most patients were overweight, with mean (SD) body mass of 26.1 (4.68) kg/m2. Tacrolimus-based triple therapy regimens were most commonly utilized (N=67, 66.3%), followed by cyclosporine-based (N=22, 21.6%). History of delayed graft function, defined as dialysis requirement within 7 days in the early Tx period, was documented for every fifth patient (N=23, 22.6%). For mortality and kidney allograft loss, the median (IQR) follow-up time was 82.6 (42.1, 85.2) months. Most patients retained adequate kidney function at baseline, with mean eGFR upon presentation at 58.4 (22.6) ml/1.73m2/min, respectively.
Most subjects were hypertensive (N=92, 90.2%), with dyslipidemia and diabetes reported for 15 (14.7%) and 25 (24.5%) individuals, respectively. Nearly all diabetics (N=29, 28.4%) were characterized with post-transplant diabetes (N=24, 23.5%). Manifest CV disease was observed in every fourth patient (N=25, 24.5%). A comparison of clinical features according to allograft status at last available follow-up is shown in Table 1.
When comparing with healthy controls, we observed higher concentrations of serum FGF-23 (mean difference [SE] 1055.39 (488.63), P<0.001) and elabela (mean difference [SE] 723.28 (252.36), P=0.003) in KTRs, alongside singificantly lower α Klotho (mean difference [SE]: -543.20 (146.59), P=0.02) levels, but apelin (mean difference [SE]: -138.26 (295.13), P=0.77) (see Figure 1).
3.2 Relationships between kidney function, vascular and metabolic parameters and circulating concentrations of apelinergic peptides and mineral-bone hormones
This cohort included patients from a broad range of graft function at inception. When examining temporal changes in eGFR, we observed that patients who maintain allograft function have stable trajectories, as compared with subjects who (at a later date) required dialysis therapy and exhibited a gradual trend of eGFR decline (Figure 2, Panel A). In an age-adjusted logistic regression model (Somer’s Dxy 0.76, optimism corrected R2=0.52), higher mean eGFR over two-year follow-up was associated with significantly lower odds of graft loss (OR 0.04, 95% CI 0.01-0.15; P<0.001). The predicted probabilities are illustrated and adjusted for age, according to the interquartile cut-offs (Figure 2, Panel B).
Baseline renal function showed no significant association with apelinergic markers (R=-0.10, P =0.36 for apelin, r=0.15, P=0.13 for elabela), in contrast to mineral bone peptides (R=-0.25, P=0.02 for FGF-23; r=0.34, P<0.001 for α Klotho, respectively). Only a weak and inverse linear trend could be observed between serum apelin and elabela (r=-0.18, P=0.08), as well as between FGF-23 and α Klotho (r=-0.15, P=0.14).
We observed significant relationships between circulating concentrations of apelinergic peptides and mineral-bone hormones with future eGFR over follow-up (Figure 3).
3.3 Developing a robust prognostic model for future kidney function
To understand the importance of baseline eGFR, a univariable random forest regression model was tuned and estimated to explain 87% of variance (R2 [SE]: 0.865 [0.01]) in future 2-year eGFR. This suggests the overarching importance of eGFR assessments over follow-up in stable KTRs. Therefore, mean eGFR over follow-up – baseline eGFR was adopted as a measure of kidney outcome agonistic of eGFR. This approach enables easier interpretation of other, clinically-relevant covariates.
For the final random forest model (agnostic of baseline eGFR), we observed stable performance with ability to explain 15% of variance (R2 [SE]: 0.146 (0.02) in future 2-year eGFR minus baseline eGFR). We also assessed the performance of this model based on aggreagated estimates over resamples, with a mean absolute error of 7.41 (0.194) and a root mean squared error of 9.01 (0.240). Global variable importance is presented in Figure 4 (left-hand panel), while prediction analysis for a sample clinical case illustrates how biochemical markers and clinical characteristics contribute to eGFR prognosis (right-hand panel).
When examining marginal effects of age and body-mass (Figure 5, Panel E & F), we observed a log-like relationship for age and non-monotonic contribution of BMI. It is probable that overall, older age and normal BMI are likely to exert a protective effect on future renal function. For APJ agonists, circulating apelin and elabela deficiency may be associated with worse prognosis for eGFR change, which contrasts with the effect of mineral-bone peptides.
3. Discussion
This was a longitudinal cohort study of long-term, stable KTRs, in whom we observed stable trajectories of graft function over two year follow-up. While post-KTx renal function is a major determinant of graft and patient survival23 and acute rejection remains a critical offender in the early stages, less is known about risk stratification in long-term KTRs. Epidemiological trends suggest an improvement in prognosis over recent years, which has been attributed to modern treatment schemes24. However, the case of stable, long-term KTRs is particularly interesting, as under the assumption of satisfactory HLA matching and compliance with immunosuppressive treatment, a number of different risk factors (including CV and metabolic) for allograft loss emerge in importance19.
There are several theoretical considerations that justify investigation of the apelinergic system as a prognostic measure in KTRs. Elabela and apelin show partially overlapping distribution throughout the vasculature15,25. Their joint receptor activation promotes vasodilation15,26 and increases cardiac contractility15,27, suggesting importance in a compensatory response to developing vascular disease. Meanwhile, within the kidney, APJ expression is more prevalent within glomeruli, indicating its hypothetical utility as a measure of glomerular, rather than tubular injury16,28,29. Kidney injury results in apelin expression30, while elabela mRNA levels are markedly reduced in ischemia-reperfusion models31, suggesting an inverse role in the reparative response. Both apelinergic agonists are attributed cardio-32,33 and renoprotective10,34,35 roles, though the activity of elabela is describes as more potent31. Thus, alterations in elabela, rather than apelin, are expected to parallel impairment within the kidney vasculature.
Among the investigated markers, significantly higher concentrations of elabela, FGF-23 and lower α Klotho concentrations were observed among KTRs, as compared with healthy controls. Both elabela and α Klotho serum concentrations showed a significant relationship with future eGFR, while a similar association for apelin was noted only among individuals with overt CV disease. Using a robust machine-learning technique, we constructed predictive models to assist in understanding the prognostic contribution of these novel markers, which may be comparable to clinically-established covariates. However, the overarching importance of eGFR at prior visits is also evident for patients with long-term, stable kidney allograft.
In patients with kidney disease, plasma apelin concentrations were first described by Malyszko et al.36. Lower concentrations of apelin were tied to CVD, with ventricular dimensions in echocardiography as its major determinants, speculated to parallel to the extent hemodynamic overload36. However, these were prevalent dialysis subjects. Another study in autosomal dominant polycystic kidney disease noted an independent association between apelin and renal function decline, speculating it may parallel dehydration status and suggesting its utility as a prognostic measure, beyond the information derived from eGFR estimates37. In patients with type 2 diabetes, serum apelin concentrations were markedly higher (compared with healthy subjects) and associated with albuminuria, an early indicator of diabetic nephropathy. Experimental evidence suggests apelin may modulate the permeability and proliferation capacity of glomerular endothelial cells, which may affect hyperfiltration of glomeruli38. In contrast, a report by Lu et al. showed no significant differences in apelin concentrations across patients with different stages of CKD when comparing age- and sex-matched controls39. Further work is necessary to clarify the relevance of circulating apelin concentrations.
Prior studies have suggested that elabela may be associated with progression of nephropathy39. Data from two study groups suggests that aside from selective expression within the kidney, it is localized tubular, rather than glomerular compartments16,40. Through its interaction with the Gi signaling pathway, elabela is a likely regulator of fluid homeostasis17, though it may also act to antagonizes aldosterone-stimulation of kidney tissue41. Numerous experimental reports identify elabela expression as an effective countermeasure to renal injury15,40, though whether this is mechanistically linked to cell cycle and/or anti-inflammatory effects6 remains unclear.
Apelinergic activity can also be ascribed a protective role in arterial calcification due to its regulation of smooth muscle cell differentiation into osteoblast lines42 and attenuation of calcification processes43. While mineral-bone peptides (FGF-23 and α Klotho) did not show any apparent association with apelinergic markers, they contributed to model predictions and shared an association with renal function over follow-up. We hypothesize that FGF-23 excess and alpha Klotho deficiency slowly develops in post-Ktx patients, in a manner akin to progressing CKD. Therefore, its supplemental prognostic value in KTx recipients may parallel recognition as risk factor of kidney dysfunction in the general population44.
Traditional CV risk factors, with underlying metabolic disorders and excessive RAAS activation, are processes that promote allograft dysfunction45. Apelin shares a close relationship with TNF α expression in adipose tissue, which has lead to its consideration as a potential bridge mediator between inflammation and insulin resistance in obesity7,46, though body-mass is likely not a major determinant of its circulating forms47. In patients with type 2 diabetes, serum elabela levels are inversely tied to proteinuria and creatinine elevation48, while experimental evidence suggests a protective role of elabela that reduces podocyte apoptosis49. Antagonism between the RAAS and APJ axes has been previously described50,51 and a simplified schematic is provided for the reader in Supplementary Figure S2. Importantly, on a systemic level, apelinergic activity may only be marked under conditions of tissue stress52,53, as demonstrated by experimental studies of myocardial contractility in failing hearts54,55. This may explain, at least in part, the conflicting findings of patient-level studies.
In model construct we utilized novel, potential markers of recently discovered biologic processes that can shape cardiorenal disease. Interpretation of circulating concentrations and causal cross-attribution, respective to their purported biologic roles, remains tentative. The apelinergic system remains incompletely understood, particularly in the post-KTx setting. We assessed both apelin and elabela in human serum using enzyme-linked immunosorbent assays, as has been performed by others. The potential variability in bioactivity and tissue affinity of specific molecular isoforms, as well as different sources of circulating peptides requires clarification in future studies. Interactions within signaling pathways tied to vascular disease (e.g., shear stress, endothelial dysfunction) or immune dysregulation may further occlude our conceptual understanding of changes in serum/plasma concentrations of specific markers. Thus, the results of the present study may be influenced by interaction with unknown clinical characteristics or incomplete understanding of biologic activity.
Many studies are limited by utilization of binary outcomes derived from repeated eGFR measures. Inter- and intra-individual variation in eGFR assessments, as well as timepoint number and temporal relation should be considered. To improve clinical interpretation, we defined the response variables as mean sequential eGFR and last available eGFR. For both outcomes, a strong association with ‘death with functioning graft’ censored dialysis was verified. To limit bias of the model predicting last available eGFR, we also utilized internal validation based on resampling techniques.
4. Conclusions
Discovering new biomarkers with prognostic potential and developing a risk model for renal function loss in KTRs holds high clinical importance, as no single parameter is sufficient to identify patients at greatest risk for allograft loss, and routine biochemical and clinical tests are often insufficient to predict graft failure accurately. The developing promise for apelinergic therapeutics in cardiac and renal diseases suggests another, supplemental role for apelin and elabela, as candidate biomarkers of a multisystem axis that counteract the effects of RAAS1,6,10.
Following the early post-KTx stage, allograft loss is difficult to predict and complex models involving demographic, clinical and biochemical data are warranted. This was an observational cohort study of long-term, stable KTRs, which examined temporal trajectories of renal function based on sequential eGFR assessments. Candidate biomarkers tied to apelinergic activity and mineral-bone imbalance were evaluated, with theoretical justification derived from modulatory activity within vasculature and renal tissue. Circulating concentrations of serum elabela, apelin, FGF-23 and α Klotho concentrations were significantly associated with mean eGFR over two year follow-up. After accounting for kidney function at initial visit, a random forest model combining basic clinical information and biomarker levels was estimated to account for 14% of variance in future eGFR.
Conflict of interest statement
The authors declare that there is no conflict of interest regarding the publication of this paper.
Data Availability
All relevant data are within the manuscript and its Supporting Information files.
Funding sources
This study was supported by a statutory grant from the Jagiellonian University Medical College (N41/DBS/000193; to K.K.). This study was performed using equipment of the Medical University of Bialystok acquired via the RPOWP 2007-2013 funding grant, Priority I, Axis 1.1, contract No. UDA-RPPD.01.01.00-20-001/15-00 dated to 26.06.2015.
Author contributions
All authors contributed significantly to the present report. All authors performed critical revision and provided acceptance of this paper. Conceptualization, K.B., K.K.; Data curation, K.B., K.N., K.S., A.S., M.B., J.M., K.K., M.K., E.K.-Ż., M.Ż., A.B-P.; Formal analysis, K.B., K.N., K.S., A.S., M.B., J.M., K.K., M.K., E.K.-Ż., M.Ż., A.B-P.; Funding acquisition, K.K. Investigation, Methodology, K.B., K.N., K.S., A.S., M.B., J.M., K.K., M.K., E.K.-Ż., M.Ż., A.B-P.; Software, K.B.; Supervision, J.M., M.K., K.K.; Writing—original draft, review and editing, K.B., K.N., K.S., A.S., M.B., J.M., K.K., M.K., E.K.-Ż., M.Ż., A.B-P. All authors have read and agreed to the published version of the manuscript.
Abbreviations
- APJ
- apelin/elabela receptor
- BMI
- Body mass index
- CV
- Cardiovascular
- CKD
- Chronic kidney disease
- CKD-EPI
- Chronic Kidney Disease-Epidemiology Collaboration
- eGFR
- Estimated glomerular filtration
- HF
- Heart failure
- KTx
- Kidney Transplant
- KTRs
- Kidney Transplant Recipients
- MI
- Myocardial infarction
- RRT
- Renal replacement therapy
- RAAS
- Renin angiotensin aldosteron
Acknowledgments
This study was conducted with the use of equipment purchased by Medical University of Bialystok as part of the RPOWP 2007–2013 funding, Priority I, Axis 1.1, contract No. UDA-RPPD.01.01.00-20-001/15-00 dated 26 June 2015.