Authors: Elie Sarraf (1Department of Anesthesiology and Perioperative Medicine, Penn State Milton S Hershey Medical Center, Hershey, PA 17036, USA), Alireza Vafaei Sadr (2Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA, USA), Vida Abedi (2Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA, USA), Anthony S Bonavia (3Division of Critical Care Medicine, Department of Anesthesiology and Perioperative Medicine, Penn State Milton S Hershey Medical Center, Hershey, PA 17036, USA)
Categories: Article, Sepsis, machine learning, Sequential Organ Failure Assessment score, MIMIC database, mortality
Source: Journal of critical care
Authors: Elie Sarraf, Alireza Vafaei Sadr, Vida Abedi, Anthony S Bonavia
The Sequential Organ Failure Assessment (SOFA) score monitors organ failure and defines sepsis but may not fully capture factors influencing sepsis mortality. Socioeconomic and demographic impacts on sepsis outcomes have been highlighted recently.
To evaluate the prognostic value of SOFA scores against demographic and social health determinants for predicting sepsis mortality in critically ill patients, and to assess if a combined model increases predictive accuracy.
The study utilized retrospective data from the MIMIC-IV database and prospective external validation from the Penn State Health cohort. A Random Forest model incorporating SOFA scores, demographic/social data, and the Charlson Comorbidity Index was trained and validated.
In the MIMIC-IV dataset of 32,970 sepsis patients, 6,824 (20.7%) died within 30 days. A model including demographic, socioeconomic, and comorbidity data with SOFA scores improved predictive accuracy beyond SOFA scores alone. Day 2 SOFA, age, weight, and comorbidities were significant predictors. External validation showed consistent performance, highlighting the importance of delta SOFA between days 1 and 3.
Adding patient-specific demographic and socioeconomic information to clinical metrics significantly improves sepsis mortality prediction. This suggests a more comprehensive, multidimensional prognostic approach is needed for accurate sepsis outcome predictions.
The Sepsis-related Organ Failure Assessment (SOFA) score, initially devised by Vincent et al. [1], was designed to understand the progression and impact of organ dysfunction in sepsis, complementing other severity scores. Later known as the ‘Sequential’ Organ Failure Assessment score, it was incorporated into the diagnostic criteria for sepsis (Sepsis-3) in 2016 [2]. The SOFA score provides an objective measure of sequential organ dysfunction or recovery during a patient’s hospitalization, contrasting with other scoring systems like the APACHE II score [3]. It relies heavily on laboratory data, which may not always reflect the current clinical condition of critically ill patients. Additionally, among the six SOFA domains, four hinge on laboratory results, potentially delaying the identification of critical changes in a patient’s condition. Although other severity-of-illness scores exist, no scoring system exclusively predicts mortality in sepsis [4, 5]. Recognizing its limitations, a review is currently underway to update and improve the SOFA score for more accurate, sepsis-specific applications [6].
Recent literature suggests that socioeconomic status and demographic background may play a more definitive role in the outcomes of sepsis than previously recognized [7–10]. In fact, artificial intelligence has facilitated data-driven approaches to uncovering existing healthcare disparities [11]. These social determinants of health (SDoH) may provide immediate, specific insight into patient risk without the need for expensive tests. Thus, our study sought to assess the predictive power of SOFA scores relative to demographic and socioeconomic factors for sepsis prognosis. We hypothesized that a combined approach using SOFA scores measured on the days following sepsis onset, together with social health determinants, could offer a more accurate prediction of mortality.
We accessed anonymized patient data from the Medical Information Mart for Intensive Care (MIMIC)-IV database (version 2.2, Jan 6 2023), which includes critical care data from Beth Deaconess Medical Center [12]. This resource was chosen for its comprehensive data and patient mortality information. Using existing code [13], we calculated SOFA scores for patients fitting Sepsis-3 criteria [2].
From the onset of suspected infection, we tracked the highest daily SOFA score and its change between each day of critical illness. Additionally, we computed daily changes in SOFA (delta SOFA) and cumulative daily SOFA scores (sum SOFA). This method ensured data completeness, with no need for imputation except when patients were discharged from the ICU. We also gathered demographic and sociological data, specifically gender, age, weight, height, insurance status (Medicare, Medicaid, private insurance), marital status and race as well as the Charlson Comorbidity Index (CCI, [14]), focusing on 30-day mortality as the primary outcome.
Using Python’s scikit-learn library (version 1.3.2)[15], we trained a Random Forest model, evaluating it through 5-fold cross-validation. We tested six model types, including (1) models that were trained on 3-day and 8-day retrospective data, (2) models with and without daily SOFA organ component measures, and (3) models with and without patient socioeconomic, demographic and comorbidity data. Metrics such as feature importance, PPV, NPV, sensitivity, specificity, and AUROC were assessed during each iteration.
We then externally validated our best predictive model with real-world, observational data from critically ill patients having sepsis and forming part of a prospective research cohort at Penn State Health (PSH) from Aug 2020 – Feb 2024. Eligible participants were adults over 18 years, identified within 48 hours of critical illness onset. Sepsis was defined according to Sepsis-3 criteria, and critical illness required continuous intravenous vasopressors or ongoing respiratory support. Patients with active hematologic malignancies or those on immune-altering therapies were excluded.
Given that the validation data set included SOFA scores measured every other day, the original model was re-tuned to only use odd-numbered days. We also compared feature importance between a streamlined dataset for a Random Forest model with both the MIMIC and PSH data. Our analysis code is available on Github, and the study adhered to TRIPOD guidelines for predictive model reporting [16].
Of 299712 unique patients included in the MIMIC-IV database, 32970 had sepsis. Among these, 57.8% were male and the mean age was 66.7 ± 16 years. Supplemental Table 1 details the SDoH factors included in the study. A total of 6824 patients (20.7%) died within 30 days of the sepsis onset.
For each of the six scenarios described above, 25 different machine learning models were trained. Figure 1A describes five metrics pertaining to the model having the highest comparative AUROC (75%, with 95% CI 73 – 77%). This model demonstrated 77% sensitivity (95% CI 75 – 77%), 74% accuracy (95% CI 72 – 76%) and a precision of 43% (95% CI 39 – 47%). The most effective model combined demographic, socioeconomic, and comorbidity data with total SOFA score and SOFA organ component measures over an 8-day period post-sepsis. However, adding just demographic, socioeconomic, and comorbidity data to the total SOFA score (i.e., without individual organ dysfunction measures) significantly enhanced 30-day mortality prediction as compared with total SOFA score alone. This improvement was evident in analyses conducted at both 3 and 8 days after sepsis onset.
With respect to feature importance, SOFA organ component scores measured on day 2 were most predictive of 30-day mortality (Figure 1B). Amongst organ systems affected by sepsis, cardiovascular dysfunction, renal dysfunction, and central nervous system dysfunction were most predictive of 30-day mortality. Equally important was the relative feature importance of age, weight, height, marital status and CCI as compared with organ-specific or overall SOFA measures (Figure 1B).
Figure 2A illustrates the results of external validation using prospective data from 105 septic, critically ill patients undergoing care at PSH. In this data set, 52.4% of patients were male, with a mean age was 66.5 ± 15 years. The microbial sources of sepsis in this cohort are enumerated in Supplementary Table 2. Nineteen patients (18.1%) died within 30 days. When individual metrics were compared between discovery and validation cohorts, there was a noted increase in sensitivity and negative predictive value, and to a minor degree AUROC, at the expense of a lower specificity and positive predictive value. The negative predictive value of the model derived from MIMIC data appeared to decrease when using data derived from every other day of critical illness. Otherwise, the model’s performance remained consistent.
Interestingly, in the streamlined model, age, weight, height and CCI continued to comprise the most important features (Figure 2B). The importance of marital status was variable as compared with results derived from the MIMIC model, while ‘delta SOFA’ between days 1 – 3 of critical illness proved to be far better correlated with mortality in the PSH data set.
Historically, a higher SOFA score, signaling severe organ failure, has been closely linked with poorer clinical outcomes [17, 18]. However, our analysis proposes a paradigm shift. Our pivotal finding is the pronounced impact of patient-specific and social risk factors on 30-day mortality, overshadowing the predictive relevance of organ dysfunction severity [10].
Our data indicate that, by day 2 of sepsis, certain organ dysfunction measures can predict mortality with some degree of reliability. This finding is consistent with previous literature that indicates enhanced prognostic value of SOFA score when measured later in a patient’s hospitalization [19, 20]. Yet, it is the integration of a patient’s age, weight, height, marital status and comorbidity profile (as encapsulated by the Charlson Comorbidity Index) that significantly amplifies prognostic precision as compared with the SOFA score alone or any of its individual organ components. By Day 8, models enriched with these variables not only sustained but also enhanced the predictive accuracy of mortality, with model performance metrics surpassing those based solely on clinical measures.
The inclusion of individual organ dysfunction parameters only modestly improved the model’s performance, reinforcing the premise that, while clinical measures of organ dysfunction are not to be overlooked, they are evidently less predictive of patient outcomes compared to SDoH. This is a critical observation, as it underscores the limitations of current clinical-only prognostic models and highlights the potential for improved risk stratification through the incorporation of other, readily available, patient-specific factors.
We externally validated the best model generated from publicly available data on a subset of patients receiving health care at our institution. We found that, while the predictive model remained stable, delta SOFA between days 1 and 3 of critical illness was a far stronger predictor of mortality in our patient population. The prognostic utility of delta SOFA has been previously reported [21], although it is not currently in widespread clinical use.
Past efforts to improve on the SOFA score have produced mixed results [22–24]. However, our research supports these endeavors, suggesting that an enhanced SOFA score incorporating SDoH could sharpen the predictive accuracy for mortality in sepsis patients. This augmented model may be particularly beneficial for patients who, due to cognitive impairments caused by their illness, are unable to provide a detailed medical history. Given the large regional variability in socioeconomic status and ethnic diversity, it may also be necessary to adjust the model to reflect this local variability. Methods such as federated learning can be utilized for this purpose [25].
Although our study offers valuable insights, it is important to acknowledge its limitations. A significant constraint was that our validation dataset recorded SOFA scores every 48 hours instead of every 24 hours. This less frequent recording could have affected the accuracy of validation, given that peak SOFA scores usually manifest within the first 24 hours of sepsis onset. Furthermore, our exclusive use of the Sepsis-3 diagnostic criteria might not encompass the full range of clinical presentations and outcomes. This is particularly relevant since other frameworks, like Sepsis-2 [26], are still widely used in clinical practice. Finally, as the validation cohort was small, with minimal numbers for many subgroups, one must be cautious in interpreting the data in a general manner; rather, it should be approached as a preliminary glance at the model performance. Further external validation, using diverse populations, a community healthcare setting and temporal separation of development/validation datasets are underway to address this limitation.
Our study’s findings have implications beyond the immediate scope of sepsis and can be applied to various models across different medical domains [27, 28]. Integrating SDoH and demographic data with clinical scoring systems can enhance predictive accuracy and provide a more holistic understanding of patient outcomes. This approach is particularly relevant in chronic diseases such as diabetes and heart failure, where socioeconomic factors profoundly impact disease progression and mortality. Moreover, incorporating SDoH into predictive models for infectious diseases like COVID-19 can improve resource allocation and patient management in pandemics, highlighting the broader applicability of our findings across diverse medical fields.
In conclusion, our investigation underscores the imperative to revisit the prognostic frameworks for sepsis. It is evident that a multidimensional approach, encompassing both clinical and non-clinical factors, is crucial for a more accurate prediction of outcomes. Our work contributes to the growing body of evidence that supports the integration of broader patient information, extending beyond the confines of physiological and laboratory measures, into prognostic models for sepsis.