Hostname: page-component-7b7567c8c-kwt9d Total loading time: 0 Render date: 2026-11-14T07:15:50.142Z Has data issue: false hasContentIssue false
  • English
  • Français

Association of COVID-19 coinfection with increased mortality among patients with Pseudomonas aeruginosa bloodstream infection in the Veterans Health Administration system

Published online by Cambridge University Press:  15 December 2023

Leila S. HojatOpen the ORCID record for Leila S. Hojat[Opens in a new window],
Brigid M. WilsonOpen the ORCID record for Brigid M. Wilson[Opens in a new window],
Federico Perez,
Maria F. MojicaOpen the ORCID record for Maria F. Mojica[Opens in a new window],
Mendel E. SingerOpen the ORCID record for Mendel E. Singer[Opens in a new window],
Robert A. Bonomo and
Lauren H. Epstein
Leila S. Hojat*
Affiliation:
Department of Medicine, Division of Infectious Diseases, Case Western Reserve University, Cleveland, OH, USADivision of Infectious Diseases and HIV Medicine, University Hospitals Cleveland Medical Center, Cleveland, OH, USA
Brigid M. Wilson
Affiliation:
Department of Medicine, Division of Infectious Diseases, Case Western Reserve University, Cleveland, OH, USAGeriatric Research Education and Clinical Center (GRECC), the VA Northeast Ohio Healthcare System, Cleveland, OH, USA
Federico Perez
Affiliation:
Department of Medicine, Division of Infectious Diseases, Case Western Reserve University, Cleveland, OH, USAGeriatric Research Education and Clinical Center (GRECC), the VA Northeast Ohio Healthcare System, Cleveland, OH, USACase Western Reserve University, Cleveland VAMC Center for Antimicrobial Resistance and Epidemiology (Case VA CARES), Cleveland, OH, USA
Maria F. Mojica
Affiliation:
Case Western Reserve University, Cleveland VAMC Center for Antimicrobial Resistance and Epidemiology (Case VA CARES), Cleveland, OH, USAResearch Service, VA Northeast Ohio Healthcare System, Cleveland, OH, USAGrupo de Resistencia Antimicrobiana y Epidemiología Hospitalaria, Universidad El Bosque, Bogotá, ColombiaDepartments of Pathology, Pharmacology, Molecular Biology and Microbiology, Biochemistry, and Proteomics and Bioinformatics, Case Western Reserve University School of Medicine, Cleveland, OH, USA
Mendel E. Singer
Affiliation:
Department of Population and Quantitative Health Sciences, Case Western Reserve University, Cleveland, OH, USA
Robert A. Bonomo
Affiliation:
Department of Medicine, Division of Infectious Diseases, Case Western Reserve University, Cleveland, OH, USAGeriatric Research Education and Clinical Center (GRECC), the VA Northeast Ohio Healthcare System, Cleveland, OH, USACase Western Reserve University, Cleveland VAMC Center for Antimicrobial Resistance and Epidemiology (Case VA CARES), Cleveland, OH, USAResearch Service, VA Northeast Ohio Healthcare System, Cleveland, OH, USADepartments of Pathology, Pharmacology, Molecular Biology and Microbiology, Biochemistry, and Proteomics and Bioinformatics, Case Western Reserve University School of Medicine, Cleveland, OH, USA
Lauren H. Epstein
Affiliation:
Infectious Diseases, US Department of Veterans Affairs Medical Center, Decatur, GA, USA
*
Corresponding author: Leila S. Hojat; Email: Leila.Hojat@UHhospitals.org

Abstract

Objective:

Pseudomonas aeruginosa bloodstream infection (PA-BSI) and COVID-19 are independently associated with high mortality. We sought to demonstrate the impact of COVID-19 coinfection on patients with PA-BSI.

Design:

Retrospective cohort study.

Setting:

Veterans Health Administration.

Patients:

Hospitalized patients with PA-BSI in pre-COVID-19 (January 2009 to December 2019) and COVID-19 (January 2020 to June 2022) periods. Patients in the COVID-19 period were further stratified by the presence or absence of concomitant COVID-19 infection.

Methods:

We characterized trends in resistance, treatment, and mortality over the study period. Multivariable logistic regression and modified Poisson analyses were used to determine the association between COVID-19 and mortality among patients with PA-BSI. Additional predictors included demographics, comorbidities, disease severity, antimicrobial susceptibility, and treatment.

Results:

A total of 6,714 patients with PA-BSI were identified. Throughout the study period, PA resistance rates decreased. Mortality decreased during the pre-COVID-19 period and increased during the COVID-19 period. Mortality was not significantly different between pre-COVID-19 (24.5%, 95% confidence interval [CI] 23.3–28.6) and COVID-19 period/COVID-negative (26.0%, 95% CI 23.5–28.6) patients, but it was significantly higher in COVID-19 period/COVID-positive patients (47.2%, 35.3–59.3). In the modified Poisson analysis, COVID-19 coinfection was associated with higher mortality (relative risk 1.44, 95% CI 1.01–2.06). Higher Charlson Comorbidity Index, higher modified Acute Physiology and Chronic Health Evaluation score, and no targeted PA-BSI treatment within 48 h were also predictors of higher mortality.

Conclusions:

Higher mortality was observed in patients with COVID-19 coinfection among patients with PA-BSI. Future studies should explore this relationship in other settings and investigate potential SARS-CoV-2 and PA synergy.

Information

Type
Original Article
Information
Antimicrobial Stewardship & Healthcare Epidemiology , Volume 3 , Issue 1 , 2023, e237
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2023. Published by Cambridge University Press on behalf of The Society for Healthcare Epidemiology of America

Introduction

Invasive infections caused by Pseudomonas aeruginosa (PA) are challenging to treat among hospitalized patients due to the organism’s high virulence and complex patterns of antimicrobial resistance, including decreased outer membrane permeability, efflux pump systems, and antibiotic-inactivating enzymes.Reference Rice1–Reference Breidenstein, de la Fuente-Núñez and Hancock3 The acquisition of resistance genes via horizontal transfer or spontaneous mutation can also induce additional mechanisms of resistance, either due to exposure to broad-spectrum antibiotics or nosocomial transmission.Reference Berrazeg, Jeannot and Ntsogo Enguéné4 Mortality associated with PA is high, particularly in the bloodstream, with an estimated 30-d mortality rate for PA bloodstream infection (PA-BSI) between 24% and 32%.Reference Gerver, Nsonwu, Thelwall, Brown and Hope5–Reference Lodise, Patel and Kwa7

The Centers for Disease Control and Prevention (CDC) published encouraging data in 2019 indicating a decrease in multidrug-resistant (MDR) PA from 46,000 estimated cases and 3,900 deaths in 2012 to 32,600 cases and 2,700 deaths in 2017.8 However in 2022, CDC reported a 32% increase in rates of MDR PA from 2019 to 2020, largely attributed to antibiotic overuse in patients with COVID-19.9 In addition to observed higher rates of MDR PA, several studies have identified a recent rise in PA-BSI incidence. A retrospective analysis of English National Health Service acute trusts between August 2020 and February 2021 observed an increase both in incidence of hospital-onset PA-BSI (4.9–6.2 cases per 100,000 bed-days) and in percentage of patients with PA-BSI and concomitant COVID-19 infection (0.7–34.1% among PA-BSI cases).Reference Sloot, Nsonwu and Chudasama10 A study performed within the Ascension health system additionally noted a 28% increase in hospital-onset PA-BSI and approximately five times higher rate of PA-BSI among patients with COVID-19 compared to patients without COVID-19.Reference Sturm, Saake, Roberts, Masoudi and Fakih11 Other published reports have demonstrated similar findings.Reference Giacobbe, Battaglini and Ball12,Reference Buetti, Tabah and Loiodice13 While these studies collectively suggest that the rise in observed PA-BSI coincides with the onset of the COVID-19 pandemic, outcomes of patients with PA-BSI and COVID-19 coinfection have not been evaluated.

In this study, we aimed to use the Veterans Health Administration (VHA) population to study the impact of concomitant COVID-19 infection on all-cause mortality in patients with PA-BSI. We hypothesized that patients with PA-BSI and COVID-19 coinfection would have a higher risk of mortality compared with patients with PA-BSI only. Further, we hypothesized that patients with PA-BSI diagnosed during the COVID-19 pandemic period would have higher mortality than those with PA-BSI prior to this period, regardless of COVID-19 coinfection status.

Materials and methods

Data source and population

Clinical data were extracted from the VHA Corporate Data Warehouse (CDW), which serves as a repository for patient data recorded in the electronic health record for all patients throughout the VHA system. The VHA is the largest integrated healthcare system in the United States and provides care to a heterogeneous population of over 9 million Veterans distributed among both urban and rural settings throughout the country.14,Reference Chang, Yoon and Esmaeili15 Patients ≥18 yr with at least one hospital admission to a VHA facility as of January 1, 2006 were eligible to be included in this analysis. We defined a case of PA-BSI as a positive blood culture collected from a patient from January 1, 2009 through June 30, 2022. We excluded eligible cases 1) among patients with unknown mortality status 30 d from the index culture; 2) without hospitalization; and 3) in which the blood culture was collected greater than 7 d before or 30 d after the date of hospital admission. Multiple blood cultures with PA collected within 30 d from the same patient were included as one case. If a patient had multiple PA-BSIs that met our inclusion criteria collected more than 30 d apart, we only included the first episode of PA-BSI; thus, each case represented a unique patient. We defined the pre-COVID-19 period as January 1, 2009 through December 31, 2019 and the COVID-19 period as January 1, 2020 through June 30, 2022. Cases during the COVID-19 period were identified as having COVID-19 coinfection if a positive test result for COVID-19 was documented during the associated hospitalization.

We collected clinical information for each case-patient including age, sex, race, ethnicity, Charlson Comorbidity Index, and modified Acute Physiology and Chronic Health Evaluation (mAPACHE) score.Reference Fortis, O’Shea and Beck16 mAPACHE score was based on clinical information collected within 48 h following the blood culture collection corresponding to the index PA-BSI occurrence for each case. Antimicrobial susceptibility testing results were obtained from laboratory reports within the CDW; individual laboratory susceptibility testing methods and standards were not available. We collected information regarding treatment against the identified PA isolate administered within 48 h of the index positive culture. Targeted treatment was defined as administration of an agent with in vitro activity against the identified PA isolate on at least one day between days 0 and 2 from the date of index positive blood culture collection.

We defined four resistance phenotype classifications: (1) difficult-to-treat resistant (DTR); (2) MDR; (3) other unclassified resistant (OUR); and (4) expected resistant (ER). DTR was defined as in vitro resistance or intermediate resistance to all tested agents within the beta-lactam and fluoroquinolone classes (and aztreonam if tested).Reference Kadri, Adjemian and Lai17 DTR determination required in vitro testing results for ≥1 expanded-spectrum cephalosporin, ≥1 carbapenem, and ≥1 fluoroquinolone. MDR was assigned for isolates not meeting criteria for DTR and with in vitro resistance or intermediate resistance to ≥1 drug in at least three categories including piperacillins, expanded-spectrum cephalosporins, fluoroquinolones, carbapenems, and aminoglycosides.18,Reference Magiorakos, Srinivasan and Carey19 OUR was any other pattern of in vitro resistance not otherwise classified as DTR or MDR. We classified cases as ER if there was no resistance detected to any tested agent within the piperacillin, expanded-spectrum cephalosporin, fluoroquinolone, carbapenem, aminoglycoside, monobactam, or polymyxin classes.

Analysis

We performed a descriptive analysis characterizing all PA-BSI cases by resistance phenotype, targeted antimicrobial treatment, and 30-d mortality rate by year. A multivariable logistic regression analysis was performed to characterize the independent association between the primary predictors of COVID-19 coinfection and COVID-19 period and the primary outcome of 30-d mortality, as defined by the odds ratio with 95% confidence interval (CI) estimates. Relative risk (RR) was calculated by modified Poisson regression using a sandwich estimator to model robust error variance.Reference Bauer, Puzniak and Yu20 Models both with and without COVID-19 coinfection as a predictor were estimated. The pre-COVID-19 period was additionally divided into early (2009–2016) and late (2017–2019) periods in the model for better characterization of possible changing mortality trends over the study period. The late-pre-COVID-19 period (2017–2019) was limited to 3 yr to create a more comparable cohort to the 3-yr COVID-19 period (2020–2022). For categorical variables that were missing at least 15% or more of the observations, the missing values were recoded to a new category, “missing data.” Data for variables with a smaller proportion of missing observations were completed with multiple imputation using the mice package (version 3.15.0) with 10 imputations. Continuous predictor variables were rescaled (e.g., age in decades) or considered as categorical terms to allow for more interpretable assessment of effect size. Sensitivity analyses were performed using a complete case analysis and uncategorized continuous data. Due to concern that the estimated effect of targeted antimicrobial treatment administered within 48 h could be biased by death prior to this timepoint, analyses were also performed with and without targeted treatment within 48 h and mAPACHE score before and after excluding patients with death within 48 h. All analyses were performed in R version 4.1.2.

Ethics

Data were obtained as part of a protocol approved by the Institutional Review Board of the Cleveland Veterans Affairs Medical Center to study MDR organisms in hospitalized patients.

Results

Characteristics of the study population

Among 7,804 PA-BSIs identified from January 1, 2009 to June 30, 2022, 6,714 PA-BSIs met our case definition (Figure 1). Among these, 5,484 (82%) were included from the pre-COVID-19 period, 1,158 (17%) were included from the COVID-19 period without COVID-19 coinfection, and 72 (1%) were included from the COVID-19 period with COVID-19 coinfection (Table 1). A majority of the case-patients were male, ≥70 yr old, white, and non-Hispanic, with multiple chronic comorbid medical conditions including uncomplicated diabetes mellitus (n = 2,892, 43.1%), cancer without metastases (n = 2,825, 42.1%), chronic pulmonary disease (n = 2,670, 39.8%), and chronic renal disease (n = 2,294, 34.2%).

Figure 1.

Flow diagram of inclusion and exclusion criteria resulting in 6,714 cases and further categorization into COVID-19 period and coinfection cohorts.

Table 1.

Demographic characteristics among patients with Pseudomonas aeruginosa bloodstream infection (n = 6,714) stratified by COVID-19 status and COVID-19 period

Note. IQR, interquartile range; DTR, difficult-to-treat resistant.

a Excludes 421 patients with missing race data.

b Excludes 278 patients with missing ethnicity data.

Antimicrobial resistance

Antimicrobial susceptibility information sufficient for classification into the defined resistance phenotypes was available for 5,537 PA-BSI cases (82%). Throughout the study period, antimicrobial resistance decreased, with a downward trend observed for both MDR and DTR phenotypes and a steady increase in the proportion of PA isolates with ER; the proportion of PA isolates with OUR remained relatively unchanged (Figure 2). Among PA-BSI cases with COVID-19 coinfection, no DTR isolates were identified, though the proportion of MDR resistance in patients with COVID-19 was higher in 2021 and 2022 than in patients without COVID-19 coinfection during the same years.

Figure 2.

Figure demonstrates that the proportion of PA-BSI isolates within the VHA system susceptible to all tested agents increased over the study period, with a similar increase in proportion of all cases in which targeted treatment was administered within 48 h. A corresponding decrease in MDR and DTR resistance was also observed over the study period. Among patients with COVID-19 in addition to PA-BSI relative to patients with PA-BSI only, more isolates had OUR in 2020–2021 and more MDR was present in 2021–2022, but no DTR was present. Bar chart excludes 1,177 patients with missing resistance phenotype data. Line plot excludes 94 patients with missing data for targeted treatment within 48 h. 2022 data represent January 1 through June 30. PA-BSI, Pseudomonas aeruginosa bloodstream infection; VHA, Veterans Health Administration; ER, expected resistant; OUR, otherwise unclassified resistant; MDR, multidrug-resistant; DTR, difficult-to-treat resistant.

Antimicrobial treatment data were available for 6,620 PA-BSI case-patients (99%). Throughout the study period, the proportion of case-patients who received targeted treatment for PA-BSI increased (Figure 2). The proportion of case-patients who received targeted treatment within 48 h was lowest among the DTR and MDR groups (20.6% and 53.4%, respectively) compared with the OUR and ER groups (78.4% and 89.8%, respectively).

COVID-19 period and coinfection

Overall, 30-d all-cause mortality decreased among case-patients during the pre-COVID-19 period; conversely, we observed an increase in 30-d mortality among all case-patients during the COVID-19 period (not stratified by COVID-19 coinfection status) (Figure 3). Case-patients identified in either period were similar in terms of age, comorbidities, and severity of illness (characterized by the mAPACHE score) (Figure 4a–c). Case-patients with COVID-19 coinfection in the COVID-19 period had a 30-d mortality rate of 47.2%, whereas case-patients without COVID-19 during the COVID-19 period and case-patients in the pre-COVID-19 period had similar, lower mortality rates of approximately 25% (Figure 4d).

Figures 3.

Figure indicates the proportion of PA-BSI cases in the VHA system meeting the primary end point of mortality at 30 d during each year of the study period. Mortality decreased across the pre-COVID-19 period to a low of 20.6% and increased during the COVID-19 period, peaking at 30.8% in 2022 through June 30. PA-BSI, Pseudomonas aeruginosa bloodstream infection; VHA, Veterans Health Administration.

Figure 4.

Distribution of clinical characteristics and mortality rate are shown stratified by COVID-19 period and COVID-19 coinfection status. Age (A), Charlson Comorbidity Index (B), and modified APACHE score (C) are noted to have a similar distribution across each cohort, despite that the mortality rate (D) is significantly higher in the COVID-19-positive/COVID-19 period patients.

The pre-COVID-19 periods were associated with lower mortality among all patients with PA-BSI in the unadjusted analyses; however, this association was no longer significant after adjusting for the presence of COVID-19 coinfection (Table 2). After adjusting for potential confounders, patients with PA-BSI and COVID-19 coinfection were associated with higher mortality compared to patients with PA-BSI without COVID-19 in both the logistic regression (odds ratio 2.15, 95% CI 1.26–3.68) and modified Poisson (RR 1.44, 95% CI 1.01–2.06) models (Table 2). Higher Charlson Comorbidity Index, higher mAPACHE score, and lack of targeted treatment for PA-BSI within 48 h were also associated with higher risk of mortality in both analyses (Figure 5).

Table 2.

Predictors of 30-d mortality among patients with Pseudomonas aeruginosa bloodstream infection (n = 6,714) with and without adjustment for COVID-19 coinfection

Note. OR, odds ratio; CI, confidence interval; RR, relative risk.

Figure 5.

Modified Poisson analysis results illustrate the positive association between COVID-19 coinfection and increased all-cause 30-d mortality risk in patients with PA-BSI in the VHA system during the study period. Additional predictors associated with 30-d mortality include Charlson Comorbidity Index, administration of targeted treatment within 48 h, and modified APACHE score.

In a sensitivity analysis, the DTR resistance phenotype was associated with higher mortality when the results were not adjusted for targeted antimicrobial treatment administered within 48 h (RR 1.32, 95% CI 1.10–1.58). Analyses including cases with no missing data, consideration of Charlson Comorbidity Index and mAPACHE as continuous variables, and exclusion of patients with death in less than 48 h did not substantially alter the relationship between the primary predictors and the primary outcome.

Discussion

In this study of hospitalized adult patients with PA-BSI identified within the VHA system, we observed approximately 50% greater adjusted risk of 30-d all-cause mortality in patients with PA-BSI and COVID-19 coinfection compared with patients with PA-BSI only. Increased mortality was additionally observed in patients with higher Charlson Comorbidity Index, higher mAPACHE score, and no targeted antibiotic treatment administered within 48 h relative to patients without these features. The Charlson Comorbidity Index and mAPACHE score remained stable between the pre-COVID-19 and COVID-19 periods; however, the proportion of patients who received targeted treatment for PA-BSI increased. Decreasing rates of PA resistance were also observed throughout the full study period (2009–2022). While the mortality rate among case-patients increased during the COVID-19 period compared to the pre-COVID-19 period, the COVID-19 period itself did not appear to be associated with increased mortality independent of COVID-19 coinfection. This remained true even when isolating the late-pre-COVID-19 period corresponding to the 3 years immediately preceding the COVID-19 period with the lowest observed mortality rate of the entire study period. These results suggest that COVID-19 coinfection may account for the increased PA-BSI mortality observed during the COVID-19 period.

To our knowledge, our study is the first description of the direct relationship between COVID-19 and mortality in patients with PA-BSI. Our findings are consistent with previous reports of poor outcomes among hospitalized patients with any BSI type and COVID-19 coinfection. A multicenter study of adults with BSI and COVID-19 coinfection in 271 US facilities showed higher rates of intensive care unit admission (73% versus 55%), longer length of stay (35 d versus 26 d), and higher rates of ventilator use (53% versus 22%) compared with patients in the prepandemic period.Reference Bauer, Puzniak and Yu20 Multiple studies have also identified higher mortality in patients with COVID-19 and BSI.Reference Zhu, Rawson and Mookerjee21–Reference Samantaray, Karan and Sharma23 Zhu et al. noted 27% higher mortality and 20 excess days of hospital stay in patients with hospital-acquired BSI during the COVID-19 pandemic period, regardless of COVID-19 coinfection status.Reference Zhu, Rawson and Mookerjee21

Several factors could account for the increased mortality among patients coinfected with BSI and COVID-19. Excessive focus on the diagnosis of COVID-19 infection may delay recognition and workup of coexisting or secondary processes, thereby delaying initiation of appropriate treatment.Reference Miyagami, Uehara and Harada24 Logistical barriers could affect management and subsequent outcomes in patients with COVID-19 related to isolation procedures and need for personal protective equipment.Reference Istanboulian, Rose, Yunusova and Dale25,Reference Petry, Pelzelmayer, Ernst, Thuerlimann and Naef26 Additionally, circumstances at the institutional level related to the pandemic period, such as nursing- or provider-to-patient ratio, hospital volume, and availability of infectious diseases expertise, could impact outcomes on an individual level. Finally, the interaction of the two infectious processes could lead to a synergistic or additive effect on infectivity and the host response. Respiratory viruses have been found to alter pulmonary tissue innate immune function and promote bacterial growth through an increased burden on macrophages and decreased dendritic cell and macrophage antigen-presenting capabilities.Reference Lucien, Canarie and Kilgore27,Reference Morgan, Casulli and Chew28 The inflammation and fluids, which accumulate within the alveoli secondary to SARS-CoV-2 virus infection in particular, are thought to serve as good substrates for PA and other bacterial respiratory pathogens.Reference Manohar, Loh and Athira29 Alteration of the respiratory tract microbiome may also play a role.Reference Hanada, Pirzadeh, Carver and Deng30

Our study has several limitations. The retrospective study design generated potential biases, particularly for mAPACHE score and resistance phenotype. Missing data may have been secondary to factors undefined within our data set, such as intensive care unit admission and differences in antimicrobial susceptibility testing protocols between facilities. Other factors not included but potentially relevant to the primary outcome were source of PA-BSI, other concomitant infections, COVID-19 treatment, indication of recurrent versus initial COVID-19 infection, ventilatory support, facility complexity, and ratio of staff support to hospital volume. Focusing on PA-BSI may limit the generalizability of the results to other types of BSI, though this must be weighed against the potential risk of diluting the results with inclusion of other BSI types involving different patient populations, settings, sources of infection, and risk factors. The mAPACHE score may not be the best measure of disease severity in context of PA-BSI, but it is more readily available relative to other measures such as the Pitt bacteremia score, which rely on more resource-intensive data collection or natural language processing tools. The small number of COVID-19 positive cases also represented a potential limitation, and it is unknown whether patients in this study with detectable SARS-CoV-2 virus had a clinical syndrome consistent with COVID-19 infection, though asymptomatic COVID-19 infection could still pose immunologic effects. Finally, this study focused on the VHA population, which may not capture the full spectrum of patients at high risk for PA infection, such as individuals with cystic fibrosis.

There are also notable strengths. This is the first study to examine the impact of the COVID-19 pandemic on outcomes for patients with PA-BSI. The VHA CDW serves as a rich database for this study, capturing data from over 150 medical centers throughout the continental United States and Puerto Rico and spanning a period of almost 20 yr.14 Moreover, the VHA data regarding resistance trends generally align with CDC data in that rates of MDR PA infection were observed to decrease in the last decade until 2020 when an increase was observed; in our VHA cohort, the increase from 12 MDR/DTR cases in 2019 (2.6%) to 21 cases in 2020 (4.2%) represents approximately a 60% increase. Another strength is that while the subset of patients with concomitant COVID-19 infection was modest in size, this study included data for nearly 7,000 unique patients, thereby serving as one of the largest studied cohorts of patients with PA-BSI.

Our findings lay the groundwork for many potential future studies. The relationship between COVID-19 and PA should be explored in more detail to determine whether characteristics of either pathogen enhance the infectivity of the other or work synergistically against the host immune response, and inclusion of observations of patients with COVID-19 with and without PA-BSI could be helpful to this end. More features related to the COVID-19 period would also be worthwhile to investigate, such as impact of staffing, hospital volume, and resources on PA-BSI-associated outcomes. Utilizing alternative measures of disease severity and complexity such as Pitt bacteremia score, Sequential Organ Failure Assessment score, neutrophil-to-lymphocyte ratio, or the Veterans Health Administration COVID-19 Index (VACO) for COVID-19 mortality could enhance these types of analyses. Further, novel disease severity scores relevant to PA-BSI with inclusion of COVID-19 infection status should be validated for inclusion in future analyses. Finally, interventions that represent innovative measures to improve care and outcomes of patients with PA-BSI, including those with concomitant infections such as COVID-19, would be hugely impactful given the prevalence and poor outcomes associated with these conditions.

Conclusions

Patients with PA-BSI had higher 30-d mortality during the COVID-19 period compared to the pre-COVID period. No targeted PA treatment within 48 h, higher mAPACHE score, and higher Charlson Comorbidity Index were also associated with higher mortality among all patients with PA-BSIs.

Acknowledgments

The authors would like to thank the Veterans Affairs Science and Health Initiative to Combat Infectious and Life-Threatening Diseases (VA SHIELD) for their support of this work.

Financial support

This work was in part supported by the National Institutes of Health (R01 AI063517-13 to RAB) and the Veterans Affairs Merit Review Program (BX001974 to RAB).

Competing interests

FP receives grant funding and research support from Accelerate, Merck, and Pfizer. RAB receives grant funding and research support from NIH, VA, VenatoRx, Merck, Wockhardt, Entasis, and the Cystic Fibrosis Foundation. The other authors report no conflicts of interest.

Footnotes

Previous Presentation of Data: A portion of the data presented in this manuscript was previously presented as an oral abstract at IDWeek 2022 in Washington D.C., and the abstract was subsequently published in Open Forum Infectious Diseases (Volume 9, Issue Supplement_2, December 2022, ofac492.027). The data and analyses have been updated and expanded for publication as a manuscript.

References

Rice, LB. Federal funding for the study of antimicrobial resistance in nosocomial pathogens: No ESKAPE. J Infect Dis 2008;197:1079–1081. doi: 10.1086/533452 CrossRefGoogle Scholar
Pang, Z, Raudonis, R, Glick, BR, Lin, TJ, Cheng, Z. Antibiotic resistance in Pseudomonas aeruginosa: Mechanisms and alternative therapeutic strategies. Biotechnol Adv 2019;37:177–192. doi: 10.1016/j.biotechadv.2018.11.013 CrossRefGoogle ScholarPubMed
Breidenstein, EBM, de la Fuente-Núñez, C, Hancock, REW. Pseudomonas aeruginosa: All roads lead to resistance. Trends Microbiol 2011;19:419–426. doi: 10.1016/j.tim.2011.04.005 CrossRefGoogle ScholarPubMed
Berrazeg, M, Jeannot, K, Ntsogo Enguéné, VY, et al. Mutations in β-Lactamase AmpC increase resistance of pseudomonas aeruginosa isolates to antipseudomonal cephalosporins. Antimicrob Agents Chemother 2015;59:6248–6255. doi: 10.1128/AAC.00825-15 CrossRefGoogle ScholarPubMed
Gerver, SM, Nsonwu, O, Thelwall, S, Brown, CS, Hope, R. Trends in rates of incidence, fatality and antimicrobial resistance among isolates of Pseudomonas spp. causing bloodstream infections in England between 2009 and 2018: results from a national voluntary surveillance scheme. J Hosp Infect 2022;120:73–80. doi: 10.1016/j.jhin.2021.11.013 CrossRefGoogle ScholarPubMed
Kang, CI, Kim, SH, Kim, HB, et al. Pseudomonas aeruginosa bacteremia: risk factors for mortality and influence of delayed receipt of effective antimicrobial therapy on clinical outcome. Clin Infect Dis 2003;37:745–751. doi: 10.1086/377200 CrossRefGoogle Scholar
Lodise, TP, Patel, N, Kwa, A, et al. Predictors of 30-day mortality among patients with Pseudomonas aeruginosa bloodstream infections: impact of delayed appropriate antibiotic selection. Antimicrob Agents Chemother 2007;51:3510–3515. doi: 10.1128/AAC.00338-07 CrossRefGoogle ScholarPubMed
Centers for Disease Control and Prevention (U.S.), National Center for Emerging Zoonotic and Infectious Diseases (U.S.), 2019. Division of Healthcare Quality Promotion. Antibiotic Resistance Coordination and Strategy Unit. Antibiotic resistance threats in the United States, 2019. https://stacks.cdc.gov/view/cdc/82532.Google Scholar
Centers for Disease Control and Prevention, National Center for Emerging and Zoonotic Infectious Diseases, Division of Healthcare Quality Promotion, 2022. COVID-19: U.S. Impact on Antimicrobial Resistance, Special Report 2022. https://stacks.cdc.gov/view/cdc/117915.Google Scholar
Sloot, R, Nsonwu, O, Chudasama, D, et al. Rising rates of hospital-onset Klebsiella spp. and Pseudomonas aeruginosa Bacteraemia in NHS acute trusts in England: a review of national surveillance data, August 2020-February 2021. J Hosp Infect 2022;119:175–181. doi: 10.1016/j.jhin.2021.08.027 CrossRefGoogle Scholar
Sturm, LK, Saake, K, Roberts, PB, Masoudi, FA, Fakih, MG. Impact of COVID-19 pandemic on hospital onset bloodstream infections (HOBSI) at a large health system. Am J Infect Control 2022;50:245–249. doi: 10.1016/j.ajic.2021.12.018 CrossRefGoogle Scholar
Giacobbe, DR, Battaglini, D, Ball, L, et al. Bloodstream infections in critically ill patients with COVID-19. Eur J Clin Invest 2020;50:e13319. doi: 10.1111/eci.13319 CrossRefGoogle ScholarPubMed
Buetti, N, Tabah, A, Loiodice, A, et al. Different epidemiology of bloodstream infections in COVID-19 compared to non-COVID-19 critically ill patients: a descriptive analysis of the Eurobact II study. Crit Care 2022;26:319. doi: 10.1186/s13054-022-04166-y CrossRefGoogle ScholarPubMed
Administration VH. VA.gov | Veterans Affairs. https://www.va.gov/health/aboutvha.asp. Accessed January 6, 2023.Google Scholar
Chang, ET, Yoon, J, Esmaeili, A, et al. Outcomes of a randomized quality improvement trial for high-risk Veterans in year two. Health Serv Res 2021;56:1045–1056. doi: 10.1111/1475-6773.13674 CrossRefGoogle ScholarPubMed
Fortis, S, O’Shea, AMJ, Beck, BF, et al. An automated computerized critical illness severity scoring system derived from APACHE III: modified APACHE. J Crit Care 2018;48:237–242. doi: 10.1016/j.jcrc.2018.09.005 CrossRefGoogle ScholarPubMed
Kadri, SS, Adjemian, J, Lai, YL, et al. Difficult-to-treat resistance in gram-negative bacteremia at 173 US hospitals: Retrospective cohort analysis of prevalence, predictors, and outcome of resistance to all first-line agents. Clin Infect Dis 2018;67:1803–1814. doi: 10.1093/cid/ciy378 Google ScholarPubMed
A.R. & Patient Safety Portal. https://arpsp.cdc.gov/. Accessed November 30, 2022.Google Scholar
Magiorakos, AP, Srinivasan, A, Carey, RB, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: An international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect 2012;18:268–281. doi: 10.1111/j.1469-0691.2011.03570.x CrossRefGoogle ScholarPubMed
Bauer, KA, Puzniak, LA, Yu, KC, et al. Epidemiology and outcomes of culture-positive bloodstream pathogens prior to and during the SARS-CoV-2 pandemic: A multicenter evaluation. BMC Infect Dis 2022;22:841. doi: 10.1186/s12879-022-07810-8 CrossRefGoogle ScholarPubMed
Zhu, NJ, Rawson, TM, Mookerjee, S, et al. Changing patterns of bloodstream infections in the community and acute care across 2 coronavirus disease 2019 epidemic waves: A retrospective analysis using data linkage. Clin Infect Dis 2022;75:e1082–e1091. doi: 10.1093/cid/ciab869 CrossRefGoogle ScholarPubMed
Rebold, N, Alosaimy, S, Morrisette, T, et al. Clinical Characteristics associated with bacterial bloodstream coinfection in COVID-19. Infect Dis Ther 2022;11:1281–1296. doi: 10.1007/s40121-022-00636-6 CrossRefGoogle ScholarPubMed
Samantaray, S, Karan, P, Sharma, A, et al. Prevalence, presentation and outcome of secondary bloodstream infections among COVID-19 patients. Infect Disord Drug Targets 2022;22:77–84. doi: 10.2174/1871526522666220418093450 CrossRefGoogle ScholarPubMed
Miyagami, T, Uehara, Y, Harada, T, et al. Delayed treatment of bacteremia during the COVID-19 pandemic. Diagnosis (Berl). 2021;8:327–332. doi: 10.1515/dx-2020-0114 CrossRefGoogle ScholarPubMed
Istanboulian, L, Rose, L, Yunusova, Y, Dale, C. Barriers to and facilitators for supporting patient communication in the adult ICU during the COVID-19 pandemic: A qualitative study. J Adv Nurs 2022;78:2548–2560. doi: 10.1111/jan.15212 CrossRefGoogle ScholarPubMed
Petry, H, Pelzelmayer, K, Ernst, J, Thuerlimann, E, Naef, R. Nurse-patient interaction during the Covid-19 pandemic: Creating and maintaining an interactive space for care. J Adv Nurs 2023;79:281–296. doi: 10.1111/jan.15486 CrossRefGoogle ScholarPubMed
Lucien, MAB, Canarie, MF, Kilgore, PE, et al. Antibiotics and antimicrobial resistance in the COVID-19 era: Perspective from resource-limited settings. Int J Infect Dis 2021;104:250–254. doi: 10.1016/j.ijid.2020.12.087 CrossRefGoogle ScholarPubMed
Morgan, DJ, Casulli, J, Chew, C, et al. Innate immune cell suppression and the link with secondary lung bacterial pneumonia. Front Immunol 2018;9:2943. doi: 10.3389/fimmu.2018.02943 CrossRefGoogle ScholarPubMed
Manohar, P, Loh, B, Athira, S, et al. Secondary bacterial infections during pulmonary viral disease: Phage therapeutics as alternatives to antibiotics? Front Microbiol 2020;11:1434. doi: 10.3389/fmicb.2020.01434 CrossRefGoogle ScholarPubMed
Hanada, S, Pirzadeh, M, Carver, KY, Deng, JC. Respiratory viral infection-induced microbiome alterations and secondary bacterial pneumonia. Front Immunol 2018;9:2640. doi: 10.3389/fimmu.2018.02640 CrossRefGoogle ScholarPubMed
Figure 0

Figure 1.Flow diagram of inclusion and exclusion criteria resulting in 6,714 cases and further categorization into COVID-19 period and coinfection cohorts.

Figure 1

Table 1.Demographic characteristics among patients with Pseudomonas aeruginosa bloodstream infection (n = 6,714) stratified by COVID-19 status and COVID-19 period

Figure 2

Figure 2.Figure demonstrates that the proportion of PA-BSI isolates within the VHA system susceptible to all tested agents increased over the study period, with a similar increase in proportion of all cases in which targeted treatment was administered within 48 h. A corresponding decrease in MDR and DTR resistance was also observed over the study period. Among patients with COVID-19 in addition to PA-BSI relative to patients with PA-BSI only, more isolates had OUR in 2020–2021 and more MDR was present in 2021–2022, but no DTR was present. Bar chart excludes 1,177 patients with missing resistance phenotype data. Line plot excludes 94 patients with missing data for targeted treatment within 48 h. 2022 data represent January 1 through June 30. PA-BSI, Pseudomonas aeruginosa bloodstream infection; VHA, Veterans Health Administration; ER, expected resistant; OUR, otherwise unclassified resistant; MDR, multidrug-resistant; DTR, difficult-to-treat resistant.

Figure 3

Figures 3.Figure indicates the proportion of PA-BSI cases in the VHA system meeting the primary end point of mortality at 30 d during each year of the study period. Mortality decreased across the pre-COVID-19 period to a low of 20.6% and increased during the COVID-19 period, peaking at 30.8% in 2022 through June 30. PA-BSI, Pseudomonas aeruginosa bloodstream infection; VHA, Veterans Health Administration.

Figure 4

Figure 4.Distribution of clinical characteristics and mortality rate are shown stratified by COVID-19 period and COVID-19 coinfection status. Age (A), Charlson Comorbidity Index (B), and modified APACHE score (C) are noted to have a similar distribution across each cohort, despite that the mortality rate (D) is significantly higher in the COVID-19-positive/COVID-19 period patients.

Figure 5

Table 2.Predictors of 30-d mortality among patients with Pseudomonas aeruginosa bloodstream infection (n = 6,714) with and without adjustment for COVID-19 coinfection

Figure 6

Figure 5.Modified Poisson analysis results illustrate the positive association between COVID-19 coinfection and increased all-cause 30-d mortality risk in patients with PA-BSI in the VHA system during the study period. Additional predictors associated with 30-d mortality include Charlson Comorbidity Index, administration of targeted treatment within 48 h, and modified APACHE score.