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

ORIGINAL RESEARCH article

Front. Endocrinol., 30 July 2025

Sec. Systems Endocrinology

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1555082

This article is part of the Research TopicIntegrated Diagnostics and Biomarker Discovery in Endocrinology and Biomedical Sciences: Volume IIView all 6 articles

Association between serum glucose potassium ratio and short- and long-term all-cause mortality in patients with sepsis admitted to the intensive care unit: a retrospective analysis based on the MIMIC-IV database

Jiaqi Lou&#x;Jiaqi Lou1†Ziyi Xiang&#x;Ziyi Xiang2†Xiaoyu Zhu&#x;Xiaoyu Zhu3†Jingyao SongJingyao Song4Shengyong CuiShengyong Cui1Jiliang LiJiliang Li1Guoying JinGuoying Jin1Neng HuangNeng Huang1Youfen Fan*Youfen Fan1*Sida Xu*Sida Xu1*
  • 1Burn Department, Ningbo No. 2 Hospital, Ningbo, Zhejiang, China
  • 2Institute of Pathology, Faculty of Medicine, University of Bonn, Bonn, Germany
  • 3Health Science Center, Ningbo University, Ningbo, Zhejiang, China
  • 4School of Mental Health, Wenzhou Medical University, Wenzhou, Zhejiang, China

Background: The glucose potassium ratio (GPR) is emerging as a biomarker for predicting clinical outcomes in various conditions. However, its value in sepsis patients admitted to the intensive care unit (ICU) remains unclear. Prior studies have shown conflicting results, with some indicating GPR’s potential as an early warning indicator of metabolic decompensation in septic patients, while others found no significant association. The current study addresses these inconsistencies by conducting the first large-scale, systematic validation of GPR in ICU sepsis patients.

Methods: This retrospective cohort study used patient records from the MIMIC-IV database to examine outcomes in sepsis patients. The primary outcomes were hospital and ICU mortality at 30, 60, and 90 days. The correlation between GPR and these outcomes was evaluated using Kaplan-Meier survival analysis, Cox regression models, and restricted cubic spline (RCS) regression analysis. Sensitivity analyses, including Propensity Score Matching (PSM) and E-value Quantification and Subgroup analyses, were performed to assess the robustness of the findings.

Results: The study included 9,108 patients with sepsis. Kaplan-Meier survival curves indicated progressively worsening survival probabilities from Q1 to Q4 for both hospital and ICU mortality across all time points. Cox analysis revealed that patients in the highest GPR quartile (Q4) had a significantly increased risk of mortality compared to those in the lowest quartile (Q1). A nonlinear relationship between GPR and mortality was identified, with a critical threshold at GPR=30. Subgroup analysis showed that the effect size and direction were consistent across different subgroups. Sensitivity analyses, including E-value quantification and propensity score matching, supported the robustness of our findings.

Conclusion: This study demonstrates that higher GPR levels strongly predict increased short- and long-term mortality risk in ICU-admitted sepsis patients. The composite nature of GPR, reflecting both hyperglycemia and hypokalemia, offers incremental prognostic value beyond single metabolic parameter. A critical threshold effect was observed at GPR=30, where risk substantially increased. This consistent association across patient subgroups positions GPR as a promising biomarker for identifying high-risk sepsis patients, warranting prospective validation.

1 Background

Sepsis, a life-threatening organ dysfunction stemming from a dysregulated host response to infection, poses a significant challenge in intensive care units (ICUs) across the globe. Despite advancements in medical care, it remains one of the leading causes of morbidity and mortality, impacting millions annually and resulting in substantial healthcare expenditures (1). The pathophysiology of sepsis is intricate, characterized by a cascade of inflammatory responses that lead to widespread cellular and metabolic abnormalities. Notably, alterations in glucose and potassium homeostasis are critical metabolic disruptions that affect cellular function and systemic homeostasis. Hyperglycemia is frequently observed in septic patients, often attributed to stress-induced hypermetabolism and insulin resistance (2). This metabolic state intensifies oxidative stress and inflammation, further compromising immune function and organ performance. Conversely, potassium imbalances, such as hypokalemia, are common due to factors like increased renal excretion and intracellular shifts caused by insulin therapy or catecholamine surges (3). These electrolyte disturbances can lead to severe complications, including cardiac arrhythmias and muscle weakness (4), thereby worsening the clinical trajectory of sepsis. In recent years, there has been a pressing need to identify reliable prognostic markers to enhance the prediction of sepsis outcomes. While markers like procalcitonin, C-reactive protein, and lactate have shown promise (5), they primarily reflect inflammatory or perfusion-related aspects. Consequently, the identification of novel prognostic biomarkers that capture the complex metabolic imbalances in sepsis remains a crucial research priority.

The serum glucose-potassium ratio (GPR) has emerged as a promising biomarker that reflects the dynamic interplay between glucose and potassium homeostasis, which is often disrupted in various pathological states. Its clinical utility has been recognized in conditions such as diabetic ketoacidosis (6), myocardial infarction (7), and heart failure (8), where it offers a composite view of metabolic derangements that singular parameters fail to capture. In these conditions, an altered GPR has been associated with increased morbidity and mortality, suggesting its potential as a prognostic tool. For instance, studies in myocardial infarction patients have demonstrated a correlation between a high GPR and adverse cardiovascular events, indicating that this biomarker could enhance risk stratification and guide treatment decisions (9, 10). However, research on the application of GPR in sepsis remains limited and has yielded mixed results. Some studies suggest that a high GPR correlates with increased mortality rates and worsened clinical outcomes in sepsis patients, positing that the ratio could serve as an early warning of metabolic decompensation (11, 12). In contrast, a study by Güler et al. (13) found no significant predictive relationship between the glucose-to-potassium ratio and mortality risk in sepsis or septic shock patients admitted to the emergency intensive care unit. These discrepancies may stem from variations in study design, patient populations, or analytical methods. Furthermore, the lack of standardized thresholds and guidelines for interpreting GPR in sepsis complicates its clinical application. Thus, the current understanding of GPR’s relevance to sepsis is limited, underscoring the need for comprehensive evaluations and validation in larger, well-characterized cohorts to establish its potential as a reliable prognostic indicator.

In this context, the MIMIC-IV database serves as a rich repository of de-identified health-related data from thousands of ICU admissions (14), offering a unique opportunity to comprehensively investigate the clinical parameters of sepsis. The database is publicly accessible via the MIMIC-IV platform and contains extensive datasets, including vital signs, laboratory results, and clinical outcomes, which facilitate large-scale retrospective analyses (15). This study aims to explore the association between the serum glucose-potassium ratio and short- and long-term all-cause mortality in ICU-admitted sepsis patients using the MIMIC-IV database. By examining this relationship, we aim to enhance the understanding of metabolic markers in sepsis and potentially identify a novel prognostic indicator that can improve risk stratification and inform treatment strategies for critically ill patients.

2 Methods

2.1 Data source and study design

We conducted a retrospective cohort study utilizing data from the MIMIC-IV database (version 2.2), which is developed and maintained by the Massachusetts Institute of Technology (MIT) and Beth Israel Deaconess Medical Center (BIDMC) (15). This database comprises two in-house systems: a customized hospital-wide electronic health record (EHR) and an ICU-specific clinical information system, encompassing data from 2008 to 2024. One of the authors (JQ L) completed the necessary authentication process and passed the Collaborative Institutional Training Initiative examination (authentication number 60691748) to access the database. Relevant variables were extracted, and patient data were de-identified to ensure privacy. Given the study’s retrospective nature and the anonymized patient data, the Human Research Ethics Committee of Ningbo No.2 Hospital waived the requirement for informed consent.

2.2 Participants

The study encompassed all sepsis patients from the MIMIC-IV v2.2 database. Sepsis was defined according to the Sepsis 3.0 criteria, which were jointly established by the American Society for Critical Care Medicine (SCCM) and the European Society for Critical Care Medicine (ESICM). Patient data were extracted using PostgreSQL. The inclusion criteria were sepsis patients aged 18 and above who were admitted to the ICU for the first time. The following exclusion criteria were applied: (1) patients under 18 years old; (2) ICU stay shorter than 48 hours; (3) patients with recurrent sepsis (only their initial ICU admission was considered); and (4) insufficient data, such as missing records for serum glucose and potassium (Figure 1).

2.3 Research procedures and definitions

Data extraction from MIMIC-IV was performed using Structured Query Language (SQL) via Navicat Premium. The extracted data encompassed a comprehensive set of variables, including patient demographics (age, height, weight, gender, insurance, race, marital status), medical history (hypertension, type 2 diabetes, heart failure, myocardial infarction, malignant tumors, chronic renal failure, cirrhosis, hepatitis, tuberculosis, pneumonia, chronic obstructive pulmonary disease, hyperlipidemia, etc.), and initial laboratory test results (white blood cell count, red blood cell count, neutrophil count, lymphocyte count, platelet count, hemoglobin, mean corpuscular volume, hematocrit, albumin, globulin, total protein, sodium, potassium, calcium, chloride, blood glucose, GPR, anion gap, blood pH, lactate, free calcium, thrombin time, fibrinogen, partial thromboplastin time, international normalized ratio, bilirubin, ALT, AST, urea nitrogen, creatinine, troponin, urine protein, creatine kinase, creatine kinase isoenzyme, N-terminal B-type natriuretic peptide precursor). Special treatments (mechanical ventilation and CRRT), clinical scores (SOFA score, APACHE III score, SAPS II, Oasis score, Charlson Comorbidity Index, SIRS score, GCS score), and clinical outcomes (length of hospital stay, in-hospital mortality, ICU stay, ICU mortality) were also recorded. The 30-day, 60-day, and 90-day mortality rates were calculated. During data cleaning, predictors with more than 30% missing data were excluded. The serum glucose-potassium ratio (GPR) was calculated using the first recorded serum glucose and potassium measurements obtained within 24 hours of ICU admission, based on the formula: GPR = serum glucose (mg/dL)/serum potassium (mmol/L) (16).

2.4 Outcomes and measures

The primary outcomes of this study were hospital mortality and ICU mortality at 30-day, 60-day and 90-day.

2.5 Statistical analysis

Continuous variables were presented as mean ± standard deviation or median (interquartile range), while categorical variables were reported as frequency and percentage. Data conforming to a normal distribution were analyzed through the t-test or analysis of variance (ANOVA).

For data not following a normal distribution, the Mann-Whitney U test or Kruskal-Wallis test was employed (17, 18). Kaplan-Meier survival analysis was utilized to assess the incidence of endpoint events across different GPR levels, with differences evaluated through the log-rank test. Kaplan-Meier curves offer a visual comparison of survival differences between groups or conditions and do not require prior assumptions about data distribution (19), so it was relatively flexible in use.

The Cox proportional hazards model was utilized to calculate the hazard ratio (HR) and 95% confidence interval (CI) between the GPR and the endpoint. This model, taking survival outcome and survival time as dependent variables, enabled simultaneous analysis of multiple factors affecting survival and analysis of the data with censored survival time, and did not necessitate the estimation of the survival distribution type (20). The GPR was analyzed both as a continuous variable and by quartiles. Cox proportional hazards models were constructed in three sequential tiers: Model 1 (univariate); Model 2 (adjusted for demographics: age, sex, height, weight, insurance, marital status, race); Model 3 (further adjusted for laboratory/clinical covariates: WBC, RBC, RDW, albumin, chloride, ALT, AST, comorbidities [hypertension, diabetes, heart failure, etc.], treatments [CRRT], and severity scores [SOFA, SAPS II, etc.]).

Restricted cubic splines (RCS) used 4 knots placed at the 5th, 35th, 65th, and 95th percentiles. Nonlinearity was tested via the significance of the second spline term. The GPR was incorporated as either a continuous or ordered variable into the model, with the first quartile of the GPR serving as the reference group. The quartile level was used for the calculation of the P-value of the trend. RCS was a non-parametric flexible fitting method that models survival curves by transforming survival times into piecewise functions at individual nodes (21) and can accommodate various types of survival time distributions without excessive assumptions.

Subgroup analyses (22) were conducted to explore potential differences across various subgroups based on age (≤ 70 years and > 70 years), sex, BMI (<27.4 kg/m², 27.4-31.2 kg/m², ≥31.2 kg/m²), age, sex, BMI, hypertension, type 2 diabetes, heart failure, CKD, stroke, AKI, CRRT, and mechanical ventilation, to evaluate the consistency of the GPR’s prognostic value for the primary outcomes. Cox models were also adopted in subgroup analyses to adjust for all variables in the patient’s baseline information.

Sensitivity analyses included: (1) E-values to quantify unmeasured confounding. To evaluate the potential impact of unmeasured confounding on the association between GPR and mortality outcomes, we also calculated E-values using the formula: E-value = RR + √(RR*(RR-1)), where RR is the hazard ratio (HR) derived from Cox regression models. This approach helped assess the robustness of our findings against unmeasured confounding.; (2) Propensity score matching (PSM) (22)). To further assess the robustness of our findings and address potential confounding factors, we conducted a propensity score matching (PSM) analysis. This method helps to reduce selection bias by balancing the distribution of observed covariates between the exposure groups (high GPR group and low GPR group). We defined the high GPR group as patients with GPR above the mean value and the low GPR group as patients with GPR below the mean value. The nearest-neighbor matching method was used to match each patient in the high GPR group with two patients in the low GPR group (1:2 matching), with a caliper width of 0.2 standard errors. Categorical variables were converted into dummy variables for the analysis. For example, marital status was categorized as divorced (1) versus others (0), married (1) versus others (0), and so on. The matching process aimed to create a more balanced comparison group by controlling for key variables such as age, sex, and SOFA score, which are known to influence outcomes in sepsis patients. In the PSM analysis, the balance assessment focuses on comparing the distribution of covariates between the treatment (high GPR) and control (low GPR) groups. The goal of balance assessment is to ensure that these groups are comparable in terms of key covariates, which is crucial for reducing selection bias and enhancing the validity of the study. It is important to note that different outcome variables do not influence the results of balance assessment, as the assessment is solely concerned with the distribution of covariates. Thus, our selection of covariates for balance assessment is based on their potential confounding effects on the relationship between GPR and hospital mortality. This approach ensures that the matched groups are balanced in terms of key covariates, providing a solid foundation for the subsequent analysis of the association between GPR and hospital mortality. After matching, we repeated the Cox regression analysis to assess the association between GPR and hospital mortality. The primary outcome was the all-cause mortality at 30-day, 60-day, and 90-day. The balance of covariates before and after matching was assessed using standardized bias and t-tests. A standardized bias of less than 10% and a p-value greater than 0.05 for the t-tests indicated successful matching. Additionally, a common support test was performed to ensure that the propensity scores of the treatment and control groups overlapped sufficiently, minimizing potential biases.

Data processing and analysis were carried out via R version 4.3.0, along with Zstats v1.0 (www.zstats.net), with statistical significance set at P<0.05 for two-tailed tests. The primary analyses utilized the following packages: Data management and transformation were conducted using dplyr and tidyr. Survival analyses including Kaplan-Meier curves, log-rank tests (via survdiff()), and univariate/multivariate Cox proportional hazards regression (via coxph()) were implemented with the survival package. Nonlinear relationships were assessed through RCS using the rms package, with knots placement and trend significance testing performed via rcs() and anova() functions. Subgroup analyses were streamlined using purrr for iterative modeling and broom for result standardization. E-value analysis was also conducted in R, utilizing packages survival for Cox regression and EValue for E-value calculation. The PSM was performed using the MatchIt package in R, which allows for various matching algorithms, including nearest neighbor, optimal, and full matching. Visualizations were generated with ggplot2 and enhanced using survminer for survival plots. For missing values in the data, the multiple imputation method of the random forest was used to interpolate the missing value data (through the R package “mice”). Features with missing values exceeding 50% were removed before interpolation.

3 Results

Among the adult patients in the MIMIC-IV database, a total of 22,517 subjects met the eligibility criteria. From the database, 148 prognostic factors were initially extracted. Following data cleaning, 80 predictors with over 30% missing data were excluded. In the end, 68 forecast factors were included in the model.

3.1 Characteristics of included patients

A total of 9,108 people were included in the study, of which 2,272 (24.95%) were in GPR quantile 1 (Q1) group (GPR ≤ 6.67), 2,282 (25.05%) people were in quantile 2 group (6.67 < GPR≤ 25.71), 2,277 people were in quantile 3 group (25.71 < GPR ≤ 40.81), and 2,277 people in quantile 4 group, accounting for 25.00% (GPR > 40.81). IQR is 15.09. The average GPR of all patients was 35.55 ± 16.49. Upon stratification into these four categories, the distribution of each variable across the groups was analyzed. All baseline data are presented in Table 1 and Supplementary Table 1.

Table 1
www.frontiersin.org

Table 1. Summary of characteristics that are statistically different of the study population.

Patients in Q1 were older and had lower body weight than those in the other groups, and there were differences in sex and marital status among the four groups. WBC, RBC, platelet, hemoglobin, hematocrit, albumin, sodium, glucose, anion gap, fibrinogen, PPT, ALT, CK, CKMB were also higher in Q4 group than in Q1 group, but RDW, potassium, hematocrit and bilirubin were lower than Q1 group. There was no significant difference in height, insurance, languages, CRRT days, ventilation, COPD, heart failure, hepatitis, tuberculosis and pneumonia (P>0.05) (Table 1).

Figure 1
Flowchart illustrating patient selection from the MIMIC-IV V2.2 database. Out of 50,920 patients, 22,517 adults diagnosed with sepsis were identified. Inclusion criteria filter for age 18-150, early hospital admission, and sepsis-3 diagnosis. Exclusion criteria remove cases with missing hospital time, ICU stay under 24 hours, and missing serum data. The resulting analysis cohort comprises 9,109 patients, divided into four groups (Q1: 2,272; Q2: 2,282; Q3: 2,277; Q4: 2,277).

Figure 1. Selection of the study population from the MIMIC-IV database.

3.2 Kaplan-Meier survival curve analysis

Kaplan-Meier curves (Figure 2) demonstrated worsening survival probabilities from Q1 to Q4 for both hospital and ICU mortality at 30-day, 60-day, and 90-day intervals (log-rank test, all P < 0.001). Specifically, a total of 9,108 people were included in the study, of which 2,272 (24.95%) were in GPR quantile 1 (Q1) group (GPR ≤ 6.67), 2,282 (25.05%) people were in quantile 2 group (6.67 < GPR≤ 25.71), 2,277 people were in quantile 3 group (25.71 < GPR ≤ 40.81), and 2,277 people in quantile 4 group, accounting for 25.00% (GPR > 40.81).

Figure 2
Six Kaplan-Meier survival curves labeled A to F show survival probabilities over time for four GPR quartile groups, distinguished by color. Each chart includes a log-rank p-value below 0.001, indicating statistical significance. Survival probabilities decrease over time across all graphs, with survivability varying among quartiles. Tables below each graph display the number at risk for specific days.

Figure 2. Kaplan-Meier survival curve of cumulative survival rate during hospitalization and ICU for groups. (A): Kaplan-Meier survival curve of cumulative survival rate during hospitalization for groups at 30-day. (B): Kaplan-Meier survival curve of cumulative survival rate during ICU for groups at 30-day. (C) Kaplan-Meier survival curve of cumulative survival rate during hospitalization for groups at 60-day. (D) Kaplan-Meier survival curve of cumulative survival rate during ICU for groups at 60-day. (E) Kaplan-Meier survival curve of cumulative survival rate during hospitalization for groups at 90-day. (F) Kaplan-Meier survival curve of cumulative survival rate during ICU for groups at 90-day. X-axis: Time (Days); Y-axis: Survival Probability. Log-rank test, all P < 0.001. Q1: dark blue; Q2: red; Q3: green; Q4: light blue.

3.3 Cox regression models for all-cause mortality (in hospital and ICU)

In the Cox regression analysis, a higher GPR was positively correlated with increased mortality rates in both the ICU and hospital settings among critically ill patients with sepsis. When the GPR was analyzed as a continuous variable, it was independently associated with a higher risk of hospital mortality both at 30-day, 60-day and 90-day (All P < 0.05). Patients in Q4 had a 15–20% higher risk of mortality compared to Q1 across all time points. At 60-day, when categorized into quartiles, Model 1 revealed that the risk of hospital mortality for Q4 were 19% higher than for Q1 (HR 1.19 [95% CI 1.08 to 1.31], P < 0.001), Model 2 revealed that the risk of hospital mortality for Q4 were 18% higher than for Q1 (HR 1.18 [95% CI 1.03 to 1.35], P < 0.001). At 90-day, when categorized into quartiles, Model 1 revealed that the risk of hospital mortality for Q4 were 20% higher than for Q1 (HR 1.20 [95% CI 1.09 to 1.32], P < 0.001), Model 2 revealed that the risk of hospital mortality for Q4 were 15% higher than for Q1 (HR 1.15 [95% CI 1.01 to 1.32], P = 0.037). The differences in Model 3 results compared to Models 1 and 2 are likely due to the additional adjustment for confounding variables such as WBC, RBC, RDW, albumin, chloride, ALT, etc.

For ICU mortality, the GPR, when used as a continuous variable, was significantly associated with an elevated risk of ICU death in Models 1, 2 and 3 (All P < 0.001). Furthermore, when the GPR was categorized into quartiles, at 30-day, Model 1 demonstrated that the risk of ICU mortality for Q4 was 1.13 times that of Q1 (HR 1.13 [95% CI 1.01 to 1.26], P < 0.001). At 60-day, Model 1 demonstrated that the risk of ICU mortality for Q4 was 1.23 times that of Q1 (HR 1.23 [95% CI 1.10 to 1.37], P < 0.001), Model 2 demonstrated that the risk of ICU mortality for Q4 was 1.21 times that of Q1 (HR 1.04 [95% CI 1.01 to 1.41], P = 0.015) (Table 2).

Table 2
www.frontiersin.org

Table 2. The association between GPR groups and in-hospital and ICU mortality at 30-day, 60-day and 90-day.

3.4 RCS regression models for all-cause mortality (in hospital and ICU)

We subsequently employed the RCS regression models to elucidate the risk and discovered a nonlinear association between the GPR and mortality. Figures 3 and 4 illustrate the results of the univariate and multivariate analyses regarding the relationship between the GPR and in-hospital, In-ICU mortality in three time points, respectively.

Figure 3
Six graphs labeled A to F depict hazard ratios with 95% confidence intervals versus GPR values. Each graph shows a curve with a red line and shaded pink area. Graphs A, C, and E display upward trends, while B, D, and F show initial decreases followed by stabilization or slight increases. Each graph provides p-values indicating overall and nonlinear significance.

Figure 3. RCS regression for GPR and in-hospital mortality. (A) Univariate analysis at 30-day (P for overall effect <0.001; P for nonlinearity <0.001). (B) Multivariate analysis at 30-day (P for overall effect <0.001; P for nonlinearity <0.001). (C) Univariate analysis at 60-day (P for overall effect <0.001; P for nonlinearity <0.001). (D) Multivariate analysis at 60-day (P for overall effect 0.007; P for nonlinearity 0.004). (E) Univariate analysis at 90-day (P for overall effect <0.001; P for nonlinearity <0.001). (F) Multivariate analysis at 90-day (P for overall effect 0.014; P for nonlinearity 0.010).

Figure 4
Six graphs labeled A to F illustrate the hazard ratios (with 95% confidence intervals) plotted against GPR values. The red lines represent hazard ratios, with surrounding shaded areas indicating confidence intervals. Significant p-values for nonlinearity are shown, suggesting differing patterns for each graph. Graphs A, C, and E display increasing hazard ratios at higher GPR, while B, D, and F show decreasing patterns.

Figure 4. RCS regression for GPR and mortality during ICU admission. (A) Univariate analysis at 30-day (P for overall effect <0.001; P for nonlinearity <0.001). (B) Multivariate analysis at 30-day (P for overall effect <0.001; P for nonlinearity <0.001). (C) Univariate analysis at 60-day (P for overall effect <0.001; P for nonlinearity <0.001). (D) Multivariate analysis at 60-day (P for overall effect 0.007; P for nonlinearity 0.004). (E) Univariate analysis at 90-day (P for overall effect <0.001; P for nonlinearity <0.001). (F) Multivariate analysis at 90-day (P for overall effect 0.014; P for nonlinearity 0.010).

Figures 3A, B present the findings of the univariate and multivariate analyses concerning the association between the GPR and hospital mortality on 30-day, respectively. Before adjusting for 30-day in-hospital mortality, the p-value for the overall effect was < 0.001, and the p-value for the nonlinear effect was also < 0.001. Following adjustment, all p-values were less than 0.05. Similarly, nonlinear associations were observed for 60-day (Figures 3C, D) and 90-day (Figures 3E, F) in-hospital mortality, both before and after adjustment for relevant factors.

For ICU mortality, on 30-day mortality (Figures 4A, B), the unadjusted p value for the overall effect was less than 0.001, the p value for the nonlinear effect was less than 0.001, and all adjusted p values were greater than 0.05. The unadjusted p value was less than 0.001 for the overall effect and less than 0.001 for the nonlinear effect on 60-day mortality (Figures 4C, D). The adjusted p value was 0.007 for the overall effect and 0.004 for the nonlinear effect. Finally, on 90-day mortality (Figures 4E, F), the unadjusted p value was less than 0.001 for the overall effect and less than 0.001 for the nonlinear effect. After adjustment, the p value of overall effect was 0.014, and the p value of nonlinear effect was 0.01. Figures 3 and 4 demonstrate that the inflection point in both multifactorial models is about 30.

3.5 Subgroup analysis

In subgroup analyses, the directionality of the effect estimates in subgroups was consistent with the overall outcomes. Subgroup analyses were stratified by age, sex, BMI, hypertension, type 2 diabetes, heart failure, CKD, stroke, AKI, CRRT, and mechanical ventilation.

The directional trends in the effect estimates for in-hospital mortality (Figure 5A) in almost subgroups were consistent with the overall outcomes before adjustment for covariates. Similarly, almost all subgroups were consistent with the overall outcome of ICU mortality (Figure 5B). In addition, there was an interaction between mechanical ventilation subgroup parameters (P < 0.01 for interaction). After adjustment for covariates, the directionality of the effect estimates in in-hospital and ICU mortality was consistent with the overall outcome in almost all subgroups except AKI and the subgroups with CRRT and no mechanical ventilation. There was no interaction between GPR and age, gender, BMI, hypertension, type 2 diabetes, heart failure, CKD, shock and mechanical ventilation (all P for interaction >0.05).

Figure 5
Two tables labeled A and B compare the hazard ratios and confidence intervals for various subgroups, including gender, age, BMI, and medical conditions like diabetes and stroke. Each table presents data, hazard ratios, and p-values for interaction. Both tables include markers indicating statistical significance.

Figure 5. Forest plots for subgroup analyses of the association between GPR and mortality. (A) Subgroup analysis of the association between GPR and in-hospital mortality after covariate adjustment. (B) Subgroup analysis of the association between GPR and ICU mortality after covariate adjustment. For both plots, hazard ratios (HRs) and 95% confidence intervals (CIs) are shown. The analysis includes subgroups based on age (≤70 years and >70 years), sex, BMI (<27.4 kg/m², 27.4–31.2 kg/m², ≥31.2 kg/m²), hypertension, type 2 diabetes, heart failure, CKD, stroke, AKI, CRRT, and mechanical ventilation. The P value for interaction is provided for each subgroup analysis.

3.6 Sensitivity analyses

The E-values for the association between GPR and mortality outcomes at different time points are as follows: For ICU mortality, the E-values are 1.60 (30-day, HR=1.13), 1.79 (60-day, HR=1.23), and 1.85 (90-day, HR=1.26). For in-hospital mortality, the E-values are 1.11 (30-day, HR=1.05), 1.66 (60-day, HR=1.19), and 1.68 (90-day, HR=1.20). An E-value of 1.60 for 30-day ICU mortality implies that an unmeasured confounder would need to be associated with both the exposure and outcome by at least 1.60-fold to fully explain the observed association. Similarly, higher E-values for other time points indicate the minimum association strength required for potential unmeasured confounders to explain the observed results.

The common support test results confirmed that the propensity scores of the high GPR and low GPR groups had sufficient overlap. The kernel density plots showed that the density lines for the two groups were closely aligned both before and after matching, indicating a large common support region. The histograms further demonstrated that most observations were within the common support range, ensuring minimal sample loss during the matching process. This confirmed the reliability of the matching process and the comparability of the matched groups. The balance assessment figures demonstrate that after PSM, the bias for all covariates was reduced to below 10%, and the t-tests showed no significant differences between the groups (p > 0.05). This indicates that the matching process successfully balanced the covariates between the high and low GPR groups. The kernel density and histogram figures show that the propensity scores of the two groups had sufficient overlap both before and after matching. After matching, the density lines and histogram bars for the two groups were closely aligned, indicating a large common support region and minimal loss of samples. This ensures that the matched groups are comparable and the results are reliable. These visualizations provide additional evidence of the effectiveness of the PSM method in reducing bias and enhancing the comparability of the groups, thereby strengthening the validity of the study findings. (Supplementary Table 2, Figures 6, 7).

Figure 6
Chart A shows density plots for propensity scores before matching, with a clear separation between treat and control groups. Chart B displays these scores after matching, showing improved overlap between groups. Chart C presents a histogram of propensity scores, categorizing treated and untreated groups into on support and off support, using different colors for each category.

Figure 6. Propensity score matching and common support assessment regarding in-hospital mortality. (A) Kernel Density Estimation Before Matching: Displays the kernel density estimates of propensity scores for the treatment group (blue line) and control group (red line) prior to matching. The overlapping regions between the two curves indicate the initial common support area. Before matching, the density curves show some overlap, but there are also areas where the propensity scores of the treatment and control groups do not align closely, suggesting a limited common support region. (B) Kernel Density Estimation After Matching: Shows the kernel density estimates of propensity scores for the treatment group (blue line) and control group (red line) following matching. After matching, the density curves of the two groups are closely aligned across a wider range of propensity scores. This close alignment demonstrates an expanded common support region, indicating that the matching process has effectively balanced the distribution of propensity scores between the treatment and control groups. (C) Histogram of Common Support: Presents a histogram displaying the distribution of propensity scores for both the treatment and control groups. The green bars represent the treated observations within the common support range, the red bars represent the untreated observations within the common support range, the blue bar represents untreated observations outside the support, and the orange bar represents treated observations outside the support. The majority of observations fall within the common support range (indicated by the green and red bars), which means that only a minimal number of samples were excluded during the matching process. This ensures that the matched groups are highly comparable and reduces the potential for bias in the subsequent analysis.

Figure 7
Panel A shows kernel density plots of propensity scores for treatment and control groups before matching, with control having a higher peak. Panel B presents similar plots after matching, showing closer alignment between groups. Panel C is a histogram of propensity scores divided into treated and untreated groups, indicating balance with bars for treated on support and untreated off support.

Figure 7. Propensity score matching and common support assessment regarding in-ICU mortality. (A) Kernel Density Estimation Before Matching: Displays the kernel density estimates of propensity scores for the treatment group (blue line) and control group (red line) prior to matching. The overlapping regions between the two curves indicate the initial common support area. Before matching, the density curves show some overlap, but there are also areas where the propensity scores of the treatment and control groups do not align closely, suggesting a limited common support region. (B) Kernel Density Estimation After Matching: Shows the kernel density estimates of propensity scores for the treatment group (blue line) and control group (red line) following matching. After matching, the density curves of the two groups are closely aligned across a wider range of propensity scores. This close alignment demonstrates an expanded common support region, indicating that the matching process has effectively balanced the distribution of propensity scores between the treatment and control groups. (C) Histogram of Common Support: Presents a histogram displaying the distribution of propensity scores for both the treatment and control groups. The green bars represent the treated observations within the common support range, the red bars represent the untreated observations within the common support range, the blue bar represents untreated observations outside the support, and the orange bar represents treated observations outside the support. The majority of observations fall within the common support range (indicated by the green and red bars), which means that only a minimal number of samples were excluded during the matching process. This ensures that the matched groups are highly comparable and reduces the potential for bias in the subsequent analysis.

4 Discussion

This study examines the association between GPR and short- and long-term all-cause mortality in ICU-admitted sepsis patients using the MIMIC-IV database. With a large sample and extensive confounder adjustment, the results show a significant link between higher GPR and increased mortality risk in both hospital and ICU settings over 90 days. The nonlinear relationship identified by restricted cubic spline regression, with an inflection point at GPR 30, adds depth to GPR’s prognostic potential. Our study is the first large-scale validation of GPR in ICU sepsis patients, addressing inconsistencies in prior literature (13, 23). The composite GPR captures synergistic metabolic dysregulation (hyperglycemia + hypokalemia), explaining its incremental prognostic value over isolated markers. The U-shaped association—lower risk in Q2/Q3 vs. Q1—may reflect protective effects of moderate metabolic stress, whereas extremes (Q1: hypokalemia; Q4: severe dysregulation) drive mortality. The former is likely to exacerbate cardiac instability, while the latter’s extreme dysregulation overrides compensatory mechanisms. This aligns with the RCS-identified threshold (GPR=30), beyond which mortality risk escalates sharply. Sensitivity analyses including E-value quantification and propensity score matching further reinforced the robustness of our primary findings. The E-values (1.60–1.85 for ICU mortality) indicate that unmeasured confounders would need strong associations to nullify our results, while PSM confirmed the mortality gradient across quartiles in matched cohorts. These findings underscore GPR’s utility as a prognostic indicator in critically ill septic patients.

This study underscores that GPR, when evaluated both as a continuous variable and within categorized quartiles, stands out as a predictive marker for mortality in septic patients requiring intensive care. In particular, patients belonging to the highest GPR quartile (Q4) consistently demonstrated notably higher mortality rates compared to those in the lowest quartile (Q1) across all measured intervals (30, 60, and 90 days) and settings (hospital and ICU), as shown by Hazard Ratios (HRs) that reflected increased risk. These findings highlight the GPR’s potential as an independent prognostic indicator beyond traditional physiological and biochemical markers often used in ICU settings. While our study offers novel insights into the prognostic role of GPR in sepsis, it builds upon a modest body of prior research investigating GPR in various medical contexts. In non-septic conditions, such as myocardial infarction (7) and heart failure (6), elevated GPRs have also demonstrated correlations with increased morbidity and mortality, signifying its broad potential as a marker of metabolic imbalance. In ischemic stroke patients, a study (24) found that GPR was positively correlated with 30-day mortality, and the relationship between them was linear. In a multicenter retrospective cohort study (25), baseline GPR serum was found to be an independent predictor of all-cause mortality within 12 months in patients with acute and subacute ischemic stroke, and the study by Zhang et al. (26) also reached a similar conclusion. Chen et al. (27) found that high GPR was an independent risk factor for in-hospital mortality in patients with Acute type A aortic dissection (ATAAD). Serum GPR was observed in 146 patients. In cases of severe traumatic brain injury is substantially associated with trauma severity and 30-day mortality (28), and a similar association has been observed in patients with traumatic brain injury undergoing emergency craniotomy (29). Similarly, another study (30) observed a significant relationship between serum GPR and admission injury severity and the 6-month prognosis acute traumatic Spinal cord injurypatients. A high GFR correlated with Hunt and Kosnik grade and was also observed in patients with aneurysmal subarachnoid hemorrhage at admission Glasgow Outcome Scale score at discharge (31, 32). The predictive value between GPR and adverse clinical outcomes was also preliminarily verified in patients with acute intracerebral hemorrhage. In a retrospective study (33), it was observed that the predictive efficacy of GRF for the diagnosis of massive pulmonary embolism and non-massive pulmonary embolism in ICU patients was higher than that of D-dimer. However, another study based on ICU patients (34) found that the mortality of patients with isolated blunt abdominal trauma was highly correlated with GFR, and the sensitivity and specificity of GRF were both higher than 70%. Such studies provide a contextual backdrop where the dysregulation of glucose and potassium levels has been similarly implicated in adverse outcomes, suggesting a possible cross-pathophysiological utility of the GPR. However, existing literature on GPR specifically within sepsis is relatively scant, and the findings have been inconclusive due to significant methodological variances and population differences.

The GPR in sepsis reflects intricate metabolic dysregulations that accompany the systemic inflammatory response characteristic of this condition. Understanding the potential pathological mechanisms that lead to changes in both glucose and potassium levels can provide valuable insights into the prognostic value and clinical significance of GPR in sepsis. In sepsis, hyperglycemia is a frequent occurrence due to a combination of increased hepatic glucose production and impaired peripheral glucose utilization. Stress-induced hormonal responses (35), including the release of cortisol, catecholamines, glucagon, and pro-inflammatory cytokines (36), like tumor necrosis factor-alpha and interleukins, stimulate hepatic gluconeogenesis and glycogenolysis. This hypermetabolic state is compounded by insulin resistance, which limits glucose uptake by peripheral tissues, further elevating blood glucose levels (37). The pathological mechanism of hyperglycemia in sepsis can exacerbate the disease’s course through a variety of pathways. Elevated glucose levels contribute to oxidative stress by generating advanced glycation end products (AGEs) (38), which promote inflammation and tissue injury. Hyperglycemia also impairs neutrophil function (39), thereby weakening the host immune response and increasing susceptibility to infections. Furthermore, it is associated with endothelial dysfunction (40, 41), facilitating microvascular thrombosis and impaired tissue perfusion, which can deteriorate organ function. Clinically, the presence of hyperglycemia in sepsis patients has been linked to worse outcomes, including increased mortality rates, prolonged ICU stay, and higher incidences of multi-organ failure (42). This underlines the importance of close glycemic control in critical care settings, although the potential benefits must be weighed against the risks of hypoglycemia. Potassium imbalances, notably hypokalemia, are also common in sepsis and can stem from several factors. These include intracellular shifts of potassium driven by insulin administration (43) (used therapeutically to control hyperglycemia), beta-adrenergic stimulation, and metabolic alkalosis, as well as increased renal losses due to activation of the renin-angiotensin-aldosterone system and nephrotoxic effects of medications or the sepsis itself. Alternatively, hyperkalemia can occur, particularly in cases of acute kidney injury or significant cellular lysis (44). The clinical consequences of potassium imbalances are profound. Hypokalemia may lead to arrhythmias, muscle weakness, and respiratory failure, while hyperkalemia can precipitate potentially fatal cardiac arrhythmias (45). Potassium levels are critical for the function of cells, particularly in excitable tissues such as nerves and muscles, including the heart, implicating disturbances in significant morbidity in septic patients (46).

The ratio of serum glucose to potassium, or GPR, synthesizes the metabolic derangements of these two crucial solutes into a single metric. While each component on its own provides insight into specific pathophysiological processes, the GPR captures the overarching metabolic stress within the body (47). A high GPR may indicate a metabolic milieu marked by severe insulin resistance, profound stress response, and possibly inadequate compensatory mechanisms for electrolyte maintenance (48). This composite biomarker might therefore reflect a higher severity of systemic physiological derangement, correlating with worse clinical outcomes. Integrating glucose and potassium levels into a single ratio could afford a fuller picture of the metabolic state in sepsis compared to evaluating each element in isolation. In clinical practice, monitoring the GPR in sepsis patients could potentially aid in identifying patients at higher risk of adverse outcomes, offering opportunities for early intervention and more tailored therapeutic strategies. However, understanding the precise interplay and optimizing clinical use of GPR necessitate further research exploring the dynamic interrelations between glucose and potassium metabolisms in the progression of sepsis.

The E-values calculated for the association between GPR and mortality outcomes provide additional insight into the robustness of our findings against unmeasured confounding. For instance, an E-value of 1.60 for 30-day ICU mortality implies that an unmeasured confounder would need to be associated with both the exposure and outcome by at least 1.60-fold to fully explain the observed association. Similarly, the propensity score matching (PSM) analysis confirmed the consistency of our findings, further strengthening the validity of the observed association between GPR and mortality in sepsis patients.

This study’s contribution to the field is highlighted by its significant dataset derived from the MIMIC-IV database, encompassing a variety of demographic and clinical variables not previously analyzed in this combination. By confirming the prognostic relevance of GPR across a diverse ICU population, our findings suggest this biomarker could play a critical role in advancing sepsis management protocols, potentially guiding therapeutic decisions to mitigate mortality risks more effectively. Future research should focus on prospective validation of GPR thresholds and exploration of GPR dynamics over the course of sepsis to better understand its prognostic implications. By identifying patients at high risk of poor outcomes early in their treatment course, clinicians could tailor more aggressive monitoring and intervention strategies, which might include tighter glucose control, more frequent electrolyte assessments, or enhanced cardiovascular monitoring. Such an approach could lead to better resource allocation in high-intensity care environments and possibly improve patient outcomes by preemptively managing predicted complications.

This study also has limitations. The MIMIC-IV database consists largely of data from patients at a single tertiary care center, potentially limiting the generalizability of findings to other settings with different demographics, socioeconomic backgrounds, or healthcare systems (49). This can result in a population that is not fully representative of broader, more diverse sepsis populations worldwide. The demographic composition within the database may not sufficiently capture the variability across different ethnic and racial groups, which can affect disease presentation and responses to treatment, potentially skewing results and interpretations. Although the study includes adjustments for factors such as age and comorbidities, the inherent diversity in these variables may not be fully comparable across different demographic groups (50), implicating variations in baseline mortality risk that might confound the association between GPR and outcomes. In addition, as a retrospective study, it is subject to inherent biases such as selection bias and information bias (51). Decisions regarding data extraction and the variables included can introduce unintended biases that might impact the overall interpretation of findings. Despite efforts to adjust for numerous confounders, it is possible that not all relevant factors were considered or measured accurately, leading to residual confounding. Factors such as medication usage, nutritional status, or patient management differences might not be fully accounted for. The timing of GPR measurement relative to the onset of sepsis or the clinical course has not been standardized (52), potentially impacting its reliability as a consistent prognostic tool. The variation in when glucose and potassium levels are recorded can introduce discrepancies in how the GPR is calculated and interpreted. What’s more, a notable limitation is the potential for missing data, as not all patients may have fully recorded laboratory measurements or clinical outcomes. The study relied on multiple imputation methods to address missing data, which may introduce bias if assumptions about missingness are incorrect (53). The dataset may lack comprehensive longitudinal data necessary to explore causal relationships over time, limiting insights into how changes in GPR might reflect disease progression or response to interventions. Certain clinical variables crucial for understanding individual patient conditions, such as specific dietary intake, detailed medication histories, and underlying genetic predispositions (54, 55), may not be captured in the database, affecting the depth of analysis. Given the nature of the database as an aggregation of EMR from clinical practice, the quality and precision of recorded data can be variable. This variability may affect the accuracy of the input data, especially laboratory measurements, and the resulting analysis (56, 57). Notably, the lack of data on treatment interventions such as insulin therapy and fluid resuscitation represent a key limitation, as these factors can significantly influence patient outcomes and may confound the relationship between GPR and mortality (58, 59).

As a path forward, prospective studies evaluating GPR longitudinally across different stages of sepsis, and within broader and more varied populations, could validate our findings. Investigations might also focus on optimal intervention strategies for patients identified as high-risk by their GPR, possibly examining the impact of targeted therapies aimed at normalizing glucose and potassium homeostasis (60). Furthermore, establishing standardized GPR thresholds and developing clinical guidelines for their use could facilitate more widespread integration of GPR into ICU protocols. Limitations of our study, such as its retrospective nature and reliance on a single database, should also be addressed in future studies to enhance generalizability (61). Additionally, detailed longitudinal data collection could enable a better understanding of the causal pathways potentially involved in the links between GPR and sepsis outcomes.

5 Conclusion

In summary, the serum glucose-potassium ratio emerges from our investigation as a promising biomarker of mortality risk in sepsis, warranting further exploration and validation in future clinical research endeavors. By enhancing our understanding and application of GPR, healthcare practitioners may improve prognostic accuracy and patient outcomes in the challenging realm of sepsis management.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics statement

This study adhered to the principles outlined in the Declaration of Helsinki. Furthermore, it is important to note that the MIMIC-IV database received approval from the Massachusetts Institute of Technology and the Beth Israel Deaconess Medical Center. The copies of the datasets used in this study are available from the MIMIC database. Due to the retrospective nature of the study, the need of informed consent was waived by the Human Research Ethics Committee of Ningbo No.2 Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.

Author contributions

JQL: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. ZX: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. XZ: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. JS: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. SC: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. JLL: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. GJ: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. NH: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. YF: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. SX: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the HwaMei Research Foundation of Ningbo No.2 Hospital (Grant No.2022HMKY48 and No.2023HMZD07), Medical Scientific Research Foundation of Zhejiang Province (Grant No.2021KY1004, No.2023RC081 and No.2025KY1395), the Project of NINGBO Leading Medical & Health Discipline (Project No.2022-F17), the Ningbo Top Medical and Health Research Program (No.2023030615) the Zhejiang Clinovation Pride (CXTD202502004), Research and development of efficient hemostatic materials (2024001), the Zhu Xiu Shan Talent Project of Ningbo No.2 Hospital (No.2023HMYQ25) and Ningbo Health Youth Technical Backbone Talent Development Program (No.2024RC-QN-02). Funders played no role in the study design, execution, or manuscript writing.

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.

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/fendo.2025.1555082/full#supplementary-material

Abbreviations

ICU, Intensive Care Units; TG, Triglyceride; HDL-C, High-Density Lipoprotein Cholesterol; BMI, Body Mass Index; MIMIC-IV, Medical Information Mart for Intensive Care IV; SOFA, Sequential Organ Failure Assessment; RCS, Restricted Cubic Splines; COPD, Chronic Obstructive Pulmonary Disease; CKD, Chronic Kidney Disease; HR, Hazard Risk; CI, Confidence Interval; HF, Heart Failure; HT, Hypertension; DM2, Diabetes Mellitus type 2; OASIS, Oxford Acute Severity of Illness Score; SAPS II, Simplified Acute Physiology Score II; WBC, White Blood Cell; Rbc, Red Blood Cell; RDW, Red blood cell Distribution Width. AKI, Acute Kidney Injury; CRRT, Continuous Renal Replacement Therapy; HLP, Hyperlipidemia; TB, Tuberculosis; ARF, acute renal failure; MT, Malignant Tumor; MI, Myocardial Infarction.

References

1. Qian H, Shang W, Zhang S, Pan X, Huang S, Li H, et al. Trends and predictions of maternal sepsis and other maternal infections among women of childbearing age: a systematic analysis for the global burden of disease study 2019. Front Public Health. (2024) 12:1428271. doi: 10.3389/fpubh.2024.1428271

PubMed Abstract | Crossref Full Text | Google Scholar

2. Cheng S, Li Y, Sun X, Liu Z, Guo L, Wu J, et al. The impact of glucose metabolism on inflammatory processes in sepsis-induced acute lung injury. Front Immunol. (2024) 15:1508985. doi: 10.3389/fimmu.2024.1508985

PubMed Abstract | Crossref Full Text | Google Scholar

3. Li J, Nie M, Lu Z, Wang Y, and Shen X. Efficacy of compound sodium acetate Ringer’s solution in early fluid resuscitation for children with septic shock: a preliminary retrospective cohort study. BMC Pediatr. (2024) 24:708. doi: 10.1186/s12887-024-05184-1

PubMed Abstract | Crossref Full Text | Google Scholar

4. Lou J, Xiang Z, Zhu X, Fan Y, Song J, Cui S, et al. A retrospective study utilized MIMIC-IV database to explore the potential association between triglyceride-glucose index and mortality in critically ill patients with sepsis. Sci Rep. (2024) 14:24081. doi: 10.1038/s41598-024-75050-8

PubMed Abstract | Crossref Full Text | Google Scholar

5. Liu S, Yang T, Jiang Q, Zhang L, Shi X, Liu X, et al. Lactate and lactylation in sepsis: A comprehensive review. J Inflammation Res. (2024) 17:4405–17. doi: 10.2147/JIR.S459185

PubMed Abstract | Crossref Full Text | Google Scholar

6. Khan AA, Ata F, Yousaf Z, Aljafar MS, Seijari MN, Matarneh A, et al. A retrospective study on comparison of clinical characteristics and outcomes of diabetic ketoacidosis patients with and without acute pancreatitis. Sci Rep. (2023) 13:4347. doi: 10.1038/s41598-023-31465-3

PubMed Abstract | Crossref Full Text | Google Scholar

7. Won P, Craig J, Choe D, Collier Z, Gillenwater TJ, and Yenikomshian HA. Blood glucose control in the burn intensive care unit: A narrative review of literature. Burns: J Int Soc Burn Injuries. (2023) 49:1788–95. doi: 10.1016/j.burns.2023.06.002

PubMed Abstract | Crossref Full Text | Google Scholar

8. Shan L, Zheng K, Dai W, Hao P, and Wang Y. J-shaped association between serum glucose potassium ratio and prognosis in heart failure with preserved ejection fraction with stronger predictive value in non-diabetic patients. Sci Rep. (2024) 14:29965. doi: 10.1038/s41598-024-81289-y

PubMed Abstract | Crossref Full Text | Google Scholar

9. Gao H, Xiang Q, Li J, Yu M, Lan Y, Ba J, et al. Clinical analysis of the serum muscle enzyme spectrum of patients with newly diagnosed Sheehan’s syndrome. Medicine. (2022) 101:e30834. doi: 10.1097/MD.0000000000030834

PubMed Abstract | Crossref Full Text | Google Scholar

10. Hagerman A, Schorer R, Putzu A, Keli-Barcelos G, and Licker M. Cardioprotective effects of glucose-insulin-potassium infusion in patients undergoing cardiac surgery: A systematic review and meta-analysis. Semin Thorac Cardiovasc Surg. (2024) 36:167–81. doi: 10.1053/j.semtcvs.2022.11.002

PubMed Abstract | Crossref Full Text | Google Scholar

11. Cheng NC, Tai HC, Chang SC, Chang CH, and Lai HS. Necrotizing fasciitis in patients with diabetes mellitus: clinical characteristics and risk factors for mortality. BMC Infect Dis. (2015) 15:417. doi: 10.1186/s12879-015-1144-0

PubMed Abstract | Crossref Full Text | Google Scholar

12. Koyuncu S, Sipahioglu H, Bol O, İlik HKZ, Dilci A, Elmaağaç M, et al. The evaluation of different treatment approaches in patients with earthquake-related crush syndrome. Cureus. (2023) 15:e47194. doi: 10.7759/cureus.47194

PubMed Abstract | Crossref Full Text | Google Scholar

13. Güler S, Demirtaş E, Üçöz Kocaşaban D, Sarıhan MB, Esen E, Ata MA, et al. Evaluation of the predictive value of the glucose-to-potassium ratio in predicting in-hospital mortality of patients with sepsis and septic shock. Untersuchung des prädiktiven Werts des Glukose-Kalium-Quotienten für die Vorhersage der Krankenhausmortalität bei Patienten mit Sepsis und septischem Schock. Medizinische Klinik Intensivmedizin und Notfallmedizin. (2025). doi: 10.1007/s00063-024-01244-7

PubMed Abstract | Crossref Full Text | Google Scholar

14. Monard C, Tebib N, Trächsel B, Kelevina T, and Schneider AG. Comparison of methods to normalize urine output in critically ill patients: a multicenter cohort study. Crit Care (Lond Engl). (2024) 28:425. doi: 10.1186/s13054-024-05200-x

PubMed Abstract | Crossref Full Text | Google Scholar

15. Seckel MA. Sepsis best practices: Definitions, guidelines, and updates. Nursing. (2024) 54:31–9. doi: 10.1097/NSG.0000000000000010

PubMed Abstract | Crossref Full Text | Google Scholar

16. Liu Y, Qiu T, Fu Z, Wang K, Zheng H, Li M, et al. Systemic immune-inflammation index and serum glucose-potassium ratio predict poor prognosis in patients with spontaneous cerebral hemorrhage: An observational study. Medicine. (2024) 103:e39041. doi: 10.1097/MD.0000000000039041

PubMed Abstract | Crossref Full Text | Google Scholar

17. Lange S, Lange R, Tabibi E, Hitschold T, Müller VI, and Naumann G. Comparison of vaginal pessaries to standard care or pelvic floor muscle training for treating postpartum urinary incontinence: a pragmatic randomized controlled trial. Geburtshilfe und Frauenheilkunde. (2024) 84:246–55. doi: 10.1055/a-2243-3784

PubMed Abstract | Crossref Full Text | Google Scholar

18. Tramonti Fantozzi MP, Ceccarelli L, Petri D, De Vita E, Agostini A, Colombatto P, et al. Hepatitis C epidemiology and treatment outcomes in Italy: Impact of the DAA era and the COVID-19 pandemic. J Viral hepatitis. (2024) 31:623–32. doi: 10.1111/jvh.13983

PubMed Abstract | Crossref Full Text | Google Scholar

19. Rebernick RJ, Bell HN, and Wakeam E. Survival analyses: A statistical review for surgeons. Semin Thorac Cardiovasc Surg. (2022) 34:1388–94. doi: 10.1053/j.semtcvs.2022.01.001

PubMed Abstract | Crossref Full Text | Google Scholar

20. Liu C, Deng L, Lin S, Liu T, Ren J, Shi J, et al. Enteral nutrition support in patients with cancer: association of short-term prognosis and medical costs with inflammation. Supportive Care Cancer. (2024) 33:50. doi: 10.1007/s00520-024-09085-y

PubMed Abstract | Crossref Full Text | Google Scholar

21. Dai T, Liu M, Bao D, Manor B, and Zhou J. Transcranial direct current stimulation alleviates the pain severity in people suffering from knee osteoarthritis: a systematic review and meta-analysis. Pain Rep. (2024) 10:e1215. doi: 10.1097/PR9.0000000000001215

PubMed Abstract | Crossref Full Text | Google Scholar

22. Pu Y, Xing N, Wang Y, Wang H, Xu J, and Li X. Differential impact of TyG and TyG-BMI indices on short- and long-term mortality in critically ill ischemic stroke patients. BMC Cardiovasc Disord. (2024) 24:754. doi: 10.1186/s12872-024-04450-5

PubMed Abstract | Crossref Full Text | Google Scholar

23. Demir FA, Ersoy İ, Yılmaz AŞ, Taylan G, Kaya EE, Aydın E, et al. Serum glucose-potassium ratio predicts inhospital mortality in patients admitted to coronary care unit. Rev da Associacao Med Bras (1992). (2024) 70:e20240508. doi: 10.1590/1806-9282.20240508

PubMed Abstract | Crossref Full Text | Google Scholar

24. Lu Y, Ma X, Zhou X, and Wang Y. The association between serum glucose to potassium ratio on admission and short-term mortality in ischemic stroke patients. Sci Rep. (2022) 12:8233. doi: 10.1038/s41598-022-12393-0

PubMed Abstract | Crossref Full Text | Google Scholar

25. Alamri FF, Almarghalani DA, Alraddadi EA, Alharbi A, Algarni HS, Mulla OM, et al. The utility of serum glucose potassium ratio as a predictive factor for haemorrhagic transformation, stroke recurrence, and mortality among ischemic stroke patients. Saudi Pharm J. (2024) 32:102082. doi: 10.1016/j.jsps.2024.102082

PubMed Abstract | Crossref Full Text | Google Scholar

26. Zhang Q, Huang Z, Chen S, Yan E, Zhang X, Su M, et al. Association between the serum glucose-to-potassium ratio and clinical outcomes in ischemic stroke patients after endovascular thrombectomy. Front Neurol. (2024) 15:1463365. doi: 10.3389/fneur.2024.1463365

PubMed Abstract | Crossref Full Text | Google Scholar

27. Chen Y, Peng Y, Zhang X, Liao X, Lin J, Chen L, et al. The blood glucose-potassium ratio at admission predicts in-hospital mortality in patients with acute type A aortic dissection. Sci Rep. (2023) 13:15707. doi: 10.1038/s41598-023-42827-2

PubMed Abstract | Crossref Full Text | Google Scholar

28. Zhou J, Yang CS, Shen LJ, Lv QW, and Xu QC. Usefulness of serum glucose and potassium ratio as a predictor for 30-day death among patients with severe traumatic brain injury. Clin Chim Acta. (2020) 506:166–71. doi: 10.1016/j.cca.2020.03.039

PubMed Abstract | Crossref Full Text | Google Scholar

29. Zhou H, Tang Y, Li Y, and Zhu T. Serum glucose-potassium ratio predicts prognosis of traumatic brain injury in patients undergoing emergency craniotomy: A retrospective study. Asian J Surg. (2023) 46:2958–9. doi: 10.1016/j.asjsur.2023.02.034

PubMed Abstract | Crossref Full Text | Google Scholar

30. Zhou W, Liu Y, Wang Z, Mao Z, and Li M. Serum glucose/potassium ratio as a clinical risk factor for predicting the severity and prognosis of acute traumatic spinal cord injury. BMC Musculoskeletal Disord. (2023) 24:870. doi: 10.1186/s12891-023-07013-5

PubMed Abstract | Crossref Full Text | Google Scholar

31. Fujiki Y, Matano F, Mizunari T, Murai Y, Tateyama K, Koketsu K, et al. Serum glucose/potassium ratio as a clinical risk factor for aneurysmal subarachnoid hemorrhage. J Neurosurg. (2018) 129:870–5. doi: 10.3171/2017.5.JNS162799

PubMed Abstract | Crossref Full Text | Google Scholar

32. Matano F, Fujiki Y, Mizunari T, Koketsu K, Tamaki T, Murai Y, et al. Serum glucose and potassium ratio as risk factors for cerebral vasospasm after aneurysmal subarachnoid hemorrhage. J Stroke Cerebrova Dis. (2019) 28:1951–7. doi: 10.1016/j.jstrokecerebrovasdis.2019.03.041

PubMed Abstract | Crossref Full Text | Google Scholar

33. Boyuk F. The predictor potential role of the glucose to potassium ratio in the diagnostic differentiation of massive and non-massive pulmonary embolism. Clin Appl thrombosis/hemostasis. (2022) 28:10760296221076146. doi: 10.1177/10760296221076146

PubMed Abstract | Crossref Full Text | Google Scholar

34. Katipoğlu B and Demirtaş E. Assessment of serum glucose potassium ratio as a predictor for morbidity and mortality of blunt abdominal trauma. Künt Abdominal Travma Hastalarında Serum Glikoz Potasyum Oranının Morbidite Ve Mortalite İçin Bir Öngörücü Olarak Değerlendirilmesi. Ulusal travma ve acil cerrahi dergisi = Turkish J Trauma Emergency Surg: TJTES. (2022) 28:134–9. doi: 10.14744/tjtes.2020.88945

PubMed Abstract | Crossref Full Text | Google Scholar

35. Wasyluk W, Wasyluk M, and Zwolak A. Sepsis as a pan-endocrine illness-endocrine disorders in septic patients. J Clin Med. (2021) 10:2075. doi: 10.3390/jcm10102075

PubMed Abstract | Crossref Full Text | Google Scholar

36. Maas MB, Lizza BD, Kim M, Abbott SM, Gendy M, Reid KJ, et al. Stress-induced behavioral quiescence and abnormal rest-activity rhythms during critical illness. Crit Care Med. (2020) 48:862–71. doi: 10.1097/CCM.0000000000004334

PubMed Abstract | Crossref Full Text | Google Scholar

37. Miyauchi H, Fujioka K, Okubo S, Nishida K, Ashina M, Ikuta T, et al. Insulin therapy for hyperglycemia in neonatal sepsis using a preterm mouse model. Pediatr Int. (2020) 62:581–6. doi: 10.1111/ped.14126

PubMed Abstract | Crossref Full Text | Google Scholar

38. Yi H, Duan Y, Song R, Zhou Y, Cui Y, Liu C, et al. Activation of glucagon-like peptide-1 receptor in microglia exerts protective effects against sepsis-induced encephalopathy via attenuating endoplasmic reticulum stress-associated inflammation and apoptosis in a mouse model of sepsis. Exp Neurol. (2023) 363:114348. doi: 10.1016/j.expneurol.2023.114348

PubMed Abstract | Crossref Full Text | Google Scholar

39. Fan W, Wang C, Xu K, Liang H, and Chi Q. Ccl5+ Macrophages drive pro-inflammatory responses and neutrophil recruitment in sepsis-associated acute kidney injury. Int Immunopharmacol. (2024) 143:113339. doi: 10.1016/j.intimp.2024.113339

PubMed Abstract | Crossref Full Text | Google Scholar

40. Yang X, Pu X, Xu Y, Zhao J, Fang X, Cui J, et al. A novel prognosis evaluation indicator of patients with sepsis created by integrating six microfluidic-based neutrophil chemotactic migration parameters. Talanta. (2025) 281:126801. doi: 10.1016/j.talanta.2024.126801

PubMed Abstract | Crossref Full Text | Google Scholar

41. Binny R, Kotsanas D, Buttery J, Korman T, and Tan K. Is neutrophil to lymphocyte ratio an accurate predictor of neonatal sepsis in premature infants? Early Hum Dev. (2025) 200:106147. doi: 10.1016/j.earlhumdev.2024.106147

PubMed Abstract | Crossref Full Text | Google Scholar

42. Guo Y, Qiu Y, Xue T, Yan P, Zhao W, Wang M, et al. Association between admission baseline blood potassium levels and all-cause mortality in patients with acute kidney injury combined with sepsis: A retrospective cohort study. PLoS One. (2024) 19:e0309764. doi: 10.1371/journal.pone.0309764

PubMed Abstract | Crossref Full Text | Google Scholar

43. Yu J, Fu Y, Zhang N, Gao J, Zhang Z, Jiang X, et al. Extracellular histones promote TWIK2-dependent potassium efflux and associated NLRP3 activation in alveolar macrophages during sepsis-induced lung injury. Inflammation Res. (2024) 73:1137–55. doi: 10.1007/s00011-024-01888-3

PubMed Abstract | Crossref Full Text | Google Scholar

44. Contrera Rolón N, Cantos J, Huespe I, Prado E, Bratti GI, Schreck C, et al. Fractional excretion of sodium and potassium and urinary strong ion difference in the evaluation of persistent AKI in sepsis. Medicina Intensiva. (2024) 49(1):1–7. doi: 10.1016/j.medine.2024.02.003

PubMed Abstract | Crossref Full Text | Google Scholar

45. Luo W, Xiong L, Wang J, Li C, and Zhang S. Development and performance evaluation of a clinical prediction model for sepsis risk in burn patients. Medicine. (2024) 103:e40709. doi: 10.1097/MD.0000000000040709

PubMed Abstract | Crossref Full Text | Google Scholar

46. Bouadma L, Mankikian S, Darmon M, Argaud L, Vinclair C, Siami S, et al. Influence of dyskalemia at admission and early dyskalemia correction on survival and cardiac events of critically ill patients. Crit Care (Lond Engl). (2019) 23:415. doi: 10.1186/s13054-019-2679-z

PubMed Abstract | Crossref Full Text | Google Scholar

47. López A, Varela JJ, Cid MM, Couñago M, and Gago N. Hydroelectrolytic and infectious complications in one year of parenteral nutrition in critical care. Complicaciones hidroelectrolíticas e infecciosas en un año de nutrición parenteral en cuidados críticos. Rev Espanola Anestesiol y Reanimacion. (2018) 65:373–9. doi: 10.1016/j.redar.2018.03.002

PubMed Abstract | Crossref Full Text | Google Scholar

48. Zhang JL, Chen YT, Chen GD, Wang T, Zhang JX, and Zeng QY. Glucose-insulin-potassium alleviates intestinal mucosal barrier injuries involving decreased expression of uncoupling protein 2 and NLR family-pyrin domain-containing 3 inflammasome in polymicrobial sepsis. BioMed Res Int. (2017) 2017:4702067. doi: 10.1155/2017/4702067

PubMed Abstract | Crossref Full Text | Google Scholar

49. Duan W, Yang F, Ling H, Li Q, and Dai X. Association between lactate to hematocrit ratio and 30-day all-cause mortality in patients with sepsis: a retrospective analysis of the Medical Information Mart for Intensive Care IV database. Front Med. (2024) 11:1422883. doi: 10.3389/fmed.2024.1422883

PubMed Abstract | Crossref Full Text | Google Scholar

50. Shi Y, Hu Y, Xu GM, and Ke Y. Development and validation of a predictive model for pulmonary infection risk in patients with traumatic brain injury in the ICU: a retrospective cohort study based on MIMIC-IV. BMJ Open Respir Res. (2024) 11:e002263. doi: 10.1136/bmjresp-2023-002263

PubMed Abstract | Crossref Full Text | Google Scholar

51. Zhang Y, Zhou C, Huang Z, and Ye X. Study of cuffless blood pressure estimation method based on multiple physiological parameters. Physiol Measurement. (2021) 42. doi: 10.1088/1361-6579/abf889

PubMed Abstract | Crossref Full Text | Google Scholar

52. Simons PIHG, Simons N, Stehouwer CDA, Schalkwijk CG, Schaper NC, and Brouwers MCGJ. Association of common gene variants in glucokinase regulatory protein with cardiorenal disease: A systematic review and meta-analysis. PLoS One. (2018) 42(5):e0206174. doi: 10.1371/journal.pone.0206174

PubMed Abstract | Crossref Full Text | Google Scholar

53. Luo J, Huang S, Lan L, Yang S, Cao T, Yin J, et al. EMR-LIP: A lightweight framework for standardizing the preprocessing of longitudinal irregular data in electronic medical records. Comput Methods Programs Biomed. (2025) 259:108521. doi: 10.1016/j.cmpb.2024.108521

PubMed Abstract | Crossref Full Text | Google Scholar

54. Russell SL, Penunuri G, and Condon C. Diverse genetic conflicts mediated by molecular mimicry and computational approaches to detect them. Semin Cell Dev Biol. (2025) 165:1–12. doi: 10.1016/j.semcdb.2024.07.001

PubMed Abstract | Crossref Full Text | Google Scholar

55. Gu D and Huang S. Letter to the editor on ‘Aspirin is associated with improved outcomes in patients with sepsis-induced myocardial injury: An analysis of the MIMIC-IV database’. J Clin Anesth. (2025) 100:111683. doi: 10.1016/j.jclinane.2024.111683

PubMed Abstract | Crossref Full Text | Google Scholar

56. Park YJ, Bae J, Yoo JK, Ahn SH, Park SY, Kim YS, et al. Effects of NF-κB inhibitor on sepsis depend on the severity and phase of the animal sepsis model. J Personal Med. (2024) 14:645. doi: 10.3390/jpm14060645

PubMed Abstract | Crossref Full Text | Google Scholar

57. Alcamo AM, Becker AE, Barren GJ, Hayes K, Pennington JW, Curley MAQ, et al. Diagnostic identification of acute brain dysfunction in pediatric sepsis and septic shock in the electronic health record: A comparison of four definitions in a reference dataset. Pediatr Crit Care Med. (2024) 25(8):740–7. doi: 10.1097/PCC.0000000000003529

PubMed Abstract | Crossref Full Text | Google Scholar

58. Pérez-Tome JC, Parrón-Carreño T, Castaño-Fernández AB, Nievas-Soriano BJ, and Castro-Luna G. Sepsis mortality prediction with Machine Learning Tecniques. Med Intensiva (Engl Ed. (2024) 48(10):584–593. doi: 10.1016/j.medine.2024.05.009

PubMed Abstract | Crossref Full Text | Google Scholar

59. Yang B, Niu K, Zhu Y, Zheng X, Li T, Wang Z, et al. Effects of ondansetron exposure during ICU stay on outcomes of critically ill patients with sepsis: a cohort study. Front Cell Infect Microbiol. (2023) 13:1256382. doi: 10.3389/fcimb.2023.1256382

PubMed Abstract | Crossref Full Text | Google Scholar

60. Yu X, Xin Q, Hao Y, Zhang J, and Ma T. An early warning model for predicting major adverse kidney events within 30 days in sepsis patients. Front Med (Lausanne). (2024) 10:1327036. doi: 10.3389/fmed.2023.1327036

PubMed Abstract | Crossref Full Text | Google Scholar

61. Fang Y, Xiong B, Shang X, Yang F, Yin Y, Sun Z, et al. Triglyceride-glucose index predicts sepsis-associated acute kidney injury and length of stay in sepsis: A MIMIC-IV cohort study. Heliyon. (2024) 10:e29257. doi: 10.1016/j.heliyon.2024.e29257

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: intensive care unit, MIMIC, mortality, sepsis, glucose potassium ratio, long term, Cox regression

Citation: Lou J, Xiang Z, Zhu X, Song J, Cui S, Li J, Jin G, Huang N, Fan Y and Xu S (2025) Association between serum glucose potassium ratio and short- and long-term all-cause mortality in patients with sepsis admitted to the intensive care unit: a retrospective analysis based on the MIMIC-IV database. Front. Endocrinol. 16:1555082. doi: 10.3389/fendo.2025.1555082

Received: 03 January 2025; Accepted: 15 July 2025;
Published: 30 July 2025.

Edited by:

Darko Stefanovski, University of Pennsylvania, United States

Reviewed by:

Xiaozhu Liu, Capital Medical University, China
Yong Sun, Xuzhou Medical University, China

Copyright © 2025 Lou, Xiang, Zhu, Song, Cui, Li, Jin, Huang, Fan and Xu. 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: Youfen Fan, MTM5MDY2ODM2MTNAMTYzLmNvbQ==; Sida Xu, bmJleXh1c2lkYUAxMjYuY29t

These authors have contributed equally to this work

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.