Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 18 August 2025

Sec. Clinical Infectious Diseases

Volume 15 - 2025 | https://doi.org/10.3389/fcimb.2025.1586087

This article is part of the Research TopicOutcome of Sepsis and Prediction of Mortality Risk - Volume IIView all 7 articles

Development and validation of web-based, interpretable predictive models for sepsis and mortality in extensive burns

Shi-Qi Wang,&#x;Shi-Qi Wang1,2†Kan Qiu&#x;Kan Qiu3†Qi-Rui ZhengQi-Rui Zheng4Bing-Jie ZhouBing-Jie Zhou5Ming-Yu LiMing-Yu Li2Hai-Yan ZhongHai-Yan Zhong1Yong ChenYong Chen1Si-Ming Yuan,*&#x;Si-Ming Yuan1,2*‡
  • 1Department of Burns and Plastic Surgery, Jinling Hospital, Jinling Clinical Medical College, Nanjing University of Chinese Medicine, Nanjing, Jiangsu, China
  • 2Department of Burns and Plastic Surgery, Jinling Hospital, Affiliated Hospital of Medical School, Nanjing University, Nanjing, Jiangsu, China
  • 3Department of Burns and Plastic Surgery, Anqing Petrochemical Hospital, Nanjing Drum Tower Hospital Group, Anqing, Anhui, China
  • 4School of Computer Science, Peking University, Beijing, China
  • 5Southeast University School of Medicine, Nanjing, Jiangsu, China

Background: Burn injuries are a common cause of trauma globally, with extensive burns (≥ 50% total body surface area burned) associated with high rates of sepsis and mortality. This study aims to identify risk factors associated with sepsis and mortality in extensively burned patients and to develop accurate, interpretable predictive models via machine learning algorithms.

Methods: A retrospective cohort study was conducted utilizing data from two Burn Critical Care Units in Eastern China from 2012-2023. A total of 237 patients with extensive burns were included. We applied ten machine learning algorithms, including random forest, gradient boosting tree (GBT), and logistic regression, to predict sepsis and mortality. The models were evaluated via AUC, precision, recall, accuracy, and F1 score, and were compared with the SOFA score performance. Model interpretability was enhanced via SHapley Additive exPlanations (SHAP).

Results: The key predictive factors for sepsis included the SOFA score, new onset shock, albumin, blood urea nitrogen (BUN), third-degree burned area, TBSA burned, white blood cell count, and inhalation injury. For mortality, the key predictive factors included alanine aminotransferase (ALT), the SOFA score, type of burn, new onset shock, third-degree burn area, TBSA burned, and sepsis. The RF model demonstrated superior performance in predicting sepsis (AUC = 0.977, accuracy = 0.945, recall = 0.964, precision = 0.930, and F1 score = 0.945). For mortality prediction, the GBT model yielded the highest AUC of 0.981 (accuracy = 0.952, recall = 0.965, precision = 0.942, and F1 score = 0.953). The sepsis prediction model outperformed the SOFA-based logistic regression model. Web-based calculators were developed to aid clinical decision-making.

Conclusion: Machine learning models, RF and GBT, demonstrate strong predictive ability for sepsis and mortality in extensive burn patients. The application of SHAP enhances model transparency, facilitating clinical interpretation and early intervention. Two web-based calculators can guide intensive care strategies and improve patient outcomes.

1 Introduction

Burns are one of the most common and severe forms of trauma worldwide, with approximately 330,000 deaths annually attributed to burns (Organization, W.H, 2002). In developing countries, burns predominantly affect young individuals and the working population, imposing a significant burden on healthcare systems, labor markets, and economic development. Extensive burns (with a burn area ≥ 50% of the total body surface area, TBSA) cause profound disruption of the skin barrier, leading to vascular leakage, hypovolemic shock, immune dysregulation, and metabolic disturbances. These systemic responses significantly increase the susceptibility to sepsis (Cato et al., 2018; Torres et al., 2021; Yang et al., 2023). Compared with patients sustaining burns of 20–50% TBSA, those with burns exceeding 50% are at particularly high risk due to the critical shortage of autologous donor skin, which delays wound closure and prolongs exposure to pathogens (Zhang et al., 2021). Reports indicate that the mortality rate associated with burns complicated by sepsis can be as high as 20.3% (Yoon et al., 2018). Early identification of sepsis risk and prompt administration of appropriate antimicrobial treatment can effectively improve prognosis and even reduce mortality (Li et al., 2022). In addition to sepsis, large burns frequently precipitate systemic inflammatory response syndrome (SIRS), further promoting organ dysfunction, including cardiovascular collapse, acute respiratory distress syndrome (ARDS), and neurogenic shock (Raith et al., 2017). Identifying risk factors associated with mortality allows burn centers to quickly integrate resources for targeted, prioritized care.

Severe burns result in profound systemic responses, including sustained inflammation and metabolic disturbances, which complicate the early recognition of sepsis. Unlike earlier criteria that emphasize systemic inflammatory markers, such as those defined by the American Burn Association (ABA) (Greenhalgh et al., 2007) and the Mann-Salinas burn-specific predictors (Mann-Salinas et al., 2013; Yan et al., 2018), the Sepsis-3 consensus definition established in 2016 prioritizes organ dysfunction over inflammatory signs (Singer et al., 2016), representing a more appropriate approach for this patient population. In a comparative study, Yan et al. (2018) evaluated the diagnostic sensitivity of these three criteria in burn patients and found that Sepsis-3 demonstrated a sensitivity of 85%, which was higher than that of the ABA criteria (60%) and the Mann-Salinas model (20%), supporting its broader clinical applicability in this context.

In addition to sepsis, accurate prediction of mortality in patients with extensive burns remains a major clinical challenge. Several classical prognostic models have been developed to estimate mortality in burn patients, including the Baux score (age + %TBSA) (Baux, 1961), the Abbreviated Burn Severity Index (ABSI) (age, gender, inhalation injury, TBSA) (Tobiasen et al., 1982), and the Flame score (age, APACHE II, TBSA, gender) (Gomez et al., 2008). However, with advances in critical care, the relevance of these models has been increasingly questioned. For example, the ABSI associates female gender with worse outcomes, though recent evidence disputes this link (Ederer et al., 2019). The Baux score, being overly simplistic, lacks predictive accuracy (Liu et al., 2015), while the Flame score is too complex for routine bedside use without computational tools.

Machine learning is a subdiscipline of artificial intelligence that uses various mathematical algorithms to identify certain data patterns, thus making predictive conclusions about selected endpoints. Unlike traditional statistical models, machine learning can independently identify complex data relationships, overcome linear limitations, and ensure stability in high-dimensional datasets (Deo, 2015). Recent studies have utilized machine learning to develop predictive models for burn patients; for instance, Tran et al. developed an automated machine-learning platform to predict sepsis in burn patients using clinical data (Tran et al., 2020), while Jeon et al. developed a model leveraging immature granulocytes as biomarkers for sepsis prediction (Jeon et al., 2021). However, these studies mainly focus on general burn patients with smaller TBSA and use single algorithms without incorporating interpretability. There is no interpretable machine learning model specifically developed for clinical application in patients with extensive burn injuries.

In this study, our objective was to develop predictive models for sepsis and mortality in patients with extensive burns, based on data obtained within the first 24 hours of admission. We conducted a retrospective study of patients with extensive burns admitted to two Burn Critical Care Units in Eastern China. Using ten commonly applied machine learning algorithms, we developed predictive models to identify the best-performing approaches. Additionally, we compared the sepsis prediction model we developed with the SOFA score prediction using receiver operating characteristic (ROC) curve analysis to further validate its discrimination. To interpret the results of the machine learning models, we combined advanced machine learning algorithms with a method based on SHapley Additive exPlanations (SHAPs). SHAP is a popular machine learning technique that can provide a deeper understanding of the complex relationships between features and predictions (Lundberg et al., 2018). Intuitive explanations help clinicians comprehensively understand how the developed models make specific predictions, thereby increasing opportunities for early intervention. Finally, we created web-based calculators for the optimal models to facilitate their use by clinical practitioners.

2 Methods

2.1 Study population

A 12-year retrospective, dual-center dataset was established, covering patients with extensive burns who were admitted to two Burn Critical Care Units, Jinling Hospital Burn Critical Care Unit and Anqing Petrochemical Hospital Burn Critical Care Unit, in Eastern China from January 2012 to December 2023. The inclusion and exclusion criteria were as follows:

www.frontiersin.org

.

This study received approval from the Institutional Review Board of Jinling Hospital (Approval No. 2023-LLSC-026). Given the retrospective study design and anonymized clinical data, informed consent from participants was waived. The flowchart of the study was presented in Figure 1. The process of data selection is illustrated in Supplementary Figure S1 (see the appendix in the electronic Supplementary Material).

Figure 1
Flowchart illustrating a data analysis pipeline from 2012 to 2023. It involves data collection from two burn critical care units in Eastern China, selecting 28 variables from 237 patients, and using SMOTE for variable selection. This leads to model training with a 70/30 train/test split. The process includes modeling, evaluation using metrics like AUC, and explaining results with SHAP. The final output is a web-based calculator for personalized predictions.

Figure 1. Flowchart of the study. SMOTE, synthetic minority over-sampling technique; SOFA, sequential organ failure assessment; AUC, area under the curve; SHAP, SHapley additive exPlanations.

2.2 Outcomes and predictors

2.2.1 Outcomes

Sepsis: Sepsis refers to an occurrence during hospitalization. The presence of sepsis and septic shock was determined according to the Sepsis 3.0 criteria, which require evidence of infection combined with an increase in the SOFA score of ≥ 2 points (Singer et al., 2016). For consistency across the study period (2013–2022), all sepsis cases were retrospectively identified using the Sepsis-3 definition by trained clinicians during data collection in 2022–2024.

Mortality: Mortality refers to in-hospital death occurring beyond 24 hours after admission. Patients were classified as either survivors or non-survivors accordingly.

2.2.2 Predictors

The following variables collected upon admission were considered as potential predictors, categorized and discussed as follows:

1. Demographic Information: age, sex.

2. Clinical Characteristics at Admission: injury to hospital time, type of burn, TBSA burned, third-degree burned area (%TBSA), open decompression, tracheotomy, presence of inhalation injury, shock status, history of hypertension or diabetes mellitus, Sequential Organ Failure Assessment (SOFA) score: emphasized as a key prognostic indicator, providing a quantitative measure of organ dysfunction severity.

3. Treatment Factors: fluid resuscitation, primary treatment.

4. Laboratory Data: White Blood Cell (WBC) count, Red Blood Cell (RBC) count, hemoglobin; platelet count, Aspartate Aminotransferase (AST), Alanine Aminotransferase (AST), Indirect Bilirubin (IBIL), Direct Bilirubin (DBIL), Total Bilirubin (TBIL), Albumin (ALB), Globulin (GLB), Blood Urea Nitrogen (BUN), Serum Creatinine (SCR), Potassium (K+), Sodium (Na+), Chloride (Cl-), Calcium (Ca2+), Procalcitonin (PCT).

Both categorical and continuous variables were included to maximize the available information and capture diverse aspects of the patient’s condition.

2.3 Data processing and handling

The data were collected using Microsoft Excel (Redmond, Washington, USA) and processed with Python 3.9. Categorical variables were presented as counts (percentages) and analyzed using the chi-square test or Fisher’s exact test when appropriate. Otherwise, variables were expressed as means ± standard deviations and analyzed using the independent-samples t-test for normally distributed data; otherwise, the Mann-Whitney U test was applied. The overall rate of missing data was less than 5%. For variables with missing values, multiple imputations were used for numerical variables and mode imputation for categorical variables, assuming data were missing at random.

The dataset included 162 patients (68.4%) who developed sepsis and 75 patients (31.6%) who did not. Additionally, 199 patients survived (83.9%), whereas 38 (16.1%) died. The significant imbalance in the outcome categories posed a challenge for model prediction accuracy. To address this, we used the synthetic minority oversampling technique (SMOTE), which generates synthetic samples by interpolating between nearest neighbors within the minority class. This approach helps balance the data distribution and improves the model’s prediction capabilities for the minority class (Chawla et al., 2002; Haibo and Garcia, 2009). Variables with a P value < 0.001 in univariate regression analysis were selected as candidate factors for constructing predictive models.

2.4 Development of ten machine learning models

The entire dataset was randomly divided into a training set and a validation set at a 7:3 ratio, where the training set was used to construct the predictive models, and the validation set was used to evaluate model performance.

Ten commonly used machine learning algorithms were employed to establish prediction models, including Random Forest (RF), Multilayer Perceptron (MLP), Adaptive Boosting (AdaBoost), Logistic Regression (LR), Extreme Gradient Boosting (XGBoost), Naive Bayes (NB), Gradient Boosting Tree (GBT), Support Vector Machine (SVM), Decision Tree (DT), and K-Nearest Neighbor (KNN) methods. Among these, RF, AdaBoost, GBT, and XGBoost are ensemble classifiers, while the rest are single classifiers. During the training process, 5-fold cross-validation combined with grid search was applied to optimize the hyperparameters for each model, ensuring robustness and generalizability.

2.5 Model evaluation

The models’ performance was evaluated using several metrics, including accuracy, precision, recall, F1 score, and AUC (area under the curve). Higher scores closer to 1 indicate better predictive performance. ROC curves were used to graphically illustrate the discriminative power of the model. Furthermore, the ROC curve of the best-performing machine learning model was compared with that of a logistic regression model based solely on the SOFA score.

2.6 Model interpretability

To improve the interpretability of the optimal model, SHapley Additive exPlanations (SHAPs) were employed. SHAP values provide localized explanations, quantifying each variable’s impact on individual predictions (Lundberg et al., 2018). Originating from cooperative game theory, SHAP values represent the average marginal contribution of a variable across all possible combinations. For each individual prediction, the SHAP value indicates whether a variable contributes positively or negatively to the predicted outcome.

2.7 Web-based calculators

Two user-friendly web-based calculators were developed using Streamlit (https://github.com/streamlit/streamlit). Streamlit is an open-source library that enables researchers to create web applications with only a few lines of Python code, without requiring expertise in front-end development.

3 Results

3.1 Data characteristics

After applying the inclusion and exclusion criteria, data from 237 patients with extensive burns were included in the final model development (Figure 1). The dataset comprised 147 male and 90 female patients. The baseline characteristics, clinical features, treatment interventions, and initial laboratory results upon admission were summarized in Table 1. Additional in-hospital information, including wound culture, blood culture, catheter culture, length of stay, and multidrug-resistant (MDR), is provided in Supplementary Table S1.

Table 1
www.frontiersin.org

Table 1. Differences in demographic and clinical characteristics between sepsis and non-sepsis groups, as well as between the survived and deceased groups.

In comparison to the non-sepsis cohort, the sepsis group exhibited a significantly higher rate of inhalation injuries (48.00% vs. 29.01%, P = 0.007). The TBSA burned was greater in the sepsis group (61.36% vs. 56.07%, P < 0.001), with a notably greater proportion of III-degree burns (18.76% vs. 9.81%, P < 0.001). In terms of laboratory findings, the sepsis group presented higher elevated white blood cell counts (25.18 ± 11.53 vs. 20.09 ± 8.55, P < 0.001). Additionally, new onset shock was observed in 49.33% of patients with sepsis (n = 37), compared to only one case (0.62%) in the non-sepsis group (P < 0.001). Abnormal liver and kidney function, reflected in elevated TBIL, BUN, and SCR, was also more prevalent in the sepsis group. Lastly, the median SOFA score was significantly higher among sepsis patients [2.00 (IQR: 2.00–3.00) vs. 0 (IQR: 0.00–1.00)], suggesting its value in early prognostic assessment.

Comparisons between the mortality and survival groups revealed that extensive burns resulting from scalds (23.68% vs. 11.56%), explosions (7.89% vs. 3.02%) and electrical injuries (5.36% vs. 3.02%) were more prevalent in the mortality group. Additionally, the TBSA burned (60.45 ± 8.79% vs. 57.23 ± 7.37%, P = 0.04) and third-degree burned area (23.63 ± 28.47% vs. 10.54 ± 11.11%, P = 0.008) were notably greater. The mortality group presented a significantly greater incidence of new onset shock (47.37% vs. 10.05%, P < 0.001) and sepsis (63.16% vs. 25.63%, P < 0.001). Among the laboratory indicators, abnormal elevations in the liver enzyme ALT (107.47 ± 149.93 vs. 28.58 ± 54.33, P = 0.003) were more pronounced in the mortality group. The SOFA score also demonstrated a marked difference, with a median score of 2.00 (IQR: 1.00, 2.00) in the mortality group versus 1.00 (IQR: 0.00, 2.00) in the survival group, emphasizing its value as a predictive tool.

3.2 Features selected in the models

The data balanced via the SMOTE algorithm is shown in Supplementary Table S2. Using sepsis occurrence as the outcome variable, univariate regression analysis was conducted to identify significant features. The regression coefficients and odds ratios (ORs) with 95% confidence intervals for these predictors are summarized in Table 2. After highly correlated features were removed, 8 key variables (P < 0.005) were selected: TBSA burned, third-degree burned area, new-onset shock, WBC, ALB, the SOFA score, and inhalation injury.

Table 2
www.frontiersin.org

Table 2. Features selected in models.

For mortality as the outcome variable, univariate regression analysis was similarly conducted. After highly correlated variables were excluded, 7 significant predictors (P < 0.001) were identified: sepsis, TBSA burned, third-degree burned area, type of burn, new onset shock, ALT, and the SOFA score.

3.3 Model evaluation

For the 10 machine learning models developed, five key metrics were used to assess their generalization performance: accuracy, precision, F1 score, recall, and AUC (area under the curve). These results are summarized in Table 3. To facilitate comparison, the five evaluation metrics are visualized as bar charts (Figure 2). Additionally, ROC curves for each model are displayed in Figures 3A, C.

Table 3
www.frontiersin.org

Table 3. Evaluation metrics of all machine learning algorithms.

Figure 2
Bar charts comparing model scores for different metrics. Chart A shows scoring of accuracy, precision, recall, F1 score, and AUC across various models like Random Forest and XGBoost. Chart B, titled “Models for Predicting Mortality,” presents similar metrics for models like GBT and Logistic Regression. Each model is color-coded, with a legend included.

Figure 2. (A) Evaluation metrics of models for predicting sepsis; (B) evaluation metrics of models for predicting mortality. XGBoost, extreme gradient boosting; KNN, k-nearest neighbors; SVM, support vector machine; GBT, gradient boosting tree; Adaboost, adaptive boosting; MLP, multi-layer perceptron; AUC, area under the curve.

Figure 3
Three Receiver Operating Characteristic (ROC) curves depict the performance of various models. Panel A shows models for predicting sepsis with Random Forest achieving the highest area under the curve (AUC) at 0.966. Panel B compares Random Forest (AUC 0.983) and SOFA Logistic Regression (AUC 0.918). Panel C illustrates models for predicting mortality, with Gradient Boosting Trees obtaining the highest AUC at 0.984. Each ROC curve displays true positive versus false positive rates, with a diagonal line representing chance.

Figure 3. ROC curves. (A) ROC curves for predicting sepsis in extensive burn patients with machine learning algorithms; (B) The ROC curve comparison between the RF model and the SOFA logistic regression model for predicting sepsis in extensive burn patients; (C) ROC curves for predicting mortality in extensive burn patients with machine learning algorithms; (B) The ROC curve comparison between the XGBoost model and the SOFA logistic regression model for predicting mortality in extensive burn patients. ROC, receiver operating characteristic; XGBoost, extreme gradient boosting; KNN, k-nearest neighbors; SVM, support vector machine; GBT, gradient boosting tree; Adaboost, adaptive boosting; MLP, multi-layer perceptron; AUC, area under the curve.

When sepsis occurrence was used as the outcome variable, the RF model emerged as the best-performing model for predicting sepsis in patients with extensive burns, as it achieved the highest AUC of 0.977. Other performance metrics for this model were similarly strong, with an accuracy of 0.945, a recall of 0.964, a precision of 0.930, and an F1 score of 0.945. As shown in Figure 3B, the ROC curve of the RF model significantly outperformed that of a logistic regression model based solely on the SOFA score (z score = 2.689, P = 0.007, DeLong’s test (DeLong et al., 1988).

When mortality was used as the outcome variable, the GBT model yielded the highest AUC (0.981), with accuracy (0.952), recall (0.965), precision (0.942), and F1 score (0.953), making it the choice for predicting mortality in patients with extensive burns.

3.4 SHAP for the established models

To understand the role of key biomarkers and other variables within the machine learning prediction models, the SHAP tree framework was applied to the best-performing models. SHAP values not only quantify the direct influence of individual biomarkers but also account for their interactions with other variables, providing an interpretable and comprehensive view of model predictions.

In the model for predicting sepsis in burn patients, positive SHAP values indicated an increased risk, whereas negative values suggested a protective effect. The greater the SHAP value deviates from 0, the stronger its impact on the prediction. As shown in Figure 4A, high albumin levels (indicated by the blue dots) were associated with a lower sepsis risk, whereas high values of features such as SOFA score, new onset shock, ALB, third-degree burned area, WBC, and inhalation injury indicated a higher risk. Figures 4B, C present individual SHAP plots for patients with (A) sepsis and (B) non-sepsis. SHAP values highlight the specific predictive features for each individual patient and illustrate the contribution of each feature to sepsis prediction. In Figure 4B, red bars indicate features that increase the risk of sepsis, driving the final prediction value f(x) = 0.98, indicating a high risk of sepsis. Notably, the SOFA score, new onset shock and TBSA burned were the most significant factors contributing to increased sepsis risk. In contrast, blue bars represent third-degree burned area, and ALB showed a negative contribution to the risk prediction; however, their effects were relatively small in magnitude compared to the dominant high-risk factors in this specific case. Figure 4C shows that while third-degree burned area and ALB slightly increased sepsis risk, the overall effect of most features, particularly those in the blue region, decreased the predicted sepsis risk. Thus, the final prediction value of 0.02 indicated a low risk of sepsis, suggesting that the patient’s overall condition was relatively stable.

Figure 4
Three visualizations illustrate a random forest model predicting sepsis. Panel A shows a SHAP summary plot with features like SOFA, new onset shock, and ALB affecting the model output, indicated by SHAP values. Panel B displays a SHAP force plot for sepsis with contributions of various features such as Inhalation Damage and WBC to the model output. Panel C shows a SHAP force plot for no sepsis, highlighting features like ALB and SOFA. The plots use color gradients to denote feature values, from low (blue) to high (red).

Figure 4. (A) SHAP beeswarm summary plot on the impact of input variables on the RF model’s predictive ability. SHAP force plot showing (B) high chance of sepsis, (C) low risk of sepsis. RF, random forest; SOFA, sequential organ failure assessment; ALB, albumin; BUN, blood urea nitrogen; III, third-degree burned area; WBC, white blood cell count; SHAP, SHapley additive exPlanations.

Figure 5A shows the contributions of biomarkers and clinical features to mortality risk in the GBT model. Elevated levels of ALT, alongside SOFA score and type of burn were the most significant factors driving mortality risk. For surviving patients (Figure 5B), certain features, including SOFA score, type of burn, the absence of new onset shock. In deceased patients (Figure 5C), a combination of high-risk biomarkers such as third-degree burned area, SOFA score, and ALT drove a high mortality prediction.

Figure 5
Three-part image showing SHAP values for a GBT model predicting mortality. Part A is a dot plot with features: ALT, SOFA, Type of burn, New onset shock, III, TBSA burned, and Sepsis, indicating SHAP value impacts. Red represents high feature value, blue low. Part B explains a survivor's lower SHAP values in blue text. Part C shows a non-survivor’s higher SHAP values in red text.

Figure 5. (A) SHAP beeswarm summary plot on the impact of input variables on the GBT model’s predictive ability. SHAP force plot showing (B) high chance of survival (C) low risk of survival. GBT, gradient boosting tree; ALT, alanine aminotransferase; SOFA, sequential organ failure assessment;III, third-degree burned area; SHAP, SHapley additive exPlanations.

3.5 Web-based calculators

To facilitate use by healthcare professionals, we developed two web-based calculators via the Streamlit platform. Users can input relevant variables to obtain the predicted probabilities of sepsis and mortality. The calculators are accessible at the two URLs: https://sepsis-prediction-model-for-severe-burn-patients.streamlit.app. And https://mortality-prediction-model-for-severe-burn-patients.streamlit.app.

Examples of how to use the calculators are illustrated in Supplementary Figure S2 (see the appendix in the electronic Supplementary Material).

4 Discussion

Managing patients with extensive burns remains clinically challenging due to their high mortality risk and the scarcity of robust clinical data and standardized treatment protocols. In patients with extensive burns, the loss of skin barrier protection, hypovolemia, dysregulated inflammatory responses, and hypermetabolism contribute to a high risk of sepsis, which remains a major cause of mortality in this population. Early identification of sepsis risk factors and the timely application of comprehensive treatments, including antimicrobial therapy, are crucial to improving patient survival. In this study, we collected data from over 12 years of extensive burn patients at two burn centers in Eastern China and used ten machine learning algorithms to develop predictive models for sepsis and mortality based on admission data. Our results indicate that random forest and GBT models performed best for predicting sepsis and mortality, respectively. Through feature importance analysis and SHAP interpretability methods, we identified key risk factors for sepsis and mortality in patients with extensive burns. Finally, the developed web-based tool will assist healthcare workers in applying these findings in clinical practice.

The SOFA score evaluates six critical organ systems, circulatory, respiratory, renal, hepatic, coagulation, and neurological, by incorporating a wide range of parameters, including vasopressor dosage, mean arterial pressure (MAP), respiratory support, serum creatinine, urine output, bilirubin levels, platelet count, and the Glasgow Coma Scale (GCS) score (Singer et al., 2016). As such, it serves as a comprehensive and robust tool for detecting organ dysfunction. In patients with extensive burns, additional factors such as TBSA burned, especially the third-degree burned area, new onset shock, WBC, ALB, BUN, and inhalation injury also play crucial roles. Our comparative analysis indicates that integrating these variables into the model significantly improves its predictive performance for sepsis (Figure 3B).

In the models we developed, the TBSA burned emerged as a highly significant biomarker in predicting both sepsis and mortality risk. Studies have shown that the mortality rate for burn patients with a TBSA burned ranging from 21% to 30% is two-thirds higher than that for those with a TBSA burned ranging from 11% to 20% (Lionelli et al., 2005). The survival likelihood decreases markedly with increasing TBSA burned: reported mortality rates reach 51.1%, 70.6%, and 82.6% for patients with TBSA of 70–79%, 80–89%, and ≥ 90%, respectively (Iqbal et al., 2013). In our cohort, where the TBSA burned was ≥ 50%, the mortality rate was 16.03%. Although lower than figures reported in some international studies, the mortality rate falls within the nationally reported range for patients with severe burns (9.79%–23.2%) (Cheng, 2017; He et al., 2022). Several factors may account for this comparatively favorable outcome. First, patients discharged within the first 24 hours were excluded from the analysis; some of these individuals may have been transferred to tertiary burn centers for advanced care or may have discontinued treatment (Cheng, 2017), potentially affecting mortality estimates. Second, a subset of patients sustained extensive burns primarily composed of superficial injuries (e.g., grade I and superficial partial-thickness burns), which inflate TBSA calculations but are associated with better prognoses Finally, data were derived from two specialized burn centers equipped with comprehensive multidisciplinary teams, including intensive care, burn surgery, nephrology, and nutritional support, which likely contributed to improved survival outcomes.

Third-degree burns, also known as full-thickness burns, involve all layers of the skin, leading to severe immune dysregulation (Greenhalgh, 2017). The complete loss of protective barriers significantly increases the risk of both local and systemic infections, thereby increasing the probability of sepsis. Severe fluid loss, electrolyte imbalances, multi-organ failure, and uncontrolled sepsis greatly increase mortality risk. These findings are consistent with previous research reports (Tan Chor Lip et al., 2019; Tsolakidis et al., 2022; Schmidt et al., 2023).

In our study cohort, scald, blast, and electrical burns were associated with higher mortality rates among patients with extensive burns. Although scald injuries typically have lower mortality than flame burns in small-area cases, the risk increases significantly as TBSA increases. For example, among patients in our cohort with scald burns and TBSA ≥ 50%, the mortality rate reached 31.03% (9 deaths out of 29 cases), which is comparable to the 34% reported by Chong et al. in patients with scald injuries over 40% TBSA (Chong et al., 2020). Similarly, Saraçoğlu et al. (Saracoglu et al., 2014) observed a 26% mortality rate in electrical burns with TBSA > 25%, which may be attributed to the high contact temperatures (2000–4000°C) and deep tissue destruction characteristic of electrical injuries.

During the course of illness in burn patients, the occurrence of new onset shock indicates an acute and life-threatening condition. Severe capillary leakage and cardiac suppression result in poor organ perfusion, exacerbating the risk of multiple organ dysfunction. Moreover, the hypermetabolic state induced by burns further depletes the patient’s energy reserves and impairs the immune system, making organ function recovery more challenging. These factors, in combination, significantly impact the prognosis of burn patients (Jeschke et al., 2020).

Among laboratory biomarkers, serum albumin (Alb) has emerged as a strong predictor of sepsis. Lower albumin levels are associated with increased risk (Aguayo-Becerra et al., 2013). Alb exerts a key regulatory role in the inflammatory response by binding to pathogen-associated molecular patterns (PAMPs), such as peptidoglycan, lipoteichoic acid, and lipopolysaccharide, forming Alb–PAMP complexes that modulate immune activation (Gioannini et al., 2002; Esparza et al., 2012). In the pathophysiological context of severe burns, increased capillary permeability and an intensified inflammatory state led to a significant decline in serum albumin levels (Fleck et al., 1985). Hypoalbuminemia, in turn, contributes to extravascular fluid accumulation, resulting in complications such as impaired wound healing and edema (Lehnhardt et al., 2005). Furthermore, decreased albumin levels reduce the body’s antioxidant defenses, rendering patients more susceptible to oxidative stress–related damage associated with sepsis (Taverna et al., 2013). Therefore, in the management of patients with extensive burns, timely supplementation with albumin and blood products may serve as an important strategy for preventing or mitigating sepsis.

WBC count, a classic and sensitive hematological marker for infection, emerged in our model as a stronger predictor of sepsis than PCT. Unfortunately, due to equipment limitations at the participating hospitals, we were only able to collect high-sensitivity CRP (hs-CRP) data for approximately half of the patients. To minimize potential bias, we therefore excluded CRP from the set of predictive variables. Nonetheless, consistency analysis of the available data indicated a high level of agreement between CRP and PCT levels. Both markers are known to be elevated in burn patients (Assicot et al., 1993; Thomas et al., 2004); however, they lack specificity for early sepsis diagnosis in this population (Lavrentieva et al., 2007; Jeschke et al., 2013). By contrast, WBC demonstrated superior predictive ability in our feature selection process. This may be attributed to its greater sensitivity to early systemic inflammatory responses, as well as its broader availability, more frequent testing, and higher temporal stability in clinical settings (Knaus et al., 1985; Rimmer et al., 2022). These findings highlight the potential of WBC as a practical and robust early indicator of sepsis in severely burned patients, especially in resource-limited settings where advanced biomarker testing may not be consistently available.

Inhalation injury has been widely recognized as a significant risk factor for both sepsis and mortality in burn patients (Temiz et al., 2020; Yazıcı et al., 2022). It may result from direct thermal damage to the airway, chemical irritation from toxic gases, and systemic toxicity due to carbon monoxide exposure. These insults can initiate a strong inflammatory response, potentially leading to complications such as pneumonia, acute respiratory distress syndrome (ARDS), and even death. In our cohort of patients with extensive burns, the incidence of inhalation injury was 35.02%, which is slightly higher than the rates reported by most hospitals in China, ranging from 19.09-32.38% (Dou and Zhang, 2021). This relatively high rate may reflect the greater injury severity in our study population. Importantly, 54 out of 83 patients with inhalation injury (65.06%) received timely endotracheal intubation. Although inhalation injury remained a strong predictive factor for sepsis, there was no statistically significant difference in its incidence between survivors and non-survivors. This finding may be partly attributed to the high rate of early airway intervention, which likely mitigated progression to life-threatening respiratory complications. According to expert consensus in China, early prophylactic tracheotomy is recommended for patients at high risk of airway obstruction. Such procedures should ideally be performed before the peak of tissue edema, in order to avoid severe neck swelling that may obscure anatomical structures and increase procedural difficulty (Ming et al., 2018).

BUN, as a marker of renal function (Jeschke et al., 2020), and ALT, as a marker of hepatic injury (Nielson et al., 2017), were both incorporated into the sepsis and mortality prediction models, respectively. Although the SOFA score incorporates key organ function parameters such as creatinine and bilirubin, it remains a categorical scoring system that may overlook subtle variations in laboratory values. The inclusion of continuous variables such as BUN and ALT in the predictive models may therefore capture early or subclinical organ dysfunction that is not adequately represented by SOFA scores alone. This highlights the limitations of SOFA when used as a standalone measure and supports the utility of integrating granular laboratory data for improved risk stratification.

Our study offers several distinctive features that contribute to advancements in predicting sepsis and mortality in critically ill burn patients:

1. Targeted population: We focused on patients with exceptionally extensive burns (TBSA burned ≥ 50%), addressing a critical research gap concerning outcomes in this high-risk population.

2. Comparative machine learning models: We employed and evaluated ten widely used traditional machine learning models to identify the most effective predictive model.

3. SOFA model comparison: We compared models with a single-variable SOFA prediction model, confirming the utility of the SOFA score in predicting sepsis in patients with extensive burns.

4. Model interpretability via SHAP: We utilized SHAP values to interpret the model both globally and on an individual level, enhancing the model’s transparency and clinical applicability.

5. User-Friendly web tools: We developed two accessible web calculators to facilitate usage by healthcare professionals.

However, there are limitations to our study. Although our dataset, which was derived from two burn centers in East China, is representative, its relatively small size may restrict the performance of the machine learning models. Moreover, the limited geographic and demographic scope could affect the generalizability of our findings across diverse regions and ethnic groups. Additionally, we acknowledge that treatment protocols may have evolved over the 12-year period covered in this study (Chai, 2013; Guo et al., 2018; Hunag, 2020). Nevertheless, by applying consistent diagnostic criteria and focusing exclusively on admission-time variables, we aimed to minimize potential confounding effects caused by temporal changes in clinical practice. Finally, the rapid advances in artificial intelligence offer opportunities for further enhancements. Multisensory data, as well as imaging data from CT and MRI, are increasingly being integrated into clinical decision-making (Burns et al., 2020; Gichoya et al., 2022). In future work, we aim to incorporate such data for more precise patient monitoring and timely interventions.

5 Conclusion

In patients with extensive burns, the random forest and GBT methods effectively predict sepsis and mortality. By integrating key biomarkers such as TBSA, albumin, total bilirubin, and serum creatinine, these models provide valuable insights into the underlying pathophysiology. The SHAP tool aids in both global and personalized model interpretation, enhancing clinical understanding. Web calculators further simplify usage for medical staff, but their clinical utility requires further validation over time and across independent studies.

Data availability statement

The data analyzed in this study is subject to the following licenses/restrictions: The datasets used and analyzed during the current study are available from the corresponding author on reasonable request. All the code used in this study has been uploaded to https://github.com/ZzZqr/sepsis_and_mortality. Requests to access these datasets should be directed to S-MY, eXVhbnNtQDE2My5jb20=.

Ethics statement

The studies involving humans were approved by the Institutional Review Board of Jinling Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

S-QW: Conceptualization, Writing – original draft, Writing – review & editing, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration. KQ: Conceptualization, Data curation, Writing – original draft. Q-RZ: Conceptualization, Data curation, Formal Analysis, Methodology, Software, Writing – original draft. B-JZ: Data curation, Methodology, Writing – original draft. M-YL: Writing – original draft. H-YZ: Data curation, Writing – original draft. YC: Data curation, Formal Analysis, Writing – original draft. S-MY: Conceptualization, Data curation, Writing – original draft.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by Natural Science Foundation of Jiangsu Province, China (No. BK20231504, awarded to S-MY), Medical Project of Jinling Hospital (No. 22LCZLXJS4, awarded to S-MY), the Foundation for Key Medical Specialty of Nanjing (SZDZK202001, awarded to S-MY) and Jiangsu Funding Program for Excellent Postdoctoral Talent (No. 2023ZB048, awarded to S-QW).

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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/fcimb.2025.1586087/full#supplementary-material

Supplementary Figure 1 | Flowchart illustrating the data selection process.

Supplementary Figure 2 | Examples of web-based calculators for prediction models in patients with extensive burns. (A) The random forest model predicts sepsis.(B) The GBT model predicts mortality.

Abbreviations

SIRS, systemic inflammatory response syndrome; SOFA, sequential organ failure assessment; TBSA, total body surface area; ABSI, abbreviated burn severity index; APACHE II, chronic health evaluation II; ROC, receiver operating characteristic; SHAPs, SHapley Additive exPlanations; WBC, white blood cells; ALB, albumin; SMOTE, synthetic minority over-sampling technique; RF, random forest; MLP, multilayer perceptron; AdaBoost, adaptive boosting; LR, logistic regression; XGBoost, extreme gradient boosting; NB, naive Bayes; GBT, gradient boosted tree; SVM, support vector machine; DT, decision tree; KNN, k-nearest neighbor; AUC, area under the curve; TBIL, bilirubin; BUN, blood urea nitrogen; SCR, serum creatinine; ALT, alanine aminotransferase; MDR, multidrug-resistant bacteria.

References

Aguayo-Becerra, O. A., Torres-Garibay, C., Macías-Amezcua, M. D., Fuentes-Orozco, C., Chávez-Tostado Mde, G., Andalón-Dueñas, E., et al. (2013). Serum albumin level as a risk factor for mortality in burn patients. Clinics (Sao Paulo) 68, 940–945. doi: 10.6061/clinics/2013(07)09

PubMed Abstract | Crossref Full Text | Google Scholar

Assicot, M., Gendrel, D., Carsin, H., Raymond, J., Guilbaud, J., and Bohuon, C. (1993). High serum procalcitonin concentrations in patients with sepsis and infection. Lancet 341, 515–518. doi: 10.1016/0140-6736(93)90277-n

PubMed Abstract | Crossref Full Text | Google Scholar

Baux, S. (1961). Contribution à l’étude du traitement local des brûlures thermiques étendues. Paris: These.

Google Scholar

Burns, J. E., Yao, J., Chalhoub, D., Chen, J. J., and Summers, R. M. (2020). A machine learning algorithm to estimate sarcopenia on abdominal CT. Acad. Radiol. 27, 311–320. doi: 10.1016/j.acra.2019.03.011

PubMed Abstract | Crossref Full Text | Google Scholar

Cato, L. D., Wearn, C. M., Bishop, J. R. B., Stone, M. J., Harrison, P., and Moiemen, N. (2018). Platelet count: A predictor of sepsis and mortality in severe burns. Burns 44, 288–297. doi: 10.1016/j.burns.2017.08.015

PubMed Abstract | Crossref Full Text | Google Scholar

Chai, J. K. (2013). Diagnosis and comprehensive management of sepsis after burn. Zhonghua Shao Shang Za Zhi 29, 105–108. doi: 10.3760/cma.j.issn.1009-2587.2013.02.003

Crossref Full Text | Google Scholar

Chawla, N. V., Bowyer, K. W., Hall, L. O., and Kegelmeyer, W. P. (2002). SMOTE: synthetic minority over-sampling technique. J. Artif. Intell. Res. 16, 321–357. doi: 10.1613/jair.953

Crossref Full Text | Google Scholar

Cheng, W. (2017). The current status of epidemiological study of burns and a multicenter study on the epidemiology and prognosis of massively burned patients in China (Master's thesis). Chinese People's Liberation Army (PLA) Medical School, Beijing, China.

Google Scholar

Chong, H. P., Quinn, L., Cooksey, R., Molony, D., Jeeves, A., Lodge, M., et al. (2020). Mortality in paediatric burns at the Women’s and Children’s Hospital (WCH), Adelaide, South Australia: 1960-2017. Burns 46, 207–212. doi: 10.1016/j.burns.2019.06.004

PubMed Abstract | Crossref Full Text | Google Scholar

DeLong, E. R., DeLong, D. M., and Clarke-Pearson, D. L. (1988). Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 44, 837–845. doi: 10.2307/2531595

Crossref Full Text | Google Scholar

Deo, R. C. (2015). Machine learning in medicine. Circulation 132, 1920–1930. doi: 10.1161/circulationaha.115.001593

PubMed Abstract | Crossref Full Text | Google Scholar

Dou, Z. and Zhang, G. A. (2021). Systematic review of the epidemiological characteristics of inhalation injury in burn patients in China. Zhonghua Shao Shang Za Zhi 37, 654–660. doi: 10.3760/cma.j.cn501120-20200306-00129

PubMed Abstract | Crossref Full Text | Google Scholar

Ederer, I. A., Hacker, S., Sternat, N., Waldmann, A., Salameh, O., Radtke, C., et al. (2019). Gender has no influence on mortality after burn injuries: A 20-year single center study with 839 patients. Burns 45, 205–212. doi: 10.1016/j.burns.2018.08.012

PubMed Abstract | Crossref Full Text | Google Scholar

Esparza, G. A., Teghanemt, A., Zhang, D., Gioannini, T. L., and Weiss, J. P. (2012). Endotoxin{middle dot}albumin complexes transfer endotoxin monomers to MD-2 resulting in activation of TLR4. Innate Immun. 18, 478–491. doi: 10.1177/1753425911422723

PubMed Abstract | Crossref Full Text | Google Scholar

Fleck, A., Raines, G., Hawker, F., Trotter, J., Wallace, P. I., Ledingham, I. M., et al. (1985). Increased vascular permeability: a major cause of hypoalbuminaemia in disease and injury. Lancet 1, 781–784. doi: 10.1016/s0140-6736(85)91447-3

PubMed Abstract | Crossref Full Text | Google Scholar

Gichoya, J. W., Banerjee, I., Bhimireddy, A. R., Burns, J. L., Celi, L. A., Chen, L. C., et al. (2022). AI recognition of patient race in medical imaging: a modelling study. Lancet Digit Health 4, e406–e414. doi: 10.1016/s2589-7500(22)00063-2

PubMed Abstract | Crossref Full Text | Google Scholar

Gioannini, T. L., Zhang, D., Teghanemt, A., and Weiss, J. P. (2002). An essential role for albumin in the interaction of endotoxin with lipopolysaccharide-binding protein and sCD14 and resultant cell activation. J. Biol. Chem. 277, 47818–47825. doi: 10.1074/jbc.M206404200

PubMed Abstract | Crossref Full Text | Google Scholar

Gomez, M., Wong, D. T., Stewart, T. E., Redelmeier, D. A., and Fish, J. S. (2008). The FLAMES score accurately predicts mortality risk in burn patients. J. Trauma 65, 636–645. doi: 10.1097/TA.0b013e3181840c6d

PubMed Abstract | Crossref Full Text | Google Scholar

Greenhalgh, D. G. (2017). Sepsis in the burn patient: a different problem than sepsis in the general population. Burns Trauma 5, 23. doi: 10.1186/s41038-017-0089-5

PubMed Abstract | Crossref Full Text | Google Scholar

Greenhalgh, D. G., Saffle, J. R., Holmes, J. H. 4., Gamelli, R. L., Palmieri, T. L., Horton, J. W., et al. (2007). American Burn Association consensus conference to define sepsis and infection in burns. J. Burn Care Res. 28, 776–790. doi: 10.1097/BCR.0b013e3181599bc9

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, F., Zhu, Y. S., Huang, J., Wu, Y. H., Sun, Z. F., Xia, X. B., et al. (2018). National experts consensus on clinical diagnosis and treatment of inhalation injury, (2018 version). Zhonghua Shao Shang Za Zhi 34, 770–775. doi: 10.3760/cma.j.issn.1009-2587.2018.11.010

PubMed Abstract | Crossref Full Text | Google Scholar

Haibo, H. and Garcia, E. A. (2009). Learning from imbalanced data. IEEE Trans. Knowledge Data Eng. 21, 1263–1284. doi: 10.1109/tkde.2008.239

Crossref Full Text | Google Scholar

He, Q. L., Gao, S. W., Qin, Y., Huang, R. C., Chen, C. Y., Zhou, F., et al. (2022). Gastrointestinal dysfunction is associated with mortality in severe burn patients: a 10-year retrospective observational study from South China. Mil Med. Res. 9, 49. doi: 10.1186/s40779-022-00403-1

PubMed Abstract | Crossref Full Text | Google Scholar

Hunag, Y. (2020). National experts consensus on prevention and treatment of burn shock, (2020 version). Zhonghua Shao Shang Za Zhi 36, 786–792. doi: 10.3760/cma.j.cn501120-20200623-00323

PubMed Abstract | Crossref Full Text | Google Scholar

Iqbal, T., Saaiq, M., and Ali, Z. (2013). Epidemiology and outcome of burns: Early experience at the country’s first national burns Centre. Burns. 39 (2), 358–362. doi: 10.1016/j.burns.2012.07.011

PubMed Abstract | Crossref Full Text | Google Scholar

Jeon, K., Lee, N., Jeong, S., Park, M.-J., and Song, W. (2021). Immature granulocyte percentage for prediction of sepsis in severe burn patients: a machine leaning-based approach. BMC Infect. Dis. 21, 1258. doi: 10.1186/s12879-021-06971-2

PubMed Abstract | Crossref Full Text | Google Scholar

Jeschke, M. G., Finnerty, C. C., Kulp, G. A., Kraft, R., and Herndon, D. N. (2013). Can we use C-reactive protein levels to predict severe infection or sepsis in severely burned patients? Int. J. Burns Trauma 3, 137–143.

Google Scholar

Jeschke, M. G., van Baar, M. E., Choudhry, M. A., Chung, K. K., Gibran, N. S., and Logsetty, S. (2020). Burn injury. Nat. Rev. Dis. Primers 6, 11. doi: 10.1038/s41572-020-0145-5

PubMed Abstract | Crossref Full Text | Google Scholar

Knaus, W. A., Draper, E. A., Wagner, D. P., and Zimmerman, J. E. (1985). APACHE II: a severity of disease classification system. Crit. Care Med. 13, 818–829. doi: 10.1097/00003246-198510000-00009

Crossref Full Text | Google Scholar

Lavrentieva, A., Kontakiotis, T., Lazaridis, L., Tsotsolis, N., Koumis, J., Kyriazis, G., et al. (2007). Inflammatory markers in patients with severe burn injury. What is the best indicator of sepsis? Burns 33, 189–194. doi: 10.1016/j.burns.2006.07.001

PubMed Abstract | Crossref Full Text | Google Scholar

Lehnhardt, M., Jafari, H. J., Druecke, D., Steinstraesser, L., Steinau, H. U., Klatte, W., et al. (2005). A qualitative and quantitative analysis of protein loss in human burn wounds. Burns 31, 159–167. doi: 10.1016/j.burns.2004.08.015

PubMed Abstract | Crossref Full Text | Google Scholar

Li, A. T., Moussa, A., Gus, E., Paul, E., Yii, E., Romero, L., et al. (2022). Biomarkers for the early diagnosis of sepsis in burns: systematic review and meta-analysis. Ann. Surg. 275, 654–662. doi: 10.1097/sla.0000000000005198

PubMed Abstract | Crossref Full Text | Google Scholar

Lionelli, G. T., Pickus, E. J., Beckum, O. K., Decoursey, R. L., and Korentager, R. A. J. B. (2005). A three decade analysis of factors affecting burn mortality in the elderly. Burns 31, 958–963. doi: 10.1016/j.burns.2005.06.006

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, T., Xie, J., Yang, F., Chen, J. J., Li, Z. F., Yi, C. L., et al. (2015). The influence of sex on outcomes in trauma patients: a meta-analysis. Am. J. Surg. 210, 911–921. doi: 10.1016/j.amjsurg.2015.03.021

PubMed Abstract | Crossref Full Text | Google Scholar

Lundberg, S. M., Nair, B., Vavilala, M. S., Horibe, M., Eisses, M. J., Adams, T., et al. (2018). Explainable machine-learning predictions for the prevention of hypoxaemia during surgery. Nat. Biomed. engineering. 10), 2. doi: 10.1038/s41551-018-0304-0

PubMed Abstract | Crossref Full Text | Google Scholar

Mann-Salinas, E. A., Baun, M. M., Meininger, J. C., Murray, C. K., Aden, J. K., Wolf, S. E., et al. (2013). Novel predictors of sepsis outperform the American Burn Association sepsis criteria in the burn intensive care unit patient. J. Burn Care Res. 34, 31–43. doi: 10.1097/BCR.0b013e31826450b5

PubMed Abstract | Crossref Full Text | Google Scholar

Ming, J., Lei, P., Duan, J. L., Tan, J. H., Lou, H. P., Di, D. Y., et al. (2018). National experts consensus on tracheotomy and intubation for burn patients, (2018 version). Zhonghua Shao Shang Za Zhi 34, E006. doi: 10.3760/cma.j.issn.1009-2587.2018.11.E006

PubMed Abstract | Crossref Full Text | Google Scholar

Nielson, C. B., Duethman, N. C., Howard, J. M., Moncure, M., and Wood, J. G. (2017). Burns: pathophysiology of systemic complications and current management. J. Burn Care Res. 38, e469–e481. doi: 10.1097/bcr.0000000000000355

PubMed Abstract | Crossref Full Text | Google Scholar

Raith, E. P., Udy, A. A., Bailey, M., McGloughlin, S., MacIsaac, C., Bellomo, R., et al. (2017). Prognostic accuracy of the SOFA score, SIRS criteria, and qSOFA score for in-hospital mortality among adults with suspected infection admitted to the intensive care unit. JAMA 317, 290–300. doi: 10.1001/jama.2016.20328

PubMed Abstract | Crossref Full Text | Google Scholar

Rimmer, E., Garland, A., Kumar, A., Doucette, S., Houston, B. L., Menard, C. E., et al. (2022). White blood cell count trajectory and mortality in septic shock: a historical cohort study. Can. J. Anesth. 69, 1230–1239. doi: 10.1007/s12630-022-02282-5

PubMed Abstract | Crossref Full Text | Google Scholar

Saracoglu, A., Kuzucuoglu, T., Yakupoglu, S., Kilavuz, O., Tuncay, E., Ersoy, B., et al. (2014). Prognostic factors in electrical burns: a review of 101 patients. Burns 40, 702–707. doi: 10.1016/j.burns.2013.08.023

PubMed Abstract | Crossref Full Text | Google Scholar

Schmidt, S. V., Drysch, M., Reinkemeier, F., Wagner, J. M., Sogorski, A., Macedo Santos, E., et al. (2023). Improvement of predictive scores in burn medicine through different machine learning approaches. Healthcare (Basel) 11 (17), 2437. doi: 10.3390/healthcare11172437

PubMed Abstract | Crossref Full Text | Google Scholar

Singer, M., Deutschman, C. S., Seymour, C. W., Shankar-Hari, M., Annane, D., Bauer, M., et al. (2016). The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA 315, 801–810. doi: 10.1001/jama.2016.0287

PubMed Abstract | Crossref Full Text | Google Scholar

Tan Chor Lip, H., Tan, J. H., Thomas, M., Imran, F.-H., and Azmah Tuan Mat, T. N. (2019). Survival analysis and mortality predictors of hospitalized severe burn victims in a Malaysian burns intensive care unit. Burns Trauma 7, 3. doi: 10.1186/s41038-018-0140-1

PubMed Abstract | Crossref Full Text | Google Scholar

Taverna, M., Marie, A. L., Mira, J. P., and Guidet, B. (2013). Specific antioxidant properties of human serum albumin. Ann. Intensive Care 3, 4. doi: 10.1186/2110-5820-3-4

PubMed Abstract | Crossref Full Text | Google Scholar

Temiz, A., Albayrak, A., Peksöz, R., Dışcı, E., Korkut, E., Tanrıkulu, Y., et al. (2020). Factors affecting the mortality at patients with burns: Single centre results. Ulus Travma Acil Cerrahi Derg 26, 777–783. doi: 10.14744/tjtes.2020.37862

PubMed Abstract | Crossref Full Text | Google Scholar

Thomas, S., Wolf, S. E., Chinkes, D. L., and Herndon, D. N. (2004). Recovery from the hepatic acute phase response in the severely burned and the effects of long-term growth hormone treatment. Burns 30, 675–679. doi: 10.1016/j.burns.2004.03.006

PubMed Abstract | Crossref Full Text | Google Scholar

Tobiasen, J., Hiebert, J. M., and Edlich, R. F. (1982). The abbreviated burn severity index. Ann. Emerg. Med. 11, 260–262. doi: 10.1016/s0196-0644(82)80096-6

PubMed Abstract | Crossref Full Text | Google Scholar

Torres, M. J. M., Peterson, J. M., and Wolf, S. E. (2021). Detection of infection and sepsis in burns. Surg. Infections 22, 20–27. doi: 10.1089/sur.2020.348

PubMed Abstract | Crossref Full Text | Google Scholar

Tran, N. K., Albahra, S., Pham, T. N., Holmes, J. H., Greenhalgh, D., Palmieri, T. L., et al. (2020). Novel application of an automated-machine learning development tool for predicting burn sepsis: proof of concept. Sci. Rep. 10, 12354. doi: 10.1038/s41598-020-69433-w

PubMed Abstract | Crossref Full Text | Google Scholar

Tsolakidis, S., Freytag, D. L., Dovern, E., Alharbi, Z., Kim, B.-S., Houschyar, K. S., et al. (2022). Infections in burn patients: A retrospective view over seven years. Medicina 58, 1066. doi: 10.3390/medicina58081066

PubMed Abstract | Crossref Full Text | Google Scholar

World Health Organization. (2002). The Injury Chart Book: A Graphical Overview of the Global Burden of Injuries. (Geneva: World Health Organization).

Google Scholar

Yan, J., Hill, W. F., Rehou, S., Pinto, R., Shahrokhi, S., and Jeschke, M. G. (2018). Sepsis criteria versus clinical diagnosis of sepsis in burn patients: A validation of current sepsis scores. Surgery 164, 1241–1245. doi: 10.1016/j.surg.2018.05.053

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, G., Tan, L., Yao, H., Xiong, Z., Wu, J., and Huang, X. (2023). Long-term effects of severe burns on the kidneys: research advances and potential therapeutic approaches. J. Inflammation Res. Volume 16, 1905–1921. doi: 10.2147/jir.s404983

PubMed Abstract | Crossref Full Text | Google Scholar

Yazıcı, H., Uçar, A. D., Namdaroğlu, O., and Yıldırım, M. (2022). Mortality prediction models for severe burn patients: Which one is the best? Ulus Travma Acil Cerrahi Derg 28, 790–795. doi: 10.14744/tjtes.2021.29540

PubMed Abstract | Crossref Full Text | Google Scholar

Yoon, J., Kym, D., Hur, J., Kim, Y., Yang, H. T., Yim, H., et al. (2018). Comparative usefulness of sepsis-3, burn sepsis, and conventional sepsis criteria in patients with major burns. Crit. Care Med. 46, e656–e662. doi: 10.1097/ccm.0000000000003144

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, P., Zou, B., Liou, Y. C., and Huang, C. (2021). The pathogenesis and diagnosis of sepsis post burn injury. Burns Trauma 9, tkaa047. doi: 10.1093/burnst/tkaa047

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: extensive burns, machine learning, sepsis, mortality, predictive model, SHAP, XGBoost, GBT

Citation: Wang S-Q, Qiu K, Zheng Q-R, Zhou B-J, Li M-Y, Zhong H-Y, Chen Y and Yuan S-M (2025) Development and validation of web-based, interpretable predictive models for sepsis and mortality in extensive burns. Front. Cell. Infect. Microbiol. 15:1586087. doi: 10.3389/fcimb.2025.1586087

Received: 02 March 2025; Accepted: 29 July 2025;
Published: 18 August 2025.

Edited by:

Stefano Busani, University Hospital of Modena, Italy

Reviewed by:

Bogdan Pavel, Carol Davila University of Medicine and Pharmacy, Romania
Martina Schiavello, University of Turin, Italy
Cristian Sorin Hariga, Carol Davila University of Medicine and Pharmacy, Romania

Copyright © 2025 Wang, Qiu, Zheng, Zhou, Li, Zhong, Chen and Yuan. 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: Si-Ming Yuan, eXVhbnNtQDE2My5jb20=

†These authors share first authorship

ORCID: Si-Ming Yuan, orcid.org/0000-0003-1380-5425

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.