Skip to main content

Prediction of acute kidney injury after cardiac surgery: model development using a Chinese electronic health record dataset

Abstract

Background

Acute kidney injury (AKI) is a major complication following cardiac surgery that substantially increases morbidity and mortality. Current diagnostic guidelines based on elevated serum creatinine and/or the presence of oliguria potentially delay its diagnosis. We presented a series of models for predicting AKI after cardiac surgery based on electronic health record data.

Methods

We enrolled 1457 adult patients who underwent cardiac surgery at Nanjing First Hospital from January 2017 to June 2019. 193 clinical features, including demographic characteristics, comorbidities and hospital evaluation, laboratory test, medication, and surgical information, were available for each patient. The number of important variables was determined using the sliding windows sequential forward feature selection technique (SWSFS). The following model development methods were introduced: extreme gradient boosting (XGBoost), random forest (RF), deep forest (DF), and logistic regression. Model performance was accessed using the area under the receiver operating characteristic curve (AUROC). We additionally applied SHapley Additive exPlanation (SHAP) values to explain the RF model. AKI was defined according to Kidney Disease Improving Global Outcomes guidelines.

Results

In the discovery set, SWSFS identified 16 important variables. The top 5 variables in the RF importance matrix plot were central venous pressure, intraoperative urine output, hemoglobin, serum potassium, and lactic dehydrogenase. In the validation set, the DF model exhibited the highest AUROC (0.881, 95% confidence interval [CI] 0.831–0.930), followed by RF (0.872, 95% CI 0.820–0.923) and XGBoost (0.857, 95% CI 0.802–0.912). A nomogram model was constructed based on intraoperative longitudinal features, achieving an AUROC of 0.824 (95% CI 0.763–0.885) in the validation set. The SHAP values successfully illustrated the positive or negative contribution of the 16 variables attributed to the output of the RF model and the individual variable’s effect on model prediction.

Conclusions

Our study identified 16 important predictors and provided a series of prediction models to enhance risk stratification of AKI after cardiac surgery. These novel predictors might aid in choosing proper preventive and therapeutic strategies in the perioperative management of AKI patients.

Background

Acute kidney injury (AKI), a common and potentially life-threatening clinical syndrome, is more and more frequent with increasing cardiac surgical volume in developed and developing countries. A meta-analysis reported an incidence of cardiac surgery-associated acute kidney injury (CSA-AKI) of 26.0–28.5% [1]. AKI not only severely affects in-hospital morbidity and mortality but also long-term prognosis. Patients who survive an episode of AKI after surgery are also at elevated risk of developing major adverse cardiovascular events, advanced chronic kidney disease, and all-cause death [2].

The current consensus definition of AKI depends on the increase of serum creatinine (Scr) and/or the presence of oliguria; however, this can lead to delayed diagnosis and treatment. Consequently, substantial efforts have been made to explore biomarkers or develop clinical prediction models in recent years. Several novel biomarkers were proposed to substitute Scr in the assessment of kidney function, such as NGAL, KIM-1, and DKK3 [3, 4]. However, an individual biomarker is inadequate in predicting AKI; that is to say, the pathophysiology of CSA-AKI is multifactorial and intricate. On the other hand, these biomarkers are costly and difficult to assay and are thus out of consideration by the majority of clinicians [5]. Clinical scoring systems (e.g., Cleveland Clinic score, Simplified Renal Index score, Mehta score, etc.) have been introduced into clinical practice for more than a decade [6,7,8], while widespread adoption of these models would be challenging. First, these models are developed following the traditional logistic regression method. Their derivation requires the statistical assumption regarding a linear relationship between covariates and outcomes, with analytical model restricting to selection of a small set of parameters that are known to be clinically relevant. Second, most of these models perform traditional feature selection methods based on a small scale of exposure variables and identify common risk factors (e.g., age, diabetes mellitus, hypertension, cardiac function, surgery type, etc.). These risk factors reflect preoperative conditions and generally present inadequate predictive power in different races or populations [9].

Recent advancements in electronic health record (EHR) systems, data accessibility, and artificial intelligence have raised great interest in developing completely electronic data-driven machine learning (ML) models for predicting specific clinical outcomes. Several studies have used ML to predict inpatient AKI using EHR data. For instance, a continuous deep learning algorithm can predict 55.8% of all inpatient episodes of AKI up to 48 h in advance and over 90% of all AKI patients requiring subsequent renal replacement therapy (RRT) [10]. Therefore, our first objective was herein to develop three tree-based ML models to predict CSA-AKI by incorporating preoperative, intraoperative, and early postoperative data from the EHR of Nanjing First Hospital. In addition, to uncover the “black-box” of ML, SHapley Additive exPlanation (SHAP) values were utilized to explain the ML model and evaluated individual variable prediction [11]. Given the generalizability of the linear model, a nomogram model was finally constructed based on the logistic regression analysis.

Method

Study design and participants

This is a retrospective, observational study. Consecutive patients who underwent cardiac surgery, admitted between January 2017 and June 2019, were recruited from Nanjing First Hospital. We enrolled patients who had received coronary artery bypass grafting (CABG), valve surgery, and a combination of both treatments. Patients were excluded if they met the following criteria: (i) aged < 18 years; (ii) preoperative AKI, end-stage renal disease, or dialysis; (iii) did not receive cardiopulmonary bypass (CPB); (iv) missing Scr data. Patient informed consent in this study was waived due to the retrospective nature of the study. This study was approved by the Ethical Committee of Nanjing First Hospital. We reported our work following TRIPOD statement guideline [12].

Data preprocessing

The study cohort was acquired from two population-based databases, consisting of patient information available from EHR in digital format: Jiangsu Province Coronary Artery Bypass Grafting Register (218.2.200.37:2356/Multicenter) and Patient Information Management Platform (218.2.200.37:2356/PatientList). The databases are specifically designed for cardiac patients and consist of perioperative clinical characteristics including patient demographics, admission assessment, comorbidities, laboratory test, medication, surgical information, and CPB data. All clinical data regarding preoperative, intraoperative and early postoperative variables were included for model derivation (preoperative laboratory biomarkers were collected at 6 a.m. the next day following hospital admission; early postoperative variables were measured within 6 h after surgery). Missing values were filled in by a second manual review of the EHR, and personal information was de-identified before delivering for analysis. The final dataset was randomly partitioned to a discovery (80% of observations) set and a validation (20%) set.

Anesthesia, CPB, and critical care

All participants received general intravenous-inhalation combined anesthesia, which was maintained intraoperatively with a continuous infusion of propofol (4–6 mg/kg/h), remifentanil (0.2–0.4 μg/kg/min), cisatracurium (0.2–0.3 mg/kg/h) as well as an intermittent addition of sufentanil and midazolam. CPB was performed with non-pulsatile perfusion, with a perfusion flow of 2.0–2.8 L/kg/min and a mean arterial pressure of 55–85 mmHg in most cases. During CPB, monitoring records included nasopharyngeal temperature, bladder temperature, rectal temperature, perfusion flow, oxygen delivery, perfusion pressure, central venous pressure (CVP), conventional ultrafiltration (CUF), and urine output. After the surgery was completed, patients were transferred to the intensive care unit (ICU) and placed on ventilators in synchronized intermittent mandatory ventilation or assist-control models set at 8–10 mL/kg tidal volume and 5 cmH2O positive end-expiratory pressure. Arterial blood gas was checked at the time of ICU admission; other laboratory measurements (e.g., blood cell analysis, liver and kidney function, coagulation function, etc.) were obtained within 6 h postoperatively.

AKI definition

Postoperative AKI was defined according to the Scr-based criteria from the Kidney Disease: Improving Global Outcomes (KDIGO) consensus definition, specifically an acute increase in Scr ≥ 50% within 7 days or ≥ 0.3 mg/dL within 48 h compared with the baseline level, or a requirement for RRT [13]. In this case, the patient’s baseline Scr was determined by the Scr level measured at hospital admission.

Model development

This study mainly comprised two stages: (1) feature selection; and (2) model development. We used a ML-based feature selection method, the sliding windows sequential forward feature selection (SWSFS), to identify the number of important variables among the 193 clinical features. Then two types of models were developed based on the selected variables including three tree-based supervised learning models and a nomogram model.

Feature selection

First, we applied a random forest (RF) classifier to calculate variable importance score (VIS) according to the Gini index via the importance function in R. To minimize random errors, we performed the RF 30 times by setting the random seeds from 1 to 30 and calculated the average Gini index of each variable. SWSFS was used to identify a set of important variables. Briefly, VIS of all the clinical features (except age and gender) was obtained from RF and ranked by the averaged Gini index in descending order. Next, the features were included one by one to the RF model based on their VIS ranks. Afterward, we plotted the model error, which measured the “out of bag (OOB)” rate of each RF model consisting of different numbers of variables. Finally, the number of features was identified based on the lowest model error rate.

Tree-based ensemble algorithms

We developed the prediction models using the following tree ensemble methods, which are the most popular and advanced ML algorithms for binary classification: extreme gradient boosting (XGBoost), RF, and deep forest (DF). Both XGBoost and RF use rules to binary split data based on decision trees. Generally, a tree with many splits will probably lead to overfitting and result in poor performance in new datasets. RF works based on the idea of the ensemble method, which collects individual decision trees, bagging and random feature selection, thus providing more accurate results and making the model more resistant to overfitting [14]. XGBoost is the engineering realization of the Gradient Boosting Decision Tree, which provides superior prediction by combining multiple decision trees in boosting ways [15]. As a novel and advanced deep learning method, DF generates a multi-layer cascade forest containing various RFs [16]. This structure has been designed to ensure the diversity of the model by including different types of forests. Each layer in the cascade forest receives the information processed by the previous stage and outputs the processing results to the next layer (Additional file 1: Fig. S1).

The area under the receiver operating characteristic curve (AUROC) was performed to compare the discrimination of various ML models. Model calibration was evaluated using a calibration plot based on the isotonic regression method [17]. Furthermore, for evaluating ML models, a series of interpretable parameters were determined: accuracy, sensitivity, specificity, positive predictive value, negative predictive value, and F score. It is worth noting that the correct interpretability of a ML model is challenging. We used SHapley Additive exPlanation (SHAP) values to explain the RF model. SHAP values are developed based on the concept of Shapley values from cooperative game theory [18]. It approximates a complex model to a linear model and evaluates variable importance to demonstrate the amount by which a given feature changes the prediction. Specificity, the SHAP values not only highlight the contributions of the individual variable to the model but also demonstrate the influence of each variable on global model effects. In addition, we used SHAP plot function to uncover the complex relationship between variables and outcomes in the RF model.

Nomogram construction

We further constructed a nomogram model based on the variables selected by SWSFS. It is a logistic regression model. Any multicollinearity of the variables was excluded by establishing the variance inflation factor (VIF); the maximum VIF was 1.41 (Additional file 1: Fig. S2). Given the high association between AKI and procedure-related factors, the important intraoperative longitudinal data were handled as group-based trajectory modelling (GBTM). This approach is designed to identify clusters of individuals following similar progressions over time. We used the Stata command, traj, a plugin based on SAS PROC TRAJ macro, to fit a semi-parameter model for longitudinal data by maximum likelihood estimation [19]. The optimal number of groups with at least 5% of patients in the smallest trajectory was determined by establishing the Bayesian information criterion. Subsequently, the trajectory pattern of longitudinal data and remaining important features were incorporated into the multivariate logistic regression analysis to generate a nomogram model. This nomogram model included dynamic intraoperative information and provided a dynamic prediction paradigm. The discrimination of the nomogram was accessed using C-statistic, an index equivalent to AUROC. Model calibration was evaluated using the Brier score and visualized with a calibration plot using the 1000 bootstrap resampling method. A lower Brier score indicates superior calibration [20]. The clinical net benefit of the nomogram was estimated by decision curve analysis [21].

Statistical analysis

For descriptive analyses, continuous variables are described as means (standard deviation) and categorical variables as frequencies with proportions. The clinical characteristics of patients who developed AKI or did not were compared using the Student’s t test, Mann Whitney U test, chi-square test, or Fisher’s exact probability method as appropriate. After a second manual review for data integrity, the dataset had missing data ranging from 0% to 1.37%. We handled the missing values by using the multiple imputation method. Statistical analyses were performed using Stata (version 13.0) with the package of traj, R (version 4.0.3) with packages of mice, rms, and rmda, and Python (version 3.8) with the packages of sklearn, deep-forest, and shap. A two-sided P value < 0.05 was considered statistically significant.

Results

Study population

From January 2017 to June 2019, 1457 consecutive patients were enrolled in the final cohort (Additional file 1: Fig. S3). The mean (standard deviation) of their age was 60 (12.3) years, and 848 (58.2%) of the participants were male. Of them, 486 (33.4%) underwent CABG, 802 (55.0%) were subjected to valve surgery, and 169 (11.6%) received concomitant CABG and valve surgery. Each individual patient had 193 clinical features. They were randomly assigned to a discovery set (n = 1170) or a validation set (n = 287). The rates of AKI were 24.3% and 24.0% in the discovery and validation sets, respectively. In both sets, differences in clinical characteristics between patients who developed AKI and those who did not are outlined in Additional file 1: Table S1.

Feature selection

The VIS of each variable was obtained from the RF algorithm and ranked in descending order. Figure 1 illustrates the importance matrix plot of the top 100 features. Afterward, all features (except age and gender) were included in SWSFS one by one in order of their VIS ranks. Based on the lowest OOB error rate (Fig. 2), SWSFS identified 14 important features including five preoperative factors (Scr, neutrophil to lymphocyte ratio [NLR], blood glucose, uric acid [UA], and high-density lipoprotein [HDL]), five intraoperative factors (urine output, ultrafiltration volume, CVP_T3, CVP_T4, and perfusion flow_T3), and four early postoperative factors (intubated PaO2/FiO2 ratio, hemoglobin, serum potassium, and lactic dehydrogenase [LDH]).

Fig. 1
figure 1

Importance matrix plot of random forest algorithm based on the discovery set

Fig. 2
figure 2

Feature selection by the sliding windows sequential forward feature selection method. First, variable importance score (VIS) of features was obtained from a random forest algorithm and ranked in descending order. The variables were then included one by one to the random forest model in order of VIS rank. Finally, the optimal number of variables (14 variables) was determined by minimum model error (red circle)

Tree-based learning models

Two covariates (age and gender) and the 14 important variables were included in the following ML models: XGBoost, RF, and DF. Figure 3A illustrates the model performance in the validation set, measured by AUROC. The DF model exhibited the largest AUROC (0.881, 95% confidence interval [CI] 0.831–0.930), followed by the RF model (0.872, 95% CI 0.820–0.923) and XGBoost model (0.857, 95% CI 0.802–0.912). Accordingly, five-fold cross-validation results of AUROC and accuracy of XGBoost, RF, and DF are summarized in Additional file 1: Table S2. The calibration plot demonstrated that all these models had well calibrated. The Brier scores were 0.109, 0.117, and 0.116 for the DF, RF, and XGBoost models, respectively (Fig. 3B). The parameters of these algorithms are presented in Additional file 1: Table S3.

Fig. 3
figure 3

Performance of machine learning models. (A) Comparison of area under the receiver operating characteristic curves among XGBoost, random forest, and deep forest models in the validation set; (B) Calibration plots of XGBoost, random forest, and deep forest models in the validation set. The Brier score is reported in the lower right legend (the smaller the value, the better the calibration). AUROC: area under the receiver operating characteristic curve

SHAP values through visualization

We used SHAP values to provide accurate attribution values for each variable within the RF model. Figure 4A describes the SHAP summary plot for each input variable in the discovery set. The y-axis displays the 16 variables ranked in order of importance with their mean absolute SHAP values. The x-axis indicates the SHAP values associated with each variable and patient, which allowed the determination of whether a feature had a negative effect on the prediction toward a non-AKI class or positive effect on the prediction toward an AKI class. (Fig. 4B). Besides, SHAP dependent plot provides a visualization of the impact for individual observations. Using zero as a dividing line, a feature’s effect can be clearly observed regarding its positive or negative contribution to the model (Fig. 5). The individual patient-level prediction was depicted using the SHAP decision plot. This plot function allows for a better understanding of the individual decision path and describes why patients A-C were predicted as AKI whereas patients D-F were not (Fig. 6).

Fig. 4
figure 4

SHAP summary plot of the random forest model. (A) Average absolute impact of variables on the final model output magnitude ordered by decreasing feature importance; (B) The plot depicts the dot estimation on the model output of random forest model. Each dot represents an individual patient from the dataset. Red represents higher SHAP value of specific features; blue represents lower SHAP value of specific features. The higher the SHAP values, the greater the risk of developing acute kidney injury development. CVP: central venous pressure; CUF: conventional ultrafiltration; LDH: lactic dehydrogenase; NLR: neutrophil to lymphocyte ratio; Scr: serum creatinine; UA: uric acid; HDL: high-density lipoprotein

Fig. 5
figure 5

SHAP dependence plot of the random forest model. Each panel demonstrates that each feature affects the output of the random forest prediction model. The x-axis represents the raw values of each feature and the y-axis indicates the SHAP values of features. When the SHAP value of a specific feature exceeds zero, it indicates an increased risk of acute renal injury. CVP: central venous pressure; CUF: conventional ultrafiltration; LDH: lactic dehydrogenase; NLR: neutrophil to lymphocyte ratio; Scr: serum creatinine; UA: uric acid; HDL: high-density lipoprotein

Fig. 6
figure 6

SHAP decision plot of individual patient-level prediction. The plots depict the decision path for predicting acute renal injury. All patients started with mean predictive values and were evaluated at each factor level to obtain the final probability of acute kidney injury. AC Delineate illustrative examples of patients predicted to be acute renal injury; DF display illustrative examples of patients predicted to be non-acute renal injury. CVP, central venous pressure; CUF: conventional ultrafiltration; LDH: lactic dehydrogenase; NLR: neutrophil to lymphocyte ratio; Scr: serum creatinine; UA: uric acid; HDL: high-density lipoprotein

Group-based trajectory modelling and nomogram model

Among the 16 important predictors, CVP_T3, CVP_T4, and perfusion flow_T3 are CPB-related factors associated with AKI development. To observe their dynamic pattern, GBTMs were created in both sets, in which the three clusters of CVP and perfusion flow following different progressions could be observed. The higher indices the trajectory groups of CVP and perfusion flow, the greater the risk of AKI development (Fig. 7). Inclusion of trajectory groups and other important variables in the multivariate logistic regression model resulted in 14 predictors (excluding UA, P > 0.1) that were statistically significant for AKI (Table 1). These independently associated risk factors were incorporated to form an AKI estimation nomogram (Fig. 8). The nomogram was internally validated using the bootstrap validation method. The nomogram achieved good discrimination in the discovery set, with a C-statistic of 0.827 (95% CI 0.800–0.854) and a bootstrap-adjusted C-statistic of 0.810. Correspondingly, in the validation set, the nomogram displayed a C-statistic of 0.824 (95% CI 0.763–0.885). In both sets, 1000 bootstrap resampling calibration plots confirmed an optimal agreement between the predicted and observed risk of AKI (Fig. 9A–D). Furthermore, the decision curve analysis revealed that the nomogram could provide clinical net benefit for most of the examined probabilities (Fig. 9E, F).

Fig. 7
figure 7

Group-based trajectory modelling of central venous pressure and perfusion flow in the discovery (A) and validation (B) sets. These variables were repeatedly measured at specific time points (T1, initial phase of cardiopulmonary bypass; T2, 20 min after aortic cross-clamping; T3, 10 min after aortic declamping; T4, before weaning off from cardiopulmonary bypass)

Table 1 Multivariate logistic regression model for predicting acute kidney injury after cardiac surgery based on discovery set
Fig. 8
figure 8

Nomogram predicting acute kidney injury after cardiac surgery. Scr: serum creatinine; NLR: neutrophil to lymphocyte ratio; HDL: high density lipoprotein; CUF: conventional ultrafiltration; CVP: central venous pressure; LDH: lactic dehydrogenase; AKI: acute kidney injury

Fig. 9
figure 9

Performance evaluation of the nomogram model. Receiver operating characteristic curve of the nomogram in the discovery (A) and validation (B) sets; 1000-resample bootstrapped calibration plot of the nomogram in the discovery (C) and validation (D) sets; decision curve analysis of the nomogram in the discovery (E) and validation (F) sets. AUROC, area under the receiver operating characteristic curve; CI: confidence interval

Discussion

In this study, we applied the SWSFS technique to screen for clinical characteristics and developed a series of models to optimize AKI prediction following cardiac surgery. Using SWSFS, we identified 14 important risk factors associated with AKI among the 193 clinical variables, thus boosting efficiency by incorporating preoperative, intraoperative, and early postoperative data from EHR. The performance of the tree ensemble ML algorithms (XGBoost, RF, and DF) was clinically satisfactory, with AUROC ranging 0.857–0.881 in the validation set. In addition, we constructed a nomogram model for AKI which also presented good performance both in terms of discrimination (AUROC 0.824) and calibration (Brier score 0.144). Our study highlights the value of EHR data in the evaluation of AKI. These important perioperative factors might be helpful in providing individualized preventive strategies and delivering proper treatments in the management of AKI after cardiac surgery.

In addition to some well-known risk factors (e.g., age, gender, hypertension, diabetes mellitus, Scr, etc.) that have been identified by previous studies, most variables are novel predictors for AKI risk prediction. Intraoperative factors reflect acute physiological responses during surgery and play pivotal roles in the development of AKI, particularly the unique physiological perturbations of CPB. In this study, we generated developmental trajectories to describe the course of factors over time. Two specific time point measurements of CVP (T3 and T4) were found to be strong predictors, highlighting the effect of intraoperative venous congestion on renal function. Traditionally, CSA-AKI is considered to be caused by renal hypoperfusion due to hypotension, inadequate perfusion flow, or renal ischemia [22, 23]. This concept has been challenged by accumulating evidence that elevated CVP is a more powerful hemodynamic determinant than mean perfusion pressure for the development of postoperative AKI [24]. More recently, Lopez et al. [25] uncovered that higher levels of CVP during cardiac surgery were independently associated with higher odds of AKI. They also demonstrated that venous congestion is more accurate than hypotension in predicting AKI. CVP-induced AKI can therefore be regarded as “congestive kidney failure”. When examining the trajectory pattern of perfusion flow, the separate clusters of perfusion flow at T3 measurement in our trajectory analysis implied that some patients (high-level cluster) might go through cardiac insufficiency, hemodynamic instability, or other unstable intraoperative conditions after aortic declamping and require additional mechanical assistance before weaning off from CPB. The positive effect of ultrafiltration volume on AKI risk may be another latent indicator of the congestive state of the body, as CUF is used for fluid removal to reduce fluid overload [26]. Taken together, these findings demonstrate that not only renal ischemia but also renal congestion plays a vital role in worsening kidney function. Our study suggested that the clinician may need to pay more attention to hemodynamic changes, in particular during the cardiac resuscitation period. Although adequate urine output does not assure normal kidney function as a result of non-pulsatile flow and cold-induced diuresis, the presence of oliguria typically indicates an acute response to renal hypoperfusion. This result is consistent with the study by Tseng et al. [27], who observed that intraoperative urine output was the most influential feature in predicting AKI.

Early postoperative laboratory biomarkers could reflect the acute pathophysiology of kidney injury. In this study, we identified four laboratory biomarkers (intubated PaO2/FiO2 ratio, serum potassium, hemoglobin, and LDH) associated with CSA-AKI. These biomarkers reflect the patient’s overall disease severity. For example, hypoxemia (low intubated PaO2/FiO2 ratio) is a severe complication after cardiac surgery and has been shown to be highly related to prolonging mechanical ventilation, respiratory complications, and in-hospital death [28]. The loss in glomerular filtration rate reserve usually leads to electrolyte disorder. It was observed that some patients experienced abrupt deterioration in renal function and decreased glomerular filtration capacity in the early postoperative period, resulting in fluid overload and elevated serum potassium levels [29]. LDH is abundant in the kidney, heart, liver, and muscle and is, therefore, most commonly measured to detect tissue damage as well as disease severity of critical patients [30]. Although LDH acts as a nonspecific biomarker for kidney injury, it demonstrates adequate predictive value for AKI risk prediction in several clinical settings [31, 32]. On the other hand, the elevated LDH in the immediate postoperative period may be an indicator of hemolysis from CPB, and CPB-induced hemolysis is associated with the development of AKI [33].

Previous studies have concluded that inflammation, oxidative stress, and endothelial dysfunction are central components of the pathogenesis of AKI. NLR, a promising marker of inflammation, has been identified as a novel predictor of AKI [34]. Moreover, it has emerged as a potential biomarker for lethal outcomes and adverse events in patients undergoing cardiac surgery [35, 36]. Interestingly, lower baseline HDL levels were independently associated with an increased risk of AKI after cardiac surgery. This relationship was also observed in Smith et al.’s study [37]. Systematically, high HDL levels inhibit systemic inflammation and reduce oxidative stress via acting as receptors of prooxidant lipids and associated antioxidant enzymes and therefore play a role in the pathogenesis of AKI [37,38,39]. However, the natural effect of HDL may be pharmacologically modified by traditional preoperative lipid-lowering treatment in cardiac patients; novel pharmacologic agents with the potential for improving HDL function are warranted. Collectively, the pathophysiology of CSA-AKI is complex and multifactorial. Our study identified a set of important factors attributed to CSA-AKI such as venous congestion, renal hypoperfusion, inflammatory response, metabolic disorder, and baseline renal function. However, many of these factors are modifiable. Using these factors may therefore provide a basis for early diagnosis, prevention, and treatment strategies in the preoperative, intraoperative, and early postoperative management of AKI patients, such as optimizing hemodynamic status, intensified early postoperative biomarkers monitoring, and individualized blood glucose and lipids management.

Before implementing ML models in clinical practice, their predictive power must be validated in different clinical settings. Previous studies developed ML models using all features as input variables [27, 40]. However, inclusion of numerous features makes the models much more complicated and difficult to validate in additional representative datasets, as most variables are irrelevant to AKI classification. Besides, many ML algorithms exhibit a decrease in predictive power when the number of variables is significantly higher than optimal [41]. For example, as outlined in Additional file 1: Fig. S4, these ML models were further evaluated with all variables as input variables, and no improvement in model discrimination was noted. Our study confirms the significance of feature selection in ML applications. We applied SWSFS technique to determine the optimal and minimum size of features, thus increasing the efficiency and usefulness of the models for further validation.

Our study has several strengths. First, a major barrier to the widespread use of ML models is their correct interpretation, as a true “black-box” can hardly be accepted by clinicians or decision-makers. As an additive model explainable approach, SHAP analysis is seldom employed in ML application. We utilized the SHAP values to demystify the ML, providing a “white-box” AKI prediction model that allowed a quick comprehension of the effect of a single feature on the model’s prediction. This explainable or “white-box” predictive technique may be helpful in ML transportability across hospitals. Second, the three tree-based ensemble learning algorithms demonstrated high calculating efficiency and may be adapted to certain medical working environments. Indeed, XGBoost and RF have the advantages of being trained quickly and providing reliable feature importance estimates and are increasingly emphasized as competitive alternatives to traditional regression methods. In addition, both XGBoost and RF algorithms are bootstrapping method applications, which can improve predictive power when available datasets are small. Notably, for the first time to our knowledge, we used DF to predict AKI after cardiac surgery. As an alternative to the deep learning framework, DF improves the robustness of the traditional deep learning method working on small-scale clinical data and provides an effective solution for binary classification. Third, given the importance of intraoperative factors on AKI development, we applied GBTM method to handle intraoperative longitudinal data in which the dynamic process of variables over time could be clearly observed. GBTM identifies distinctive clusters of individual trajectories within the population, determining the subgroup patients at different risk levels; it provides a critical time point for clinical decision making. By monitoring intraoperative hemodynamic parameters, it is expected to build more accurate dynamic early warning systems that can help clinicians timely identify patients at risk of postoperative AKI. Fourth, the incidence of AKI in our study was generally in line with the report from a meta-analysis [1], indicating that our patient cohort is representative of cardiac patients in general. Finally, the dataset contained only a few missing data because most of the missing values were filled in during a second manual review of EHR. Therefore, the impact of missing data on the prediction models is negligible.

However, several limitations should also be considered. First, the models were developed based on the dataset derived from a single center with patients undergoing on-pump cardiac surgery. Consequently, before the models can be implemented in clinical practice for new patients, their predictive performance would require training and evaluation on other races, nationalities, or additional datasets. Second, postoperative urine output (less than 0.5 mL/kg/h for 6 h) was not used to define AKI due to its unavailability in the majority of patients. However, we were unlikely to miss vital clinical patients as urine output may have been maintained by diuretics. Moreover, due to intensified monitoring in ICU, persistent oliguria is uncommon and transient oliguria may simply imply insufficient volume resuscitation. Third, although SHAP values can explain most traditional ML models, they cannot explain the DF algorithm. In other words, the correct interpretation of DF remains challenging. We are working on developing a more advanced algorithm for DF that could provide variable importance estimation. Fourth, we did not include traditional scoring systems to make model comparisons, because the most extensively used and robust models for AKI are those designed for AKI requiring RRT. Given the high incidence of subclinical AKI and its strong association with adverse outcomes, more advanced models should be developed to predict any-stage AKI after cardiac surgery.

Taken together, our study highlights the potential of tree-based ensemble methods in generating robust AKI prediction tools. By contrast to clinical scoring systems or biomarkers, ML algorithm is a completely data-driven prediction tool. As for clinical practice, a web-based online tool can be developed to facilitate the application of the ML model. By inputting the values of features of a patient, the online tool would estimate the probability of developing AKI. With advances in EHR, it can be easily transferred to the EHR system to calculate AKI risk by automatically reading features.

ML algorithms for clinical data analysis have revolutionized the traditional way of conducting cardiovascular research. As clinicians continue to gather significant amounts of patient data through EHR, more novel associations between specific features and AKI will be identified. Future work is ongoing on the development of more advanced ML algorithms and EHR systems for real-time adjustment of the AKI risk after cardiac surgery, which in return will optimize treatment and enhance prognosis.

Conclusions

In this study, based on the 16 important perioperative predictors, we successfully established three tree-based ML models and a nomogram model to optimize CSA-AKI risk prediction.

Availability of data and materials

The original contributions presented in the study are included in the article and additional files. Further data that support the findings of this study are available from the corresponding author on reasonable request.

Abbreviations

AKI:

Acute kidney injury

CSA-AKI:

Cardiac surgery-associated acute kidney injury

EHR:

Electronic health record

SWSFS:

Sliding windows sequential forward feature selection

XGBoost:

Extreme gradient boosting

RF:

Random forest

DF:

Deep forest

AUROC:

Area under the receiver operating characteristic curve

SHAP:

SHapley Additive explanation

Scr:

Serum creatinine

ML:

Machine learning

RRT:

Renal replacement therapy

CABG:

Coronary artery bypass grafting

CPB:

Cardiopulmonary bypass

ICU:

Intensive care unit

VIS:

Variable importance score

OOB:

Out of bag

GBTM:

Group-based trajectory modelling

CVP:

Central venous pressure

CUF:

Conventional ultrafiltration

NLR:

Neutrophil to lymphocyte ratio

UA:

Uric acid

HDL:

High-density lipoprotein

LDH:

Lactic dehydrogenase

References

  1. Neugarten J, Sandilya S, Singh B, Golestaneh L. Sex and the risk of AKI following cardio-thoracic surgery: a meta-analysis. Clin J Am Soc Nephrol. 2016;11:2113–22.

    Article  Google Scholar 

  2. James MT, Bhatt M, Pannu N, Tonelli M. Long-term outcomes of acute kidney injury and strategies for improved care. Nat Rev Nephrol. 2020;16:193–205.

    Article  Google Scholar 

  3. Pickkers P, Darmon M, Hoste E, Joannidis M, Legrand M, Ostermann M, et al. Acute kidney injury in the critically ill: an updated review on pathophysiology and management. Intensive Care Med. 2021;47:835–50.

    Article  Google Scholar 

  4. Schunk SJ, Zarbock A, Meersch M, Küllmar M, Kellum JA, Schmit D, et al. Association between urinary dickkopf-3, acute kidney injury, and subsequent loss of kidney function in patients undergoing cardiac surgery: an observational cohort study. Lancet. 2019;394:488–96.

    Article  CAS  Google Scholar 

  5. Ugwuowo U, Yamamoto Y, Arora T, Saran I, Partridge C, Biswas A, et al. Real-time prediction of acute kidney injury in hospitalized adults: implementation and proof of concept. Am J Kidney Dis. 2020;76:806–14.

    Article  CAS  Google Scholar 

  6. Thakar CV, Arrigain S, Worley S, Yared JP, Paganini EP. A clinical score to predict acute renal failure after cardiac surgery. J Am Soc Nephrol. 2005;16:162–8.

    Article  Google Scholar 

  7. Wijeysundera DN, Karkouti K, Dupuis JY, Rao V, Chan CT, Granton JT, et al. Derivation and validation of a simplified predictive index for renal replacement therapy after cardiac surgery. JAMA. 2007;297:1801–9.

    Article  CAS  Google Scholar 

  8. Mehta RH, Grab JD, O’Brien SM, Bridges CR, Gammie JS, Haan CK, et al. Bedside tool for predicting the risk of postoperative dialysis in patients undergoing cardiac surgery. Circulation. 2006;114:2208–16.

    Article  Google Scholar 

  9. Jiang W, Xu J, Shen B, Wang C, Teng J, Ding X. Validation of four prediction scores for cardiac surgery-associated acute kidney injury in chinese patients. Braz J Cardiovasc Surg. 2017;32:481–6.

    PubMed  PubMed Central  Google Scholar 

  10. Tomašev N, Glorot X, Rae JW, Zielinski M, Askham H, Saraiva A, et al. A clinically applicable approach to continuous prediction of future acute kidney injury. Nature. 2019;572:116–9.

    Article  Google Scholar 

  11. Lundberg SM, Erion GG, Lee S-I. Consistent individualized feature attribution for tree ensembles. arXiv [Preprint]. 2018. https://arxiv.org/abs/1802.03888. Accessed 12 Feb 2018.

  12. Collins GS, Reitsma JB, Altman DG, Moons KG. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMJ. 2015;350:g7594.

    Article  Google Scholar 

  13. Stevens PE, Levin A. Evaluation and management of chronic kidney disease: synopsis of the kidney disease: improving global outcomes 2012 clinical practice guideline. Ann Intern Med. 2013;158:825–30.

    Article  Google Scholar 

  14. Breiman L. Random forests. Mach Learn. 2001;45:5–32.

    Article  Google Scholar 

  15. Chen T, Guestrin C. XGBoost: A Scalable Tree Boosting System. arXiv [Preprint]. 2016. https://arxiv.org/abs/1603.02754. Accessed 9 Mar 2016.

  16. Zhou ZH, Feng J. Deep Forest: towards an alternative to deep neural networks. arXiv [Preprint]. 2017. https://arxiv.org/abs/1702.08835v1. Accessed 28 Feb 2017.

  17. Niculescu-Mizil A, Caruana R. Predicting good probabilities with supervised learning. In: Machine learning: Proceedings of the 22nd international conference. 2005. See section 4 (Qualitative Analysis of Predictions).

  18. Shapley LS. A value for n-person games. Princeton: Princeton University Press; 1953.

    Google Scholar 

  19. Nagin DS, Jones BL, Passos VL, Tremblay RE. Group-based multi-trajectory modeling. Stat Methods Med Res. 2018;27:2015–23.

    Article  Google Scholar 

  20. Cook NR. Statistical evaluation of prognostic versus diagnostic models: beyond the ROC curve. Clin Chem. 2008;54:17–23.

    Article  CAS  Google Scholar 

  21. Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Making. 2006;26:565–74.

    Article  Google Scholar 

  22. Kanji HD, Schulze CJ, Hervas-Malo M, Wang P, Ross DB, Zibdawi M, et al. Difference between pre-operative and cardiopulmonary bypass mean arterial pressure is independently associated with early cardiac surgery-associated acute kidney injury. J Cardiothorac Surg. 2010;5:71.

    Article  Google Scholar 

  23. O’Neal JB, Shaw AD, Billings FT. Acute kidney injury following cardiac surgery: current understanding and future directions. Crit Care. 2016;20:187.

    Article  Google Scholar 

  24. Gambardella I, Gaudino M, Ronco C, Lau C, Ivascu N, Girardi LN. Congestive kidney failure in cardiac surgery: the relationship between central venous pressure and acute kidney injury. Interact Cardiovasc Thorac Surg. 2016;23:800–5.

    Article  Google Scholar 

  25. Lopez MG, Shotwell MS, Morse J, Liang Y, Wanderer JP, Absi TS, et al. Intraoperative venous congestion and acute kidney injury in cardiac surgery: an observational cohort study. Br J Anaesth. 2021;126:599–607.

    Article  CAS  Google Scholar 

  26. Manning MW, Li YJ, Linder D, Haney JC, Wu YH, Podgoreanu MV, et al. Conventional ultrafiltration during elective cardiac surgery and postoperative acute kidney injury. J Cardiothorac Vasc Anesth. 2021;35:1310–8.

    Article  Google Scholar 

  27. Tseng PY, Chen YT, Wang CH, Chiu KM, Peng YS, Hsu SP, et al. Prediction of the development of acute kidney injury following cardiac surgery by machine learning. Crit Care. 2020;24:478.

    Article  Google Scholar 

  28. Esteve F, Lopez-Delgado JC, Javierre C, Skaltsa K, Carrio ML, Rodríguez-Castro D, et al. Evaluation of the PaO2/FiO2 ratio after cardiac surgery as a predictor of outcome during hospital stay. BMC Anesthesiol. 2014;14:83.

    Article  Google Scholar 

  29. Gao XP, Zheng CF, Liao MQ, He H, Liu YH, Jing CX, et al. Admission serum sodium and potassium levels predict survival among critically ill patients with acute kidney injury: a cohort study. BMC Nephrol. 2019;20:311.

    Article  Google Scholar 

  30. Zhang D, Shi L. Serum lactate dehydrogenase level is associated with in-hospital mortality in critically Ill patients with acute kidney injury. Int Urol Nephrol. 2021;53:2341–8.

    Article  CAS  Google Scholar 

  31. Kym D, Cho YS, Yoon J, Yim H, Yang HT. Evaluation of diagnostic biomarkers for acute kidney injury in major burn patients. Ann Surg Treat Res. 2015;88:281–8.

    Article  Google Scholar 

  32. Guan C, Li C, Xu L, Zhen L, Zhang Y, Zhao L, et al. Risk factors of cardiac surgery-associated acute kidney injury: development and validation of a perioperative predictive nomogram. J Nephrol. 2019;32:937–45.

    Article  CAS  Google Scholar 

  33. Vermeulen Windsant IC, Snoeijs MG, Hanssen SJ, Altintas S, Heijmans JH, Koeppel TA, et al. Hemolysis is associated with acute kidney injury during major aortic surgery. Kidney Int. 2010;77:913–20.

    Article  Google Scholar 

  34. Weedle RC, Da Costa M, Veerasingam D, Soo AWS. The use of neutrophil lymphocyte ratio to predict complications post cardiac surgery. Ann Transl Med. 2019;7:778.

    Article  CAS  Google Scholar 

  35. Wang Q, Li J, Wang X. The neutrophil-lymphocyte ratio is associated with postoperative mortality of cardiac surgery. J Thorac Dis. 2021;13:67–75.

    Article  Google Scholar 

  36. Haran C, Gimpel D, Clark H, McCormack DJ. Preoperative neutrophil and lymphocyte ratio as a predictor of mortality and morbidity after cardiac surgery. Heart Lung Circ. 2021;30:414–8.

    Article  Google Scholar 

  37. Smith LE, Smith DK, Blume JD, Linton MF, Billings FT 4th. High-density lipoprotein cholesterol concentration and acute kidney injury after cardiac surgery. J Am Heart Assoc. 2017;6:e006975.

    Article  Google Scholar 

  38. Navab M, Hama SY, Cooke CJ, Anantharamaiah GM, Chaddha M, Jin L, et al. Normal high density lipoprotein inhibits three steps in the formation of mildly oxidized low density lipoprotein: step 1. J Lipid Res. 2000;41:1481–94.

    Article  CAS  Google Scholar 

  39. Navab M, Hama SY, Anantharamaiah GM, Hassan K, Hough GP, Watson AD, et al. Normal high density lipoprotein inhibits three steps in the formation of mildly oxidized low density lipoprotein: steps 2 and 3. J Lipid Res. 2000;41:1495–508.

    Article  CAS  Google Scholar 

  40. Lee HC, Yoon HK, Nam K, Cho YJ, Kim TK, Kim WH, et al. Derivation and validation of machine learning approaches to predict acute kidney injury after cardiac surgery. J Clin Med. 2018;7:322.

    Article  Google Scholar 

  41. Saeys Y, Inza I, Larrañaga P. A review of feature selection techniques in bioinformatics. Bioinformatics. 2007;23:2507–17.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

None.

Funding

This work was supported by the National Natural Science Foundation of China (No. 8217021245 to Prof. Xin Chen, and No. 82173620 to Prof. Yang Zhao).

Author information

Authors and Affiliations

Authors

Contributions

HZ and ZW conceived and designed the study. HZ and WC obtained, organized, and cleaned the dataset. HZ, ZW, YT, XC, DY, YW, and ZY performed the data analysis. HZ, ZW, MY, and ZY wrote the manuscript. HZ and ZW revised the manuscript. YZ and XC supervised the whole process. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yang Zhao or Xin Chen.

Ethics declarations

Ethics approval and consent to participate

The study was reviewed and approved by the Ethics Committee of Nanjing First Hospital, with the need for individual patient consent waived.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Figure S1. Illustration of cascade forest structure. Figure S2. Variance expansion factor assessment of 14 variables included in the logistic regression model. Figure S3. Formation of the discovery and validation sets with acute kidney injury after cardiac surgery. Figure S4. Comparison of area under the receiver operating characteristic curves among the XGBoost, random forest, and deep forest models in the validation set. Table S1. Clinical characteristics of patients in the discovery and validation sets who did or did not develop acute kidney injury after cardiac surgery. Table S2. Five-fold cross-validation results of AUROC and accuracy of XGBoost, random forest, and deep forest model. Table S3. Classifiers’ predictive performance in the validation set.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, H., Wang, Z., Tang, Y. et al. Prediction of acute kidney injury after cardiac surgery: model development using a Chinese electronic health record dataset. J Transl Med 20, 166 (2022). https://doi.org/10.1186/s12967-022-03351-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-022-03351-5

Keywords