Article Text

COVID-19 hospital prevalence as a risk factor for mortality: an observational study of a multistate cohort of 62 hospitals
  1. Mohamad Ghazi Fakih1,
  2. Allison Ottenbacher2,
  3. Baligh Yehia1,
  4. Richard Fogel1,
  5. Collin Miller2,
  6. Angela Winegar2,
  7. Christine Jesser3,
  8. Joseph Cacchione1
  1. 1Clinical & Network Services, Ascension Health, Saint Louis, Missouri, USA
  2. 2Ascension Data Science Institute, Ascension Health, Saint Louis, Missouri, USA
  3. 3Ascension Clinical Research Institute, Ascension Health, Saint Louis, Missouri, USA
  1. Correspondence to Dr Mohamad Ghazi Fakih, Clinical & Network Services, Ascension Health, Saint Louis, MO, 63134, USA; Mohamad.Fakih{at}ascension.org

Abstract

Background The associated mortality with COVID-19 has improved compared with the early pandemic period. The effect of hospital COVID-19 patient prevalence on COVID-19 mortality has not been well studied.

Methods We analysed data for adults with confirmed SARS-CoV-2 infection admitted to 62 hospitals within a multistate health system over 12 months. Mortality was evaluated based on patient demographic and clinical risk factors, COVID-19 hospital prevalence and calendar time period of the admission, using a generalised linear mixed model with site of care as the random effect.

Results 38 104 patients with COVID-19 were hospitalised, and during their encounters, the prevalence of COVID-19 averaged 16% of the total hospitalised population. Between March–April 2020 and January–February 2021, COVID-19 mortality declined from 19% to 12% (p<0.001). In the adjusted multivariable analysis, mid and high COVID-19 inpatient prevalence were associated with a 25% and 41% increase in the odds (absolute contribution to probability of death of 2%–3%) of COVID-19 mortality compared with patients with COVID-19 in facilities with low prevalence (<10%), respectively (high prevalence >25%: adjusted OR (AOR) 1.41, 95% CI 1.23 to 1.61; mid-prevalence (10%–25%): AOR 1.25, 95% CI 1.13 to 1.38). Mid and high COVID-19 prevalence accounted for 76% of patient encounters.

Conclusions Although inpatient mortality for patients with COVID-19 has sharply declined compared with earlier in the pandemic, higher COVID-19 hospital prevalence remained a common risk factor for COVID-19 mortality. Hospital leaders need to reconsider how we provide support to care for patients in times of increased volume and complexity, such as those experienced during COVID-19 surges.

  • COVID-19
  • mortality (standardised mortality ratios)
  • hospital medicine

Data availability statement

All data relevant to the study are included in the article or uploaded as supplemental information. Deidentified participant data are available from the Ascension Data Science Institute team.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Introduction

In January 2020, the first case of COVID-19 was reported in the USA.1 By 25 August 2021, the USA has reported over 38.1 million confirmed cases and 631 600 deaths.2 Over the 12 months of the pandemic, we have witnessed an improvement in in-hospital COVID-19 mortality especially compared with the initial few months.3 4 In addition, our understanding of the clinical course and management of the disease has rapidly evolved. Some reports attributed the improvement in COVID-19 mortality partly to a change in hospitalised patient population characteristics, combined with an enhancement of medical management.5 Others focused on the relationship between COVID-19-infected patient load on the intensive care environment and the associated adverse outcomes.6

Many therapeutic and management strategies have been deployed to optimise patient outcomes with varying results.7–9 The use of dexamethasone has been associated with improved mortality, while the outcomes for other therapies were mixed. Recognising that COVID-19 can increase the risk of thrombosis, clinical guidelines recommend the use of pharmacological prophylaxis for venous thromboembolism in those hospitalised with COVID-19.10 Access to new antiviral and immune-based therapies also increased.11 In addition, the availability of COVID-19 vaccines in 2021 has provided essential protection against infection, deterioration and hospitalisation. On the other hand, the care of patients infected with COVID-19 requires substantial human, operational and clinical resources. Sudden large increases in the number of patients cared for with COVID-19 have led to increased strain on hospital capacity, especially within the critical care environment.6 12

To assess the impact of hospital strain on COVID-19 mortality, we evaluated changes in hospitalised patient population characteristics and hospital COVID-19 patient prevalence over the first 12 months of the pandemic in a large multistate health system in 11 US states.

Methods

Study design, setting and population

We conducted an observational study of a multistate health system in the USA. Adults (≥18 years old) with confirmed SARS-CoV-2 infection (positive result by PCR testing of a nasopharyngeal sample, or a COVID-19 encounter diagnosis of International Classification of Diseases, 10th Revision (ICD-10) code U07.1) admitted to 62 Ascension hospitals between 1 March 2020 and 28 February 2021 were included. Hospitals are in 11 states: Alabama (n=5), Florida (n=7), Indiana (n=7), Kansas (n=3), Maryland (n=1), Michigan (n=12), New York (n=2), Oklahoma (n=4), Tennessee (n=5), Texas (n=6) and Wisconsin (n=10). Overall, 16 (26%) facilities had fewer than 100 beds, 22 (35%) facilities had between 100 and 300 beds, and 24 (39%) facilities had more than 300 beds. As a system, standard protocols were established early in the pandemic, regularly updated and implemented across the sites of care for COVID-19 screening, testing and management. We followed the Strengthening the Reporting of Observational Studies in Epidemiology reporting guidelines.13

Data collection

Patient sociodemographic characteristics, laboratory results and healthcare utilisation, sourced from electronic health records and administrative data, were abstracted, quality assured and assembled into a uniform database. Race was self-reported and categorised as white, black/African American or other. Patient age was standardised by subtracting the mean and dividing by the SD. Based on the date of hospital admission, patients were grouped into six consecutive 2-month time periods from March 2020 to February 2021 to evaluate trends.

The Agency for Healthcare Research and Quality (AHRQ) Elixhauser Comorbidity Index (ECI), a method of assessing comorbidities based on 29 groups of ICD-10 diagnosis codes, was calculated for each patient.14–16 The index ranged from −23 to 67 in our sample, with higher numbers indicating greater comorbidity. The ECI variable used in the model was standardised by subtracting the mean and dividing by the SD. The county-level overall Social Vulnerability Index (SVI)17 percentile rank was mapped to each encounter based on county name. However, because county was often missing in the source data, the zip code associated with the encounter was used to identify the county, which was then used to map SVI data. The categorical variable for SVI used in the model is computed based on overall percentile rank for all counties in the USA (higher percentile=higher vulnerability): low (<0.25), mid-low (0.25–<0.50), mid-high (0.50–<0.75), high (≥0.75).

To account for the environment where patients infected with COVID-19 may require more extensive utilisation of resources and frontline worker support, we created a variable to represent the strain on the hospital during different levels of COVID-19 hospital prevalence. The percent of inpatient hospital census that was COVID-19 positive was calculated daily for each facility by dividing the number of patients with COVID-19 by the total number of admitted patients. For each encounter, a patient-specific measure of the hospital COVID-19 prevalence was calculated by averaging (taking the mean) the daily facility per cent COVID-19 positive (referenced above) for each day between the admit date and discharge date. Each patient was assigned the corresponding categorical variable indicator for per cent COVID-19 in-hospital prevalence as follows: low (<10%), mid (10%–25%), high (>25%). The categories for low, mid and high prevalence were assigned based on the distribution of encounters with approximately 25% of encounters experiencing low prevalence; 50% of encounters experiencing mid-prevalence and 25% of encounters experiencing high prevalence.

Statistical analysis

Statistical comparisons of patient characteristics were assessed across (1) the six 2-month time periods and (2) levels of COVID-19 hospital prevalence. Unadjusted associations between characteristics and COVID-19 hospital prevalence and time periods were evaluated using the Χ2 test of independence for categorical variables and Kruskal-Wallis tests for continuous variables.

A generalised linear mixed model with site of care as the random effect was used to evaluate the relationship between COVID-19 prevalence group and all-cause in-hospital COVID-19 mortality. Variables that were independently associated with mortality, or identified through a literature review of previous studies that evaluated the impact of patient and clinical characteristics on COVID-19 mortality, or expert clinical recommendation, were included in the final model. Model diagnostics and fit were evaluated on the entire dataset with assessments of multicollinearity and linearity between the independent variables and the logit of the binary outcome variable. Age was considered a categorical variable for descriptive analysis but used as a continuous variable in the full model, which calculated the increased odds of mortality associated with a 1 SD increase in age (17 years). This analysis was conducted for the full sample, and separately for patients admitted to the intensive care unit (ICU). A sensitivity analysis that excluded the first 2 months of the pandemic was conducted. Two-sided testing was used and p<0.05 was considered statistically significant. The margins package in R was used to calculate the average marginal effect for all independent variables. The average marginal effect can be interpreted as the absolute per cent difference in the probability of mortality for each level of a categorical variable relative to the reference group or a one unit change for continuous variables. All statistical analyses were performed using R V.4.0.2.

Results

Characteristics of patients hospitalised with COVID-19

The final sample included 38 104 patients with confirmed COVID-19, representing 7% of all admissions during 12 months of the study, and 16% of the average census during their encounters. A total of 19 270 patients (51%) were male and 8941 (23%) were black/African American. Median age was 67 years (IQR: 54–78); 21 692 (57%) were covered by Medicare insurance; median ECI was 9 (IQR: 0–17) and 8180 (21%) had high social vulnerability. Common (>30% of patients) comorbidities included diabetes (42%), hypertension (42%) and obesity (35%). On admission, 33% of patients had an oxygen saturation <94%, 21% had a respiratory rate ≥24 and 13% had a temperature ≥38°C (table 1). Overall, 24% of patient encounters occurred when there was a low COVID-19 hospital prevalence (<10%), 53% of encounters occurred during times of mid-COVID-19 hospital prevalence (10%–25%), and 23% occurred at times of high COVID-19 hospital prevalence (>25%). The distribution of COVID-19 hospital prevalence changed across the 12-month period, ranging from 51% high COVID-19 prevalence in March–April to 2% high COVID-19 census in September–October (figure 1). The number of patients infected with COVID-19 who were hospitalised during the 12 months of the study varied, with 54% admitted between November 2020 and February 2021.

Figure 1

The number of patients and COVID-19 hospital prevalence per facility, between March 2020 and February 2021.

Table 1

Characteristics of inpatients with COVID-19 admitted between March 2020 and February 2021, by COVID-19 hospital prevalence

Distributions of race, age, sex, insurance, comorbidities (excluding asthma and transplant), ECI, social vulnerability, vital signs, ICU admission and hospital COVID-19 prevalence all significantly varied over this time period (p<0.05) (online supplemental appendix 1).

Supplemental material

Outcomes of patients hospitalised with COVID-19

Overall, 16 756 (44%) patients infected with COVID-19 were admitted to the ICU and 5009 (13%) patients with COVID-19 died in the hospital (table 1). In-hospital COVID-19 mortality initially decreased from 19% in March–April, to 10% in May–June, then climbed back to 14% in November–December and dropped again in January–February 2021 to 12% (p<0.01). After excluding March–April, there was a significant increase in mortality between September–October and November–December (12% to 14%, p<0.001), and a significant decrease between November–December and January–February (14% to 12%, p<0.001). The percent of patients with COVID-19 admitted to the ICU ranged from 42% to 46% during the 12-month time frame. Mortality was highest for patients older than 85 years, ranging from 31% in March–April to 18% in January–February (figure 2).

Figure 2

Mean COVID-19 hospital prevalence and mortality of inpatients infected with COVID-19 between March 2020 and February 2021, by age.

COVID-19 prevalence and in-hospital COVID-19 mortality

Adjusting for time period, demographics, comorbidities and baseline vital signs, there were increased odds of COVID-19 mortality for patients admitted and cared for during periods of mid (adjusted OR (AOR) 1.25, 95% CI 1.13 to 1.38) and high (AOR 1.41, 95% CI 1.23 to 1.61) COVID-19 prevalence, compared with times with low COVID-19 hospital prevalence. Higher adjusted odds of mortality for mid (AOR 1.32, 95% CI 1.18 to 1.49) and high (AOR 1.53, 95% CI 1.30 to 1.80) COVID-19 prevalence, compared with low, were also observed for the subset of patients who required ICU stay. After excluding the early pandemic period and examining the data in the 10 months from May 2020 to February 2021, an independent association of increased risk of mortality with higher COVID-19 prevalence remained (mid AOR 1.15, 95% CI 1.03 to 1.29 and high AOR 1.21, 95% CI 1.03 to 1.42) (online supplemental appendix 2). During the 10-month period, the odds of mortality remained significantly higher for ICU patients with COVID-19 in the mid-prevalence group, but was not significant for the high COVID-19 prevalence group (mid AOR 1.19, 95% CI 1.05 to 1.36; high AOR 1.20, 95% CI 0.99 to 1.45).

Other factors associated with in-hospital COVID-19 mortality

Adjusting for other factors, compared with patients admitted in March–April 2020, patients admitted in subsequent months had lower odds of COVID-19 in-hospital mortality, May–June (AOR 0.66, 95% CI 0.55 to 0.79), July–August (AOR 0.67, 95% CI 0.57 to 0.78), September–October (AOR 0.71, 95% CI 0.61 to 0.84), November–December (AOR 0.72, 95% CI 0.63 to 0.82), and January–February (AOR 0.58, 95% CI 0.51 to 0.68) (table 2). Older age (AOR 1.88, 95% CI 1.78 to 1.99, reported race as ‘other’ (AOR 1.33, 95% CI 1.16 to 1.52), male sex (AOR 1.27, 95% CI 1.19 to 1.36) and being a Medicare recipient (AOR 1.17, 95% CI 1.05 to 1.29) were significantly associated with increased odds of COVID-19 mortality compared with their respective reference groups. Chronic conditions (obesity, diabetes, liver disease, ischaemic heart disease) and less favourable vital signs on admission (respiratory rate ≥24 breaths/min, oxygen saturation <94%) were also significantly associated with increased odds of COVID-19 mortality.

Table 2

Adjusted ORs (AORs) of in-hospital mortality among patients with COVID-19 admitted between March 2020 and February 2021

Among the factors that were associated with in-hospital COVID-19 mortality, ECI contributed the highest average marginal effect (8%) for in-hospital mortality between March 2020 and February 2021 (figure 3). Other patient characteristics that had 5% or greater contribution include respiratory rate on admission ≥24 breaths/min (7%), obesity (6%), age (6%), oxygen saturation on admission <94% (6%) and liver disease (5%). The average marginal effect of diabetes was 3.6% and close to the effect of high COVID-19 prevalence (3%).

Figure 3

Characteristics with increased (absolute) per cent contribution to probability of death (average marginal effects) among patients with COVID-19. *COVID-19 prevalence compared with low COVID-19 prevalence category (<10%).

Additionally, factors found to be significant contributors of increased odds of COVID-19 mortality during the 12-month study were compared across each 2-month time period with results showing that while overall mortality decreased among patients hospitalised with COVID-19, the contribution of each risk factor to in-hospital COVID-19 mortality did not significantly differ across time periods (data not shown).

Discussion

In this large multistate cohort of patients hospitalised with COVID-19 over a 12-month period, we observed an improvement in COVID-19 hospital mortality between the first 2 months of the pandemic compared with later time periods. We also identified a significant association between hospital COVID-19 prevalence during the patient stay and corresponding COVID-19 mortality. Hospitals assumed the care for a large proportion of patients who rapidly decompensated and necessitated considerable levels of support. Adjusting for patient, clinical and hospital characteristics, the odds of dying for patients with COVID-19 infection increased between 25% and 41% when COVID-19 hospital prevalence was greater than 10% compared with patients with COVID-19 hospital prevalence less than 10%. These associations were slightly higher for patients with COVID-19 requiring intensive care support. Increased COVID-19 hospital prevalence strains health system resources via increased volumes of sicker patients who require more intense support,18 coupled with increased shortages of healthcare workers either being infected with COVID-1919 20 or quarantined due to exposure.

Hospitals experienced surge periods where rapid escalation of admissions occurred. More than three-quarters of patients were cared for in periods of either mid or high hospital COVID-19 prevalence. The risk of higher mortality with increased hospital COVID-19 prevalence persisted on sensitivity analysis that excluded the first two pandemic months. Higher COVID-19 hospital prevalence remained an independent risk factor for COVID-19 mortality, even as treatment knowledge improved over time. For example, the management of COVID-19 respiratory failure, treatment with steroids and antivirals, and thrombosis prophylaxis were standardised across our hospitals by May 2020. Our health system multidisciplinary team of clinical experts continues to update our clinical guidance to manage COVID-19-infected patients. The updates are regularly disseminated to clinicians and implemented locally. Acute increases in the number of inpatients infected with COVID-19 may overwhelm clinicians and disrupt standard operations, especially in high-demand areas, such as the critical care environment.6 12 ICU strain21 and bed shortages22 have been historically associated with increased mortality, but our findings show a more pronounced effect associated with COVID-19 infection. In addition, nurse staffing shortages have been perceived to adversely impact the care of patients during the pandemic.23 24 A recent survey of nurses from more than 250 US hospitals reported burnout and understaffing, suggesting that the nursing workforce capacity is not prepared to withstand substantial increases in resource utilisation.25 Hospitals have attempted to mitigate shortages by recruiting travelling nurses or retraining others from different disciplines. It is uncertain that such accommodations optimally support the care of patients with COVID-19. Medical specialties, like infectious diseases and critical care medicine, heavily vested in the management of patients with COVID-19, have also been disproportionately impacted. A recent survey of more than 2000 critical care physicians revealed increased physical and emotional exhaustion, combined with staffing shortages.26 Proactive planning to adjust workforce staffing based on pandemic levels may provide more capacity and less risk for burnout.27

Some processes may have been altered during the pandemic to reduce potential exposure of frontline workers to patients. For example, the physical care of patients was often limited to essential nurses and physicians, while other clinicians performed ‘electronic medical record rounding’, leading to a reduction in the number of evaluations of patients per day.28 Intravenous pumps were placed in hallways in some hospitals to avoid frequent entry of nurses to rooms of patients with COVID-19.29 Process changes and increased staff workload may result in a delay in recognition of a change in patient condition. We have previously reported that hospitals with increased COVID-19 prevalence had significantly higher central line-associated bloodstream infection (CLABSI) events.30 In addition, patients infected with COVID-19 were greater than five times likely to develop CLABSIs compared with other patients during the pandemic.30 Our clinical and quality teams are re-evaluating our processes to ensure the standards of care are hardwired, patient decompensation is recognised early and prompt management is instituted for patients who are clinically deteriorating.

Our study confirms previously described associations linking intrinsic patient factors such as older age and greater comorbidity with higher in-hospital COVID-19 mortality.31 32 We observed that age, obesity, liver disease and diabetes were independently associated with increased COVID-19 mortality risk. Similar to what we have reported previously, black/African Americans did not have higher risk of COVID-19 mortality compared with white individuals over the study period.33 Another important element was the patient’s clinical status at hospital presentation, where respiratory compromise was associated with twice the odds of COVID-19 mortality.

Our findings have important implications on quality of care and outcomes. First, we need to reconsider how we best support the care for patients in times of increased volume and complexity, such as those experienced during COVID-19 surges. Traditionally, clinical resources have been allocated based on the physical location of patient care (intensive care vs general medical units). The COVID-19 pandemic stressed resources even for those cared for in ICUs.6 Second, a re-evaluation of the healthcare workforce structure and its needs is essential.34 Many hospitals have limited capacity and have shown inadequate flexibility in buffering acute demands such as those seen during the pandemic. A national plan may be necessary to deploy resources—physically and virtually—and weather large disruptive events. There is a substantial opportunity to build partnerships among government, local entities and health systems. Third, the value of investing in quality and infection prevention, coupled with a robust analytics infrastructure, have surfaced as essential elements to identify acute changes in processes and associated adverse outcomes. Finally, national physician, nursing and specialty societies should re-examine our current state and whether we are ready to support workforce challenges in the future.

Our study has limitations. We have not specifically evaluated the therapeutic interventions that patients received during the study. We used hospital COVID-19 prevalence as a surrogate for hospital strain and healthcare workforce stress, such as staffing shortage or availability of specialty services. Results were dependent on the categorisation of COVID-19 prevalence and may vary if alternate cut-off points were specified. In addition, we have not included in the analysis the frequency of testing for patients hospitalised during the 12-month period or the potential effect of the circulating SARS-CoV-2 strain over time. Lastly, while our hospitals are located in diverse settings across multiple states, our experience may not be fully representative of other health systems in the USA or generalisable to other countries. Moreover, we have had robust system support to ensure supplies are always available to provide proper patient care. Our findings may be accentuated in other settings outside of the USA where resources may not be as available. Limitations notwithstanding, our results are based on a cohort of almost 40 000 patients with COVID-19 and representative of 62 hospitals in 11 states, making our findings generalisable to the US populations.

We conclude that COVID-19 hospital prevalence is an important measure that is significantly associated with higher inpatient mortality with COVID-19 infection, independent of intrinsic patient factors. Hospital leaders should further evaluate capacity and support to address COVID-19 surges while collaborating with policymakers to support enhanced hospital preparedness to care for patients with more intense needs.

Data availability statement

All data relevant to the study are included in the article or uploaded as supplemental information. Deidentified participant data are available from the Ascension Data Science Institute team.

Ethics statements

Patient consent for publication

Ethics approval

The study was approved with a waiver of informed consent by the Ascension Seton Institutional Review Board.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Twitter @TheQualityDoc

  • Contributors All authors have had substantial contributions to the work: MGF, JC, BY, RF and AW for the conception and design of the work; AW, AO and CM for the acquisition of data; AW, AO, CM, MGF and CJ for data analysis. All authors were involved in drafting and revising the work critically for important intellectual content and have given final approval of the version to be published. MGF is responsible for the overall content as guarantor. MGF accepts full responsibility for the finished work and/or the conduct of the study, had access to the data and controlled the decision to publish.

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.