Next Article in Journal
De Novo Large Deletions in the PHEX Gene Caused X-Linked Hypophosphataemic Rickets in Two Italian Female Infants Successfully Treated with Burosumab
Next Article in Special Issue
Video-Based versus On-Site Neonatal Pain Assessment in Neonatal Intensive Care Units: The Impact of Video-Based Neonatal Pain Assessment in Real-World Scenario on Pain Diagnosis and Its Artificial Intelligence Application
Previous Article in Journal
Fuzzy Logic-Based System for Identifying the Severity of Diabetic Macular Edema from OCT B-Scan Images Using DRIL, HRF, and Cystoids
Previous Article in Special Issue
Applicability of American College of Radiology Appropriateness Criteria Decision-Making Model for Acute Appendicitis Diagnosis in Children
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploiting Machine Learning Technologies to Study the Compound Effects of Serum Creatinine and Electrolytes on the Risk of Acute Kidney Injury in Intensive Care Units

1
Graduate Institute of Biomedical Electronics and Bioinformatics, National Taiwan University, Taipei City 10617, Taiwan
2
Graduate Institute of Networking and Multimedia, National Taiwan University, Taipei City 10617, Taiwan
3
Department of Computer Science and Information Engineering, National Kaohsiung University of Science and Technology, Kaohsiung City 807618, Taiwan
4
Department of Emergency Medicine, National Taiwan University Hospital, Taipei City 10002, Taiwan
5
Department of Computer Science and Information Engineering, National Taiwan University, Taipei City 10617, Taiwan
*
Authors to whom correspondence should be addressed.
Submission received: 25 June 2023 / Revised: 27 July 2023 / Accepted: 29 July 2023 / Published: 31 July 2023

Abstract

:
Assessing the risk of acute kidney injury (AKI) has been a challenging issue for clinicians in intensive care units (ICUs). In recent years, a number of studies have been conducted to investigate the associations between several serum electrolytes and AKI. Nevertheless, the compound effects of serum creatinine, blood urea nitrogen (BUN), and clinically relevant serum electrolytes have yet to be comprehensively investigated. Accordingly, we initiated this study aiming to develop machine learning models that illustrate how these factors interact with each other. In particular, we focused on ICU patients without a prior history of AKI or AKI-related comorbidities. With this practice, we were able to examine the associations between the levels of serum electrolytes and renal function in a more controlled manner. Our analyses revealed that the levels of serum creatinine, chloride, and magnesium were the three major factors to be monitored for this group of patients. In summary, our results can provide valuable insights for developing early intervention and effective management strategies as well as crucial clues for future investigations of the pathophysiological mechanisms that are involved. In future studies, subgroup analyses based on different causes of AKI should be conducted to further enhance our understanding of AKI.

1. Introduction

Acute kidney injury (AKI) is a condition frequently encountered in medical care [1]. The underlying pathophysiological processes of AKI ultimately lead to a decline in renal function. As a result, the patients suffer from the accumulation of waste products, an imbalance of electrolytes, and a widespread inflammatory response that affects organs beyond the kidneys [2]. According to a recent study, 20% to 50% of the patients in an intensive care unit (ICU) suffered from AKI [3]. Therefore, how to assess the risk of AKI is a critical issue for clinicians in an ICU [4]. However, several early signs of AKI, including edema, hypertension, and oliguria, are non-specific. Therefore, the current practice only monitors the level of serum creatinine and the volume of urine output in order to assess the risk of AKI [5,6].
Due to the observation above, scientists have been investigating the physiological signs that may be associated with the development of AKI. Leaf et al. conducted a review of the pathophysiology of dysregulated mineral metabolism, specifically focusing on calcium, phosphate, parathyroid hormone, and vitamin D metabolites in the context of AKI [7]. A review conducted by Yokota et al. found that the most common comorbidities associated with AKI in elderly patients included respiratory failure, cardiovascular disease, hypertension, diabetes, surgical complications, and liver disease [8].
As the kidneys play a crucial role in regulating the balance of calcium, phosphate, and magnesium, it is conceivable that an imbalance of serum electrolytes may be associated with the development of AKI. In this respect, a previous study reported that acute phosphate nephropathy was an early condition of AKI and might subsequently progress to chronic renal failure [9]. Furthermore, a number of studies were conducted to investigate how the levels of serum electrolytes, including chloride, phosphorus, magnesium, potassium, sodium, and calcium, were associated with the development of AKI [10,11,12]. Suetrong et al. observed a linear correlation between the concentration of serum chloride and the development of AKI among sepsis/septic shock patients [13]. Marttinen et al. reported a similar result and showed that the temporal chloride level was associated with an increased risk of AKI [14]. The work by Moon et al. revealed that a high level of serum phosphorus increased the risk of AKI [15]. Cheungpasitporn et al. showed that both hypomagnesemia and hypermagnesemia led to an increased risk of in-hospital AKI [16]. Thongprayoon et al. observed a U-shaped association between the level of serum ionized calcium and in-hospital AKI. Furthermore, both hypocalcemia and hypercalcemia were reported to be associated with an increased risk of hospital-acquired AKI [17,18], and Chen et al. discovered that abnormal levels of serum sodium or potassium before an AKI diagnosis were more likely to lead to AKI progression and a poor prognosis [19]. Nevertheless, Yessayan et al. reported that the concentration of hyperchloremia and the onset of AKI within 72 h of admission were not correlated [20]. Finally, Morooka et al. divided pediatric patients into three groups based on their serum magnesium values and investigated the association between magnesium levels and outcomes [21].
In addition to the studies addressed above, the latest trend is to exploit various machine learning algorithms, including artificial neural networks [22], support vector machines (SVMs) [23], Bayesian networks [24], random forests (RFs) [25], etc., to predict incidences of AKI, and Song et al. reviewed how the conventional logistic regression (LR) and various machine learning methods performed in this respect [26]. A representative study was conducted by Tomasev et al. [27]. In their study, the authors employed a recurrent neural network to build their prediction models based on a cohort of 703,782 cases collected from the medical facilities of the U.S. Department of Veterans Affairs.
Though the effects of several serum electrolytes on the development of AKI have been well reported, a comprehensive investigation into how these serum electrolytes interact in the context of the development of AKI has not been conducted [28,29,30,31,32,33,34]. It is conceivable that such studies can provide crucial clues for developing new clinical guidelines to assess the risk of AKI. Accordingly, we initiated this study aiming not only to illustrate how these factors interact with each other but also to provide new insights for developing new clinical practices. Our analyses focused on ICU patients who had no prior history of AKI and were free of AKI-related comorbidities such as diabetes and hypertension as well as common causes of AKI such as hypovolemia and heart failure. By focusing on this group of patients, we were able to eliminate the confounding influences of these conditions and examine the associations between the levels of serum electrolytes and renal function in a more controlled manner.
In this study, we exploited decision tree (DT) models [28,29,30] and RF models [31,32]. Compared to the other commonly exploited machine learning models, such as SVMs [23] and deep neural networks (DNNs) [22], DT and RF models are favorites in many applications due to the interpretable decision rules exhibited by these models. Figure 1 shows a DT structure that summarizes the main results of this study. A user can figure out the decision rules by traversing the tree structure from the root node, which is at the top of the structure and colored yellow. They can proceed by following the branch originating from the root node that matches the condition of the case. The path ends at one of the leaf nodes at the bottom level of the tree. The “n+” and “n” symbols at each node denote the number of positive cases and the number of negative cases, respectively, in our study cohort that met the criteria specified along the path from the root node to this particular node. If a path ends at a red node, the prediction is positive. On the other hand, if a path ends at a green node, the prediction is negative. Based on these interpretable decision rules, physicians can have a comprehensive understanding of how these key factors interact with each other and develop new clinical guidelines accordingly. On the other hand, due to the non-linear transformations and the large number of coefficients involved in the prediction process, it is essentially impossible for a user to interpret the mathematics equations that an SVM or DNN model follows to make a prediction.

2. Materials and Methods

2.1. Study Cohort

Our study cohort was extracted from the Medical Information Mart for Intensive Care (MIMIC)-IV database, version 1.0, published in March 2021 [33,34]. The MIMIC database has been carefully de-identified to protect patient privacy. Its use for research purposes has been approved by the institutional review boards of the Massachusetts Institute of Technology (Protocol No. 0403000206) and Beth Israel Deaconess Medical Center (Protocol No. 2001-P-001699/14). These approvals indicate that the appropriate ethical considerations have been taken into account to ensure the responsible and lawful use of the database for research purposes.
Figure 2 shows the flow that we followed to generate our study cohort. Initially, the dataset contained 256,878 clinical records collected at the emergency department and the intensive care unit between 2008 and 2019. According to the 2012 Kidney Disease: Improving Global Outcomes (KDIGO) recommendation statements [35,36,37,38], AKI is defined by any of the following criteria: (1) an increase in the level of serum creatinine by 0.3 mg/dL (26.5 µmol/L) or more within 48 h or (2) an increase in the level of serum creatinine to 1.5 times the baseline level within 7 days. As the guideline requires two readings of the serum creatinine level and our study focused on patients in ICUs, 205,482 records in the database were excluded due to a lack of required information after admission into ICUs. As a result, only 51,396 records, all of which corresponded to the first available data after ICU admission, were included for subsequent analyses.
Since one patient could be admitted into the ICU more than one time, for a patient who had suffered from AKI, we included only the record corresponding to his/her stay in the ICU during which the patient suffered from AKI the first time. On the other hand, for a patient who had never suffered from AKI, we included only the record corresponding to his/her first stay in the ICU. As a result, only 41,878 records corresponding to 41,878 individual cases remained. In the next step, we employed the criteria provided in Table 1 to exclude those patients whose medical records showed AKI-related comorbidities [39] so that the interferences from other factors such as renal impairment, cardiac failure, diabetes, and electrolyte imbalances would be avoided. After this step, only 17,085 cases remained in the dataset. Finally, we employed the following excluding criteria to further screen the dataset: (1) the record of the case did not include all the readings listed in Table 2; (2) one or more readings in the record were in the highest 0.1% or the lowest 0.1% of the distributions; and (3) one or more readings for the case were not made within 168 h of admission. In the end, our study cohort contained 550 AKI-positive cases and 12,152 AKI-negative cases. A demographic analysis of the study cohort is presented in Table 2.
Etiologically, the causes of AKI can be classified into three broad categories: pre-renal azotemia, intrinsic renal parenchymal damage, and post-renal obstruction. Tailoring treatment plans according to the specific causes of renal injury are crucial for improving patient outcomes. For instance, hypovolemia, often diagnosed by assessing a fluid status imbalance, insufficient renal perfusion, or inferior vena cava collapse, is a common clinical presentation associated with pre-renal azotemia. On the other hand, post-renal injury occurs when the urinary tract is partially or completely blocked due to functional or structural derangements anywhere from the renal pelvis to the tip of the urethra. Since the treatment plans for post-renal AKI patients are significantly different from the plans for non-post-renal AKI patients, we classified the AKI patients in our study cohort into two categories: post-renal AKI and non-post-renal AKI. According to several previous studies, the incidences of post-renal AKI accounted for less than 5% of all AKI cases [1,40,41]. In our study cohort, 24 out of 550 AKI cases, i.e., 4.4%, were post-renal, and the percentage was in line with the previous studies. Supplementary Table S1 shows the ICD-9 and ICD-10 codes employed to identify post-renal AKI cases. Table 3 shows the statistics of the post-renal AKI patents and non-post-renal AKI patients with respect to the features listed in Table 2.
Table 2. Demographic analysis of the study cohort.
Table 2. Demographic analysis of the study cohort.
Feature550 Cases with AKI
(Mean ± SD)
12,152 Cases without AKI
(Mean ± SD)
p-Value
Age (years)65.68 ± 14.6960.34 ± 17.67p < 0.001 *
Gender p < 0.001 *
Male (%)349 (63.45%)6757 (55.60%)
Female (%)201 (36.55%)5395 (44.40%)
Serum
BUN (mg/dL)26.74 ± 15.3918.06 ± 8.90p < 0.001 *
Creatinine (mg/dL)1.36 ± 0.640.86 ± 0.26p < 0.001 *
Chloride (mEq/L)110.37 ± 6.60107.39 ± 5.28p < 0.001 *
Potassium (mEq/L)4.79 ± 0.754.47 ± 0.63p < 0.001 *
Sodium (mEq/L)142.81 ± 5.77141.23 ± 4.59p < 0.001 *
Magnesium (mg/dL)2.53 ± 0.522.28 ± 0.44p < 0.001 *
Phosphorus (mg/dL)4.40 ± 1.343.80 ± 0.93p < 0.001 *
Non-ionized calcium (mg/dL)8.76 ± 0.738.73 ± 0.710.346
The symbol * indicates statistical significance. For categorical variables, the p-values were calculated based on the χ2 test [42,43]. For continuous variables, the p-values were calculated based on the t-test [42,43]. SD represents standard deviation.

2.2. Machine Learning Models

As mentioned earlier, we used DT and RF models in order to investigate the compound impacts of two or more factors and provide a clear picture of how these factors interact with each other. In particular, we focused on the compound effects of serum creatinine, BUN, and the 6 serum electrolytes listed in Table 2. The serum creatinine and BUN were included because in medical practice the levels of serum creatinine and BUN as well as the BUN-to-creatinine ratio are measured to clarify different types of renal function impairment, including pre-renal azotemia, intrinsic renal parenchymal disease, and post-renal obstruction. The 6 serum electrolytes listed in Table 2 were included because previous studies had reported their associations with the development of AKI.
In order to address the needs in different clinical scenarios, we generated prediction models with varying levels of sensitivity and examined the prediction rules embedded in these models. In this respect, we set the parameters of the machine learning packages to various combinations and then employed a 5-fold cross-validation [22] to evaluate the levels of sensitivity delivered by the prediction models generated with these alternative parameter settings. In the 5-fold cross-validation process, the study cohort was randomly and evenly partitioned into 5 subsets. For each combination of parameter settings, every subset was employed to evaluate the prediction models generated with the other 4 subsets. Then, the evaluation results of these 5 subsets were collected to calculate the performance data, i.e., the sensitivity, specificity, positive predictive value (PPV), etc., corresponding to this particular parameter combination. Supplementary Table S2 shows the software packages employed to generate the DT and RF models as well as the alternative parameter settings employed to generate the prediction models in the 5-fold cross-validation process. In this respect, we tried a large number of possible parameter combinations in order to generate prediction models that delivered sensitivity at the levels of 0.95 and 0.80. Furthermore, as we had only 550 positive cases in our study cohort, we employed the 5-fold cross-validation process instead of the 10-fold cross-validation process, which may be more commonly used in machine learning research, so that each partition would contain a good number of positive cases.

3. Results

As mentioned above, in order to address the needs in different clinical scenarios, we generated prediction models with varying levels of sensitivity. In the subsequent discussions, we will focus on the prediction models with sensitivity at the levels of 0.95 and 0.80. Table 4 summarizes the performances of the DT, RF, and LR models observed during the 5-fold cross-validation procedure. The performances of the LR models were included to provide a reference because LR models are widely employed in biomedical research communities. Detailed performance data are presented in Supplementary Table S3.
The performance data in Table 4 reveal that with respect to the specificity, the positive predictive value (PPV), the relative risk, and the area under the receiver operating characteristic curve (AUC), the DT model that delivered sensitivity at the level of 0.95 performed significantly superior to the RF model that delivered the same level of sensitivity. It was also observed that the RF model that delivered sensitivity at the level of 0.80 performed marginally superior to the rival DT model in terms of specificity, PPV, and relative risk but performed inferior to the rival DT model in terms of AUC. Based on these observations, we concluded that the overall performance of the DT models was superior to that of the RF models. Therefore, in the subsequent discussions, we will focus on the DT models and the decision rules embedded in the models.
Figure 3a,b show the DT models generated by feeding the entire study cohort into the decision tree package with the combinations of parameters cp and prior set to (0.005 and 0.5835) and (0.01 and 0.744), respectively. According to the 5-fold cross-validation addressed above, with cp and prior set to these two combinations, the generated DT models should deliver sensitivity at the levels of 0.80 and 0.95, respectively. One interesting observation regarding the DT model shown in Figure 3a is that the model predicts a patient with a serum creatinine level higher than 1.25 mg/dL to be at high risk. This prediction rule comes very close to the serum creatinine level of 1.3 mg/dL commonly used by physicians to determine whether a patient is at high risk of progression to AKI. It is also observed that the DT model shown in Figure 3b predicts a patient with a serum creatinine level higher than 0.95 mg/dL to be at high risk. This observation implies that 0.95 mg/dL can be employed as an alternative threshold if the physician wants to increase the sensitivity of his/her medical judgment.
The DT model shown in Figure 3a further reveals that for a patient with a serum creatinine level between 0.95 and 1.25 mg/dL, his/her level of serum magnesium can be used as a warning sign. If the reading is higher than 2.45 mg/dL, the patient is at high risk. If not, we should further examine his/her level of serum chloride. If the patient’s level of serum chloride is over 106.5 mEq/L, the patient is at high risk.
The blue polygons in Figure 3a,b encircle the structure shared by these two DT models. According to the shared structure, for a patient with a serum creatinine level between 0.75 and 0.95 mg/dL, we should further examine his/her levels of serum magnesium and chloride. A patient is at high risk if (1) his/her level of serum chloride is higher than 113.5 mEq/L or (2) his/her level of serum chloride is between 105.5 and 113.5 mEq/L and his/her level of serum magnesium is higher than 2.35 mg/dL. Finally, since only a very limited number of positive cases in our study cohort met the criteria defined by the lower right parts of the tree structures in Figure 3a,b, we should be able to ignore the corresponding decision rules. In summary, the structures of the two DT models shown in Figure 3 illustrate that the levels of serum creatinine, chloride, and magnesium are the three major factors associated with the development of AKI. Though the level of serum phosphorus is present in these DT models, the nodes corresponding to the level of serum phosphorus are located in the lower parts of the structures, which implies that these nodes play less significant roles in the decision rules.

4. Discussion

As of today, the clinical practice to assess the risk of AKI is based on the 2012 KDIGO Clinical Practice Guideline for Acute Kidney Injury, which monitors only the level of serum creatinine and the volume of urine output. Since AKI could lead to many complications and even fatality, identifying the risk factors of AKI and exploiting machine learning technologies to predict AKI incidences have attracted a lot of attention in biomedical research communities. In this respect, several serum electrolytes have been reported to be associated with the development of AKI. Nevertheless, the compound effects of serum creatinine, BUN, and clinically relevant serum electrolytes have yet to be thoroughly investigated. With this observation, we initiated this study aiming not only to illustrate how these factors interact with each other but also to provide new insights for developing new clinical practices for assessing AKI risk. In particular, we focused on ICU patients who had no prior history of AKI and were free of AKI-related comorbidities. By focusing on this specific group of patients, we were able to eliminate the confounding influences of these conditions and examine the associations between the levels of serum electrolytes and renal function in a more controlled manner. Furthermore, our results can provide valuable insights for developing early intervention and effective management strategies as well as for investigating the pathophysiology of AKI.
The performance data in Table 4 show that for those patients without a prior history of AKI or AKI-related comorbidities, the relative risks with these alternative prediction models were fairly high, ranging from 9.84 to 16.89. This implies that the group of patients predicted to be positive suffered significantly higher risk than the groups of patients predicted to be negative. However, the low PPVs suggest that there would be a large number of false positives if these prediction models were put into practical use. Nevertheless, according to the numbers shown in Figure 3a, this particular DT model, if put into practical use, should predict around 57% of the patients to be negative and deliver a sensitivity around 80%. Meanwhile, according to the numbers shown in Figure 3b, this particular DT model, if put into practical use, should predict around 51% of the patients to be negative and deliver a sensitivity around 95%. Therefore, a physician who employs the DT models developed in this study to assess the risks of AKI for his/her patients only needs to focus on about 50% of the patients, while the physician can expect this group of patients to suffer about 10 times the risk of the group of patients predicted to be at low risk.
Among the 10 variables listed in Table 2, only serum creatinine, chloride, magnesium, and phosphorus are present in the DT models shown in Figure 3a,b. It must be noted that this observation does not imply that serum potassium, sodium, and non-ionized calcium are not associated with the development of AKI. In fact, as mentioned earlier, previous studies have reported that serum potassium, sodium, and non-ionized calcium are all associated with the development of AKI. What happened must be that when building the prediction model, the DT algorithm figured out that the levels of serum chloride, magnesium, and phosphorus provided more information than the levels of serum potassium, sodium, and non-ionized calcium. The DT algorithm further figured out that the additional information provided by the levels of serum potassium, sodium, and non-ionized calcium after the levels of serum chloride, magnesium, and phosphorus had been incorporated was insignificant.
The DT models shown in Figure 3a,b identify the levels of serum creatinine, chloride, and magnesium as the three major factors associated with the development of AKI. Though the level of serum phosphorus is present in these two figures, all three nodes corresponding to the level of serum phosphorus are located in the lower levels of the structures. Furthermore, only a very limited number of positive cases in our study cohort met the criteria defined by these low-level structures. Therefore, in practice, we can ignore the role of serum phosphorus.
Since the level of serum creatinine is one of the major factors monitored in the current clinical practice, our study suggests that for those patients without a prior history of AKI or AKI-related comorbidities, the levels of serum chloride and magnesium should be taken into consideration in order to enhance the clinical guidelines. In this respect, the current clinical guideline, which monitors only the level of serum creatinine and the volume of urine output, may lead to misdiagnoses and/or delayed treatments in some cases because the level of serum creatinine generally reflects the degree of renal damage and should be considered as a delayed indicator of AKI. Furthermore, decreased urine output is a non-specific symptom and may only be evident once the AKI has progressed. Therefore, by incorporating the assessments of the serum chloride and magnesium levels into the enhanced clinical guideline, healthcare professionals can obtain a more comprehensive understanding of a patient’s renal function and the risk of AKI. Furthermore, the numbers shown in Table 2 reveal that the distributions of the levels of serum creatinine for patients with AKI and patients without AKI must overlap to a large degree because the standard deviation of the level of serum creatinine for patients with AKI, which is 0.64, is larger than the difference between the means of these two groups of patients, which is 0.5. This implies that additional assessments must be incorporated if we would like to evaluate the risk of AKI of a patient more accurately. Finally, with respect to the decrease in urine output among AKI patients, it is a non-specific symptom and may only be evident once the AKI has progressed. Together, these observations imply that for an ICU patient without a prior history of AKI or AKI-related comorbidities, healthcare professionals can obtain a more comprehensive understanding of the patient’s renal function and risk of AKI by incorporating assessments of serum chloride and magnesium levels into the enhanced clinical guideline. Accordingly, healthcare professionals will be able to evaluate and manage treatments more precisely and ultimately prevent disease progression and deterioration.
It must be noted that our results can only be immediately applied to ICU patients without a prior history of AKI or AKI-related comorbidities. For ICU patients with AKI-related comorbidities, further studies are needed. In this respect, we can partition the patients into several groups depending on the types of comorbidities that they suffer from so that patients in the same group have similar pathophysiological mechanisms. Then, we can apply the procedure presented in this article to each group of patients in order to develop a specific prediction model for each group and identify the critical factors accordingly.
One of the major limitations of our study is due to the different causes of AKI. As the causes of AKI are essential for physicians to develop effective treatment plans, in-depth subgroup analyses based on different categories of renal injury should be conducted to gain valuable insights into the different pathophysiological mechanisms involved and guide appropriate treatment strategies tailored to each subgroup. In this study, based on the information available in the MIMIC-IV dataset, we classified the AKI patients into two categories: post-renal and non-post-renal. The statistics in Table 3 reveal that there were no statistical differences between the levels of the eight serum ingredients for the post-renal and non-post-renal AKI patients. Therefore, our prediction models should be generally applicable to both post-renal and non-post-renal AKI patients. Nevertheless, in-depth subgroup analyses should be conducted in the future.
In addition to the limitation addressed above, this is a retrospective study based on data extracted from the MIMIC-IV database. Therefore, the results derived from this study should not be extensively applied in the decision process without taking into consideration the ethnic composition of the patients and the medical interventions that these patients may have received. Furthermore, our study was based on clinical records collected in ICUs. This implies that the patients involved had serious health conditions. The data in Table 2 also show that these patients were relatively old. Therefore, the results observed in our analyses should not be generalized to patients with different health conditions and in different age groups. Finally, our results only illustrate the associations between the investigated risk factors and the incidences of AKI. In other words, causal inferences have yet to be studied.

5. Conclusions

This study has led to an in-depth understanding of the compound effects of serum creatinine, chloride, and magnesium with respect to the development of AKI in ICUs. As we focused on patients who had no prior history of AKI and were free of AKI-related comorbidities, our study provides valuable insights for developing early intervention and effective management strategies. Furthermore, this understanding provides crucial clues not only for future enhancement of clinical practices but also for future investigation of the pathophysiological mechanisms that are involved.

Supplementary Materials

The following supporting information can be downloaded at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/diagnostics13152551/s1, Table S1: The ICD-9 and ICD-10 codes employed to identify post-renal AKI cases. Table S2: Summary of the software packages employed and parameter ranges. Table S3: The detailed performance data observed in the 5-fold cross-validation process.

Author Contributions

H.-H.L. initiated the study, and all authors contributed to the development of the methods. Software development and data analysis were performed by Y.-T.W., W.-S.K.L. and Y.-J.O. supervised the study and wrote the first version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by National Taiwan University grant #FD107016.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Boards of the Massachusetts Institute of Technology (Protocol No. 0403000206) and Beth Israel Deaconess Medical Center (Protocol No. 2001-P-001699/14).

Informed Consent Statement

Not applicable because the MIMIC-IV database has been carefully de-identified to protect patient privacy.

Data Availability Statement

MIMIC-IV ver. 1.0 was provided by physionet.

Acknowledgments

The authors would like to acknowledge Albert Li for his professional comments and advice.

Conflicts of Interest

The authors declare no conflict of interest in this work.

References

  1. Loscalzo, J.; Fauci, A.S.; Kasper, D.L.; Hauser, S.L.; Longo, D.L.; Jameson, J.L.; Harrison, T.R. Harrison’s Principles of Internal Medicine, 21st ed.; International edition; Loscalzo, J., Fauci, A.S., Kasper, D.L., Hauser, S.L., Longo, D.L., Jameson, J.L., Eds.; McGraw-Hill: New York, NY, USA, 2022. [Google Scholar]
  2. Hoste, E.A.J.; Kellum, J.A.; Selby, N.M.; Zarbock, A.; Palevsky, P.M.; Bagshaw, S.M.; Goldstein, S.L.; Cerda, J.; Chawla, L.S. Global epidemiology and outcomes of acute kidney injury. Nat. Rev. Nephrol. 2018, 14, 607–625. [Google Scholar] [CrossRef]
  3. Case, J.; Khan, S.; Khalid, R.; Khan, A. Epidemiology of acute kidney injury in the intensive care unit. Crit. Care Res. Pract. 2013, 2013, 479730. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Huang, C.T.; Liu, K.D. Exciting developments in the field of acute kidney injury. Nat. Rev. Nephrol. 2020, 16, 69–70. [Google Scholar] [CrossRef]
  5. Mercado, M.G.; Smith, D.K.; Guard, E.L. Acute Kidney Injury: Diagnosis and Management. Am. Fam. Physician 2019, 100, 687–694. [Google Scholar] [PubMed]
  6. Koyner, J.L. Subclinical Acute Kidney Injury Is Acute Kidney Injury and Should Not Be Ignored. Am. J. Respir. Crit. Care Med. 2020, 202, 786–787. [Google Scholar] [CrossRef]
  7. Leaf, D.E.; Christov, M. Dysregulated Mineral Metabolism in AKI. Semin. Nephrol. 2019, 39, 41–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Yokota, L.G.; Sampaio, B.M.; Rocha, E.P.; Balbi, A.L.; Sousa Prado, I.R.; Ponce, D. Acute kidney injury in elderly patients: Narrative review on incidence, risk factors, and mortality. Int. J. Nephrol. Renovasc Dis. 2018, 11, 217–224. [Google Scholar] [CrossRef] [Green Version]
  9. Blaine, J.; Chonchol, M.; Levi, M. Renal control of calcium, phosphate, and magnesium homeostasis. Clin. J. Am. Soc. Nephrol. 2015, 10, 1257–1272. [Google Scholar] [CrossRef] [Green Version]
  10. Bellomo, R.; Kellum, J.A.; Ronco, C. Acute kidney injury. Lancet 2012, 380, 756–766. [Google Scholar] [CrossRef]
  11. Kellum, J.A.; Romagnani, P.; Ashuntantang, G.; Ronco, C.; Zarbock, A.; Anders, H.J. Acute kidney injury. Nat. Rev. Dis. Primers 2021, 7, 52. [Google Scholar] [CrossRef]
  12. Lv, Q.; Li, D.; Wang, Y.; Yu, P.; Zhao, L.; Chen, S.; Wang, M.; Fu, G.; Zhang, W. Admission electrolyte and osmotic pressure levels are associated with the incidence of contrast-associated acute kidney injury. Sci. Rep. 2022, 12, 4714. [Google Scholar] [CrossRef]
  13. Suetrong, B.; Pisitsak, C.; Boyd, J.H.; Russell, J.A.; Walley, K.R. Hyperchloremia and moderate increase in serum chloride are associated with acute kidney injury in severe sepsis and septic shock patients. Crit. Care 2016, 20, 315. [Google Scholar] [CrossRef] [Green Version]
  14. Marttinen, M.; Wilkman, E.; Petaja, L.; Suojaranta-Ylinen, R.; Pettila, V.; Vaara, S.T. Association of plasma chloride values with acute kidney injury in the critically ill—A prospective observational study. Acta Anaesthesiol. Scand. 2016, 60, 790–799. [Google Scholar] [CrossRef]
  15. Moon, H.; Chin, H.J.; Na, K.Y.; Joo, K.W.; Kim, Y.S.; Kim, S.; Han, S.S. Hyperphosphatemia and risks of acute kidney injury, end-stage renal disease, and mortality in hospitalized patients. BMC Nephrol. 2019, 20, 362. [Google Scholar] [CrossRef] [Green Version]
  16. Cheungpasitporn, W.; Thongprayoon, C.; Erickson, S.B. Admission hypomagnesemia and hypermagnesemia increase the risk of acute kidney injury. Ren. Fail. 2015, 37, 1175–1179. [Google Scholar] [CrossRef] [Green Version]
  17. Thongprayoon, C.; Cheungpasitporn, W.; Mao, M.A.; Sakhuja, A.; Erickson, S.B. Admission calcium levels and risk of acute kidney injury in hospitalised patients. Int. J. Clin. Pract. 2018, 72, e13057. [Google Scholar] [CrossRef]
  18. Thongprayoon, C.; Cheungpasitporn, W.; Chewcharat, A.; Mao, M.A.; Bathini, T.; Vallabhajosyula, S.; Thirunavukkarasu, S.; Kashani, K.B. Impact of admission serum ionized calcium levels on risk of acute kidney injury in hospitalized patients. Sci. Rep. 2020, 10, 12316. [Google Scholar] [CrossRef] [PubMed]
  19. Chen, D.N.; Du, J.; Xie, Y.; Li, M.; Wang, R.L.; Tian, R. Relationship between early serum sodium and potassium levels and AKI severity and prognosis in oliguric AKI patients. Int. Urol. Nephrol. 2021, 53, 1171–1187. [Google Scholar] [CrossRef] [PubMed]
  20. Yessayan, L.; Neyra, J.A.; Canepa-Escaro, F.; Vasquez-Rios, G.; Heung, M.; Yee, J.; Acute Kidney Injury in Critical Illness Study, G. Effect of hyperchloremia on acute kidney injury in critically ill septic patients: A retrospective cohort study. BMC Nephrol. 2017, 18, 346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Morooka, H.; Tanaka, A.; Kasugai, D.; Ozaki, M.; Numaguchi, A.; Maruyama, S. Abnormal magnesium levels and their impact on death and acute kidney injury in critically ill children. Pediatr. Nephrol. 2022, 37, 1157–1165. [Google Scholar] [CrossRef] [PubMed]
  22. James, G.; Witten, D.; Hastie, T.; Tibshirani, R. An Introduction to Statistical Learning: With Applications in R; Springer: New York, NY, USA, 2013. [Google Scholar]
  23. Cortes, C.; Vapnik, V. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  24. Koller, D.; Friedman, N. Probabilistic Graphical Models: Principles and Techniques/Daphne Koller, Nir Friedman; MIT Press: Cambridge, MA, USA, 2009. [Google Scholar]
  25. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  26. Song, X.; Liu, X.; Liu, F.; Wang, C. Comparison of machine learning and logistic regression models in predicting acute kidney injury: A systematic review and meta-analysis. Int. J. Med. Inform. 2021, 151, 104484. [Google Scholar] [CrossRef]
  27. Tomasev, N.; Glorot, X.; Rae, J.W.; Zielinski, M.; Askham, H.; Saraiva, A.; Mottram, A.; Meyer, C.; Ravuri, S.; Protsyuk, I.; et al. A clinically applicable approach to continuous prediction of future acute kidney injury. Nature 2019, 572, 116–119. [Google Scholar] [CrossRef]
  28. Therneau, T.; Atkinson, B. Rpart: Recursive Partitioning and Regression Trees. Available online: https://CRAN.R-project.org/package=rpart (accessed on 21 March 2022).
  29. Rokach, L.; Maimon, O. Data Mining with Decision Trees Theory and Applications/Lior Rokach, Oded Maimon, 2nd ed.; World Scientific: Hackensack, NJ, USA, 2015. [Google Scholar]
  30. Cho, S.; Hong, H.; Ha, B.C. A hybrid approach based on the combination of variable selection using decision trees and case-based reasoning using the Mahalanobis distance: For bankruptcy prediction. Expert. Syst. Appl. 2010, 37, 3482–3488. [Google Scholar] [CrossRef]
  31. Gao, W.; Wang, J.; Zhou, L.; Luo, Q.; Lao, Y.; Lyu, H.; Guo, S. Prediction of acute kidney injury in ICU with gradient boosting decision tree algorithms. Comput. Biol. Med. 2021, 140, 105097. [Google Scholar] [CrossRef]
  32. Chen, T.; Guestrin, C. XGBoost. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13 August 2016; pp. 785–794. [Google Scholar]
  33. Johnson, A.; Bulgarelli, L.; Pollard, T.; Horng, S.; Celi, L.A.; Mark, R. MIMIC-IV, Version 1.0; PhysioNet: Cambridge, MA, USA, 2021. [CrossRef]
  34. Goldberger, A.L.; Amaral, L.A.; Glass, L.; Hausdorff, J.M.; Ivanov, P.C.; Mark, R.G.; Mietus, J.E.; Moody, G.B.; Peng, C.K.; Stanley, H.E. PhysioBank, PhysioToolkit, and PhysioNet: Components of a new research resource for complex physiologic signals. Circulation 2000, 101, E215–E220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Khwaja, A. KDIGO Clinical Practice Guidelines for Acute Kidney Injury. Nephron Clin. Pract. 2012, 120, c179–c184. [Google Scholar] [CrossRef]
  36. Acute Kidney Injury Work Group. Section 2: AKI Definition. Kidney Int. Suppl. 2012, 2, 19–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kidney International Work Group. Summary of Recommendation Statements. Kidney Int. Suppl. 2012, 2, 8–12. [Google Scholar] [CrossRef] [Green Version]
  38. Palevsky, P.M.; Liu, K.D.; Brophy, P.D.; Chawla, L.S.; Parikh, C.R.; Thakar, C.V.; Tolwani, A.J.; Waikar, S.S.; Weisbord, S.D. KDOQI US commentary on the 2012 KDIGO clinical practice guideline for acute kidney injury. Am. J. Kidney Dis. 2013, 61, 649–672. [Google Scholar] [CrossRef]
  39. Quan, H.; Sundararajan, V.; Halfon, P.; Fong, A.; Burnand, B.; Luthi, J.C.; Saunders, L.D.; Beck, C.A.; Feasby, T.E.; Ghali, W.A. Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data. Med. Care 2005, 43, 1130–1139. [Google Scholar] [CrossRef] [PubMed]
  40. Liano, F.; Pascual, J. Epidemiology of acute renal failure: A prospective, multicenter, community-based study. Madrid Acute Renal Failure Study Group. Kidney Int. 1996, 50, 811–818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Basile, D.P.; Anderson, M.D.; Sutton, T.A. Pathophysiology of acute kidney injury. Compr. Physiol. 2012, 2, 1303–1353. [Google Scholar] [CrossRef] [Green Version]
  42. Dowdy, S.; Wearden, S.; Chilko, D. Statistics for Research; John Wiley & Sons, Incorporated: Hoboken, NJ, USA, 2004. [Google Scholar]
  43. MacFarland, T.W.; Yates, J.M. Introduction to Nonparametric Statistics for the Biological Sciences Using R; Springer International Publishing AG: Cham, Switzerland, 2016. [Google Scholar]
Figure 1. A DT structure that summarizes the main results of this study. The root node is colored yellow.
Figure 1. A DT structure that summarizes the main results of this study. The root node is colored yellow.
Diagnostics 13 02551 g001
Figure 2. The flow for generating the study cohort. Table 1 lists the ICD-9 and ICD-10 codes employed to exclude the cases with AKI-related comorbidities/diseases. Criterion (1): (i) For a patient who had suffered from AKI, we included only the record corresponding to his/her stay in the ICU during which the patient suffered from AKI the first time. (ii) For a patient who had never suffered from AKI, we included only the record corresponding to his/her first stay in the ICU. Criterion (2): (i) the record of the case did not include all the readings listed in Table 2; (ii) one or more readings in the record were in the highest 0.1% or the lowest 0.1% of the distributions; or (iii) one or more readings in the record were not measured within 168 h of admission.
Figure 2. The flow for generating the study cohort. Table 1 lists the ICD-9 and ICD-10 codes employed to exclude the cases with AKI-related comorbidities/diseases. Criterion (1): (i) For a patient who had suffered from AKI, we included only the record corresponding to his/her stay in the ICU during which the patient suffered from AKI the first time. (ii) For a patient who had never suffered from AKI, we included only the record corresponding to his/her first stay in the ICU. Criterion (2): (i) the record of the case did not include all the readings listed in Table 2; (ii) one or more readings in the record were in the highest 0.1% or the lowest 0.1% of the distributions; or (iii) one or more readings in the record were not measured within 168 h of admission.
Diagnostics 13 02551 g002
Figure 3. DT models with two different levels of sensitivity. (a) The DT model with sensitivity at the level of 0.80. (b) The DT model with sensitivity at the level of 0.95. The blue polygons encircle the structure shared by these 2 DT models. The root node is colored yellow.
Figure 3. DT models with two different levels of sensitivity. (a) The DT model with sensitivity at the level of 0.80. (b) The DT model with sensitivity at the level of 0.95. The blue polygons encircle the structure shared by these 2 DT models. The root node is colored yellow.
Diagnostics 13 02551 g003
Table 1. Excluding criteria for the cases with AKI-related comorbidities/diseases.
Table 1. Excluding criteria for the cases with AKI-related comorbidities/diseases.
Comorbidities/DiseasesICD-9ICD-10
Renal failure 1403.11, 403.91, 404.12, 404.92, 584.5–584.9, 585.1–585.9, 586, V42.0, V45.1, V56.0, V56.8I12.0, I13.1, N17.0–N17.2, N17.8, N17.9, N18.1–N18.9, N19, N25.0, Z49.0–Z49.2, Z94.0, Z99.2
Congestive heart failure398.91, 402.11, 402.91, 404.11, 404.13, 404.91, 404.93, 428.0–428.9I09.9, I11.0, I13.0, I13.2, I25.5, I42.0, I42.5–I42.9, I50.0–I50.9, P29.0
Diabetes250.0–250.7, 250.9E10.0–E10.9, E11.0–E11.9, E12.0–E12.9, E13.0–E13.9, E14.0–E14.9
Fluid and electrolyte disorders276.0–276.9E22.2, E86.0, E86.1, E86.9, E87.0–E87.8
1 Including end-stage renal disease, AKI, and chronic kidney disease.
Table 3. Statistical analysis of the characteristics of the post-renal AKI patients and the non-post-renal AKI patients in our study cohort.
Table 3. Statistical analysis of the characteristics of the post-renal AKI patients and the non-post-renal AKI patients in our study cohort.
Feature24 Cases with
Post-Renal AKI
(Mean ± SD)
526 Cases with
Non-Post-Renal AKI
(Mean ± SD)
p-Value
Age (years)74.16 ± 12.5465.30 ± 14.660.0037 *
Gender 0.0007 *
Male (%)23 (95.66%)326 (61.98%)
Female (%)1 (4.34%)200 (38.02%)
Serum
BUN (mg/dL)27.54 ± 10.7226.70 ± 15.550.7944
Creatinine (mg/dL)1.40 ± 0.571.36 ± 0.640.7459
Chloride (mEq/L)110.08 ± 7.30110.38 ± 6.570.5228
Potassium (mEq/L)4.58 ± 0.564.80 ± 0.760.2258
Sodium (mEq/L)142.81 ± 5.28142.84 ± 5.790.9599
Magnesium (mg/dL)2.54 ± 0.382.53 ± 0.530.1666
Phosphorus (mg/dL)4.07 ± 0.934.41 ± 1.360.5254
Non-ionized calcium (mg/dL)8.66 ± 0.548.76 ± 0.740.8278
The symbol * denotes statistical significance. The p-values were calculated based on the two-sample t-test. SD stands for standard deviation.
Table 4. Summary of the performances observed during the 5-fold cross-validation process.
Table 4. Summary of the performances observed during the 5-fold cross-validation process.
Level of SensitivityModelSensitivitySpecificityPPVAUCRelative Risk
0.95DT0.9490.4790.0760.76716.893
LR0.9490.4140.0680.85513.872
RF0.9490.3820.0650.66613.012
0.80DT0.7980.7210.1160.8239.84
LR0.7990.7730.1370.85711.982
RF0.7990.7320.1190.76610.141
PPV stands for positive predictive value, also known as precision. AUC stands for the area under the receiver operating characteristic curve.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, H.-H.; Wang, Y.-T.; Yang, M.-H.; Lin, W.-S.K.; Oyang, Y.-J. Exploiting Machine Learning Technologies to Study the Compound Effects of Serum Creatinine and Electrolytes on the Risk of Acute Kidney Injury in Intensive Care Units. Diagnostics 2023, 13, 2551. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics13152551

AMA Style

Liu H-H, Wang Y-T, Yang M-H, Lin W-SK, Oyang Y-J. Exploiting Machine Learning Technologies to Study the Compound Effects of Serum Creatinine and Electrolytes on the Risk of Acute Kidney Injury in Intensive Care Units. Diagnostics. 2023; 13(15):2551. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics13152551

Chicago/Turabian Style

Liu, Hsin-Hung, Yu-Tseng Wang, Meng-Han Yang, Wei-Shu Kevin Lin, and Yen-Jen Oyang. 2023. "Exploiting Machine Learning Technologies to Study the Compound Effects of Serum Creatinine and Electrolytes on the Risk of Acute Kidney Injury in Intensive Care Units" Diagnostics 13, no. 15: 2551. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics13152551

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop