Hospital load and increased COVID-19 related mortality - a nationwide study in Israel

The rapid spread of Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) worldwide and the disease caused by the virus, coronavirus disease 19 (COVID-19), has caused a global pandemic (1) with devastating social and economic consequences. Throughout the pandemic, several health systems were overwhelmed in light of the rapid emergence of new cases within a short period of time. Notable examples include Lombardy in Italy (2) in which ICU capacity has reached its limit in March 2020, and New York city in the USA (3). The heavy workload imposed on hospital services might have negatively affected patients’ outcomes and exacerbated mortality rates.

Here, we assessed excess mortality using a model developed for predicting patient mortality based on data of day-by-day patient disease course (4). We show that during a peak of hospitalizations in September and October 2020, patient deaths significantly exceeded the model’s mortality predictions, while reverting to match the predictions as patient load subsided in October, and showing signs of renewed excess mortality as hospital load is increasing again since late December.

We analysed data on COVID-19 related mortality in Israel from 15/07/2020 to 06/01/2021. The first individual with COVID-19 in Israel was identified on February 21st 2020. In response, the Israeli Ministry of Health (MOH) has gradually employed a series of social distancing measures in order to mitigate the spread of the virus (5). Following the relaxation of these measures in May 2020, the number of new patients has substantially increased and on the 10th of September 2020, Israel became the country with the highest rate of COVID-19 infections per capita worldwide (6). This period was referred to as the second wave of COVID-19 disease and as a result, Israel was the first country worldwide to impose a second lockdown in mid September 2020. The restrictions were gradually loosened during October 2020, followed by an additional increase in the number of cases leading the government to impose a third lockdown in January 2021. In-hospital mortality of individuals with COVID-19 throughout the study period is described in Figure 1A.

• Download figure
• Open in new tab

Figure 1:

A. In hospital COVID-19 related mortality cases by date. Daily mortality is marked by gray circles, and a 7-day average is plotted as a solid line. B. Number of COVID-19-related hospitalizations by day in all hospitals in Israel. Colors depict different clinical states. A threshold of 500 severe and critical patients per day is presented by a dashed horizontal line. Dashed vertical lines separate weeks which are over/under this threshold. C. COVID-19 related deaths for sets of patients that were first hospitalized in different weeks. Observed deaths for each week are marked by a solid line. Expected deaths as predicted from the model are marked in dashed black line. Monte-Carlo-based pointwise 10%-90% confidence predictions are marked in gray. Weeks in which excess deaths were observed (positive difference between true and expected curves) are filled in red; and weeks in which deaths were overestimated (negative difference between true and expected curves) are filled in green. X-axis is marked by the week number.

Abbreviations: SE - Standard Error

Mortality predictions were made using Monte-Carlo methods based on a multistate survival analysis (Fig S1) and a set of Cox regression models, first constructed and validated on a nationwide cohort during the first stages of the pandemic in Israel (4). The model predicts the clinical course of individual patients and adjusts for right censoring, recurrent events, competing events, left truncation, and time-dependent covariates. The original aim of the model was to allow timely allocation of sufficient healthcare resources and skilled medical professionals by medical centers. Weekly predictions of mortality and number of severe cases based on the model were presented and utilized by policy-makers in Israel.

Patient data included information on age, sex, date of positive SARS-CoV-2 polymerase chain reaction (PCR) test, date of hospitalization, and clinical outcome (death or discharge from hospitalization) for each individual. In addition, daily information on disease severity during hospitalization was available. Classification of disease severity was based on the following clinical criteria, applied on 13 of July 2020 by the Israeli NOH: Mild illness - individuals who have any of the various signs and symptoms of COVID 19 (e.g., fever, cough, malaise, loss of taste and smell); Moderate illness - individuals who have evidence of pneumonia by a clinical assessment or imaging; Severe illness - individuals who have respiratory rate >30 breaths per minute, SpO2 <93% on room air at sea level, or ratio of arterial partial pressure of oxygen to fraction of inspired oxygen (PaO2/FiO2) <300 mmHg and Ventilated/Critical (denoted in this paper as Critical) - individuals with respiratory failure who require ventilation (invasive or non-invasive), multiorgan dysfunction or shock (7). These criteria were determined based on NIH (8) and WHO (1) definitions.

We applied the prediction model on nationwide hospitalization data originated from the Israeli MOH to predict mortality of hospitalized patients with COVID-19 based on their age, sex, and clinical state on their first day of hospitalization. Overall, from 15/07/2020 to 06/01/2021, 19,336 individuals were hospitalized with COVID-19 infection in Israel and were included in the analysis. Mean age was 59.33± 21.81 years old (median age was 63 years old), and 9,421 (48.72%) were females. We first divided the data into 25 weeks, and assigned each patient according to week of hospitalization. Patients’ characteristics are presented in Table 1. Patients who died on the same day of hospital admission were not included in the analysis.

The model, already trained and validated on a cohort of 2,703 hospitalized patients early in the course of the pandemic(4), was modified, re-trained and validated on data of 5,966 individuals, hospitalized from 15/07/2020 to 08/09/2020 (time-period I). The model was then applied to data of individuals hospitalized from 09/09/2020 to 06/01/2021. We divided this interval into three time periods (denoted II, III, IV) according to hospital load: Periods II (09/09/2020 to 27/10/2020) and IV (15/12/2020 to 06/01/2021) are those where the daily number of severe+critical patients exceeded 500. Excess mortality was calculated by the observed versus predicted mortality. Strikingly, during time periods II and IV, the observed 14-day mortality was respectively 22.1% (Standard Error (SE) 3.1%) and 27.9% (SE 5.3%) higher than predicted by the model (Table 2), whereas our model accurately predicted per-day COVID-19 related mortality during periods I and III (Fig 1A), where the observed 14-day mortality was 0.1% (SE 5.1%) and 10.3% (SE 5.7%) higher than predicted by the model respectively. Similar trends were observed for 28-day mortality (Table 2). Cumulative expected and actual death curves for each hospitalization week are presented in Fig S2; Fig S3 presents calibration plots by hospitalization week.

Of note, the threshold defined by policy makers in Israel as an upper bound in which the healthcare system won’t be able to adequately treat patients was 800 severe+critical patients (9). In addition, we note that the increase in observed mortality is despite the fact that throughout the pandemic, clinical experience in treatment of COVID-19 patients increased, along with a better understanding in pharmacologic (such as corticosteroids (10) and remdesivir (11)) and non-pharmacologic (such as proning (12)) treatment modalities that may be beneficial for the patients.

We postulate that the excess mortality is most likely due to the rapid escalation in the number of hospitalized patients with COVID-19 during these time-periods in Israel, which may have resulted in an insufficiency of health-care resources, thereby negatively affecting patient outcomes. Theoretically, several other explanations may account for the increased mortality observed. First, it is possible that the excess deaths may be driven by a more vulnerable patient population, with a higher risk for mortality, that was admitted for hospitalization around time-periods II and IV. However, the model adjusts for age, sex and clinical state upon first day of hospitalization, and the predictions therefore take these differences into account. Although our data did not include information on patient comorbidities that may also influence the mortality rate from COVID-19, it was previously shown that this information is not necessarily essential for accurate mortality prediction for hospitalized patients when utilizing a multistate survival model (4). Accounting for clinical state also somewhat reduces the likelihood of increased mortality being due to deferred hospitalizations. Second, it is possible that during these time periods a more virulent strain was circulating in Israel. However, there is no evidence for the existence of such a strain and the fact that the time period lasted for only 7 weeks, and was followed by a time-period in which the model achieved accurate predictions makes it unlikely. The fact that increased mortality is once again observed in time period IV, in which hospital load is once again high, further strengthens our hypothesis.

Our study has several limitations. First, the increased mortality observed during periods of high rate of COVID-19 related morbidity may be due to factors specific to the Israeli healthcare system. The Israeli healthcare system is universal and mandates all citizens to join one of the official non-profit health insurance organizations. Regional (13) or financial (14) disparities affecting the availability of health-care resources in other countries as well as racial variation (15) may affect COVID-19-related mortality. Future studies should be conducted in order to determine if this effect is also observed in other countries with different healthcare systems, and the specific threshold of patients representing healthcare capacity in each healthcare system. Full model code is available for use at https://github.com/JonathanSomer/covid-19-multi-state-model. Second, our findings may be influenced by the availability of diagnostic testing for COVID-19 as well as accurate and complete documentation of the disease severity state by clinicians. However, testing policy and physicians documentation practices did not change significantly in Israel throughout the study period. Moreover, we only included data after 13 of July 2020, in which uniform criteria for COVID-19 disease severity were applied by the MOH in all hospitals in Israel.

In conclusion, here we have shown that the mortality of hospitalized patients with COVID-19 in Israel was associated with health-care burden, reflected by the simultaneous number of hospitalized patients in a severe condition. Our work emphasizes that even in countries in which the healthcare system did not reach a specific point defined as insufficiency, the increase in hospital workload was associated with quality of care and patient mortality, ruling out factors related to change in the hospitalized population. In addition, our study highlights the importance of quantifying excess mortality in order to assess quality of care, and define an appropriate carrying capacity of severe patients in order to guide timely healthcare policies and allocate appropriate resources.

Supplementary appendix

• Download figure
• Open in new tab

Fig S1.

Multi-state model. Patients disease course transitions between 5 possible clinical states: mild or moderate, severe, critical, discharged and deceased. Each transition was modelled using a set of Cox regression models, adjusting for right censoring, recurrent events, competing events, left truncation, and time-dependent covariates.

• Download figure
• Open in new tab

Fig S2.

Cumulative expected and observed COVID-19 in-hospital deaths by week of first hospitalization. Each sub-figure presents curves for the set of patients first hospitalized in a given week. Black curves are expected deaths with orange bands representing Monte-Carlo-based pointwise 10%-90% confidence predictions. Blue and green curves are actual deaths for weekly sets of patients. Blue curves are drawn for patients which the model was trained on (weeks -8 to -1), and green curves are drawn for patients which the model was not trained on (weeks 0 to 16).

• Download figure
• Open in new tab

Fig S3.

Calibration curves - true in-hospital deaths of patients with COVID-19 versus expected deaths by week of first hospitalization. Points are drawn for time bins of 5 days in the range of 5-100, days from hospitalization to death. Blue curves are drawn for patients which the model was trained on (weeks -8 to -1), and green curves are drawn for patients which the model was not trained on (weeks 0 to 16). Diagonal dashed line represents perfect calibration.