Association between the stress hyperglycemia ratio and 28-day all-cause mortality in critically ill patients with sepsis: a retrospective cohort study and predictive model establishment based on machine learning | Cardiovascular Diabetology | Full Text
Skip to main content

Association between the stress hyperglycemia ratio and 28-day all-cause mortality in critically ill patients with sepsis: a retrospective cohort study and predictive model establishment based on machine learning

Abstract

Background

Sepsis is a severe form of systemic inflammatory response syndrome that is caused by infection. Sepsis is characterized by a marked state of stress, which manifests as nonspecific physiological and metabolic changes in response to the disease. Previous studies have indicated that the stress hyperglycemia ratio (SHR) can serve as a reliable predictor of adverse outcomes in various cardiovascular and cerebrovascular diseases. However, there is limited research on the relationship between the SHR and adverse outcomes in patients with infectious diseases, particularly in critically ill patients with sepsis. Therefore, this study aimed to explore the association between the SHR and adverse outcomes in critically ill patients with sepsis.

Methods

Clinical data from 2312 critically ill patients with sepsis were extracted from the MIMIC-IV (2.2) database. Based on the quartiles of the SHR, the study population was divided into four groups. The primary outcome was 28-day all-cause mortality, and the secondary outcome was in-hospital mortality. The relationship between the SHR and adverse outcomes was explored using restricted cubic splines, Cox proportional hazard regression, and Kaplan‒Meier curves. The predictive ability of the SHR was assessed using the Boruta algorithm, and a prediction model was established using machine learning algorithms.

Results

Data from 2312 patients who were diagnosed with sepsis were analyzed. Restricted cubic splines demonstrated a "U-shaped" association between the SHR and survival rate, indicating that an increase in the SHR is related to an increased risk of adverse events. A higher SHR was significantly associated with an increased risk of 28-day mortality and in-hospital mortality in patients with sepsis (HR > 1, P < 0.05) compared to a lower SHR. Boruta feature selection showed that SHR had a higher Z score, and the model built using the rsf algorithm showed the best performance (AUC = 0.8322).

Conclusion

The SHR exhibited a U-shaped relationship with 28-day all-cause mortality and in-hospital mortality in critically ill patients with sepsis. A high SHR is significantly correlated with an increased risk of adverse events, thus indicating that is a potential predictor of adverse outcomes in patients with sepsis.

Introduction

Infectious diseases have long been at the forefront of global health challenges. Sepsis, which is a severe complication of infections, is closely associated with abnormal immune responses. Sepsis often leads to multiorgan dysfunction, thereby posing a threat to the lives of affected individuals [1]. Despite significant progress in the management of infectious diseases due to widespread antibiotic use and advancements in medical technology, the incidence of sepsis remains high, and its mortality rate is a cause for concern [2].

Stress hyperglycemia is a physiological phenomenon characterized by a significant increase in blood glucose levels in response to physiological or pathological stress. This phenomenon is often associated with the regulation of hormone secretion, immune responses, and neural system activity. Stress hyperglycemia is typically linked to the release of hormones such as adrenaline and cortisol; furthermore, it is associated with insulin resistance and increased hepatic gluconeogenesis [3]. In some cases, stress hyperglycemia is an adaptive physiological response that helps to provide energy for the body to cope with physiological or pathological stress. However, for patients with chronic diseases or those facing prolonged stress, this hyperglycemic state may be associated with unfavorable clinical outcomes. A previous study indicated an independent correlation between acute hyperglycemia upon admission and adverse outcomes in the early and late stages of acute myocardial infarction (AMI) [4]. To explain and evaluate the impact of background blood glucose, researchers introduced HbA1c as a baseline glucose level when assessing stress hyperglycemia, and they proposed the stress hyperglycemia ratio (SHR) as a new indicator to assess a more accurate representation of acute hyperglycemia [5]. Previous studies have shown that in patients with acute coronary syndrome, critically ill patients, and patients with critical myocardial infarction, a higher SHR is associated with an increased risk of adverse outcomes [6,7,8].

Due to the abnormal activation of the immune system and the massive release of various cytokines in patients with sepsis, they often experience a state of stress [9]. However, it is currently unclear whether the SHR is correlated with adverse outcomes in critically ill patients with sepsis. Therefore, this study aimed to assess the association between the SHR and adverse outcomes in critically ill patients with sepsis.

Methods

Data source

The data utilized in this study are derived from MIMIC-IV2.2, which is an electronic health record dataset encompassing over 50,000 patients admitted to the intensive care unit (ICU) at Beth Israel Deaconess Medical Center (BIDMC) in Boston, Massachusetts, from 2008 to 2019. The Institutional Review Board of BIDMC granted a waiver of informed consent and approved the sharing of research resources. The author (JLZ) obtained access to the database (certificate number: 45848364).

Inclusion and exclusion criteria

Inclusion Criteria:

  1. 1.

    Patients aged between 18 and 90 years.

  2. 2.

    Patients diagnosed with sepsis according to the recommendations of the Third International Consensus Definitions for Sepsis and Septic Shock (Sepsis-3).

Exclusion Criteria:

  1. 1.

    Patients with an ICU stay of less than 24 h.

  2. 2.

    Missing serum glucose and glycated hemoglobin in the first laboratory test.

  3. 3.

    For patients with multiple ICU admissions, only data from the first hospitalization were included.

Outcome

The primary outcome was 28-day all-cause mortality, and the secondary outcome was in-hospital mortality.

Data extraction

Data extraction was performed using pgAdmin software. Patient characteristics, including age, sex, and weight, were collected. Information on comorbidities, such as hypertension, type 2 diabetes, type 1 diabetes, heart failure, malignancy, chronic kidney disease (CKD), stroke, pneumonia and septic shock, was extracted based on the International Classification of Diseases coding system. Vital signs (heart rate, systolic blood pressure (SBP), respiratory rate, arterial oxygen saturation (SaO2) and laboratory tests [red blood cell count (RBC), white blood cell count (WBC), platelet count (PLT), serum sodium, serum potassium, serum calcium, anion gap, pH, carbon dioxide pressure (PaCO2), arterial oxygen pressure (PaO2), lactate, prothrombin time international normalized ratio (INR), total bilirubin, aspartate aminotransferase (AST), blood urea nitrogen, serum creatinine, serum glucose, and glycosylated hemoglobin (HbA1c)] were extracted. The use of steroids (glucocorticoids) was also extracted. Additionally, the severity of illness was assessed using the sequential organ failure assessment (SOFA). The stress hyperglycemia ratio (SHR) is defined as the index calculated using the following formula: SHR = (admission blood glucose (mg/dl))/(28.7 × HbA1c(%) − 46.7)] [5].

Statistical analysis

As the current study is a retrospective analysis, no sample size calculations were conducted. Variables with a missing data rate exceeding 20% were excluded, and for those with less than a missing data rate below 20%, multiple imputation was employed. The variance inflation factor (VIF) was used to assess multicollinearity among variables. Variables with a VIF exceeding 5 were removed due to multicollinearity concerns. Patients were categorized into four groups based on quartiles of the SHR. Normally distributed continuous variables are presented as the means (standard deviations [SDs]) and were analyzed using analysis of variance (ANOVA). Nonnormally distributed variables were analyzed using the Mann‒Whitney U test or the Kruskal‒Wallis test. Categorical variables are expressed as numbers and percentages and were analyzed using the χ2 test or Fisher’s exact test. Kaplan‒Meier survival curves were used to compare the 28-day survival rates among the four groups according to the log-rank test. Proportional hazard regression models (Cox regression models) were used to assess the hazard ratio (HR) and 95% confidence interval (95% CI) for event occurrence. Model I was unadjusted for covariates. To adjust for the impacts of patient general condition and vital signs, SOFA scores, and steroid use on outcomes, Model II was adjusted for age, weight, sex, heart rate, respiratory rate, systolic blood pressure, SOFA scores, and steroid use (glucocorticoids). The time-dependent receiver operating characteristic (ROC) curve was used to compare the SHR and laboratory indicators of each continuous variable. A two-tailed P value < 0.05 was considered to indicate statistical significance. Statistical analyses were conducted using R software (version 4.3.1).

Restricted cubic splines

In this study, we collected data on survival (the outcome variable); the SHR (the continuous predictor variable); and age, weight, heart rate, SBP, respiratory rate, sex, SOFA scores and steroid use (the covariates). The potential nonlinear relationships between the change in the SHR and survival were examined by a Cox regression model with restricted cubic spline (RCS). Knots between 3 and 7 were tested, and the model with the lowest Akaike information criterion value was selected for the RCS.

Subgroup analysis

A subgroup analysis was conducted based on prespecified criteria, including age, sex, and type 2 diabetes status, and univariate analysis and multivariate analysis were performed. The multivariate analysis was adjusted for age, weight, sex, heart rate, respiratory rate, systolic blood pressure, SOFA scores, and steroids (glucocorticoids). The multivariable analysis of the male and female subgroups were not adjusted for sex. Patients were stratified into two groups based on age (< 65 years and ≥ 65 years). Cox proportional hazards regression analysis was performed for each subgroup, and the results were visually presented using forest plots, illustrating hazard ratios (HRs) and 95% confidence intervals (CIs).

Establishment and validation of the prediction models

Boruta’s algorithm is a method used to determine the most important features in a dataset. It identifies importance by comparing the Z value of each feature with the Z value of the corresponding “shadow feature”. In the algorithm, all real features are copied and shuffled, and then the Z value of each feature is obtained through the random forest model. Additionally, the Z values of the ‘shadow features’ are generated by randomly shuffling the real features [10]. If the Z value of a true feature is significantly higher than the maximum Z value of the shadow feature in multiple independent tests, the true feature is marked as “important” (green area), also known as an acceptable variable. Otherwise, it is marked as “unimportant” (red area), also known as unacceptable variables. Acceptable variables are variables that are retained during the feature selection process and are considered to contribute to the performance of the model. Unacceptable variables are excluded from the final feature selection by the algorithm because they fail to show predictive power for the target variable during the feature selection process. In addition, the Boruta algorithm was also used to explore the importance of the SHR as a predictor variable.

Acceptable variables are incorporated into the machine learning algorithm. The dataset was divided into training and validation sets at a 7:3 ratio. The filtered variables were individually analyzed using the Cox Proportional Hazards Survival Learner (coxph) algorithm, Rpart Survival Trees Survival Learner (dt) algorithm, Survival DeepSurv Learner (deepsurv) algorithm, Survival Random Forest SRC Learner (rsf), and Extreme Gradient Boosting Survival Learner (xgboost) algorithm to predict the 28-day mortality risk in critically ill sepsis patients. Hyperparameter tuning is performed during the establishment of machine learning models. The training and validation sets were utilized for model establishment and evaluation, respectively. The ROC curve and its corresponding area under the curve (AUC) were used to determine model performance. Decision curve analysis (DCA) was employed for assessing clinical effectiveness, while calibration curves were used to evaluate the accuracy of the model in predicting absolute risk.

Results

Baseline characteristics

Data from 2,312 patients diagnosed with sepsis were extracted from the MIMIC-IV (Fig. 1). Table S1 displays the variance inflation factors, indicating that there was no multicollinearity among the variables. Figure S1 illustrates the proportions of missing data for each variable, and Table 1 presents the baseline characteristics of the study subjects. There were 1396 males (60.38%), 1046 (45.24%) patients with hypertension, 906 (39.19%) patients with type 2 diabetes, 74 (3.2%) patients with type 1 diabetes, 805 (34.82%) patients with heart failure, 248 (10.73%) patients with malignancy, 450 (19.46%) patients with CKD, 327 (14.14%) patients with stroke, 372 (16.09%) patients with septic shock and 925 (40.01%) patients with pneumonia. Patients were divided into quartiles: quartile 1 (0.21 ≤ SHR < 0.915), quartile 2 (0.915 ≤ SHR < 1.14), quartile 3 (1.14 ≤ SHR < 1.45), and quartile 4 (1.45 ≤ SHR ≤ 7.41), with each group consisting of 578 individuals. Patients in Quartile 4 exhibited a higher heart rate, white blood cell count, serum potassium, anion gap, total bilirubin, aspartate aminotransferase, blood urea nitrogen, serum creatinine, lactate, SOFA score, and steroid use ratio as well as lower serum sodium, pH, and PaO2 levels.

Fig. 1
figure 1

Selection of the study population from the MIMIC-IV database

Table 1 Patient demographics and baseline characteristics

Clinical outcomes

Regarding 28-day mortality and in-hospital mortality, Quartile 4 had higher mortality (Table 2). In the Cox regression analysis, the results of Models I and II showed that with Quartile 1 as the reference, the risk of mortality in Quartile 3 and Quartile 4 increased significantly (Table 3, Table 4). Kaplan–Meier curves revealed that patients in Quartile 4 had the lowest 28-day survival probability, and the difference was significant (Fig. 2).

Table 2 28-day all-cause mortality and in-hospital mortality
Table 3 Cox regression model (28-day all-cause mortality)
Table 4 Cox regression model (in-hospital mortality)
Fig. 2
figure 2

28-day KM survival curve. KM curves showing the survival rates at 28 days for each quartile. SHR: Quartile 1 (0.21–0.915), Quartile 2 (0.915–1.14), Quartile 3 (1.14–1.45), and Quartile 4 (1.45–7.41)

Restricted cubic spline

The RCS analysis was adjusted for the effects of age, weight, heart rate, systolic blood pressure, respiratory rate, sex, SOFA scores and steroid use. The RCS analysis for 28-day all-cause mortality (Fig. 3) and in-hospital mortality (Fig. 4) both indicated a U-shaped association between the SHR and mortality risk. The turning point of the RCS curve is approximately SHR = 0.85, representing the inflection point in the relationship between the SHR and mortality risk. According to the Cox proportional hazards regression model, the turning point lies in Quartile 1, where the mortality risk is minimized, thus demonstrating consistency between the two results.

Fig. 3
figure 3

RCS analysis of 28-day all-cause mortality. Curves represent estimated adjusted hazard ratios, and shaded ribbons represent 95% confidence intervals. The vertical dotted line represents the lowest point of the curve (SHR = 0.85), which represents the lowest hazard ratio. The horizontal dashed line represents a hazard ratio of 1.0. HR hazard ratio, CI confidence interval

Fig. 4
figure 4

RCS results for in-hospital mortality. Curves represent estimated adjusted hazard ratios, and shaded ribbons represent 95% confidence intervals. The vertical dotted line represents the lowest point of the curve, which represents the lowest hazard ratio. The horizontal dashed line represents a hazard ratio of 1.0. HR hazard ratio; CI, confidence interval

Subgroup analysis

The results presented a subgroup analysis of 28-day all-cause mortality (Fig. 5). In subgroups defined by age < 65, age ≥ 65, male sex, female sex, and type 2 diabetes status, Quartile 4 consistently demonstrated a greater risk of mortality, regardless of whether covariates were adjusted. This finding indicates that, irrespective of baseline levels, a higher SHR is associated with an increased 28-day mortality risk in critically ill sepsis patients (HR > 1 in each subgroup, P < 0.05).

Fig. 5
figure 5

Subgroup forest plot for 28-day all-cause mortality. Adjusted for age, weight, sex, heart rate, respiratory rate, systolic blood pressure, SOFA scores, and use of steroids (glucocorticoids). SHR: Quartile 1 (0.21–0.915), Quartile 2 (0.915–1.14), Quartile 3 (1.14–1.45), and Quartile 4 (1.45–7.41). Diabetes I type 1 diabetes, Diabetes II type 2 diabetes, HR hazard ratio, CI confidence interval

Boruta Algorithm

Figure 6 shows the feature selection results based on the Boruta algorithm. Variables in the green area are identified as important features, and variables in the red area are unimportant features in the Boruta algorithm.

Fig. 6
figure 6

Feature selection based on the Boruta algorithm. The horizontal axis is the name of each variable, and the vertical axis is the Z value of each variable. The box plot shows the Z value of each variable during model calculation. The green boxes represent important variables, and the red boxes represent unimportant variables

Establishment and validation of the prediction model

Table S2 shows the hyperparameter, tuning scope and optimal hyperparameter of the four models. Figure 7 displays the ROC curves of various models, and model performance is represented by AUC values. The AUC of coxph was 0.8242, the AUC of dt was 0.7277, the AUC of deepsurv was 0.7393, the AUC of rsf was 0.8322, and the AUC of xgboost was 0.8068. Figure S2 shows the calibration curve of each model. The calibration curves of the coxph, dt, deepsurv, rsf, and xgboost models agree well with the reference line, indicating that they have excellent prediction performance. According to the DCA curve (Figure S3), each model showed a large net benefit, indicating that the established model has robust clinical validity. Figure 8 demonstrates the performance of the SHR and continuous laboratory data. The AUC of the SHR was 0.7081, which was higher than that of blood glucose (AUC = 0.6302) and HbA1c (AUC = 0.5149).

Fig. 7
figure 7

ROC curves of the machine learning algorithms. coxph Cox proportional hazards survival learner, dt Rpart Survival Trees Survival Learner, deepsurv Survival DeepSurv Learner, rsf Survival Random Forest SRC Learner, xgboost extreme gradient boosting survival learner, T days, AUC area under the curve

Fig. 8
figure 8

Performance of the SHR and laboratory data. WBC white blood cell count, RBC red blood cell count, INR prothrombin time international normalized ratio, AST aspartate aminotransferase, HbA1c glycosylated hemoglobin, SHR stress hyperglycemia ratio, AUC area under the curve

Discussion

The results of this study indicate a significant association between higher SHR and increased 28-day mortality and in-hospital mortality in critically ill patients with sepsis. These findings remained consistent across age, sex, and subgroup with type 2 diabetes after adjusting for covariates, demonstrating the robustness of the study results. To date, this is the first study of the relationship between the SHR and adverse outcomes in critically ill patients with sepsis.

The Boruta algorithm has become a widely used method in feature selection, determining which features are most important for predicting the target variable by simulating randomness [10]. The feature selection results of the Boruta algorithm in this study suggest that SHR significantly occupies the green area, exhibiting a high Z score in feature selection. This indicates that the SHR may play a crucial role in this study, showing a significant association with the study objectives.

The Boruta algorithm results indicate that the SHR plays an important role in predicting 28-day all-cause mortality in patients with sepsis, but we also recognize that this does not mean that it is a decisive factor. First, the Boruta algorithm is a powerful feature selection method, but it may also be affected by correlations between data features. Therefore, even though the SHR is highly important in the model, this does not necessarily mean that it is the most decisive factor. Second, we found through Cox regression analysis that a greater SHR is associated with an increased risk of 28-day mortality in patients with sepsis, which is consistent with the Boruta algorithm suggesting that the SHR is an important feature. This finding provides evidence that the SHR can be used as a predictor of 28-day all-cause mortality in patients with sepsis. Therefore, we believe that the SHR can be used as a predictor of 28-day all-cause mortality in patients with sepsis.

By incorporating acceptable variables into various machine learning algorithms, the results suggest that these predictive models all demonstrate good performance. It can be reasonably inferred that the stress hyperglycemia ratio may be a significant predictive indicator or influencing factor in this study.

Relation to previous research

Previous studies have predominantly focused on vascular diseases. Regarding cerebrovascular diseases, research suggests a significant positive correlation between higher SHR quartiles and an increased risk of ischemic stroke transformation to hemorrhage in patients [11]. Additionally, the SHR is independently associated with the severity of post-acute ischemic stroke brain edema, unfavorable functional outcomes, and post-acute ischemic stroke mortality [12]. In cardiovascular diseases, the results of a multicenter prospective study in China demonstrated a correlation between high SHR and increased long-term mortality [13]. In patients without preexisting diabetes, stress hyperglycemia has been shown to be associated with adverse outcomes in stroke [14] and acute myocardial infarction patients [15]. However, stress hyperglycemia does not necessarily represent acute-onset elevated blood glucose; it may reflect poor blood glucose control in patients’ medical history. Therefore, solely exploring the relationship between blood glucose at admission and disease status is not rigorous. Investigating the relationship between the SHR and disease by introducing HbA1c for the correction of past blood glucose levels is a more reasonable approach.

SHR and sepsis

An elevated SHR indicates a hyperglycemic stress state regardless of past blood glucose levels. Our study confirmed this finding by demonstrating a significant association between a greater SHR and adverse outcomes in critically ill patients with sepsis, irrespective of whether they had type 2 diabetes. This study emphasizes the relationship between the SHR and sepsis, a highlight of our research. First, we explored the potential mechanisms underlying the increased risk of mortality in sepsis patients with high SHRs. Animal experiments suggest that a three-hour infusion of glucose leading to hyperglycemia significantly impairs immune function [16] and activates cytokines while promoting oxidative stress in the liver [17]. High glucose levels can stimulate monocytes to enhance the synthesis and release of interleukin-6 (IL-6) [18], and elevated serum IL-6 exacerbates insulin resistance [19], prompting the liver to release glucose and contributing to hyperglycemia [20]. IL-6 is a pleiotropic cytokine involved in various pathophysiological processes, including inflammation and tissue damage [21], and has been shown to correlate with adverse outcomes in critically ill patients [22]. Additionally, tumor necrosis factor-alpha (TNF-α), a major proinflammatory factor in the pathogenesis of sepsis, is associated with adverse outcomes [23] and can induce insulin resistance [24]. Second, high blood glucose itself may have proinflammatory effects, as lowering blood glucose to 110 mg/dL has been demonstrated to have anti-inflammatory effects in critically ill patients. Our results also show that the risk of adverse events in sepsis patients is minimized when the SHR is approximately 0.85, indicating a potential anti-inflammatory effect of a lower SHR, which plays a role in reducing the occurrence of adverse events.

Furthermore, another critical factor leading to increased mortality in sepsis patients is the occurrence of disseminated intravascular coagulation (DIC). A large, multicenter, prospective study in Japan revealed that the mortality rate of sepsis patients with DIC was twice that of those without DIC [25]. In vitro studies suggest that acute elevation of glucose can cause endothelial damage, promoting abnormal activation of intravascular coagulation, which may lead to sepsis-associated DIC [26]. This could be one of the reasons why a high SHR increases the occurrence of adverse events. Elevated blood glucose promotes the aggregation of monocytes and macrophages [27], leading to the production and release of biologically active molecules such as cytokines, IL-6, and interleukin-8. These factors contribute to inflammation and thrombus formation, thereby promoting DIC.

Additionally, patients with sepsis may develop a hyperglycemic state through multiple mechanisms. First, stress associated with critical illness is characterized by activation of the hypothalamic‒pituitary‒adrenal (HPA) axis and increased release of glucocorticoids, epinephrine, glucagon, and growth hormone [28,29,30,31], which can lead to the occurrence of stress hyperglycemia. Second, the release of multiple cytokines, such as interleukin 6 (IL-6) and tumor necrosis factor-α (TNF-α), independently and synergistically with catecholamines promotes hepatic glucose production, whereas interleukin 1 (IL-1) and TNF-α also inhibit insulin release, seemingly in a concentration-dependent manner [32], jointly promoting the occurrence of hyperglycemia. Third, glucose is usually absorbed across cell membranes via carrier-mediated facilitated transport systems [33], and glucocorticoids inhibit glucose transporter 4 (GLUT4), thereby impairing insulin-mediated glucose uptake in skeletal muscle and indirectly promoting the occurrence of hyperglycemia. The above mechanisms indicate that the occurrence of stress hyperglycemia is promoted in patients with sepsis via multiple mechanisms.

The higher SHR and increased risk of all-cause mortality in patients with sepsis may be attributed to the inflammatory response caused by blood glucose fluctuations. First, a higher SHR indicates stress hyperglycemia, which is the result of the complex interaction of hormones such as catecholamines, glucocorticoids, and cytokines [34, 35]. Stress hyperglycemia may lead to increased mitochondrial reactive oxygen species production in endothelial cells, which may cause endothelial dysfunction [36]. Second, stress hyperglycemia contributes to the nonenzymatic glycation of platelet glycoproteins and may be one of the causes of platelet activation, which may lead to an increased risk of thrombosis [37]. The above mechanisms may also contribute to the increased risk of death in patients with sepsis.

Impact on clinical practice

In a previous study on the association between SHR and critically ill patients with acute myocardial infarction, a higher SHR was only correlated with an increased risk of adverse events in nondiabetic patients, and similar results were not found in diabetic patients [8]. Our study revealed that a greater SHR is associated with an increased risk of 28-day all-cause mortality and in-hospital mortality in critically ill patients with sepsis, regardless of whether they have type 2 diabetes. This suggests that the SHR can be used to predict the risk of adverse events in a broader range of sepsis patients, providing a reliable indicator for clinicians in the diagnosis and treatment of critically ill patients with sepsis.

Limitations of the study

This study has several limitations. First, this was a retrospective study relying on past records, which may be subject to information bias. Second, not all potential confounding factors could be controlled, limiting causal inferences. Third, sample selection may be influenced by known or unknown factors, resulting in samples that are not representative in certain aspects, potentially affecting the external validity of the study.

Conclusion

In conclusion, irrespective of the presence of type 2 diabetes, the SHR exhibited a U-shaped relationship with 28-day all-cause mortality and in-hospital mortality in critically ill patients with sepsis. A higher SHR is significantly correlated with an increased risk of adverse events. The SHR can be used to predict adverse outcomes in critically ill patients with sepsis. However, multicenter, prospective studies are still needed to validate these results.

Availability of data and material

Publicly available datasets were analyzed in this study. These data can be found at https://mimic.mit.edu/.

References

  1. Singer M, Deutschman CS, Seymour CW, Shankar-Hari M, Annane D, Bauer M, Bellomo R, Bernard GR, Chiche JD, Coopersmith CM, Hotchkiss RS, Levy MM, Marshall JC, Martin GS, Opal SM, Rubenfeld GD, van der Poll T, Vincent JL, Angus DC. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA. 2016;315(8):801–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Rudd KE, Johnson SC, Agesa KM, Shackelford KA, Tsoi D, Kievlan DR, Colombara DV, Ikuta KS, Kissoon N, Finfer S, Fleischmann-Struzek C, Machado FR, Reinhart KK, Rowan K, Seymour CW, Watson RS, West TE, Marinho F, Hay SI, Lozano R, Lopez AD, Angus DC, Murray CJL, Naghavi M. Global, regional, and national sepsis incidence and mortality, 1990–2017: analysis for the Global Burden of Disease Study. Lancet. 2020;395(10219):200–11.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Marik PE, Bellomo R. Stress hyperglycemia: an essential survival response! Crit Care Med. 2013;41(6):e93–4.

    Article  PubMed  Google Scholar 

  4. Ferreira JA, Baptista RM, Monteiro SR, Gonçalves FM, Monteiro PF, Gonçalves LM. Admission hyperglycemia and all-cause mortality in diabetic and non-diabetic patients with acute myocardial infarction: a tertiary center analysis. Intern Emerg Med. 2021;16(8):2109–19.

    Article  PubMed  Google Scholar 

  5. Roberts GW, Quinn SJ, Valentine N, et al. Relative hyperglycemia, a marker of critical illness: introducing the stress hyperglycemia ratio. J Clin Endocrinol Metab. 2015;100(12):4490–7.

    Article  CAS  PubMed  Google Scholar 

  6. Yang J, Zheng Y, Li C, Gao J, Meng X, Zhang K, Wang W, Shao C, Tang YD. The impact of the stress hyperglycemia ratio on short-term and long-term poor prognosis in patients with acute coronary syndrome: insight from a large cohort study in Asia. Diabetes Care. 2022;45(4):947–56.

    Article  CAS  PubMed  Google Scholar 

  7. Zhang C, Shen HC, Liang WR, Ning M, Wang ZX, Chen Y, Su W, Guo TT, Hu K, Liu YW. Relationship between stress hyperglycemia ratio and allcause mortality in critically ill patients: Results from the MIMIC-IV database. Front Endocrinol (Lausanne). 2023;14:1111026.

    Article  PubMed  Google Scholar 

  8. Liu J, Zhou Y, Huang H, Liu R, Kang Y, Zhu T, Wu J, Gao Y, Li Y, Wang C, Chen S, Xie N, Zheng X, Meng R, Liu Y, Tan N, Gao F. Impact of stress hyperglycemia ratio on mortality in patients with critical acute myocardial infarction: insight from American MIMIC-IV and the Chinese CIN-II study. Cardiovasc Diabetol. 2023;22(1):281.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Hotchkiss RS, Monneret G, Payen D. Immunosuppression in sepsis: a novel understanding of the disorder and a new therapeutic approach. Lancet Infect Dis. 2013;13(3):260–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Yue S, Li S, Huang X, Liu J, Hou X, Zhao Y, Niu D, Wang Y, Tan W, Wu J. Machine learning for the prediction of acute kidney injury in patients with sepsis. J Transl Med. 2022;20(1):215.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Yuan C, Chen S, Ruan Y, Liu Y, Cheng H, Zeng Y, Chen Y, Cheng Q, Huang G, He W, He J. The stress hyperglycemia ratio is associated with hemorrhagic transformation in patients with acute ischemic stroke. Clin Interv Aging. 2021;16:431–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Deng Y, Wu S, Liu J, Liu M, Wang L, Wan J, Zhang S, Liu M. The stress hyperglycemia ratio is associated with the development of cerebral edema and poor functional outcome in patients with acute cerebral infarction. Front Aging Neurosci. 2022;14:936862.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Cui K, Fu R, Yang J, Xu H, Yin D, Song W, Wang H, Zhu C, Feng L, Wang Z, Wang Q, Lu Y, Dou K, Yang Y, CAMI Registry Investigators. Stress hyperglycemia ratio and long-term mortality after acute myocardial infarction in patients with and without diabetes: a prospective, nationwide, and multicentre registry. Diabetes Metab Res Rev. 2022;38(7):e3562.

    Article  CAS  PubMed  Google Scholar 

  14. Capes SE, Hunt D, Malmberg K, Pathak P, Gerstein HC. Stress hyperglycemia and prognosis of stroke in nondiabetic and diabetic patients: a systematic overview. Stroke. 2001;32(10):2426–32.

    Article  CAS  PubMed  Google Scholar 

  15. Malmberg K. Prospective randomised study of intensive insulin treatment on long term survival after acute myocardial infarction in patients with diabetes mellitus. DIGAMI (Diabetes Mellitus, Insulin Glucose Infusion in Acute Myocardial Infarction) Study Group. BMJ. 1997;314(7093):1512–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Kwoun MO, Ling PR, Lydon E, Imrich A, Qu Z, Palombo J, Bistrian BR. Immunologic effects of acute hyperglycemia in nondiabetic rats. JPEN J Parenter Enteral Nutr. 1997;21(2):91–5.

    Article  CAS  PubMed  Google Scholar 

  17. Ling PR, Mueller C, Smith RJ, Bistrian BR. Hyperglycemia induced by glucose infusion causes hepatic oxidative stress and systemic inflammation, but not STAT3 or MAP kinase activation in liver in rats. Metabolism. 2003;52(7):868–74.

    Article  CAS  PubMed  Google Scholar 

  18. Morohoshi M, Fujisawa K, Uchimura I, Numano F. Glucose-dependent interleukin 6 and tumor necrosis factor production by human peripheral blood monocytes in vitro. Diabetes. 1996;45(7):954–9.

    Article  PubMed  Google Scholar 

  19. Heesen M, Bloemeke B, Heussen N, Kunz D. Can the interleukin-6 response to endotoxin be predicted? Studies of the influence of a promoter polymorphism of the interleukin-6 gene, gender, the density of the endotoxin receptor CD14, and inflammatory cytokines. Crit Care Med. 2002;30(3):664–9.

    Article  CAS  PubMed  Google Scholar 

  20. Ritchie DG. Interleukin 6 stimulates hepatic glucose release from prelabeled glycogen pools. Am J Physiol. 1990;258(1 Pt 1):E57-64.

    CAS  PubMed  Google Scholar 

  21. Papanicolaou DA, Wilder RL, Manolagas SC, Chrousos GP. The pathophysiologic roles of interleukin-6 in human disease. Ann Intern Med. 1998;128(2):127–37.

    Article  CAS  PubMed  Google Scholar 

  22. Harbarth S, Holeckova K, Froidevaux C, Pittet D, Ricou B, Grau GE, Vadas L, Pugin J, Geneva Sepsis Network. Diagnostic value of procalcitonin, interleukin-6, and interleukin-8 in critically ill patients admitted with suspected sepsis. Am J Respir Crit Care Med. 2001;164(3):396–402.

    Article  CAS  PubMed  Google Scholar 

  23. Damas P, Canivet JL, de Groote D, Vrindts Y, Albert A, Franchimont P, Lamy M. Sepsis and serum cytokine concentrations. Crit Care Med. 1997;25(3):405–12.

    Article  CAS  PubMed  Google Scholar 

  24. Hotamisligil GS. The role of TNFalpha and TNF receptors in obesity and insulin resistance. J Intern Med. 1999;245(6):621–5.

    Article  CAS  PubMed  Google Scholar 

  25. Gando S, Saitoh D, Ogura H, Fujishima S, Mayumi T, Araki T, Ikeda H, Kotani J, Kushimoto S, Miki Y, Shiraishi S, Suzuki K, Suzuki Y, Takeyama N, Takuma K, Tsuruta R, Yamaguchi Y, Yamashita N, Aikawa N, Japanese Association for Acute Medicine Sepsis Registry Study Group. A multicenter, prospective validation study of the Japanese Association for Acute Medicine disseminated intravascular coagulation scoring system in patients with severe sepsis. Crit Care. 2013;17(3):R111.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Higgins SJ, De Ceunynck K, Kellum JA, Chen X, Gu X, Chaudhry SA, Schulman S, Libermann TA, Lu S, Shapiro NI, Christiani DC, Flaumenhaft R, Parikh SM. Tie2 protects the vasculature against thrombus formation in systemic inflammation. J Clin Invest. 2018;128(4):1471–84.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Vaidyula VR, Boden G, Rao AK. Platelet and monocyte activation by hyperglycemia and hyperinsulinemia in healthy subjects. Platelets. 2006;17(8):577–85.

    Article  CAS  PubMed  Google Scholar 

  28. Marik PE, Zaloga GP. Adrenal insufficiency in the critically ill: a new look at an old problem. Chest. 2002;122(5):1784–96.

    Article  PubMed  Google Scholar 

  29. Hart BB, Stanford GG, Ziegler MG, Lake CR, Chernow B. Catecholamines: study of interspecies variation. Crit Care Med. 1989;17(11):1203–22.

    Article  CAS  PubMed  Google Scholar 

  30. Van den Berghe G. Neuroendocrine pathobiology of chronic critical illness. Crit Care Clin. 2002;18(3):509–28.

    Article  PubMed  Google Scholar 

  31. Siegel JH, Cerra FB, Coleman B, Giovannini I, Shetye M, Border JR, McMenamy RH. Physiological and metabolic correlations in human sepsis. Invit Comment Surg. 1979;86(2):163–93.

    CAS  Google Scholar 

  32. Mehta VK, Hao W, Brooks-Worrell BM, Palmer JP. Low-dose interleukin 1 and tumor necrosis factor individually stimulate insulin release but in combination cause suppression. Eur J Endocrinol. 1994;130(2):208–14.

    Article  CAS  PubMed  Google Scholar 

  33. Shepherd PR, Kahn BB. Glucose transporters and insulin action–implications for insulin resistance and diabetes mellitus. N Engl J Med. 1999;341(4):248–57.

    Article  CAS  PubMed  Google Scholar 

  34. Barth E, Albuszies G, Baumgart K, Matejovic M, Wachter U, Vogt J, Radermacher P, Calzia E. Glucose metabolism and catecholamines. Crit Care Med. 2007;35(9 Suppl):S508–18.

    Article  CAS  PubMed  Google Scholar 

  35. Andrews RC, Walker BR. Glucocorticoids and insulin resistance: old hormones, new targets. Clin Sci (Lond). 1999;96(5):513–23.

    Article  CAS  PubMed  Google Scholar 

  36. Giacco F, Brownlee M. Oxidative stress and diabetic complications. Circ Res. 2010;107(9):1058–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Ferroni P, Basili S, Falco A, Davì G. Platelet activation in type 2 diabetes mellitus. J Thromb Haemost. 2004;2(8):1282–91.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

None.

Funding

This study was supported by the Scientific Research Projects of Medical and Health Institutions of Longhua District, Shenzhen (2021057).

Author information

Authors and Affiliations

Authors

Contributions

JLZ contributed to the study concept and study design. FJY performed the statistical analysis and data interpretation. FJY, XHC and XQQ were responsible for the quality control of the data and algorithms. LLW and XYW performed the literature research and data extraction. All authors contributed to the writing of the manuscript and approved the final version.

Corresponding author

Correspondence to Jialiang Zhu.

Ethics declarations

Consent for publication

Not applicable.

Competing interests

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

Additional information

Publisher's Note

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

Supplementary Information

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yan, F., Chen, X., Quan, X. et al. Association between the stress hyperglycemia ratio and 28-day all-cause mortality in critically ill patients with sepsis: a retrospective cohort study and predictive model establishment based on machine learning. Cardiovasc Diabetol 23, 163 (2024). https://doi.org/10.1186/s12933-024-02265-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12933-024-02265-4

Keywords