Skip to main content

ORIGINAL RESEARCH article

Front. Surg., 07 February 2023
Sec. Heart Surgery

Machine learning in the prediction of cardiac surgery associated acute kidney injury with early postoperative biomarkers

\r\nRui Fan,,&#x;Rui Fan1,2,†Wei Qin,&#x;Wei Qin2,†Hao ZhangHao Zhang3Lichun GuanLichun Guan4Wuwei WangWuwei Wang2Jian LiJian Li2Wen ChenWen Chen2Fuhua Huang
Fuhua Huang2*Hang Zhang
Hang Zhang4*Xin Chen,
\r\nXin Chen1,2*
  • 1School of Medicine, Southeast University, Nanjing, China
  • 2Department of Thoracic and Cardiovascular Surgery, Nanjing First Hospital, Nanjing Medical University, Nanjing, China
  • 3Department of Nephrology, Nanjing First Hospital, Nanjing Medical University, Nanjing, China
  • 4Department of Thoracic Surgery, Shanghai General Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

Purpose: To establish novel prediction models for predicting acute kidney injury (AKI) after cardiac surgery based on early postoperative biomarkers.

Patients and methods: This study enrolled patients who underwent cardiac surgery in a Chinese tertiary cardiac center and consisted of a discovery cohort (n = 452, from November 2018 to June 2019) and a validation cohort (n = 326, from December 2019 to May 2020). 43 biomarkers were screened using the least absolute shrinkage and selection operator and logistic regression to construct a nomogram model. Three tree-based machine learning models were also established: eXtreme Gradient Boosting (XGBoost), random forest (RF) and deep forest (DF). Model performance was accessed using area under the receiver operating characteristic curve (AUC). AKI was defined according to the Kidney Disease Improving Global Outcomes criteria.

Results: Five biomarkers were identified as independent predictors of AKI and were included in the nomogram: soluble ST2 (sST2), N terminal pro-brain natriuretic peptide (NT-proBNP), heart-type fatty acid binding protein (H-FABP), lactic dehydrogenase (LDH), and uric acid (UA). In the validation cohort, the nomogram achieved good discrimination, with AUC of 0.834. The machine learning models also exhibited adequate discrimination, with AUC of 0.856, 0.850, and 0.836 for DF, RF, and XGBoost, respectively. Both nomogram and machine learning models had well calibrated. The AUC of sST2, NT-proBNP, H-FABP, LDH, and UA to discriminate AKI were 0.670, 0.713, 0.725, 0.704, and 0.749, respectively. In addition, all of these biomarkers were significantly correlated with AKI after adjusting clinical confounders (odds ratio and 95% confidence interval of the third vs. the first tertile: sST2, 3.55 [2.34–5.49], NT-proBNP, 5.50 [3.54–8.71], H-FABP, 6.64 [4.11–11.06], LDH, 7.47 [4.54–12.64], and UA, 8.93 [5.46–15.06]).

Conclusion: Our study provides a series of novel predictive models and five biomarkers for enhancing the risk stratification of AKI after cardiac surgery.

Introduction

Over the past few decades, advances in surgical techniques, anesthesia, and perioperative care have all contribute to enhancing in-hospital survival for cardiac patients. However, surgical morbidity is still a leading health burden in both developed and developing countries. As one of the most serious complications, cardiac surgery-associated acute kidney injury (CSA-AKI) was reported with the rate ranging 26.0%–28.5% (1). CSA-AKI not only severely affects acute morbidity and mortality, but also affects long-term prognosis.

Currently, because of a general lack of effective treatments for AKI, early identification of the condition is an important step to minimize the incidences of AKI and its mediated adverse events. Several attempts have been made to develop AKI prediction scores such as Cleveland Clinic score, Simplified Renal Index score, and Mehta score (24). However, these models were mainly derived based on traditional clinical factors from a patient's history (e.g., age, hypertension, diabetes mellitus, baseline kidney function), which provided limited information for AKI classification. On the other hand, due to differences in demographic characteristics and comorbidities, widely adoption of these models to other races or populations is of great challenge (5, 6).

Biomarkers have been proposed to clinical practice for more than 10 years. Unfortunately, serum creatinine (Scr), the main biomarker of AKI, does not increase until 50% of the kidney function is lost, potentially leading to diagnosis and treatment delay. Many novel biomarkers were proposed to substitute Scr in the assessment of kidney function, such as NGAL, KIM-1, L-FABP, and IGFBP-7 (7). However, in practice, a single biomarker is insensitive in its prediction for AKI. Therefore, it is necessary to consider a more comprehensive biomarker-based model to improve the accuracy and robustness of AKI prediction. The early postoperative period has important clinical implications as this time period is the optimal actionable window to optimize the management of postoperative AKI. Some of biomarkers arose during this period may reflect acute physiological responses for kidney function. In this study, we aimed to develop two types of models for CSA-AKI based on early postoperative biomarkers, applying multivariate logistic regression method and machine learning (ML) algorithms. The models were developed using a dataset of 452 patients from Nanjing First Hospital.

Material and methods

Study design and participants

This study contains two separate patient cohorts from our cardiac center. For the discovery cohort, participants admitted between November 2018 and June 2019 were retrospectively obtained from the Patient Information Management Platform of Nanjing First Hospital (218.2.200.37:2356/PatientList). This is a population-based database that consisted of patient information available from electronic health record in digital format. The clinical characteristics of patients were recorded in real time by medical personnel. The prediction models generated from the discovery stage were further validated in an independent patient cohort, which were prospectively enrolled between December 2019 and May 2020. This study was approved by the Ethical Committee of Nanjing First Hospital and informed consent was obtained from patients or their legal representatives. This study followed the TRIPOD statement guidelines for reporting (8).

We enrolled patients who underwent cardiac surgery with cardiopulmonary bypass (CPB) including coronary artery bypass grafting (CABG), valve surgery, and concomitant CABG and valve surgery. Patients were excluded if they met any of the following criteria: (i) aged <18 years old; (ii) AKI, dialysis, or end-stage renal disease at or before admission; (iii) emergency surgery; (iv) preoperative acute heart failure or hemodynamic instability; (v) missing data of Scr.

Sample test

Laboratory biomarkers included cardiac biomarkers, arterial blood gas, blood biochemistry and blood cell analysis, and coagulation function. Arterial blood gases were collected immediately after intensive care unit (ICU) admission. In our center, most elective surgery can be done before 18.00 pm. Fasting venous blood samples were obtained on the morning (6:00 am) of the first postoperative day to assay cardiac biomarkers, blood biochemistry and blood cells, and coagulation function. When patients were in ICU, blood samples were collected for assessment of renal function around 6.00–7.00 am every day. When patients were transformed into normal ward, blood samples were collected for assessment of renal function around 8.00–9.00 am. All clinical samples were tested in the department of medical laboratory of Nanjing First Hospital.

End point definition

The primary outcome was postoperative any-stage AKI, which was defined according to the Kidney Disease Improving Global Outcomes (KDIGO) clinical guideline (9), specifically an acute increase in Scr over 50% within 7 days, or 26.5 μmol/l elevation within 48 h compared with the reference Scr, or presence of oliguria (urine output less than 0.5 ml/kg/hr for 6 h), or a requirement of renal replacement therapy (RRT). The Scr levels measured before surgery were used as the reference value. Estimated glomerular filtration rate (eGFR) was calculated using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) creatinine equation (10).

Model development

We developed a series of prediction models for CSA-AKI including a nomogram model and three tree-based ML models. Before developing the nomogram, we applied least absolute shrinkage and selection operator (LASSO) regression to identify a set of important biomarkers. The LASSO is a compressed estimation method based on the idea of data dimension reduction. It achieves the goal of variable selection by constructing a penalized function (λ) that compresses coefficients of irrelevant prediction variables towards zero (11). A particular advantage of this technique is that it avoids both overfitting and overestimation during model derivation. The tuning parameter (λ) was determined in the LASSO using 10-fold cross-validation on the basis of minimum criteria and a standard error. The identified biomarkers were then incorporated in a multivariate logistic regression analysis to generate a nomogram model.

We performed the following tree-based ML algorithms to develop the prediction models, which were the most popular and advanced ML methods used for the problem of classification, including eXtreme Gradient Boosting (XGBoost), RF, and deep forest (DF). XGBoost is a learning framework mainly consisting of two parts: simple decision tree and gradient boosting algorithm. XGBoost applies a second Taylor expansion on the loss function and simultaneously implements the first derivative and the second derivative (12). RF is a classic tree-based algorithm which combines multiple decision trees through majority voting (13). Both XGBoost and RF are ensemble learning methods characterized by good accuracy, robustness, and high calculating efficiency. As a more advanced deep learning approach, DF generates a multi-layer cascade forest, a structure contains many different RFs. The purpose of this design is to include different types of forests to ensure the diversity of the model. Each layer in the cascade forest receives the information processed by the previous stage and then concatenates with the original vector to be input to the next layer (14). (Supplementary Figure S1) To correctly interpret a ML model, we used Shapley Additive exPlanation (SHAP) values to explain the complex relationship between variables and outcomes in the RF model. The SHAP algorithm works based on the concept of Shapley values used in game theory (15). It provides a method to estimate the positive or negative contribution of individual feature to the overall model prediction via the GradientSHAP approximation (16).

Statistical analysis

Continuous variables are presented as means ± standard deviations or medians (interquartile range), and categorized variables as frequencies with proportions. Comparisons between non-AKI group and AKI group were carried out using the t-test, Mann-Whitney U-test, chi-square test, or Fisher's exact probability method as appropriate. All demographic characteristics and clinical outcomes were available after a second manual review of the medical records. For laboratory data, our dataset has missing values ranging 0–3.8%. The missing data were handled using multiple imputation method. All laboratory biomarkers were loge transformed for subsequent analysis.

The discrimination of prediction model was evaluated by area under the receiver operating characteristic (AUC). Model calibration was accessed using Brier score, and visualized with a 1000-resample bootstrapped calibration plot. A lower Brier score indicates superior model calibration (17). For interpreting ML model, we used SHAP summary plot to show feature importance and depicted the effect of specific feature on model output.

To further explore the association between single biomarker and AKI, we used the violin plots to show the distribution of the single biomarker among non-AKI and AKI groups. Receiver operating characteristic (ROC) curves were also formulated to determine the discriminative ability and optimal cut-off values of each biomarker (identified by the maximum Youden index). In addition, patients were divided into three group based on the tertiles of biomarker concentration distribution. Multivariate logistic regression analysis was performed to identify whether these biomarkers were predictors of AKI, independent of clinical confounders. Odds ratios (ORs) and corresponding two-sided 95% confidence intervals (CIs) were reported. Statistical analyses were performed using R (version 4.0.3) with the packages of mice, glmnet, and rms, and Python (version 3.8) with the packages of sklearn, deep-forest, and shap. A two-sided P value less than 0.05 indicated statistically significant.

Results

Characteristics of the cohorts

Overall, 452 participants, admitted between November 2018 and June 2019, comprised the discovery cohort; 326 participants, admitted between December 2019 and May 2020, comprised the validation cohort (Figure 1). The rates of AKI were 25.4% and 29.1% in the discovery and validation cohorts, respectively. In the discovery cohort, patients who developed AKI showed a significant association on the univariate analysis included advanced age (P = 0.009), more diabetes mellitus (P = 0.038), hypertension (P = 0.047), and combined surgery (P = 0.008), longer duration of CPB and aortic cross-clamping (all P < 0.001). Patients who developed AKI had more death or on RRT (P = 0.044), more pulmonary infection (P = 0.004), longer ventilation time, ICU stay and hospital stay (all P < 0.001). Similar differences were observed in the validation cohort except diabetes mellitus (P = 0.440), hypertension (P = 0.403), and death or on RRT (P = 0.496) (Table 1).

FIGURE 1
www.frontiersin.org

Figure 1. Formation of the discovery and validation cohorts of acute kidney injury after cardiac surgery.

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of patients in the discovery and validation cohorts who did or did not develop acute kidney injury.

43 early postoperative biomarkers were summarized in the discovery cohort (Table 2). Patients who developed AKI had higher indices of soluble ST2 (sST2), N terminal pro-brain natriuretic peptide (NT-proBNP), heart-type fatty acid binding protein (H-FABP), aspartate transaminase, lactic dehydrogenase (LDH), uric acid (UA), total bile acid, serum potassium, serum magnesium, or had lower indices of high-density lipoprotein (HDL), PO2, hemoglobin, hematocrit (all P < 0.05).

TABLE 2
www.frontiersin.org

Table 2. Early postoperative biomarkers among patients in the discovery cohort who did or did not develop acute kidney injury.

Variable selection

In the discovery cohort, 43 biomarkers were included in the variable selection procedure. The LASSO identified six biomarkers predisposing to AKI: sST2, NT-proBNP, H-FABP, LDH, UA, and HDL (Supplementary Figure S2). Inclusion of these six variables in a logistic regression model resulted in five variables (excluded HDL, P > 0.05) that were independently statistically significant predictors of AKI and were included in the final model (Table 3).

TABLE 3
www.frontiersin.org

Table 3. Multivariate logistic regression model for predicting acute kidney injury based on discovery cohort.

Nomogram and model performance

An AKI prediction nomogram model was constructed (Figure 2A). The nomogram demonstrated good discrimination in the discovery cohort, with a C-statistic of 0.871 (95%CI 0.834–0.908). Correspondingly, in the validation cohort, the nomogram achieved a C-statistic of 0.834 (95%CI 0.783–0.885) (Figures 2B,C). The calibration plots revealed good calibration in both cohorts (Figures 2D,E).

FIGURE 2
www.frontiersin.org

Figure 2. Biomarker-based nomogram model and model performance. (A) Nomogram to predict CSA-AKI; Receiver operating characteristic curve of the nomogram in the discovery (B) and validation (C) cohorts; 1000-resample bootstrapped calibration plot of the nomogram in the discovery (D) and validation (E) cohorts. sST2, soluble ST2; NT-proBNP, N terminal pro-brain natriuretic peptide; H-FABP, heart-type fatty acid-binding protein; LDH, lactic dehydrogenase; UA, uric acid; AUC, area under the receiver operating characteristic curve.

Ml models

We constructed ML models using XGBoost, RF, and DF algorithms with all the biomarkers as input variables. In the validation cohort, the AUCs were 0.856 (95%CI 0.813–0.899) for RF model, 0.850 (95%CI 0.805–0.894) for DF model, and 0.836 (95%CI 0.790–0.883) for XGBoost model (Figure 3A). Different accuracy metrics were displayed in Table 4. The DF model exhibited the best calibration (Brier score: 0.143), followed by RF (Brier score: 0.157) and XGBoost (Brier score: 0.182) (Figure 3B). The SHAP values were used to highlight individual contributions of the variables in the RF model. Figure 3C describes the SHAP summary plot, showing the SHAP values in order of the important variables that contribute to AKI. According to the summary plot, each dot represents one patient and the horizontal location of each dot indicates whether the effect of a feature is associated with a higher or lower risk of AKI (Figure 3D).

FIGURE 3
www.frontiersin.org

Figure 3. Model performance and SHAP analysis of machine learning models. (A) Receiver operating characteristic curve to assess discrimination of XGBoost, random forest, and deep forest models in the validation cohort; (B) Calibration plots for XGBoost, random forest, and deep forest models in the validation cohort (5 bins); (C) SHAP summary plots of top 20 features in the random forest. Feature importance ranked by SHAP values; (D) Dot estimation for each feature's attribution value to model output. The higher SHAP value of features, the higher risk of acute kidney injury. SHAP, Shapely Additive exPlanations; H-FABP, heart-type fatty acid-binding protein; NT-proBNP, N terminal pro-brain natriuretic peptide; sST2, soluble ST2; UA, uric acid; LDH, lactic dehydrogenase; HDL, high density lipoprotein; AST, aspartate transaminase; ALT, alanine aminotransferase; TBA, total bile acid.

TABLE 4
www.frontiersin.org

Table 4. Accuracy metrics of three machine models.

Association between single biomarker and AKI

We further investigated the association of sST2, NT-proBNP, H-FABP, LDH, and UA with AKI after pooling the two cohorts together. The violin plots showed that patients in the AKI group had significantly higher levels of sST2, NT-proBNP, H-FABP, LDH, and UA than patients without AKI (Supplementary Figure S3). ROC curves were used to test the overall discriminative ability of these five biomarkers for AKI. We observed that the AUCs of sST2, NT-proBNP, H-FABP, LDH, and UA to discriminate AKI were 0.670 (95%CI, 0.627–0.713), 0.713 (95%CI, 0.669–0.758), 0.725 (95%CI, 0.686–0.765), 0.704 (95%CI, 0.664–0.745), and 0.749 (95%CI, 0.710–0.788). Meanwhile, the optimal cut-off values and corresponding specificity and sensitivity were determined (Figure 4). Moreover, we divided patients into three groups based on the tertiles of biomarker concentration distribution. With the first tertile used as the reference category, increased levels of these five biomarkers were associated with a higher risk of AKI. These associations persisted after multifactorial adjustment, including baseline eGFR (Table 5).

FIGURE 4
www.frontiersin.org

Figure 4. Receiver operating characteristic curve for sST2, NT-proBNP, H-FABP, LDH, and UA to predict CSA-AKI. AKI, acute kidney injury; AUC, area under the receiver operating characteristic curve; CI, confidence interval; sST2, soluble ST2; NT-proBNP, N terminal pro-brain natriuretic peptide; H-FABP, heart-type fatty acid-binding protein; LDH, lactic dehydrogenase; UA, uric acid.

TABLE 5
www.frontiersin.org

Table 5. Multivariate regression analysis on the association between sST2, NT-proBNP, H-FABP, LDH, and UA with acute kidney injury in patients of the total cohort.

Discussion

In this study, using early postoperative biomarkers, we constructed a series of models for enhancing risk stratification of AKI after cardiac surgery. The discrimination of these models was clinically satisfactory, with AUCs ranging 0.834–0.856 in the validation cohort. In addition, both LASSO and SHAP analysis identified sST2, NT-proBNP, H-FABP, LDH, and UA as the most influential predictors of AKI.

To the best of our knowledge, the biomarker-based nomogram model in this study consists of four cardiac biomarkers (sST2, NT-proBNP, H-FABP, and LDH) and a metabolism related biomarker (UA), which have not been included in previous studies. Overall, they are novel biomarkers for predicting AKI. sST2 is initially found to be a member of the interleukin-1 (IL-1) receptor family that serves as a decoy for IL-33, regulating immune and inflammatory responses. Emerging evidence from epidemiological studies supported that elevated serum sST2 levels were associated with mortality and adverse clinical outcomes in patients suffering from heart failure, coronary artery disease, arrhythmia, and stroke (1821). Lobdell et al. demonstrated a significant association of high preoperative concentrations of sST2 as a prognostic indicator of AKI among patients undergoing CABG (22). We observed a similar statistical difference in early postoperative phase. These findings may promote the extension of the predictive ability of sST2 from cardiovascular disease to kidney injury.

As an intracellular transport protein responsible for transporting free fatty acid in cardiomyocytes, H-FABP has been recognized as a diagnostic and prognostic marker for acute coronary syndrome (2324). Schaub et al. explored the predictive value of H-FABP for AKI at four time points perioperatively (25). They found that patients who developed AKI had higher H-FABP than those who did not. Published literatures gave plausible explanations of the relationship between H-FABP and AKI in cardiac patients (2527). First, high H-FABP levels were only observed in patients who had received on-pump cardiac surgery, but not in off-pump procedures. This indicated that ischemia-reperfusion injury and inflammatory might play a role in the elevation of H-FABP. Second, H-FABP was an indicator of hemodynamic instability after cardiac surgery. Third, elevated H-FABP levels were also found in patients with elevated venous pressure or venous congestion, which in return would affect kidney circulation and lead to “congestive kidney failure”.

Like sST2 and H-FABP, NT-proBNP and LDH are also cardiac biomarkers associated with development of AKI. NT-proBNP is mainly cleared by renal excretion, and its change has a significant modifying effect on kidney function decline (28). Besides, postoperative cardiac insufficiency, which is predicted by NT-proBNP, promotes AKI via hypotension. LDH anomaly is known to often occur in certain diseases such as renal failure, cardiovascular events, and hepatic damage. Therefore, LDH was usually measured to detect tissue damage as well as patient's overall disease severity (29). Moreover, in the CPB setting, the elevated LDH in the immediate postoperative period may be an indicator of CPB-induced hemolysis, which is associated with the development of AKI (30). UA is the final product of endogenous and exogenous purine metabolism. A population-based cohort study has shown that high serum UA levels were positively associated with elevated levels of pro-inflammatory cytokines (e.g., interleukin-6, high sensitivity C-reactive protein, tumor necrosis factor-α), which was considered as the central components of the pathogenesis of AKI (31). Kidney function may benefit from UA lowering therapy; therefore, UA may serve as a novel potential target for AKI prevention (32).

The advantage of our study included the use of SHAP values to uncover the black box of MLthe DF algorithm to predict CSA-AKI. As one of the advanced tree-based learning methods, the DF model demonstrated better predictive power than the conventional ML model. The DF boosts predictive information by integrating multiple RFs, providing an effective option to investigate binary problem and improving the robustness of the standard deep learning methods working on small-scale data. Compared with logistic regression, machine learning can include much more variables and resolve nonlinear interactions. This can improve the effectiveness of prediction. Even though AUC of machine learning we used was 0.856, which was similar to the logistic regression, it was higher than AUC of logistic regression (0.834). This means machine learning showed better performance than traditional analysis. To feasibility, these biomarkers can be collected easily from venous blood and tested in labs easily.

Several limitations in this study should also be noted. First, the models were derived from a single-center dataset. Therefore, before models can be implemented in clinical practice, their predictive power needed to be validated in external datasets. Second, we did not consider innovative biomarkers such as renal tubule-associated biomarkers (e.g., NGAL, KM-1, MMP-7) or biomarkers of inflammation (e.g., IL-6, IL-10, and TNF-α), which may also be ideal predictors of AKI. But they remained poorly studied as a biomarker-based prediction tool in general or cardiac surgery. Third, we did not include traditional scoring systems (e.g., Cleveland Clinic score, Simplified Renal Index score, Mehta score) to make model comparisons because most of these models were specially designed to predict AKI requiring RRT. Given the high incidence of mild or moderate AKI and its strong association with adverse outcomes, more efforts should be made to predict any-stage AKI after cardiac surgery.

Conclusion

We successfully constructed a nomogram and three tree-based ML models to predict CSA-AKI based on early postoperative biomarkers. Our study identified five important biomarkers (sST2, NT-proBNP, H-FABP, LDH, and UA) associated with CSA-AKI.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

This study was approved by the Ethical Committee of Nanjing First Hospital and informed consent was obtained from patients or their legal representatives.

Author contributions

RF: collected data, performed the study, wrote the manuscript. WQ, LG and HZ: designed, first drafting of the manuscript. WW, JL and WC: analyzed and interpreted data. FH: reviewed the manuscript. HZ and XC: critical review of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (Number, 8217021245) and the Science and Technology Plan of Changzhou (Number, CE20205047).

Conflict of interest

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

Publisher's note

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

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fsurg.2023.1048431/full#supplementary-material.

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(12):2113–22. doi: 10.2215/CJN.03340316

PubMed Abstract | CrossRef Full Text | Google Scholar

2. 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(1):162–8. doi: 10.1681/ASN.2004040331

PubMed Abstract | CrossRef Full Text | Google Scholar

3. 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(16):1801–9. doi: 10.1001/jama.297.16.1801

PubMed Abstract | CrossRef Full Text | Google Scholar

4. 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(21):2208–16. quiz 2208. doi: 10.1161/CIRCULATIONAHA.106.635573

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Che M, Wang X, Liu S, Xie B, Xue S, Yan YC, et al. A clinical score to predict severe acute kidney injury in Chinese patients after cardiac surgery. Nephron. (2019) 142(4):291–300. doi: 10.1159/000499345

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ranucci M, Aloisio T, Cazzaniga A, Dedda UD, Gallazzi C, Pistuddi V. Validation of renal-risk models for the prediction of non-renal replacement therapy cardiac surgery-associated acute kidney injury. Int J Cardiol. (2018) 272:49–53. doi: 10.1016/j.ijcard.2018.07.114

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Beker BM, Corleto MG, Fieiras C, Musso CG. Novel acute kidney injury biomarkers: their characteristics, utility and concerns. Int Urol Nephrol. (2018) 50(4):705–13. doi: 10.1007/s11255-017-1781-x

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Collins GS, Reitsma JB, Altman DG, Moons KG. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. Br Med J. (2015) 350:g7594. doi: 10.1136/bmj.g7594

CrossRef Full Text | Google Scholar

9. Kellum JA, Lameire N. Diagnosis, evaluation, and management of acute kidney injury: a KDIGO summary (part 1). Crit Care. (2013) 17(1):204. doi: 10.1186/cc11454

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF, Harold IF, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. (2009) 150(9):604–12. doi: 10.7326/0003-4819-150-9-200905050-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Tibshirani R. Regression shrinkage and selection via the Lasso. J R Statist Soc B. (1996) 73(1):273–82. http://www.jstor.org/stable/2346178

Google Scholar

12. Ismail BM, Faisal S. Bioactive molecule prediction using extreme gradient boosting. Molecules. (2016) 21(8):983. doi: 10.3390/molecules2108983

CrossRef Full Text | Google Scholar

13. Jones FC, Plewes R, Murison L, MacDougall MJ, Sinclair S, Davies C, et al. Random forests as cumulative effects models: a case study of lakes and rivers in Muskoka, Canada. J Environ Manage. (2017) 201:407–24. doi: 10.1016/j.jenvman.2017.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Guo Y, Liu S, Li Z, Shang X. BCDForest: a boosting cascade deep forest model towards the classification of cancer subtypes based on gene expression data. BMC Bioinformatics. (2018) 19(Suppl 5):118. doi: 10.1186/s12859-018-2095-4

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Eisenman RL. A profit-sharing interpretation of Shapley value for N-person games. Behav Sci. (1967) 12(5):396–8. doi: 10.1002/bs.3830120506

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hathaway QA, Roth SM, Pinti MV, Sprando DC, Kunovac A, Durr AJ, et al. Machine-learning to stratify diabetic patients using novel cardiac biomarkers and integrative genomics. Cardiovasc Diabetol. (2019) 18(1):78. doi: 10.1186/s12933-019-0879-0

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Cook NR. Statistical evaluation of prognostic versus diagnostic models: beyond the ROC curve. Clin Chem. (2008) 54(1):17–23. doi: 10.1373/clinchem.2007.096529

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Broch K, Ueland T, Nymo SH, Kjekshus J, Hulthe J, Muntendam P, et al. Soluble ST2 is associated with adverse outcome in patients with heart failure of ischaemic aetiology. Eur J Heart Fail. (2012) 14(3):268–77. doi: 10.1093/eurjhf/hfs006

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Richards AM, Di Somma S, Mueller T. ST2 In stable and unstable ischemic heart diseases. Am J Cardiol. (2015) 115(7 Suppl):48b–58b. doi: 10.1016/j.amjcard.2015.01.041

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Vílchez JA, Pérez-Cuellar M, Marín F, Gallego P, Manzano-Fernandez S, Valdés M, et al. sST2 levels are associated with all-cause mortality in anticoagulated patients with atrial fibrillation. Eur J Clin Invest. (2015) 45(9):899–905. doi: 10.1111/eci.12482

CrossRef Full Text | Google Scholar

21. Wolcott Z, Batra A, Bevers MB, Sastre C, Khoury J, Sperling M, et al. Soluble ST2 predicts outcome and hemorrhagic transformation after acute stroke. Ann Clin Transl Neurol. (2017) 4(8):553–63. doi: 10.1002/acn3.435

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lobdell KW, Parker DM, Likosky DS, Rezaee M, Ballmoons MW, Alam SS, et al. Preoperative serum ST2 level predicts acute kidney injury after adult cardiac surgery. J Thorac Cardiovasc Surg. (2018) 156(3):1114–23.e1112. doi: 10.1016/j.jtcvs.2018.03.149

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Roffi M, Patrono C, Collet JP, Mueller C, Valgimigli M, Andreotti F, et al. 2015 ESC guidelines for the management of acute coronary syndromes in patients presenting without persistent ST-segment elevation: task force for the management of acute coronary syndromes in patients presenting without persistent ST-segment elevation of the European society of cardiology (ESC). Eur Heart J. (2016) 37(3):267–315. doi: 10.1093/eurheartj/ehv320

PubMed Abstract | CrossRef Full Text | Google Scholar

24. O'Donoghue M, de Lemos JA, Morrow DA, Murphy SA, Buros JL, Cannon CP, et al. Prognostic utility of heart-type fatty acid binding protein in patients with acute coronary syndromes. Circulation. (2006) 114(6):550–7. doi: 10.1161/CIRCULATIONAHA.106.641936

CrossRef Full Text | Google Scholar

25. Schaub JA, Garg AX, Coca SG, Testani JM, Shlipak MG, Eikelboom J, et al. Perioperative heart-type fatty acid binding protein is associated with acute kidney injury after cardiac surgery. Kidney Int. (2015) 88(3):576–83. doi: 10.1038/ki.2015.104

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chowdhury UK, Malik V, Yadav R, Seth S, Ramakrishnan L, Kalaivani M, et al. Myocardial injury in coronary artery bypass grafting: on-pump versus off-pump comparison by measuring high-sensitivity C-reactive protein, cardiac troponin I, heart-type fatty acid-binding protein, creatine kinase-MB, and myoglobin release. J Thorac Cardiovasc Surg. (2008) 135(5):1110–9, 1119.e1111-1110. doi: 10.1016/j.jtcvs.2007.12.029

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Jiang D, Shen M, Yuan X, Wang M, Li S, Jiang W, et al. Serum heart-type fatty acid-binding protein as a predictor for the development of sepsis-associated acute kidney injury. Expert Rev Mol Diagn. (2019) 19(8):757–65. doi: 10.1080/14737159.2019.1642750

PubMed Abstract | CrossRef Full Text | Google Scholar

28. McCallum W, Tighiouart H, Kiernan MS, Huggins GS, Sarnak MJ. Relation of kidney function decline and NT-proBNP with risk of mortality and readmission in acute decompensated heart failure. Am J Med. (2020) 133(1):115–22.e112. doi: 10.1016/j.amjmed.2019.05.047

PubMed Abstract | CrossRef Full Text | Google Scholar

29. 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(11):2341–8. doi: 10.1007/s11255-021-02792-z

PubMed Abstract | CrossRef Full Text | Google Scholar

30. 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(10):913–20. doi: 10.1038/ki.2010.24

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Lyngdoh T, Marques-Vidal P, Paccaud F, Preisig M, Waeber G, Bochud M, et al. Elevated serum uric acid is associated with high circulating inflammatory cytokines in the population-based colaus study. PLoS One. (2011) 6(5):e19901. doi: 10.1371/journal.pone.0019901

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Ejaz AA, Dass B, Lingegowda V, Shimada M, Beaver TM, Ejaz NI, et al. Effect of uric acid lowering therapy on the prevention of acute kidney injury in cardiovascular surgery. Int Urol Nephrol. (2013) 45(2):449–58. doi: 10.1007/s11255-012-0192-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: acute kidney injury, cardiac surgery, biomarker, nomogram, machine learning, random forest

Citation: Fan R, Qin W, Zhang H, Guan L, Wang W, Li J, Chen W, Huang F, Zhang H and Chen X (2023) Machine learning in the prediction of cardiac surgery associated acute kidney injury with early postoperative biomarkers. Front. Surg. 10:1048431. doi: 10.3389/fsurg.2023.1048431

Received: 19 September 2022; Accepted: 11 January 2023;
Published: 7 February 2023.

Edited by:

Robert Jeenchen Chen, Stanford University, United States

Reviewed by:

Haiqing Li, Shanghai Jiao Tong Universtity School of Medicine, China
Luis Gustavo Modelli De Andrade, Sao Paulo State University, Brazil

© 2023 Fan, Qin, Zhang, Guan, Wang, Li, Chen, Huang, Zhang and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Fuhua Huang huangfuhua@sina.com Hang Zhang zhanghangnmu@163.com Xin Chen stevecx@njmu.edu.cn

These authors have contributed equally to this work and share first authorship

Specialty Section: This article was submitted to Heart Surgery, a section of the journal Frontiers in Surgery

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