The impact of deceased donor maintenance on delayed kidney allograft function: A machine learning analysis

Abstract

Background

This study evaluated the risk factors for delayed graft function (DGF) in a country where its incidence is high, detailing donor maintenance-related (DMR) variables and using machine learning (ML) methods beyond the traditional regression-based models.

Methods

A total of 443 brain dead deceased donor kidney transplants (KT) from two Brazilian centers were retrospectively analyzed and the following DMR were evaluated using predictive modeling: arterial blood gas pH, serum sodium, blood glucose, urine output, mean arterial pressure, vasopressors use, and reversed cardiac arrest.

Results

Most patients (95.7%) received kidneys from standard criteria donors. The incidence of DGF was 53%. In multivariable logistic regression analysis, DMR variables did not impact on DGF occurrence. In post-hoc analysis including only KT with cold ischemia time<21h (n = 220), urine output in 24h prior to recovery surgery (OR = 0.639, 95%CI 0.444–0.919) and serum sodium (OR = 1.030, 95%CI 1.052–1.379) were risk factors for DGF. Using elastic net regularized regression model and ML analysis (decision tree, neural network and support vector machine), urine output and other DMR variables emerged as DGF predictors: mean arterial pressure, ≥ 1 or high dose vasopressors and blood glucose.

Conclusions

Some DMR variables were associated with DGF, suggesting a potential impact of variables reflecting poor clinical and hemodynamic status on the incidence of DGF.

About the Authors:

Silvana Daher Costa

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliations Department of Clinical Medicine, Faculty of Medicine, Federal University of Ceará, Fortaleza, Ceará, Brazil, Walter Cantídio University Hospital, Fortaleza, Ceará, Brazil, Hospital Geral de Fortaleza, Fortaleza, Ceará, Brazil

Luis Gustavo Modelli de Andrade

Roles Formal analysis, Methodology, Resources, Software, Writing – original draft, Writing – review & editing

Affiliation: Department of Internal Medicine, Universidade Estadual Paulista—UNESP, Botucatu, São Paulo, Brazil

Francisco Victor Carvalho Barroso

Roles Investigation, Resources, Writing – review & editing

Affiliation: Department of Clinical Medicine, Faculty of Medicine, Federal University of Ceará, Fortaleza, Ceará, Brazil

Cláudia Maria Costa de Oliveira

Roles Data curation, Resources, Writing – review & editing

Affiliations Walter Cantídio University Hospital, Fortaleza, Ceará, Brazil, Hospital Geral de Fortaleza, Fortaleza, Ceará, Brazil

Elizabeth De Francesco Daher

Roles Conceptualization, Resources, Writing – review & editing

Affiliation: Department of Clinical Medicine, Faculty of Medicine, Federal University of Ceará, Fortaleza, Ceará, Brazil

Paula Frassinetti Castelo Branco Camurça Fernandes

Roles Data curation, Resources, Writing – review & editing

Affiliation: Walter Cantídio University Hospital, Fortaleza, Ceará, Brazil

Ronaldo de Matos Esmeraldo

Roles Data curation, Resources, Writing – review & editing

Affiliation: Hospital Geral de Fortaleza, Fortaleza, Ceará, Brazil

Tainá Veras de Sandes-Freitas

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliations Department of Clinical Medicine, Faculty of Medicine, Federal University of Ceará, Fortaleza, Ceará, Brazil, Hospital Geral de Fortaleza, Fortaleza, Ceará, Brazil

ORCID logo http://orcid.org/0000-0002-4435-0614

Introduction

Brazilian studies have reported incidences of delayed graft function (DGF) between 50 and 70%, 2 to 3-fold higher than the rates described by American and European cohorts, despite similar or more favorable recipient and donor demographics [1–6]. A Brazilian study reported 22.7% incidence of delayed kidney function in a cohort of simultaneous pancreas-kidney transplants, despite a short mean cold ischemia time of 14h and the use of ideal donors [7]. With similar demographics, international cohorts reported incidences of 4–5% [8, 9]. Notable, similar to demonstrated in American and European cohorts, DGF in Brazilian transplant recipients has negative impact on short and long-term outcomes [4, 6, 10].

There is no robust evidence explaining the high DGF incidence in our country, but it is likely that the suboptimal maintenance care of potential donors before organ recovery has an important role. Of note, recent studies have shown that achieving optimal donor maintenance parameters is associated with significant decrease in DGF occurrence [11, 12].

Traditionally, studies evaluating risk factors for DGF adopt standard statistical approaches, such as logistic regression. These models are useful in analysis using few independent variables, mainly when the effect of the predictor on the outcome is linear and homogeneous. The assumptions required to regression-based models are often not reached in clinical research and important predictor variables may be obscured. Machine learning (ML) methods can improve precision and accuracy in predicting events, by using more sensitive statistical methods, with data mining techniques and complex data interactions modeling non-linear interactions [13–15]. As an example, Decruyenaere et al. demonstrated that logistic regression was not the ideal method for DGF prediction in a Belgian cohort and ML methods performed better discriminative capacity [16]. Additionally, other regression-based models are useful and present better performance, depending on the number of events, number of predictors, variables characteristic and distribution [17].

This study aimed to evaluate the risk factors for DGF, including in the analysis donor maintenance-related (DMR) variables, which were thoroughly investigated from multidisciplinary records. To increase the analysis accuracy and properly investigate the impact of donor maintenance on DGF occurrence, we selected a cohort of brain dead donor (DBD) kidney transplants (KT) performed in a Brazilian region where DGF incidence is high despite the predominance of ideal donors. In addition, we used ML methods for data analysis beyond regression models.

Materials and methods

Study design

This study is a retrospective analysis from all deceased donor KT recipients older than 16 years of age, performed between January 1st 2015 and December 31st 2017 at two Brazilian transplant centers, located in a region with locally predominant use of standard criteria donors [18]. Preemptive, multiorgan transplants, recipients of machine perfused grafts and those who lost their grafts or died within 7 days after KT were excluded. In compliance with Brazilian law, all donors were brain dead.

Data were retrospectively collected by systematic review of medical charts and electronic database. Patient records and information was anonymized and de-identified prior to analysis. Due to the observational and retrospective nature of the study, with data anonymously analyzed, informed consent was not obtained.

The study was performed in accordance with ethical standards of National Health Council Resolution 466/12 and Declaration of Helsinki, and was approved by Institutional Review Board (IRB) of the Federal University of Ceará (Ethics Committee approval number: 2.004.286) and by the IRBs of all hospitals involved in the donation and transplantation processes: Walter Cantídio University Hospital (2.183.661), Instituto José Frota (2.183.661) and Hospital Geral de Fortaleza (2.059.876).

Definitions

Delayed graft function was defined as the requirement for at least one dialysis session during the first week after KT, regardless of the clinical indication [19]. DGF duration was assessed by the time until the last dialysis session, and by the number of sessions performed in this period.

To better test the hypothesis that variables not included in traditional DGF predictive models could explain the high incidence in our country, we opted to calculate the expected incidence of DGF using the nomogram described by Irish et al [1] as a starting point of the study. Among all available DGF predictive models, Irish nomogram has demonstrated the best predictive power in validation studies including non-American patients [20, 21]. This nomogram was developed using United Network for Organ Sharing / Organ Procurement and Transplantation Network (UNOS/OPTN) database and include the following variables: recipient ethnicity, gender, body mass index (BMI), and history of previous KT, diabetes or blood transfusion; time on dialysis, peak panel reactive antibodies (PRA), human leucocyte antigens mismatches (HLA MM); donor age, weight, history of hypertension and terminal serum creatinine (sCr), cause of death, donation after cardiac death, cold ischemia time (CIT), and warm ischemia time [1]. Since Irish nomogram does not allow the inclusion of machine-perfused kidneys, we opted to exclude them.

Donor maintenance parameters were evaluated using the Donor Management Goals (DMG) previously described by US Department of Health and Human Services, Health Resources and Services Administration (HRSA) [11, 22], with adaptations considering local peculiarities and the retrospective nature of the study: the lowest arterial blood gas pH during hospitalization was captured and patients who presented values between 7.3 and 7.45 were considered in the goal; the highest serum sodium (Na+) during hospital stay was recorded and patients who presented values between 135 and 155 mEq/L were considered in the goal; the target for the highest blood glucose was ≤ 150 mg/dL; diuresis in the last 24h prior the recovery surgery was considered adequate when between 0.5 and 3 mL/Kg/h; the lowest mean arterial pressure was in the goal when between 60 and 110 mmHg; and target for vasopressors was the use of ≤ 1 vasoactive drug with norepinephrine <0.5 μg/Kg/min. The highest creatine phosphokinase (CPK), history of reversed cardiac arrest and acute kidney injury during hospitalization were also evaluated.

Statistical analysis

Categorical variables were presented as frequency and percentage and compared using Chi-square or Fisher tests. Normally distributed continuous variables were summarized as mean and standard deviation and compared using Student’s t-test. Median was included in description of non-parametric continuous variables and comparison was performed using Mann Whitney-test. A multivariable logistic regression model was fitted to compute covariate-adjusted odds ratios (OR) for DGF. Twenty-seven variables were included in the model. Collinear variables (“final sCr” and “difference between final and initial sCr”) and those with more than 10% of missing values (“CPK”) were excluded. Diabetic donors were also excluded, since this was a near-zero variance predictor (all patients who received grafts from diabetic donors developed DGF). A p-value of <0.15 in univariable analysis was considered statistically significant for including variables in multivariable analysis. For all other analysis, a p-value of <0.05 was considered statistically significant. Statistical analysis was performed using SPSS v.23.0 software (SPSS, Inc., Chicago, IL, USA).

Predictive models for DGF were constructed using a supervised model, according to the following steps: data acquisition, data processing, model construction, and model evaluation. In data processing step, the outcome category was balanced to achieve equal proportions. In pre-process we perform recursive feature elimination to select the predictable variables with better correlation with the outcome. Variables with moderate / strong correlation with DGF (r > 0.80) were included in the model.

Population was randomly divided into training and testing set with a stratified 70:30 split. Six supervised ML algorithms were developed with this subset of variables in the training set: neural network (NN), support vector machine (SVM), C5 decision tree (DT), CHAID DT, k-nearest neighbors (KNN) and logistic regression with stepwise selection. Area under the receiver operating curve (AUC) was calculated to test the ability of each model to distinguish patients in testing set. SPSS Modeler v.18.1 (IBM, Armonk, NY, USA) was used to construct predictive models.

As an additional sensitivity analysis, we performed Elastic Net regression, a regularization model that mix Lasso and Ridge regressions. While Lasso regression enhances the prediction accuracy and interpretability of the statistical model by variable selection and regularization, Ridge regression improves prediction by shrinking large regression coefficients to reduce overfitting. Elastic net regression also provides a more interpretable model when compared to black-box results of ML methods. For this analysis, numerical predictors were normalized by a Box and Cox transformation and after center and scale. Median were imputed for missing values. Categorical predictors were dummy encoded. We used an alpha of 0.1 and 10-fold cross-validation to search lambda. After an optimal search of a lambda, we perform an elastic net with R package glmnet [23]. Those coefficients inferior or greater than zero are considered relevant predictors. Variables whose coefficient was zero were considered not important for DGF prediction.

Results

Recipient and donor demographics

From 954 KT performed in the period, 30 were living donor transplants, 70 were allocated to patients younger than 16 years, 16 were multiorgan transplants, 8 were preemptive KT, 365 recipients received machine perfused grafts, and 22 lost the graft or died in the first week. The final analysis included 443 DD KT. The incidence of DGF was 53%, the mean time of DGF was 11.8 ± 15.0 days (median 7 days) and mean number of dialysis sessions was 5.0 ± 5.1 (median 4). According to Irish nomogram, the expected incidence of DGF was 19% and this tool showed poor predictive accuracy (AUC 0.685).

Recipient and donor demographic characteristics are detailed in Tables 1 and 2. Patients in DGF group were older (45.6 ± 14.4 vs. 42.6 ± 15.0 years old, p = 0.030), presented higher prevalence of pretransplant diabetes (21.3 vs. 13.5%, p = 0.034), and longer time on dialysis (36 vs. 27 months, p = 0.001). Donors in DGF group showed higher mean age (33.2 ± 12.6 vs. 28.5 ± 12.3 years old, p<0.001), higher body mass index (BMI) (25.7 ± 3.5 vs. 24.9 ± 4.0 Kg/m2, p = 0.038), higher prevalence of hypertension (7.2 vs. 3.4%, p = 0.002), higher terminal sCr (1.2 ± 0.7 vs. 1.0 ± 0.5 mg/dL, p = 0.001), and higher Kidney Donor Profile Index (KDPI) (35.1 ± 23.0 vs. 28.4 ± 19.7%, p = 0.001). Of note, 433 donors (97.7%) presented KDPI ≤ 85%.

[Figure omitted. See PDF.]

Table 1. Recipient demographic and clinical characteristics.

https://doi.org/10.1371/journal.pone.0228597.t001

[Figure omitted. See PDF.]

Table 2. Donor demographic and clinical characteristics.

https://doi.org/10.1371/journal.pone.0228597.t002

The main perfusion solution was histidine-tryptophan-ketoglutarate (HTK) (83.1%), followed by University of Wisconsin (UW) (13.1%) and Institute Georges Lopez-1 (IGL-1) (3.8%) and this variable did not impact on DGF incidence. There were no differences between groups regarding vascular anastomosis time (36.5 ± 12.3 vs. 35.6 ± 11.1 min, p = 0.440) and anti-thymoglobulin induction therapy (99.1 vs. 96.6%, p = 0.090). However, DGF group presented longer CIT (21.7 ± 3.8 vs. 20.1 ± 4.1h, p<0.001).

Donor maintenance

Donor maintenance parameters data (Table 3) showed the poor clinical and hemodynamic conditions experienced by donors: high need for vasopressors (95.9%), substantial incidence of reversed cardiac arrest episodes (12.2%), high CPK (median 951 UI/L), increased serum Na+ (160.0 ± 13.8 mEq/L) and poor blood glucose control (193.3 ± 77.7 mg/dL). DGF group had lower percentage of donors reaching blood glucose target (26 vs. 37.5%, p = 0.010) and lower urine output in the last 24h prior to recovery surgery (median 0.9 vs. 1.2 mL/Kg/h, p = 0.005).

[Figure omitted. See PDF.]

Table 3. Donor maintenance parameters.

https://doi.org/10.1371/journal.pone.0228597.t003

Risk factors for DGF

In logistic regression analysis, variables independently associated with DGF were: recipient history of diabetes (OR 1.922, 95% CI 1.119–3.302, p = 0.018), time on dialysis (OR 1.009, 95% CI 1.004–1.014, p<0.001), donor hypertension (OR 2.331, 95% CI 1.247–4.355, p = 0.008), final sCr (OR 1.947, 95% CI 1.320–2.872, p = 0.001), and CIT (OR 1.115 95% CI 1.058–1.175, p<0.001). However, logistic regression presented poor predictive performance, both in training (AUC 0.686) and testing set (AUC 0.695). To best explore the impact of donor maintenance variables on DGF incidence, logistic regression multivariable analysis was repeated in two subgroups of CIT, based on the median value of the total cohort (21h). The incidences of DGF were 48.2% and 57.8% in CIT<21h and CIT≥21h subgroups, respectively. Of note, in CIT<21h subgroup, in addition to the previously demonstrated variables (time on dialysis, donor hypertension and CIT), urine output (OR 0.639, 95% CI 0.444–0.919, p = 0.016) and serum Na+ (OR 1.030, 95% CI 1.007–1.053, p = 0.010) were also risk factors for DGF (Table 4).

[Figure omitted. See PDF.]

Table 4. Logistic regression analysis of risk factors for DGF, according to CIT groups.

https://doi.org/10.1371/journal.pone.0228597.t004

The results of Elastic Net regression are illustrated in Fig 1. Variables in red were risk factors for DGF and variables in blue were associated with reduced risk of DGF. This model presented better predictive performance when compared to logistic regression (AUC 0.749).

[Figure omitted. See PDF.]

Fig 1. Elastic net regression.

Variable in red (coefficient > zero) are risk factors for DGF and variables in blue (coefficient < zero) are protective. Variables are disposed by importance. DGF: delayed graft function; CIT: cold ischemia time; yo: years old; BMI: body mass index, PRA: panel reactive antibodies; HLA MM: human leucocyte antigen mismatches; sCr: final serum creatinine; Na+: sodium; VAT: vascular anastomosis time.

https://doi.org/10.1371/journal.pone.0228597.g001

The three statistical methods of better performance in analyzes using ML techniques were: boosted DT using C5.0 algorithm (AUC 0.791), boosting NN (AUC 0.886), and SVM with polynomial kernel (AUC 0.784). Fig 2 illustrates ML results and predictive performance. Variable are presented using the feature importance graph, in which higher the value of the feature/variable, more important it was to predict DGF. The feature importance was normalized between 0 and 1 by dividing by the sum of all feature importance values. As demonstrated in Fig 2, urine output, mean arterial pressure, and the use of more than one vasoactive drug or high dose vasopressor were predictive of DGF in all models. Blood glucose was predictive of DGF in DT and SVM models.

[Figure omitted. See PDF.]

Fig 2. Figs on the left demonstrate the area under the receiver-operating curve (AUC) for each method in the testing set.

Figs on the right show global variable importance for prediction of delayed graft function. (A) Decision Tree; (B) Support Vector Machine; (C) Neural Network. BMI: body mass index; sCr: final serum creatinine; CIT: cold ischemia time.

https://doi.org/10.1371/journal.pone.0228597.g002

Discussion

This study suggest that poor donor clinical and hemodynamic status may impact on DGF occurrence, and this might explain the high incidence of DGF in Brazil, where the incidence is significantly higher than that the predicted by available formulas.

In our cohort, DGF incidence was almost 3-fold higher than the predicted by the nomogram described by Irish et al, suggesting an important role of variables not included in the prediction model. In fact, none of the available predictive models includes in final formula variables reflecting donor maintenance. Except for terminal serum creatinine (that could reflect renal consequences of hypovolemia, shock and other causes of acute kidney injury), only the score developed by Chapal et al analyzed some variable related to donor care (type of vasopressor) [24].

Only 4.3% of patients received kidneys from expanded criteria donors and 97.7% had KDPI below 85%, suggesting good structural quality kidneys. On the other hand, 12.2% had a reversed cardiac arrest episode before organ recovery surgery, 47.4% presented acid-base disorders, 61.6% had hypo- or hypernatremia, and 68.6% showed inadequate glycemic control, reflecting poor clinical and hemodynamic conditions. The impact of poor donor maintenance in our country has previously suggested in a study including simultaneous pancreas-kidney transplants. In this cohort, despite favorable demographics, the incidence of delayed kidney graft function was 22.7% and donor hypernatremia was an independent risk factor for DGF [7].

The challenge of properly maintaining potential deceased donors seems not to be exclusive of our population. Previous American and Canadian studies have reported low adherence to donor-care bundles at the time of consent for donation. Importantly, the early achievement of donor management goals was associated with a reduced risk of DGF in these studies [11, 12].

The high CIT was notable and its contribution to DGF is unequivocal. All statistical analysis demonstrated that each additional hour matters, even in KT with CIT < 21h. The large territorial extension, the allocation model predominantly based on HLA compatibility, and the absence of specific allocation policies for “marginal” donors contribute to the long CIT in our country [25]. However, it is noteworthy that almost half of KT with CIT <21h presented DGF, suggesting the contribution of other factors beyond the CIT.

Due to the high negative impact of CIT on DGF incidence, potentially masking other predictors, we performed a post-hoc analysis including a regression model on a sub-sample of patients with CIT inferior to 21h. As hypothesized, variables reflecting donor maintenance now emerged. Additionally, in ML analysis, donor maintenance related variables, such as blood pressure, use of high dose vasopressors, urine output and blood glucose, were also associated with DGF. To further explore the contribution of other than the traditional variables to DGF incidence using a more interpretable model, we performed a sensitivity analysis using elastic net regression. Again, beyond the traditional conditions associated to DGF occurrence, variables reflecting donor management were risk factors (serum Na+, blood glucose) or protective (high dose vasopressors and diuresis).

Recently published studies demonstrated the impact of donor hemodynamics as a predictor of DGF in transplantation from donors after cardiac death [26, 27]. However, evidences are scarce on transplantation of brain dead donors.

Our study has limitations that should be pointed out. First, it was a retrospective cohort; therefore, the capture of variables was limited to those available on medical records. We cannot assure that our results are generalizable to other transplant centers worldwide. Donor maintenance related variables (and DMGs) are not static and we could not evaluate them dynamically. Since we could not follow variables over time and assess adherence to care bundles, it is possible that some clinical and hemodynamic data reflected the severity of the disease that led to brain death and not the patient care. Besides, it’s possible that some clinical and hemodynamic parameters reflected patient situation before they became a consented donor, and thus may not reflect donor maintenance, but general intensive care. Finally, variables reflecting perioperative care were not available.

In conclusion, DGF incidence in Brazil is significantly higher than that predicted by available models. Although our data do not allow us to draw definitive conclusions, this study suggests that donor illness severity and hemodynamic instability might contribute to this scenario. Prospective studies are needed to robustly conclude how donor management impacts on DGF incidence. Additionally, a cohort including machine-perfused grafts may be useful to explore if pumping might mitigate kidney damage secondary to poor donor clinical and hemodynamic status.

We believe bringing this issue up is crucial in our setting. Scant donor care is probably a reflection of the poor economic conditions of our country. However, some educational actions could be taken, focusing on early recognition of potential donors and training staff who care for brain death patients.

Supporting information

[Figure omitted. See PDF.]

S1 Database.

https://doi.org/10.1371/journal.pone.0228597.s001

(XLSX)

Citation: Costa SD, de Andrade LGM, Barroso FVC, Oliveira CMCd, Daher EDF, Fernandes PFCBC, et al. (2020) The impact of deceased donor maintenance on delayed kidney allograft function: A machine learning analysis. PLoS ONE 15(2): e0228597. https://doi.org/10.1371/journal.pone.0228597

References

1. Irish WD, Ilsley JN, Schnitzler MA, Feng S, Brennan DC. A risk prediction model for delayed graft function in the current era of deceased donor renal transplantation. Am J Transplant. 2010;10(10):2279–86. pmid:20883559.

2. Fonseca I, Teixeira L, Malheiro J, Martins LS, Dias L, Castro Henriques A, et al. The effect of delayed graft function on graft and patient survival in kidney transplantation: an approach using competing events analysis. Transpl Int. 2015;28(6):738–50. pmid:25689397.

3. Gavela Martinez E, Pallardo Mateu LM, Sancho Calabuig A, Beltran Catalan S, Kanter Berga J, Avila Bernabeu AI, et al. Delayed graft function after renal transplantation: an unresolved problem. Transplant Proc. 2011;43(6):2171–3. pmid:21839224.

4. de Sandes-Freitas TV, Felipe CR, Aguiar WF, Cristelli MP, Tedesco-Silva H, Medina-Pestana JO. Prolonged Delayed Graft Function Is Associated with Inferior Patient and Kidney Allograft Survivals. PLoS One. 2015;10(12):e0144188. pmid:26679933; PubMed Central PMCID: PMC4683001.

5. Azevedo LS, Castro MC, Monteiro de Carvalho DB, d'Avila DO, Contieri F, Gonçalves RT, et al. Incidence of delayed graft function in cadaveric kidney transplants in Brazil: a multicenter analysis. Transplant Proc. 2005;37(6):2746–7. pmid:16182798.

6. Helfer MS, Vicari AR, Spuldaro F, Gonçalves LF, Manfro RC. Incidence, risk factors, and outcomes of delayed graft function in deceased donor kidney transplantation in a Brazilian center. Transplant Proc. 2014;46(6):1727–9. pmid:25131022.

7. Rangel EB, Melaragno CS, Gonzalez AM, Linhares MM, de Sá JR, Salzedas A, et al. Delayed kidney allograft function after simultaneous pancreas-kidney transplantation. Transplant Proc. 2010;42(9):3655–9. doi: S0041-1345(10)01308-4 [pii] pmid:21094834.

8. Gutierrez P, Marrero D, Hernandez D, Vivancos S, Perez-Tamajon L, Rodriguez de Vera JM, et al. Surgical complications and renal function after kidney alone or simultaneous pancreas-kidney transplantation: a matched comparative study. Nephrol Dial Transplant. 2007;22(5):1451–5. pmid:17327283.

9. Bunnapradist S, Cho YW, Cecka JM, Wilkinson A, Danovitch GM. Kidney allograft and patient survival in type I diabetic recipients of cadaveric kidney alone versus simultaneous pancreas kidney transplants: a multivariate analysis of the UNOS database. J Am Soc Nephrol. 2003;14(1):208–13. pmid:12506153.

10. Butala NM, Reese PP, Doshi MD, Parikh CR. Is Delayed Graft Function Causally Associated With Long-Term Outcomes After Kidney Transplantation? Instrumental Variable Analysis. Transplantation. 2013;95(8):1008–14. pmid:23591726; PubMed Central PMCID: PMC3629374.

11. Malinoski DJ, Patel MS, Ahmed O, Daly MC, Mooney S, Graybill CO, et al. The impact of meeting donor management goals on the development of delayed graft function in kidney transplant recipients. Am J Transplant. 2013;13(4):993–1000. pmid:23406284.

12. Cardinal H, Lamarche F, Grondin S, Marsolais P, Lagace AM, Duca A, et al. Organ donor management and delayed graft function in kidney transplant recipients: A multicenter retrospective cohort study. Am J Transplant. 2019;19(1):277–84. pmid:30253052.

13. Hannun PGC, Andrade LGM. The future is coming: promising perspectives regarding the use of machine learning in renal transplantation. J Bras Nefrol. 2018. pmid:30353909.

14. Goldstein BA, Navar AM, Carter RE. Moving beyond regression techniques in cardiovascular risk prediction: applying machine learning to address analytic challenges. Eur Heart J. 2017;38(23):1805–14. pmid:27436868; PubMed Central PMCID: PMC5837244.

15. Yoo KD, Noh J, Lee H, Kim DK, Lim CS, Kim YH, et al. A Machine Learning Approach Using Survival Statistics to Predict Graft Survival in Kidney Transplant Recipients: A Multicenter Cohort Study. Sci Rep. 2017;7(1):8904. pmid:28827646; PubMed Central PMCID: PMC5567098.

16. Decruyenaere A, Decruyenaere P, Peeters P, Vermassen F, Dhaene T, Couckuyt I. Prediction of delayed graft function after kidney transplantation: comparison between logistic regression and machine learning methods. BMC Med Inform Decis Mak. 2015;15:83. pmid:26466993; PubMed Central PMCID: PMC4607098.

17. Zou HH, Trevor. Regularization and variable selection via the elastic net. J R Statist Soc B. 2005;part 2(67):301–20.

18. Silva SF, Silva SL, Nascimento AC, Parente MM, Albuquerque CA, Rodrigues AA, et al. Profile of organ donors in Ceara, northeastern Brazil, from 1998 to 2012. Transplant Proc. 2014;46(6):1692–4. pmid:25131014.

19. Mallon DH, Summers DM, Bradley JA, Pettigrew GJ. Defining delayed graft function after renal transplantation: simplest is best. Transplantation. 2013;96(10):885–9. pmid:24056620.

20. Rodrigo E, Minambres E, Ruiz JC, Ballesteros A, Pinera C, Quintanar J, et al. Prediction of delayed graft function by means of a novel web-based calculator: a single-center experience. Am J Transplant. 2012;12(1):240–4. pmid:22026730.

21. Kers J, Peters-Sengers H, Heemskerk MBA, Berger SP, Betjes MGH, van Zuilen AD, et al. Predicition models for delayed graft function: external validation on The Dutch Prospective Renal Transplantation Registry. Nephrol Dial Transplant. 2018;33(7):1259–68. pmid:29462353.

22. Franklin GA, Santos AP, Smith JW, Galbraith S, Harbrecht BG, Garrison RN. Optimization of donor management goals yields increased organ use. Am Surg. 2010;76(6):587–94. pmid:20583513.

23. Friedman J, Hastie T, Tibshirani R. Regularization Paths for Generalized Linear Models via Coordinate Descent. J Stat Softw. 2010;33(1):1–22. pmid:20808728; PubMed Central PMCID: PMC2929880.

24. Chapal M, Le Borgne F, Legendre C, Kreis H, Mourad G, Garrigue V, et al. A useful scoring system for the prediction and management of delayed graft function following kidney transplantation from cadaveric donors. Kidney Int. 2014;86(6):1130–9. pmid:24897036.

25. Sandes-Freitas TV. The scenario of delayed graft function in Brazil. J Bras Nefrol. 2019. pmid:30806446.

26. Allen MB, Billig E, Reese PP, Shults J, Hasz R, West S, et al. Donor Hemodynamics as a Predictor of Outcomes After Kidney Transplantation From Donors After Cardiac Death. Am J Transplant. 2016;16(1):181–93. pmid:26361242.

27. Sun Q, Huang Z, Zhou H, Lin M, Hua X, Hong L, et al. New Factors Predicting Delayed Graft Function: a Multi-Center Cohort Study of Kidney Donation After Brain Death Followed by Circulatory Death. Kidney Blood Press Res. 2018;43(3):893–903. pmid:29871001.