Abstract

Objective. The hemoglobin-to-red cell distribution width ratio (HRR) is associated with the prognosis of sepsis-associated encephalopathy (SAE). This study aimed to determine the relationship between HRR and SAE and to clarify the possible mechanism of HRR as a prognostic factor for SAE. Methods. A multivariate Cox proportional-hazards regression model was used to assess the correlation between HRR and all-cause mortality. Piecewise linear regression and smooth-curve Cox proportional-hazards regression models were used to observe whether there was a nonlinear relationship between HRR and all-cause mortality in SAE. Results. This study included 8853 patients with SAE. A nonlinear relationship between HRR and SAE was observed through a two-segment regression model. The left inflection point for the HRR threshold was calculated to be 15.54, which was negatively correlated with all-cause mortality (HR = 0.83, 95% CI = 0.76–0.91, ). Subgroup analyses revealed significant interactions between white blood cell count, glucose, and patients who received dialysis and HRR. The inverse correlation between HRR and SAE was more pronounced in patients who did not receive vasopressin (HR = 0.91, 95% CI = 0.87–0.96, ) than in those who did receive vasopressin (HR = 0.94, 95% CI = 0.88–1.02, ) and was significantly more pronounced in patients without myocardial infarction (HR = 0.91, 95% CI = 0.88–0.96, ) than in those with myocardial infarction (HR = 0.94, 95% CI = 0.87–1.02, ). Conclusion. This large retrospective study found a nonlinear relationship between all-cause mortality and HRR in patients with SAE in intensive care units, with low HRR being inversely associated with increased all-cause mortality in patients with SAE.

1. Introduction

Sepsis-associated encephalopathy (SAE), which is a syndrome of diffuse neurological dysfunction, with no associated infectious foci in the central nervous system (CNS). Since host dysregulation does not occur in CNS infections, clear diagnostic criteria are still lacking for clinical practice [1, 2]. SAE occurrence not only affects the autonomic nervous system and neuroendocrine function of patients with general sepsis but also leads to different degrees of cognitive decline, and even accelerates the occurrence and progression of multiple organ failure in sepsis, thus somewhat prolonging the length of stay in intensive care units (ICU). Statistics suggest that the mortality rate caused by SAE in patients with sepsis in ICU has been well over 60% [3]. However, for SAE diagnosis, in addition to the basic screening of the nervous system, EEG, CT, and MRI examinations can only be used to exclude organic CNS lesions. Identifying more risk factors is therefore important for early assessments of survival and prognosis of patients with SAE [4].

Various clinical biomarkers are currently in the diagnosis and prognosis of SAE. Hemoglobin (Hb) is widely applied to the overall anemic condition of patients with sepsis, including nutrition, inflammation, and general condition. There is a very strong correlation between sepsis and encephalopathy caused by shock and prognosis [5]. The red-blood-cell distribution width (RDW) was used as a biochemical index to measure red blood cell (RBC) variation, and increased RDW is not only indicative of malnutrition caused by oversized or undersized RBC, stress response inflammation, dyslipidemia, hypertension, and other abnormal reactions but can also predict early coronary heart disease, cerebral infarction, pulmonary embolism, other cardiovascular and cerebrovascular diseases, and cancer. Previous studies have investigated the effect of RDW on patients with sepsis [68]. Hb and RDW are two important clinical parameters. The Hb/RDW ratio (HRR) can especially reflect the early abnormal manifestations of SAE. Meanwhile, no specific SAE biomarker has been found for use in clinical practice [9, 10].

This study therefore utilized data from the Medical Information Mart for Intensive Care-IV (MIMIC-IV) database, which is a generalized and publicly available database that includes baseline data from many patients with critical diseases, to analyze the connection between the HRR and in-hospital deaths due to SAE. The aim was to determine whether HRR can be used as an important specific indicator to distinguish SAE from other inflammatory diseases and to provide new insights for early and effective SAE interventions.

2. Materials and Methods

2.1. Data Source

This retrospective study was based on the MIMIC-IV database [11, 12]. The MIMIC database is a multiparameter, free-access, open-source, critical care database established by the Massachusetts Institute of Technology. This database is determined to be under the Safe Harbor Provision of the Health Insurance Portability and Accountability Act [13, 14]. It contains detailed clinical information on more than 380,000 patients in ICU during 2008–2019. We completed the collaborating agency training program (certification number: 45351934). The requirement to obtain individual patient consent was waived since the program does not affect clinical care, and all protected health information was anonymized.

2.2. Study Population and Data Extraction

In this study, we used structured query language (SQL) via Navicat Premium (version 15.0) to extract meaningful clinical baseline data based on the subject_id of each ICU patient. According to the third internationally recognized definition of sepsis (pyosis-3), sepsis was diagnosed using the latest sepsis 3 criteria, defined as life-threatening infection combined with sequential organ failure assessment (SOFA) score ≥ 2 [15]. SOFA was included according to the abovementioned criteria. Patients were screened according to the following criteria [16]: (1) patients dying or leaving within 24 h since ICU admission; (2) patients with primary brain injury (traumatic brain injury, ischemic stroke, hemorrhagic stroke, epilepsy, or intracranial infection); (3) patients with severe burn and trauma; (4) patients with chronic alcohol or drug abuse; (5) patients with severe liver disease or encephalopathy caused by kidney disease; and (6) patients without a laboratory indicator of HRR.

2.3. Sepsis-Associated Encephalopathy

As far as we know, Glasgow Coma Scale (GCS) has been proved to be an important tool to assess the degree of consciousness disturbance in SAE patients. Sepsis patients with GCS ≤ 14 and excluding other unrelated factors can be initially diagnosed as sepsis-related encephalopathy. Therefore, SAE was further defined as sepsis with GCS score ≤14 or delirium diagnosed according to ICD9 code on the first day of admission, and GCS scores assessed before surgery or before sedation and anesthesia will be included as far as possible, and excluded psychiatric disorders, neurological disorders, alcoholism, or trauma from the initial diagnosis of admission [17].

2.4. Data Extraction

We used SQL and the PostgreSQL tool (version 15.0) to extract meaningful clinical data. We followed the principle that if a measurement was made multiple times for a patient, we only used the first one. Thus, the first measurement of vital signs and laboratory data within the first 24 hours of ICU admission was included in this study. The initially selected laboratory measurements included ethnicity, age, sex, pregnancy, abrosia, dialysis, heart rate (beats/minute), mean blood pressure (MBP, mmHg), respiratory rate (breaths/minute), temperature (°C), peripheral oxygen saturation (SpO2, %), Oxford Acute Severity of Illness Score (OASIS), Simplified Acute Physiology Score II (SAPSII), Logistic Organ Dysfunction System (LODS) score, SOFA score, GCS score, myocardial infarction, congestive heart failure, peripheral vascular disease, chronic pulmonary disease, diabetes, paraplegia, RDW (×109/L), calcium (mg/dL), potassium (mmol/L), Hb (g/L), sodium (mmol/L), glucose (mg/dL), hematocrit (%), platelets (×109/L), the Hb/RDW ratio (HRR) [9, 18, 19], chloride (mmol/L), prothrombin time (PT, seconds), white blood cell (WBC) count (×109/L), anion gap (×109/L), bicarbonate (mmol/L), creatinine (mEq/L), international normalized ratio (INR), norepinephrine, and vasopressin. ICU all-cause mortality referred to in-hospital survival and was used as the outcome of this study.

3. Statistical Analysis

All baseline values of continuous variables are expressed as mean ± standard deviation or median values, and the values of categorical variables are expressed as counts and proportions. All continuous variables in this study were compared using t-tests or Wilcoxon rank-sum tests, and the resulting proportions were compared using chi-square or Fisher’s exact tests. Differences among HRR tertiles were determined using one-way ANOVA or Kruskal–Wallis H tests. We used a Cox proportional-hazards model to determine the association between HRR and in-hospital mortality in patients with SAE, and developed three different models based on this: an unadjusted model, with no adjustment for covariates; model 1, partially adjusted model (for age and sex only); and model 2, fully adjusted for potential confounders. The results are expressed as hazard ratios (HRs) with 95% confidence intervals (CIs). Results are provided for unadjusted, microadjusted, and fully adjusted analyses.

Furthermore, a generalized additive model was applied to verify the nonlinear relationship between HRR and SAE after it was confirmed using a logarithmic likelihood ratio test (). We used cubic-spline function fitting and the smooth curves of two linear regression models of all-cause mortality in the smooth figure calculation for the HRR threshold. Moreover, a recursive method was used to automatically calculate the inflection points in the fitted graph and infer the likelihood ratio of the optimal model. We analyzed the relationship between HRR and all-cause mortality in patients with SAE between different subgroups to determine whether there were differences in the effect of HRR. The statistical software package R (the R Foundation; https://www.r-project.org) was used to perform data analyses in this study.

4. Results

4.1. Baseline Results

After the above screening criteria, a total of 8853 SAE patients were included in the final cohort, and the inclusion process is shown in Figure 1. In all included population, the median age was 67 years (IQR=55–78 years) and comprised 4007 males (45.30%) and 4845 females (54.70%).. Final outcomes of in-hospital death and survival occurred in 693 (7.80%) and 8159 (92.2%) patients, respectively, with a significant difference (). Quartiles 1 (≥3.98 to <14.00), 2 (≥14.00 to <17.10), and 3 (≥17.10 to ≤28.66) of HRR values comprised 2937, 2967, and 2948 cases, respectively. Differences between groups were assessed using Pearson’s chi-square test for categorical variables and t-test for continuous variables. All the variables are listed in Table 1.

4.2. HRR and ICU Mortality

We used multivariate Cox proportional-hazards regression models to determine the relationship between HRR and ICU mortality in patients with SAE, with the unadjusted and adjusted models highlighted in Table 2. We observed that low HRR (≥3.98 to <14.00) was significantly associated with an increased all-cause mortality risk in patients with SAE (unadjusted model: HR = 0.87, 95% CI = 0.82–0.92, ; partially adjusted model: HR = 0.85, 95% CI = 0.80–0.90, ; fully adjusted model: HR = 0.86, 95% CI = 0.79–0.93, ), while high HRR (≥14.00, <17.10 and ≥17.10, ≤28.66) was not significantly associated with SAE mortality ().

4.3. Subgroup Analyses

As shown in Figure 2, we included sex, age, myocardial infarction, chronic lung disease, WBC count, PT, glucose, vasopressin, and dialysis in the subgroup analyses. Except for WBC count, the interaction between glucose and dialysis treatment was significant (interaction ), and most of the other formation interactions were not significant (). In the subgroup analyses, age <60 years (HR = 0.92, 95% CI = 0.88–0.97, ), WBC ≥ 13 × 109/L (HR = 0.91, 95% CI = 0.67–0.96, ), potassium <4.6 mmol/L (HR = 0.90, 95% CI = 0.86–0.94, ), glucose ≥165 mg/dL (HR = 0.89, 95% CI = 0.84–0.94, ), no myocardial infarction (HR = 0.91, 95% CI = 0.88–0.96, ), and chronic lung disease (HR = 0.90, 95% CI = 0.86–0.95, ), as well as these strata of previous dialysis (HR = 0.95, 95% CI = 0.90–1.00, ) and vasopressin therapy (HR = 0.91, 95% CI = 0.87–0.96, ) had significant interactions. The association between HRR and all-cause mortality in patients with SAE was similar.

4.4. Analysis of the Nonlinear Relationship between the HRR and All-Cause Mortality

This study found a significant nonlinear relationship between HRR and patients with SAE. We included the following confounding variables that were meaningful to the model: sex, age, sodium, diabetes, calcium, LODS score, SOFA score, GCS score, chloride, myocardial infarction, respiratory rate, temperature, SpO2, PT, WBC count, peripheral vascular disease, MBP, paraplegia, norepinephrine, vasopressin, anion gap, bicarbonate, blood urea nitrogen, creatinine, INR, Hb, OASIS, SAPSII, norepinephrine, chronic pulmonary disease, and congestive heart failure. Figure 3 and Table 3 display the nonlinearity of the relationship between HRR and SAE all-cause mortality. Notably, there was an inverse association between the left inflection-point threshold and all-cause mortality (HR = 0.83, 95% CI 0.76–0.91, ), while there was no significant correlation at the right inflection-point threshold (HR = 0.95, 95% CI = 0.91–1.01, ).

5. Discussion

No previous study that we were aware of investigated the connection between HRR and SAE mortality, and this was the first time that a certain relationship was assumed between the ratio of patients with SAE and HRR based on the large MIMIC-IV database. This retrospective study demonstrated a significant nonlinear association between HRR and all-cause mortality in patients with SAE. In addition, even in the presence of certain traditional factors for life-threatening cerebrovascular disease, patients with low HRR were still at a higher risk of death. This study therefore found that lower HRR remained an important predictor of mortality in patients with SAE.

SAE is most common in patients with sepsis in ICU and one of its most serious complications. Its pathogenesis has not yet been fully elucidated, but it may be related to autoimmune injury of the central body inflammatory reaction [20, 21], a blood flow regulating mechanism for the damage caused by cerebral perfusion abnormalities and microvascular dysfunction [22, 23], or a large release of cytokines in sepsis that leads to neurotransmitter imbalance [24]. Polyneuropathy is often more severe in the later stage of SAE development, which not only affects peripheral nerves but also rapidly causes cognitive deficits, personality changes, generalized epilepsy, and other symptoms. Many patients who survive SAE often have a depressive tendency, posttraumatic stress disorder, and other consciousness disorders. The clinical diagnosis of SAE is still performed as an exclusion diagnosis, and laboratory tests have not identified specific SAE biomarkers. It is therefore highly significant to identify specific early warning and diagnostic indicators among the risk factors for SAE in critically ill patients with sepsis for facilitating its early intervention and treatment [25, 26]. HRR has been previously explored as a novel inflammatory marker in many studies on the relationship between HRR and the clinical outcomes of various diseases. There is increasing evidence that HRR is a new prognostic indicator in critically ill patients with malignant diseases. For example, Qin et al. explored the nonlinear relationship between all-cause mortality and HRR value in patients with ischemic stroke [27], and Rahamim et al. found that HRR was an important prognostic tool for predicting mortality in heart failure and the hospitalization of patients with cardiovascular disease [9]. Chi et al. found that low HRR were associated with a dual risk of disease progression or cancer recurrence in patients and used them to predict poor prognosis in patients with cancer [28]. The present single-center, retrospective study confirmed the hypothesis that HRR affects the prognosis of patients with SAE. We found that SAE had an obvious nonlinear relationship with HRR, with a low HRR greatly increasing the risk of a poor prognosis for SAE. A piecewise linear regression model was used in an attempt to calculate the critical inflection point for HRR. It was found that HRR ≤ 15.57 was negatively correlated with all-cause mortality.

The underlying mechanism of HRR predicting the risk of mortality in SAE remains unclear. We attempted to explore the link from a pathophysiological point of view. First, SAE caused by prolonged septic infection often manifests as the nervous system involvement, cerebral hemorrhage, cerebral ischemia, significant hypotension, organ dysfunction, and shock, leading to massive red blood cell loss and failure to timely replantation of bone marrow hematopoiesis, resulting in anemia [29, 30]. Hb somewhat reflects the severity and prognosis of patients with SAE [31]. On the other hand, as another parameter of HRR, RDW is a variation coefficient for measuring the heterogeneity of RBC sizes [32]. Generally, a larger RDW means greater heterogeneity of RBC sizes, or larger differences in shape and size [33, 34]. It can be better explained by the decreased RBC deformation ability under the influence of inflammatory cytokines in patients with sepsis aggravating the abnormal microcirculation in blood perfusion, or the CNS inflammation induced by the oxidative stress response, which is an important indicator of the potential accelerated development of inflammation in SAE [3537]. The present study was the first to use the HRR as an independent predictive biochemical index, which was found to be significantly more stable than RDW alone, and excluded the influence of traditional factors such as cerebrovascular disease, malignancy, infectious liver disease, severe liver and kidney dysfunction, and trauma. A clear nonlinear relationship was observed between low HRR and poor prognosis in patients with SAE. Previous studies have confirmed the importance of GCS score as a neuroscientific assessment tool to identify patients with SAE at a high risk of short-term mortality early [38, 39]. Therefore, in order to strengthen early interventions for patients with SAE, HRR, and GCS assessments can be used as further prognostic indicators for patients with SAE. We look forward to providing supplementary prognostic information for patients with SAE at earlier stages.

The study had some limitations. First, this large retrospective study had a single-center retrospective design, and the selection of population and covariate factors was somewhat biased, which needs to be addressed by designing prospective multicenter studies in the future. Secondly, the laboratory indicators extracted in this study are relatively few, and the dynamic changes of biochemical indicators are not analyzed in detail. Other undiscovered confounding factors may exist, and it is difficult to further prove the causal association between HRR and SAE all-cause mortality. Therefore, further longitudinal data collection and more advanced statistical methods are needed to conduct causal inference. Thirdly, regarding the cutoff range of HRR values, there has been heterogeneity in previous studies on different critical diseases, and there is no evident standard cutoff value for this ratio in clinical practice. Finally, long-term follow-up events were not analyzed, instead the focus being on in-hospital all-cause mortality among patients with SAE. We look forward to further research with more rigor and depth, especially to address these issues through studies with a multicenter, follow-up, and prospective design.

6. Conclusion

In summary, this large retrospective study has identified a nonlinear association between all-cause mortality and HRR in patients with SAE, with a low HRR being negatively associated with increased ICU mortality in patients with SAE.

Data Availability

Publicly available datasets are analyzed in this study. The data are available on the MIMIC-IV website at https://mimic-iv.mit.edu/.

Ethical Approval

All procedures performed in studies involving human participants were in accordance with the Ethical Standards of the Institutional and National Research Committee and with the 1964 Helsinki Declaration and its later amendments or comparable Ethical Standards.

Consent was obtained for the original data collection. The need for informed consent was waived for studies using this database.

Disclosure

Xiaxuan Huang and Shiqi Yuan are the co-first authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest to disclose.

Authors’ Contributions

Jun Lyu conceptualized the research aims and planned the analyses. Xiaxuan Huang and Shiqi Yuan guided the literature review, wrote the first draft of the paper, and contributed equally to this study. Shanyuan Tan and Yitong Ling assisted with the study design and performed data collection. Shiqi Yuan, Xiaxuan Huang, and Tao Huang participated in data analysis and interpretation. Xiaxuan Huang, Shiqi Yuan and Hongtao Cheng provided comments and approved the final manuscript. Xiaxuan Huang and Shiqi Yuan have contributed equally to this work.

Acknowledgments

The authors would like to thank the Massachusetts Institute of Technology and the Beth Israel Deaconess Medical Center for the MIMIC project. This study was supported by the Clinical Frontier Technology Program of the First Affiliated Hospital of Jinan University, China (No. JNU1AF-CFTP-2022-a01235) and the Science and Technology Projects in Guangzhou, China (No. 202201020054).