ORIGINAL RESEARCH article

Front. Pharmacol., 26 August 2022

Sec. Drugs Outcomes Research and Policies

Volume 13 - 2022 | https://doi.org/10.3389/fphar.2022.938028

Development and internal validation of simplified predictive scoring (ICU-SEPSA score) for mortality in patients with multidrug resistant infection

  • 1. Division of Drug Information, Department of Pharmacy, Phrae Hospital, Phrae, Thailand

  • 2. Division of Clinical Pharmacy, Department of Pharmaceutical Care, School of Pharmaceutical Sciences, University of Phayao, Phayao, Thailand

  • 3. Department of Pharmacy Practice, Faculty of Pharmaceutical Sciences, Chulalongkorn University, Bangkok, Thailand

  • 4. College of Pharmacy, Rangsit University, Phathum Thani, Thailand

  • 5. Department of Pharmacy Practice, University of Illinois at Chicago, Chicago, IL, United States

  • 6. Division of Social and Administration Pharmacy, Department of Pharmaceutical Care, School of Pharmaceutical Sciences, University of Phayao, Phayao, Thailand

  • 7. Center of Health Outcomes Research and Therapeutic Safety (Cohorts), School of Pharmaceutical Sciences, University of Phayao, Phayao, Thailand

  • 8. Unit of Excellence on Clinical Outcomes Research and IntegratioN (UNICORN), School of Pharmaceutical Sciences, University of Phayao, Phayao, Thailand

  • 9. Novel Bacteria and Drug Discovery Research Group, Microbiome and Bioresource Research Strength, Jeffrey Cheah School of Medicine and Health Sciences, Monash University Malaysia, Bander Sunway, Malaysia

  • 10. Biofunctional Molecule Exploratory Research Group, Biomedicine Research Advancement Centre, School of Pharmacy, Monash University Malaysia, Bander Sunway, Malaysia

  • 11. Division of Ambulatory Medicine, Department of Medicine, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok, Thailand

  • 12. Division of Pharmaceutical Care, Department of Pharmacy, Phrae Hospital, Phrae, Thailand

Article metrics

View details

4

Citations

2,7k

Views

1,2k

Downloads

Abstract

Background: Mortality from multidrug-resistant (MDR) pathogens is an urgent healthcare crisis worldwide. At present we do not have any simplified screening tools to predict the risk of mortality associated with MDR infections. The aim of this study was to develop a screening tool to predict mortality in patients with multidrug-resistant organisms.

Methods: A retrospective cohort study to evaluate mortality risks in patients with MDR infections was conducted at Phrae Hospital. Univariable and multivariable analyses were used to classify possible risk factors. The model performance was internally validated utilizing the mean of three measures of discrimination corrected by the optimism using a 1000-bootstrap procedure. The coefficients were transformed into item scores by dividing each coefficient with the lowest coefficient and then rounding to the most adjacent number. The area under the receiver operating characteristic curve (AuROC) was used to determine the performance of the model.

Results: Between 1 October 2018 and 30 September 2020, a total of 504 patients with MDR infections were enrolled. The ICU-SEPSA score composed of eight clinical risk factors: 1) immunocompromised host, 2) chronic obstructive pulmonary disease, 3) urinary tract infection, 4) sepsis, 5) placement of endotracheal tube, 6) pneumonia, 7) septic shock, and 8) use of antibiotics within the past 3 months. The model showed good calibration (Hosmer-Lemeshow χ2 = 19.27; p-value = 0.50) and good discrimination after optimism correction (AuROC 84.6%, 95% confidence interval [Cl]: 81.0%–88.0%). The positive likelihood ratio of low risk (score ≤ 5) and high risk (score ≥ 8) were 2.07 (95% CI: 1.74–2.46) and 12.35 (95% CI: 4.90–31.13), respectively.

Conclusion: A simplified predictive scoring tool wad developed to predict mortality in patients with MDR infections. Due to a single-study design of this study, external validation of the results before applying in other clinical practice settings is warranted.

Introduction

The burden of infections caused by multidrug-resistant organisms (MDROs) is increasing worldwide, while the available antibiotics become less effective over time. (Morrison and Zembower, 2020) Fifteen new antibiotics have been developed and launched in the market in a recent year, but only two of them can be effectively used to treat multidrug resistant gram-negative bacterial infections. (Carlet et al., 2012) Consequently, bacterial antimicrobial resistance is a critical and urgent problem that needs to be addressed globally. In 2019, it was estimated that the highest rate of antimicrobial-resistant burden was found in sub-Saharan Africa. Approximately 27.3 deaths per 100,000 population had been reported in the region. Other affected areas with high incidences included South Asia, Europe, Latin America, and Oceania. Most common organisms leading to deaths associated with resistance were 1) Pseudomonas aeruginosa, 2) Acinetobacter baumannii, 3) Escherichia coli, 4) Staphylococcus aureus, 5) Streptococcus pneumoniae, and 6) Klebsiella pneumoniae. Predominantly, respiratory infections, bloodstream infections, and intra-abdominal infections were among the most common infectious diseases associated with antimicrobial-resistance. (Murray et al., 2022) According to collective evidence from 1,023 hospitals in Thailand in 2010, almost one-third of hospitalized patients diagnosed with infectious diseases had multidrug-resistant (MDR) infections. The previous report revealed that A. baumannii (57.6%) and P. aeruginosa (20%–30%) were the leading organisms of MDR infections, respectively. (Anudit et al., 2016) Additionally, infections with gram-negative bacteria (i.e., Enterobacteriaceae) are crucial challenges in any health-care settings. It is widely recognized that they have become difficult to treat and often require a complex treatment plan. Despite optimal and aggressive therapy, they frequently lead to increased morbidity and mortality, resulted in heightening anxiety worldwide. (Kanchanasurakit et al., 2020).

MDR infections are associated with many problems including prolonged length of hospital stays, increased cost of treatment and mortality rates. (Cosgrove, 2006; Sunenshine et al., 2007; Tanwar et al., 2014; Lim et al., 2016; Founou et al., 2017) The previous studies have shown a relationship between higher rates of MDROs and mortality. (Lim et al., 2016; Mave et al., 2017; Teerawattanapong et al., 2018) The actual causes for this association are not well defined. Comorbidities and altered immune response of patients may also contribute to the increased risk of mortality. Nevertheless, other factors including the patient’s clinical status and source or severity of infection may have played a pivotal role. (Gandra et al., 2019) At present, we do not have any tools to predict the risk of mortality in these patients.

To the extent of our knowledge, we believe that a tool used for prediction of mortality in patients with MDR infections has not been established. The previous study has identified a risk scoring system to predict risk of acquiring MDR infection, however mortality outcome was not addressed. (Tseng et al., 2017) While there were tools to predict mortality in patients diagnosed with bacteremia, those tools did not specifically focus on patients with MDROs. (Shapiro et al., 2003, 2007; Taniguchi et al., 2017) The objective of this study was to develop a risk scoring model to predict mortality for patients with MDR infections.

Materials and methods

Setting

This study was carried out at Phrae Hospital, a 500-bed medical center, located in the northern part of Thailand.This work followed the Transparent Reporting of a multivariable prediction model for Individual Prognosis Or Diagnosis (TRIPOD) Statement (Supplementary Table S1). (Moons et al., 2015).

Study design and study population

We conducted a single-center, retrospective cohort study of hospitalized adults (aged ≥ 18 years) with at least one positive culture of MDROs during the admission at Phrae Hospital between 1 October 2018, and 30 September 2020. A total of 504 patients were enrolled in the study. For individuals with multiple positive culture results within the same admission, we only used the data from their first result for analysis. Patients with negative culture results as well as those with positive culture results with mono resistance or susceptible to antibiotics were excluded. In addition, patients with missing data in the medical record were excluded.

Multidrug resistant bacteria was defined as acquired non-susceptibility to at least one agent in three or more antimicrobial categories. (Falagas and Karageorgopoulos, 2008; Magiorakos et al., 2012) Antibiotics are classified into one of the following classes: aminoglycosides, penicillins, penicillins plus beta-lactamase inhibitors, cephalosporins, fluoroquinolones, tetracyclines, macrolides, carbapenems, folate pathway inhibitors, glycopeptides, polymyxins, and lincosamides. (Magiorakos et al., 2012) In this study, the infectious disease physician diagnosed and registered patients with MDR infections depending on susceptibility tests and clinical manifestations.

Data collection and definition

Demographic characteristics including age, sex, status of patients, receipt of intensive care unit (ICU) admission, duration of hospitalization, comorbidities, source of infection, risk factors for MDR, disease severity, and use of medical devices were collected. Comorbidities were identified according to the International Classification of Diseases, 10th Revision, Clinical Modification (ICD-10) code version 2016. These were: essential hypertension (I10), diabetes mellitus type 2 (E11.0-E11.9), dyslipidemia (E78.5), heart failure (I50.0-I50.9), chronic kidney disease (N18.3-N18.9), atrial fibrillation (I48.0-I48.9), chronic obstructive pulmonary disease (J44.9), gout (M10.9), liver cirrhosis (K74.6), coronary artery disease (I20.0, I21.3-I21.4, I25.2), epilepsy (G40.9), Alzheimer’s disease (G30.9), Parkinson’s disease (G20), cerebrovascular disease including stroke (I63.0-I63.9), depression (F32.9), anxiety (F41.1-41.2, 41.9), and schizophrenia (F20.9).

Immunocompromised hosts were defined as patients with active cancer, current use of chemotherapy or immunosuppressant drugs, Acquired Immune Deficiency Syndrome (AIDS), chronic alcohol use, dialysis including hemodialysis and continuous ambulatory peritoneal dialysis, and uncontrolled diabetes mellitus with HbA1C >7%. (Delves and Roitt, 2000) Sepsis and septic shock were identified according to the Sepsis-3 criteria. (Singer et al., 2016) Sepsis was defined as life-threatening organ dysfunction caused by a dysregulated host response to infection. For clinical operationalization, organ dysfunction can increase the Sequential Organ Failure Assessment (SOFA) score of 2 points or more. Septic shock can be clinically identified if patients were in need of vasopressor to maintain a mean arterial pressure of 65 mmHg or greater and serum lactate level greater than 2 mmol/L (>18 mg/dl) in the absence of hypovolemia (Kanchanasurakit et al., 2021). In this study, mortality was defined as an in-hospital mortality secondary to infections. The cause of death was confirmed by physicians specialized in infectious disease.

Statistical analysis

Baseline characteristics are reported using frequency (percentage), mean (SD), or median (interquartile range) as appropriate. To analyze mortality associated with MDROs-infected patients, we conducted univariable analysis using chi-square or Fisher’s exact test. Odds ratios (ORs) and 95% confidence intervals (95%CI) were evaluated for each factor. Variables with p-value ≤ 0.1 were then selected into the multivariable model using logistic regression. Any risk factors with p-value < 0.05 were integrated in the process of developing the risk prediction model for mortality. Additionally, Akaike information criterion (AIC) and Bayesian information criteria (BIC) were used to help identify variables that would be incorporated in the predictive model. The model with the lowest AIC and BIC values was chosen. (Chakrabarti and Ghosh, 2011) Model performance in the original sample was evaluated using the C-statistic, Hosmer-Lemeshow statistic, calibration slope, and Brier score (Huang et al., 2020).

Sample size

We calculated sample size by applying 20 outcome events per predictor variable (EPV). (Austin and Steyerberg, 2017) Based on the previous study that used five binary variables (Zhou et al., 2019), at least 100 patients with an outcome are required.

Model internal validation

Internal validation was performed using a bootstrap procedure with 1,000 bootstrapped samples. We calculated the bootstrap performance (i.e., calibration slope, calibration-in-the-large [CITL] and C-statistic) of the final model to derive an optimism adjusting factor to correct the model for overfitting. To develop a scoring system, each variable was assessed and turned into a point using their weighted coefficients. These were then rounded to the most adjacent number to be used as a score.

Score development

To derive the point-based system, optimism adjusted coefficient and the intercept from the final multivariable model were used. We weighted and transformed each coefficient by dividing it by the lowest coefficient in the model and rounding the results to the most adjacent number. (Harrell, 1996; Moons et al., 2002) To examine the discriminative power of the predictive score we developed, a receiver operating characteristic (ROC) curve was plotted and a Hosmer-Lemeshow goodness-of-fit test was performed. (Steyerberg et al., 2010) The cut-off score was chosen to classify patients into three groups based on the risk of mortality: low (≤ 49.54%), moderate (56.98%–81.42%), and high (≥ 85.53%), respectively. The sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), accuracy, positive likelihood ratio (LR+), and negative likelihood ratio (LR-) were assessed. (Leisenring and Pepe, 1998; Janssens et al., 2005) The multi-collinearity was obtained using the variance inflation factor (VIF).

Score validation

The score was internally validated using bootstrapping. A total of 1,000 bootstrap samples were drawn and the final model containing all included variables was fitted in each bootstrap sample. This process was followed by re-simplification of the score and point assignment using the logistic regression model in each bootstrap sample. Performance of the re-fitted, simplified score was appraised in each bootstrap sample and in the original sample using the performance measures described above. Mean bootstrap and bootstrap test performances were assessed. Subsequently, optimism and optimism-corrected performance estimates were calculated.

Results

Characteristics of patients

From, 1 October 2018, through 30 September 2020, a total of 4,217 patients with culture results were identified during the study period. After excluding 3,713 patients with negative culture results, positive culture results with mono-resistance to antibiotics, missing data, or susceptible of antibiotics, 504 patients were included in the final analysis. There were 172 deaths (35%) during the follow-up (Figure 1). The mean age of the patients were 65 years, and the majority of patients (58%) were male. The mean duration of hospitalization was 18.39 days. Compared between the survival and non-survivor groups, there was no difference in bedridden status and comorbid conditions, except for heart failure, atrial fibrillation, and chronic obstructive pulmonary disease. Patients that died had recent admission at other facilities and later transferred to Phrae Hospital, admitted in the ICU, known immunocompromised status, and had recent antibiotic use with the past 3 months. Both groups had similar reports of wound infection. The survival group had less severe infections as well as fewer uses of medical devices. Baseline characteristics of both groups are shown in Table 1.

FIGURE 1

FIGURE 1

Enrollment of study patients.

TABLE 1

CharacteristicsSurvival (n = 332)Death (n = 172)Total (n = 504)p-value
Age (years), mean (±SD)64.6 (15.1)66.5 (14.5)65.2 (15.0)0.56
Male sex, No (%)187 (56.3)106 (61.6)293 (58.13)0.25
Duration of hospitalization, median (interquartile range), days9 (5–20.5)15 (7–28.5)11 (5–23)<0.001
Status of patients, No (%)
 ICU admission75 (22.59)57 (33.14)132 (26.19)0.011
 Referral from other hospital136 (40.96)96 (55.81)232 (46.03)<0.001
 Bedridden patients42 (12.65)28 (16.28)70 (13.89)0.26
 Immunocompromised host115 (34.64)83 (48.26)198 (39.29)<0.001
Inter-current conditions, No (%)
 Hypertension161 (48.49)89 (51.74)250 (49.60)0.49
 Type 2 Diabetes mellitus89 (26.81)47 (27.33)136 (26.98)0.90
 Dyslipidemia68 (20.48)29 (16.86)97 (19.25)0.33
 Heart failure7 (2.11)12 (6.98)19 (3.77)<0.001
 Chronic kidney disease88 (26.51)50 (29.07)138 (27.38)0.54
 Atrial fibrillation12 (3.61)19 (11.05)31 (6.15)0.001
 Chronic obstructive pulmonary disease20 (6.02)26 (15.12)46 (9.13)0.001
 Gout33 (9.94)10 (5.81)43 (8.53)0.13
 Liver cirrhosis9 (2.71)7 (4.07)16 (3.17)0.41
 Coronary artery disease20 (6.02)9 (5.23)29 (5.75)0.72
 Neurological disease7 (2.11)7 (4.07)14 (2.78)0.20
 Cerebrovascular disease30 (9.04)14 (8.14)44 (8.73)0.74
 Mental disorder6 (5.3)2 (2.7)8 (1.59)0.58
Source of infection, No (%)
 Bacteremia241 (72.59)79 (45.93)320 (63.49)<0.001
 Urinary tract48 (14.46)46 (26.74)94 (18.65)0.001
 Pneumonia14 (4.22)41 (23.84)55 (10.91)<0.001
 Intra-abdominal0 (0)5 (2.91)5 (0.99)0.002
 Wound29 (8.73)13 (7.56)42 (8.33)0.65
Treatment history, No (%)
 Hospitalization within the past 3 months142 (42.77)86 (50.00)228 (45.24)0.12
 Surgery within the past 3 months26 (7.83)9 (5.23)35 (6.94)0.28
 Use of antibiotics within the past 3 months120 (36.14)86 (50.00)206 (40.87)0.003
Severity of infection, No (%)
 Sepsis130 (39.16)113 (65.70)243 (48.21)<0.001
 Septic shock42 (12.65)68 (39.53)110 (21.83)<0.001
Medical devices, No (%)
 Foley catheter148 (44.58)133 (77.33)281 (55.75)<0.001
 Endotracheal tube76 (22.89)116 (67.44)192 (38.10)<0.001

Baseline characteristics of patients with MDROs.

Abbreviations: ICU, intensive care unit; SD, standard deviation.

Predictors for mortality in patients with multidrug-resistant organisms

Univariate analysis for mortality showed a significant association with hospitalization more than 14 days (OR 1.88, 95% CI: 1.27–2.78, p-value = 0.001), admission at ICU (OR 1.70, 95% CI: 1.10–2.61 p-value = 0.014), immunocompromised host (OR 1.76, 95% CI: 1.19–2.60, p-value = 0.004), heart failure (OR 3.48, 95% CI: 1.23–10.62, p-value = 0.012), atrial fibrillation (OR 3.31, 95% CI: 1.48–7.67, p-value = 0.002), chronic obstructive pulmonary disease (OR 2.78, 95% CI: 1.44–5.42, p-value = 0.002), bacteremia (OR 0.32, 95% CI: 0.21–0.48, p-value< 0.001), urinary tract infection (OR 2.16, 95% CI: 1.33–3.49, p-value = 0.001), pneumonia (OR 7.11, 95% CI: 3.63–14.56, p-value< 0.001), use of antibiotics within the last 3 months (OR 1.77, 95% CI: 1.19–2.61, p-value = 0.003), sepsis (OR 2.98, 95% CI: 1.99–4.46, p-value< 0.001), septic shock (OR 4.51, 95% CI: 2.83–7.23, p-value< 0.001), placement of foley catheter (OR 4.24, 95% CI: 2.75–6.61, p-value< 0.001), and placement of endotracheal tube (OR 6.98, 95% CI: 4.54–10.73, p-value< 0.001) (Table 2).

TABLE 2

FactorsTotal (n = 504)
Odds ratio95% CIp-value
Male1.250.84–1.850.29
Aged ≥ 65 years1.120.76–1.670.568
Hospitalization > 14 days1.881.27–2.780.001*
Admission in ICU1.701.10–2.610.014*
Bedridden1.340.77–2.320.279
Immunocompromised host1.761.19–2.600.004*
Hypertension1.140.77–1.670.512
Type 2 Diabetes mellitus1.030.66–1.580.916
Dyslipidemia0.790.47–1.300.343
Heart failure3.481.23–10.620.012*
Chronic kidney disease1.140.74–1.740.600
Atrial fibrillation3.311.48–7.670.002*
Chronic obstructive pulmonary disease2.781.44–5.420.002*
Gout0.560.24–1.200.132
Liver cirrhosis1.520.47–4.680.429
Coronary artery disease0.860.34–2.030.841
Neurological disease1.970.58–6.690.254
Cerebrovascular disease0.890.42–1.790.868
Mental disorder0.640.06–3.620.722
Bacteremia0.320.21–0.48<0.001*
Urinary tract2.161.33–3.490.001*
Pneumonia7.113.63–14.56<0.001*
Wound0.850.40–1.750.735
Hospitalization within 3 months1.340.91–1.970.132
Surgery within 3 months0.650.26–1.470.356
Use of antibiotics within 3 months1.771.19–2.610.003*
Sepsis2.981.99–4.46<0.001*
Septic shock4.512.83–7.23<0.001*
Foley catheter4.242.75–6.61<0.001*
Endotracheal tube6.984.54–10.73<0.001*

| Association between factors and mortality using univariate analysis.

*Factor was significant. Abbreviations: Cl, confidence interval; ICU, intensive care unit.

Model development

In the multivariable analysis, some variables were independent predictors of mortality. These were chronic obstructive pulmonary disease (OR 3.59, 95% CI:1.61–7.99, p-value = 0.002), pneumonia (4.86, 95% CI: 2.29–10.27, p-value<0.001), urinary tract infection (OR 3.24, 95% CI: 1.85–5.66, p-value<0.001), use of antibiotics within the past 3 months (OR 1.89, 95% CI: 1.18–3.01, p-value = 0.008), immunocompromised host (OR 2.24, 95% CI: 1.39–3.62, p-value = 0.001), sepsis (OR 2.24, 95% CI: 1.31–3.82, p-value = 0.003), septic shock (OR 2.02, 95% CI: 1.09–3.74, p-value = 0.026), and placement of endotracheal tube (OR 5.30, 95% CI: 3.27–8.59, p-value<0.001) (Table 3 and Figure 2). Since the VIF values of each variable were less than 5 (averaged value was 1.17), multi-collinearity was not a major identified in our study. We have analyzed the AIC and BIC scores to help identify variables that would be incorporated in the predictive model. The model with the lowest AIC and BIC values was chosen as the best model for predicting mortality in patients with MDR infection (Supplementary Table S2).

TABLE 3

PredictorsCoefficientp-valueTransformed coefficientsAssigned score
Immunocompromised host
 No0
 Yes0.7601520.0011.271.5
Chronic obstructive pulmonary disease
 No0
 Yes1.2042860.0022.012
Urinary tract infection
 No0
 Yes1.108349<0.0011.852
Sepsis
 No0
 Yes0.7585340.0031.271.5
On the Endotracheal tube
 No0
 Yes1.571675<0.0012.632.5
Pneumonia
 No0
 Yes1.489380<0.0012.492.5
Septic shock
 No0
 Yes0.6616890.0261.111
Use of Antibiotics within 3 months
 No0
 Yes0.5984580.00811

| Significant factors of mortality in patients with MDROs using multivariable analysis (ICU-SEPSA score).

The coefficient of model intercept is -3.010593 (p-value <0.001).

FIGURE 2

FIGURE 2

Visualize used for predicting mortality in patients with MDR infection.

This scoring system was created by setting the cut-off point according to the differentiation plot and performance of diagnostic parameters to identify patients at risk of increased mortality. Cut-off scores of 5 and 8 were chosen to classify patient into 3 groups (Table 4). Patients with a score ≥ 8 were placed in the high-risk group. Patients in this group were correctly predicted for the mortality outcome 86.49% (32/37) of the time. Patients with a score ≤ 5 were placed in the low-risk group and correctly predicted for the mortality outcome 79.95% (299/374) of the time. Correct prediction of survival or death occurred 80.54% [(299 + 32)/(374 + 37)] of the time. The incorrect prediction rate was 19.46% [(75 + 5)/(374 + 37)]. Based on this scoring system and one cut-off point, the original performance model showed good discrimination (AUC = 84.60%, 95%Cl: 81%–88%; Hosmer-Lemeshow χ2 = 19.27; p-value = 0.50) and calibration (calibration slope, 1.000; Brier score, 0.149; product moment correlation between observed and predicted probability, 1.000). Internally validated performance that gave an AUC similar original performance model and calibration (calibration slope, 1.061; Brier score, 0.149; product moment correlation between observed and predicted probability, 0.983) (Figures 3, 4A,B). (Marcin and Romano, 2007) The probability of mortality is presented in Figure 5.

TABLE 4

ParametersLow risk (score ≤ 5)Moderate risk (score 5.5–7.5)High risk (score ≥ 8)Total (n = 504)
Total3749337504
 Survival group299285332
 Death group756532172
Diagnostic performance (95% Confidence Interval)
 Sensitivity (%)90.06 (86.32–93.06)18.60 (13.09–25.24)
 Specificity (%)56.40 (48.64–63.93)98.49 (96.52–99.51)
 Positive predictive value (%)79.95 (77.02–82.59)86.49 (71.75–94.16)
 Negative predictive value (%)74.62 (67.45–80.65)70.02 (68.47–71.52)
 Accuracy78.57 (74.73–82.08)71.23 (67.06–75.15)
 Likelihood ratio (+)2.07 (1.74–2.46)12.35 (4.90–31.13)
 Likelihood ratio (-)0.18 (0.12–0.25)0.83 (0.77–0.89)
 InterpretationSurvive (74.62% certainty)Death (86.49% certainty)

Distribution of the risk of mortality in patients with MDROs and diagnostic performance.

FIGURE 3

FIGURE 3

Receiver operating characteristic (ROC) curve of risk scoring system.

FIGURE 4

FIGURE 4

The calibration curve of original performance expressed in (A); Internally validated performance expressed in (B). EO: expected-observed ratio; CITL: calibration-in-the-large; AUC: area under the receiver operating characteristic curve; CI: confidence interval.

FIGURE 5

FIGURE 5

Probability of mortality in patients with MDROs stratified by the score.

Discussion

Due to lack of predictive scoring for mortality in patients with multidrug resistant infection, we aimed to develop a risk scoring model to predict mortality in this population. Our study revealed that a clinical scoring system could be used to screen and predict the mortality in patients with MDROs. Our newly formulated tool, named “ICU-SEPSA score”, was able to classify patients at risk of mortality.

During our study, we collected data of 504 patients from October 2018 to September 2020. Using newly identified eight factors, we were able to categorize patients into three groups in relation to mortality risk: low, moderate, and high. These factors were 1) Immunocompromised host, 2) chronic obstructive pulmonary disease, 3) urinary tract infection, 4) sepsis, 5) placement of endotracheal tube, 6) pneumonia, 7) septic shock, and 8) use of antibiotics within the past 3 months. We observed similar findings compared to previous studies. (Miravitlles et al., 1999; Johnson et al., 2011; Lee et al., 2016; Moreau et al., 2018; Capsoni et al., 2019; Cillóniz et al., 2019; Du et al., 2019) While the exact relationship between immunocompromised hosts and mortality is unknown, the study by Moreau et al. (2018) found that immunocompromised hosts were linked to recent use of antibiotics. This is a well-known risk factor for developing MDROs and increased risk of mortality. In our study, we found that sepsis and septic shock were associated with higher mortality rates due to multiple organ dysfunction. (Capsoni et al., 2019) Use of any medical devices that posed an increased risk of bacteria growth either through lower bacteria clearance or epithelial damage that allowed formation of biofilm, can result in higher risk of mortality. In our study, we found that endotracheal tube was associated with mortality in MDROs infections. (Nseir et al., 2006) In addition, a previous systematic review and meta-analysis showed that intubation was considered as an invasive techniques associated with adverse outcomes. (Du et al., 2019) It has been established that some comorbidities are associated with mortality, particularly chronic obstructive pulmonary disease, but the mechanism is unclear. A previous study suggested that patients with acute exacerbation of chronic obstructive pulmonary disease required intubation or ventilators more often. (Nseir et al., 2006) The sputum culture of patients with severe pulmonary function impairment, indicated by forced expiratory volume in one second (FEV-1) <50%, has shown to contain more potentially pathogenic microorganisms, such as Haemophilus influenzae or P. aeruginosa. (Miravitlles et al., 1999) In our study, we found that certain types of infections (urinary tract infection and pneumonia) led to increased risks of mortality. This may be the result of receiving insufficient empiric treatment leading to the development of MDROs infections. (Lee et al., 2016; Cillóniz et al., 2019) The last factor correlated to mortality in patients with MDR pathogens was the use of antibiotics within the past 3 months. We hypothesized that prior antibiotic use led to higher likelihood of antimicrobial resistance. (Johnson et al., 2011)

Our study has a number of strengths. To the extent of our knowledge, this was the first time that risk scoring tool to predict mortality secondary to any MDROs infections was identified. Our data analysis was not restricted to a particular type of infections. We evaluated various sources of cultures, including blood, urine, sputum, catheter, ascitic fluid and pus. Our tool was user friendly. It was simple and did not require any specific or expensive laboratory investigation. The newly established tool, named “ICU-SEPSA score” consisted of eight predictors, is easy to remember for clinicians. We foresee a great impact of this tool to help assess mortality risks in MDRO-infected patients. Lastly, the power of this model was sufficient to detect a difference between the two groups.

There are several limitations to this study. First, the data was collected retrospectively from medical records and electronic database. There was a potential absence of information that could impact data analysis and outcome such as type of organism, smoking status, history of alcohol use, and pertinent past medical history. In addition, during medical record review, we were unable to assess the severity of some variables such as urinary tract infection and pneumonia. Third, this study was designed as a single-center study, therefore, the ability to generalize our results may be limited. Lastly, it is worth noting that our model has not been conducted in the external validation cohort.

Our newly formulated tool may be incorporated into clinical practice to help healthcare providers estimate the risk of mortality in this population. A high score corresponds to a greater risk of mortality in patients with MDROs infections that should assist clinicians to promptly identify the appropriate therapy. To implement our findings in clinical setting, we suggested that patients who classified into the high-risk group (score ≥ 8) should be considered for admission in the ICU for close monitoring. These patients should receive combination antibiotic therapy for synergistic effect and reduced risk of modified risk factors. Frequent vital sign monitoring is warranted for moderate-risk group (score 5.5–7.5) as well as the high-risk group. At Phrae Hospital, vital sign monitoring as frequently as every 2 h had been suggested. In addition, every effort to minimize prolonged use of any catheters should be considered. Patients in the low-risk group (score ≤ 5) should still have vital signs monitored regularly according to the standard of care at that clinical setting. Regardless of which group patients are categorized into, MDROs infections should be immediately assessed. Aggressive therapy with close monitoring is key to lower risks of mortality and optimize overall outcome.

Conclusion

A simplified predictive scoring tool wad developed to predict mortality in patients with MDROs infections. The score was composed of eight factors including 1) immunocompromised host, 2) chronic obstructive pulmonary disease, 3) urinary tract infection, 4) sepsis, 5) placement of endotracheal tube, 6) pneumonia, 7) septic shock, and 8) use of antibiotics within the past 3 months. Due to a single-study design of this study, external validation of the results before applying in other clinical practice settings is warranted.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by The Ethical Committee for Clinical Research in Phrae Hospital (22/2563). The patients/participants provided their written informed consent to participate in this study.

Author contributions

Study concept and design (SK, TS, SS, PP, AP, and CS); acquisition of data (SK, CS, and CT); analysis and interpretation of data (SK, CS); drafting of the manuscript (SK, TS, CS, and CT); critical revision of the manuscript (SK, NH, SS, AP, and PP).

Funding

This work was supported by a grant from the Unit of Excellence on Clinical Outcomes Research and IntegratioN (UNICORN) [Grant number: FF65-UoE005], School of Pharmaceutical Sciences, University of Phayao. The funding source had no role in the study design, collection, analysis, and interpretation of data.

Acknowledgments

We deeply thank the officials of Phrae Hospital for their support in data collection. We thank Dillon Dickow (The University of Illinois at Chicago College of Pharmacy) for proofreading and assistance with English language review of the manuscript prior to submission.

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.

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/fphar.2022.938028/full#supplementary-material

Abbreviations

AuROC, The area under the receiver operating characteristic curve; CITL, calibration-in-the-large; EPV, events per predictor variable; FEV1, forced expiratory volume in one second; ICD-10, the International Classification of Diseases, 10th Revision, Clinical Modification; ICU, intensive care unit; LR+, positive likelihood ratio; LR-, negative likelihood ratio; MDR, multidrug-resistant; MDROs, multidrug-resistant organisms; NPV, negative predictive value; ORs, odds ratios; PPV, positive predictive value; ROC, receiver operating characteristic; TRIPOD, the Transparent Reporting of a multivariable prediction model for Individual Prognosis Or Diagnosis; 95%CI, 95% confidence intervals.

References

  • 1

    AnuditC.KooltheatN.PotupP.Pankla SranujitR.UsuwanthimK. (2016). Nosocomial infection of multidrug-resistant Acinetobacter baumannii in Thailand. Am. J. Infect. Control44, 11611163. 10.1016/j.ajic.2016.03.068

  • 2

    AustinP. C.SteyerbergE. W. (2017). Events per variable (EPV) and the relative performance of different strategies for estimating the out-of-sample validity of logistic regression models. Stat. Methods Med. Res.26, 796808. 10.1177/0962280214558972

  • 3

    CapsoniN.BelloneP.AlibertiS.SotgiuG.PavanelloD.VisintinB.et al (2019). Prevalence, risk factors and outcomes of patients coming from the community with sepsis due to multidrug resistant bacteria. Multidiscip. Respir. Med.14, 23. 10.1186/s40248-019-0185-4

  • 4

    CarletJ.JarlierV.HarbarthS.VossA.GoossensH.PittetD.et al (2012). Ready for a world without antibiotics? The pensières antibiotic resistance call to action. Antimicrob. Resist. Infect. Control1, 11. 10.1186/2047-2994-1-11

  • 5

    ChakrabartiA.GhoshJ. K. (2011). “AIC, BIC and recent advances in model selection,” in Handbook of the philosophy of science. Editors BandyopadhyayP. S.ForsterM. R. B. T.-P. S. (Amsterdam: North-Holland), 583605. 10.1016/B978-0-444-51862-0.50018-6

  • 6

    CillónizC.DominedòC.TorresA. (2019). Multidrug resistant gram-negative bacteria in community-acquired pneumonia. Annu. Updat. Intensive Care Emerg. Med.2019, 459475. 10.1007/978-3-030-06067-1_36

  • 7

    CosgroveS. E. (2006). The relationship between antimicrobial resistance and patient outcomes: mortality, length of hospital stay, and health care costs. Clin. Infect. Dis.42, S82S89. 10.1086/499406

  • 8

    DelvesP. J.RoittI. M. (2000). The immune system. First of two parts. N. Engl. J. Med.343, 3749. 10.1056/NEJM200007063430107

  • 9

    DuX.XuX.YaoJ.DengK.ChenS.ShenZ.et al (2019). Predictors of mortality in patients infected with carbapenem-resistant acinetobacter baumannii: A systematic review and meta-analysis. Am. J. Infect. Control47, 11401145. 10.1016/j.ajic.2019.03.003

  • 10

    FalagasM. E.KarageorgopoulosD. E. (2008). Pandrug resistance (PDR), extensive drug resistance (XDR), and multidrug resistance (MDR) among gram-negative bacilli: need for international harmonization in terminology. Clin. Infect. Dis.46, 11211122. 10.1086/528867

  • 11

    FounouR. C.FounouL. L.EssackS. Y. (2017). Clinical and economic impact of antibiotic resistance in developing countries: A systematic review and meta-analysis. PLoS One12, e0189621. 10.1371/journal.pone.0189621

  • 12

    GandraS.TsengK. K.AroraA.BhowmikB.RobinsonM. L.PanigrahiB.et al (2019). The mortality burden of multidrug-resistant pathogens in India: A retrospective, observational study. Clin. Infect. Dis.69, 563570. 10.1093/cid/ciy955

  • 13

    HarrellF. (1996). Regression coefficients and scoring rules. J. Clin. Epidemiol.49, 819. 10.1016/0895-4356(95)00068-2

  • 14

    HuangY.LiW.MacheretF.GabrielR. A.Ohno-MachadoL. (2020). A tutorial on calibration measurements and calibration models for clinical prediction models. J. Am. Med. Inf. Assoc.27, 621633. 10.1093/jamia/ocz228

  • 15

    JanssensA. C. J. W.DengY.BorsboomG. J. J. M.EijkemansM. J. C.HabbemaJ. D. F.SteyerbergE. W. (2005). A new logistic regression approach for the evaluation of diagnostic test results. Med. Decis. Mak.25, 168177. 10.1177/0272989X05275154

  • 16

    JohnsonM. T.ReichleyR.Hoppe-BauerJ.DunneW. M.MicekS.KollefM. (2011). Impact of previous antibiotic therapy on outcome of Gram-negative severe sepsis. Crit. Care Med.39, 18591865. 10.1097/CCM.0b013e31821b85f4

  • 17

    KanchanasurakitS.SantimaleeworagunW.McPhersonC. E.PiriyachananusornN.BoonsongB.KatwilatP.et al (2020). Fosfomycin dosing regimens based on Monte Carlo simulation for treated carbapenem-resistant enterobacteriaceae infection. Infect. Chemother.5, 516529. 10.3947/ic.2020.52.4.516

  • 18

    KanchanasurakitS.SuthumpoungP.SantimaleeworagunW.NakaranurackC.HuynhN. S.SrisawatC.et al (2021). Effectiveness of thiamine therapy in mortality rate in patients with septic shock: A systematic review and meta-analysis. Int. J. Crit. Illn. Inj. Sci.11, 8694. 10.4103/IJCIIS.IJCIIS_159_20

  • 19

    LeeY.-C.HsiaoC.-Y.HungM.-C.HungS.-C.WangH.-P.HuangY.-J.et al (2016). Bacteremic urinary tract infection caused by multidrug-resistant enterobacteriaceae are associated with severe sepsis at admission: Implication for empirical therapy. Med. Baltim.95, e3694. 10.1097/MD.0000000000003694

  • 20

    LeisenringW.PepeM. S. (1998). Regression modelling of diagnostic likelihood ratios for the evaluation of medical diagnostic tests. Biometrics54, 444452. 10.2307/3109754

  • 21

    LimC.TakahashiE.HongsuwanM.WuthiekanunV.ThamlikitkulV.HinjoyS.et al (2016). Epidemiology and burden of multidrug-resistant bacterial infection in a developing country. Elife5, e18082. 10.7554/eLife.18082

  • 22

    MagiorakosA.-P.SrinivasanA.CareyR. B.CarmeliY.FalagasM. E.GiskeC. G.et al (2012). Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect.18, 268281. 10.1111/j.1469-0691.2011.03570.x

  • 23

    MarcinJ. P.RomanoP. S. (2007). Size matters to a model’s fit. Crit. Care Med.35, 22122213. 10.1097/01.CCM.0000281522.70992.EF

  • 24

    MaveV.ChandanwaleA.KagalA.KhadseS.KadamD.BharadwajR.et al (2017). High burden of antimicrobial resistance and mortality among adults and children with community-onset bacterial infections in India. J. Infect. Dis.215, 13121320. 10.1093/infdis/jix114

  • 25

    MiravitllesM.EspinosaC.Fernández-LasoE.MartosJ. A.MaldonadoJ. A.GallegoM. (1999). Relationship between bacterial flora in sputum and functional impairment in patients with acute exacerbations of COPD. Study Group of Bacterial Infection in COPD. Chest116, 4046. 10.1378/chest.116.1.40

  • 26

    MoonsK. G. M.AltmanD. G.ReitsmaJ. B.IoannidisJ. P. A.MacaskillP.SteyerbergE. W.et al (2015). Transparent reporting of a multivariable prediction model for individual Prognosis or Diagnosis (TRIPOD): explanation and elaboration. Ann. Intern. Med.162, W1W73. 10.7326/M14-0698

  • 27

    MoonsK. G. M.HarrellF. E.SteyerbergE. W. (2002). Should scoring rules be based on odds ratios or regression coefficients?J. Clin. Epidemiol.55, 10541055. 10.1016/s0895-4356(02)00453-5

  • 28

    MoreauA.-S.Martin-LoechesI.PovoaP.SalluhJ.RodriguezA.ThilleA. W.et al (2018). Impact of immunosuppression on incidence, aetiology and outcome of ventilator-associated lower respiratory tract infections. Eur. Respir. J.51, 1701656. 10.1183/13993003.01656-2017

  • 29

    MorrisonL.ZembowerT. R. (2020). Antimicrobial resistance. Gastrointest. Endosc. Clin. N. Am.30, 619635. 10.1016/j.giec.2020.06.004

  • 30

    MurrayC. J.IkutaK. S.ShararaF.SwetschinskiL.Robles AguilarG.GrayA.et al (2022). Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet399, 629655. 10.1016/S0140-6736(21)02724-0

  • 31

    NseirS.Di PompeoC.CavestriB.JozefowiczE.NyungaM.SoubrierS.et al (2006). Multiple-drug-resistant bacteria in patients with severe acute exacerbation of chronic obstructive pulmonary disease: Prevalence, risk factors, and outcome. Crit. Care Med.34, 29592966. 10.1097/01.CCM.0000245666.28867.C6

  • 32

    ShapiroN. I.HowellM. D.TalmorD.DonninoM.NgoL.BatesD. W. (2007). Mortality in emergency department sepsis (MEDS) score predicts 1-year mortality. Crit. Care Med.35, 192198. 10.1097/01.CCM.0000251508.12555.3E

  • 33

    ShapiroN. I.WolfeR. E.MooreR. B.SmithE.BurdickE.BatesD. W. (2003). Mortality in emergency department sepsis (MEDS) score: a prospectively derived and validated clinical prediction rule. Crit. Care Med.31, 670675. 10.1097/01.CCM.0000054867.01688.D1

  • 34

    SingerM.DeutschmanC. S.SeymourC. W.Shankar-HariM.AnnaneD.BauerM.et al (2016). The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA315, 801810. 10.1001/jama.2016.0287

  • 35

    SteyerbergE. W.VickersA. J.CookN. R.GerdsT.GonenM.ObuchowskiN.et al (2010). Assessing the performance of prediction models: a framework for traditional and novel measures. Epidemiology21, 128138. 10.1097/EDE.0b013e3181c30fb2

  • 36

    SunenshineR. H.WrightM.-O.MaragakisL. L.HarrisA. D.SongX.HebdenJ.et al (2007). Multidrug-resistant Acinetobacter infection mortality rate and length of hospitalization. Emerg. Infect. Dis.13, 97103. 10.3201/eid1301.060716

  • 37

    TaniguchiL. U.PiresE. M. C.VieiraJ. M. J.AzevedoL. C. P. de (2017). Systemic inflammatory response syndrome criteria and the prediction of hospital mortality in critically ill patients: a retrospective cohort study. Rev. Bras. Ter. intensiva29, 317324. 10.5935/0103-507X.20170047

  • 38

    TanwarJ.DasS.FatimaZ.HameedS. (2014). Multidrug resistance: an emerging crisis. Interdiscip. Perspect. Infect. Dis.2014, 541340. 10.1155/2014/541340

  • 39

    TeerawattanapongN.PanichP.KulpokinD.Na RanongS.KongpakwattanaK.SaksinanonA.et al (2018). A systematic review of the burden of multidrug-resistant healthcare-associated infections among intensive care unit patients in southeast Asia: The rise of multidrug-resistant acinetobacter baumannii. Infect. Control Hosp. Epidemiol.39, 525533. 10.1017/ice.2018.58

  • 40

    TsengW.-P.ChenY.-C.YangB.-J.ChenS.-Y.LinJ.-J.HuangY.-H.et al (2017). Predicting multidrug-resistant gram-negative bacterial colonization and associated infection on hospital admission. Infect. Control Hosp. Epidemiol.38, 12161225. 10.1017/ice.2017.178

  • 41

    ZhouH.YaoY.ZhuB.RenD.YangQ.FuY.et al (2019). Risk factors for acquisition and mortality of multidrug-resistant acinetobacter baumannii bacteremia: A retrospective study from a Chinese hospital. Med. Baltim.98, e14937. 10.1097/MD.0000000000014937

Summary

Keywords

mortality, multidrug-resistant infection, screening tool, antibiotic resistant, predictive scoring

Citation

Sirichayanugul T, Srisawat C, Thummakomut C, Prawang A, Huynh NS, Saokaew S, Phisalprapa P and Kanchanasurakit S (2022) Development and internal validation of simplified predictive scoring (ICU-SEPSA score) for mortality in patients with multidrug resistant infection. Front. Pharmacol. 13:938028. doi: 10.3389/fphar.2022.938028

Received

09 May 2022

Accepted

02 August 2022

Published

26 August 2022

Volume

13 - 2022

Edited by

Alvaro Francisco Lopes Sousa, University of São Paulo, Brazil

Reviewed by

Haoquan Huang, Department of Anesthesiology, Sun Yat-sen University, China

Zikria Saleem, University of Lahore, Pakistan

Godfrey Mutashambara Rwegerera, Sir Ketumile Masire Teaching Hospital, Botswana

Updates

Copyright

*Correspondence: Pochamana Phisalprapa, ; Sukrit Kanchanasurakit, ;

†ORCID: Taranee Sirichayanugul, orcid.org/0000-0002-3187-9623; Chansinee Srisawat, orcid.org/0000-0002-7166-6075; Chawin Thummakomut, orcid.org/0000-0001-8773-1639; Abhisit Prawang, orcid.org/0000-0002-5809-8886; Nina S Huynh, orcid.org/0000-0003-1120-4669; Surasak Saokaew, orcid.org/0000-0002-1382-0660; Pochamana Phisalprapa, orcid.org/0000-0003-1995-4405; Sukrit Kanchanasurakit, orcid.org/0000-0002-1268-2665

This article was submitted to Drugs Outcomes Research and Policies, a section of the journal Frontiers in Pharmacology

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics