What pushed Israel out of herd immunity? Modeling COVID-19 spread of Delta and Waning immunity

INTRODUCTION

First discovered in India in October 2020, the Delta variant became the dominant variant by March 2021 [1], and caused the deadliest surge to date, with over 200,000 documented COVID-19 associated deaths in two months [2]. Delta quickly became the dominant variant in most countries after introduction, replacing all other variants of concern (VOC) [1]. The rapid dominance of Delta is attributed to its high infectiousness, and is predicted to have a basic reproduction number, R₀, of 3.2 to 8, with a mean of 5.08, as opposed to Alpha, which was estimated to be approximately 4 [3]. A recent study confirmed that the viral load of individuals infected with Delta is significantly higher than Alpha [4].

Vaccination roll-out in Israel of the COVID-19 vaccines started on December 20, 2020. By the end of February 2021, over a third of the whole population and over 80% of adults over the age of sixty, were already fully vaccinated by the BNT162b2 vaccine developed by BioNTech and Pfizer [5]. This successfully conducted vaccination campaign resulted in a drop in cases in Israel from over 1,000 per million per day to less than 3 cases per million per day in May 2021. By May, Israel was thought to have reached the herd immunity threshold (HIT), which prompted the lifting of non-pharmaceutical interventions (NPIs), including halt of mandated masking by the end of that month.

Nevertheless, a new outbreak began in mid-June, first, with dozens of infections coming from abroad, triggering some outbreaks within local communities that went on to spread across the country. By the end of August, Israel went back to over 1,000 confirmed cases per million per day. This outbreak, from its very beginning, was completely dominated by infections with the Delta variant. Additionally, since the beginning of the outbreak, a high increase in breakthrough infections in vaccinated individuals was reported. Later, the waning of immunity was confirmed in multiple studies worldwide, showing that the BNT162b2 vaccine effectiveness (VE) is reduced to below 50% after five months of the second dose [6–9]. Finally, the Israeli government reinstated some NPIs in late June, including indoor masking, green passports and mandatory isolations.

What caused this surge to start in Israel in June 2021 or what propelled Israel above the HIT is a hot topic of debate. In the HIT formula, both the R₀ and the fraction of protected individuals are taken into account [10]. Consequently, both the higher infectiousness of Delta and the reduced effectiveness of the vaccine in those vaccinated in January and February may have pushed Israel above the HIT. Using a novel spatial-dynamic model of disease spread, which has shown high success in modeling the spread in Israel throughout the COVID-19 pandemic [11, 12], we compared models of infections that combine an R₀ of Alpha or Delta and decline and no decline in the effectiveness. Our analysis suggests that the current outbreak in Israel is caused by a combination of both the Delta variant being more contagious than the Alpha variant and the decline of VE as time passes. However, we find that even in the presence of the Alpha variant, the decline in VE would have pushed Israel above the HIT. Moreover, we found that the Delta variant itself would push Israel above the HIT even without waning immunity, but with a different age mixture of confirmed cases compared to the scenario in which there is waning in the efficacy but the dominant variant is Alpha.

METHODS

In this work, we use a diffusive Monte-Carlo (MC) model introduced by De-Leon and Pederiva [11, 12] to model the spread of the disease. We used this model to examine the expected effects of vaccinations on both confirmed cases (CC) and severe hospitalizations, as well as to estimate the vaccine’s current efficacy. As opposed to other infection models, such as the Susceptible Infected Removed (SIR) [13–15], this particle model enables us to distinguish between different age groups and treat each one separately, assuming that the infection occurs throughout the population simultaneously. Furthermore, a particle model can be adjusted to adapt the rate of population immunization to accurately examine the different effects of the vaccine on subgroups of the vaccinated population as well as the entire population.

It is necessary to differentiate between R_t (the theoretical replication number of the virus) and R_e (the effective replication number of the virus) in order to estimate the spread of a pandemic in the presence of vaccines. R_t estimates the number of encounters between carriers and healthy individuals that would have resulted in an infection if the vaccines had not been available and is defined by how contagious the current variant is, whereas, R_e is affected by vaccination rates and protection offered by the vaccines. In this work, based on the easing of social restrictions in Israel in February 2021 and Delta’s increasing dominance since June, we estimate the theoretical R_t in Israel in July 2021 at 3 for the Delta variant and 2.2 for the Alpha variant (we perform the modelling on additional R_t levels). It is important to note, that in contrast to R₀, the basic reproduction rate, R_t takes into account the NPIs that are still customary in Israel, especially the duty of isolation for verified patients and those who have been exposed to a verified patient and the duty of wearing masks indoors.

Recent studies from Israel, Qatar, UK and [6, 9, 16–18] have shown reduction in effectiveness in preventing COVID-19 infection and transmission blocking in individuals vaccinated more than 5 months before exposure. While the exact effectiveness after this waning humoral immunity, reports from Israel have estimated the effectiveness of the vaccine against infection to be between 20%-50% for those vaccinated six months ago, compared to 80%-90% for those vaccinated recently.

Data on deaths and vaccination rates was download on September 1st from OurWorldInData.org [?]. Only countries with a population over 1 million and with available information are presented.

RESULTS

After a wide vaccination campaign in Israel, the number of new confirmed and hospitalized cases declined from mid-January 2021, first among adults over 60 years of age, and later in the entire population [19, 20]. After almost complete elimination of cases, a new outbreak began in mid-June 2021. We applied our spatial-dynamic MC algorithm with information on current vaccination rates, COVID-19 recoveries, and the number of confirmed cases in Israel from July 2020 to June 2021 to model difference scenarios of the expected daily confirmed cases in the outbreak that started in June 2021 (Figure 1). We first ran our model to estimate cases in the presence of the Delta variant and in time-dependent decline of VE. This model performed well in predicting the real-world counts of CC (normalized mean absolute difference = 2%). We next ran a model with a theoretical reproduction rate (R_t) that corresponds to the Alpha variant, and also without a time-dependent decline in VE. This model suggested that if both factors have not happened, Israel would still be below the HIT, and the June outbreak could have been contained by the vaccinations alone. We next tested two additional models: a model with Delta but no time-dependent decline in VE (Delta model), and a model with Alpha and with a time-dependent decline in effectiveness (Waning model). These models show that both factors - more infectious variant and waning immunity - were able to push Israel above the HIT independently. In both models, the number of CC would have grown steadily up to approximately 2000 daily cases in the Waning model and 3,500 daily cases in the Delta model by mid-August. According to our analyses, with decline in VE, an R_t below 2 could have allowed to contain this outbreak without additional NPIs, or at least delay it (Supplementary Figure 1).

• Download figure
• Open in new tab

FIG. 1: Predicted daily cases in four scenarios.

The daily new cases (a, total population, b 60+). Solid blue band: the delta variant with time-decline of the efficacy against infection to 20%-50% 150 days after the second dose; Dashed green band: alpha variant with time-decline of the VE against infection to 20%-50% 150 days after the second dose; Dotted red band: the delta variant without a time-decline of the VE against infection; Dotted-dashed magenta band: the alpha variant without a time-decline of the VE against infection. Dots: Observed data (7-day moving average)

Interestingly, due to the fact that in the Waning model the reduced VE is expected first in the elderly population who were vaccinated first, the number of CC in 60 year old and over individuals is predicted to be higher in the Waning model compared to the Delta model (250 cases vs. 150 cases in mid-August) (Figure 1b). In the Delta+Waning model, this number grows to 700 cases (in the real-world there were 716 cases). By examining the age distribution of cases at the beginning of the outbreak, we found that our prediction matched what actually occurred. While during the winter outbreak the fraction of cases among adults 60 and over was significantly lower than their fraction in the total population (10% of cases vs. 15% of the total population, One-sample proportional test p-value = 0.009) due to the higher vaccination rates in older populations, in the June outbreak we observed an increased percentage in the adult population that were affected (17% vs 15%, one-sample proportional test p-value = 0.7), more or less as expected from random spread in the population (Supplementary Figure 2).

Since the Delta variant became the dominant variant in most countries by July 2021, and waning immunity was observed in other vaccinations, it is anticipated that countries with low vaccination rates or that have vaccinated their population early, will experience a spike in infections in August. On the other hand, countries that have only recently vaccinated a high fraction of their population will be able to contain a Delta outbreak. By analyzing the fraction of newly vaccinated individuals and the number of COVID-19 deaths in 50 countries around the world, we observed a significant negative correlation between recent vaccination rate and COVID-19 deaths (R= -0.47, p-value = 0.0005) (Figure 2).

• Download figure
• Open in new tab

FIG. 2: Effect of vaccinations on controlling Delta.

The scatter plot shows the percent of population vaccinated in April to July 2021, which are expected to be effective in August vs. the number of deaths per 1 million individuals in each country. Only countries with at population above 1 million people are presented. Size of dots is the percent of population vaccinated before April, expected to be less effective in August. D.ata was collected by [2].

DISCUSSION

In this study, we used a spatio-dynamic Monte Carlo algorithm to estimate the effect of the Delta variant on the COVID-19 outbreak that started in Israel in June 2021. We compared four different models with a combination of varying VE and theoretical reproduction rate. Based on the analysis, both the introduction of Delta and its higher reproduction rate as well as a time-decline in VE can explain this outbreak. Interestingly, we observe that the time-decline VE is responsible for the fact that the first to be infected are those vaccinated first, which are most prominently older populations, who are also at higher risk for severe disease. It is therefore recommended that counties that have vaccinated their elderly population in recent months monitor not just the number of confirmed cases, but also their age distribution and immunization status, and be aware of when this population will lose some of its immune defenses. Such monitoring will allow to detect in real-time the protective attenuation of the vaccine and provide this population another dose of the vaccine on time, which according to recent studies, has been found to be very effective against the delta variant [21, 22].

It is important to keep in mind that our models attempt to predict events that did not appear in reality. In the real world all the parameters of the model are inter-dependent. Increased reproduction rates might cause the government to impose restrictions sooner and cause a reduction in social interactions. A high reproduction rate may also result in more people getting vaccinated. As such, our analysis focuses on factors contributing to the June outbreak in Israel and its global implications, not on the prediction of cases in alternative scenarios.

Interestingly, when we examine the effect of R_t on the spread of the pandemic (Supplementary Figure 1) we find that in case of decline in the effectiveness of the vaccine against infection (which is the case today) the theoretical reproduction number, R_t should be relatively low (R_t < 2), in order not to push Israel above the HIT, which could be the case with a less contagious variant and/or presence of effective NPIs.

In summary, we argue that both the Delta variant higher infectiousness and the reduction of effectiveness of the vaccines contributed to pushing Israel above the HIT. If this is the case, the booster shots provided to all Israelis since August, are expected to reduce severe cases and the burden on hospitals, but not to bring Israel back to herd immunity without reinstating effective NPIs.

SUPPLEMENTARY FIGURES

• Download figure
• Open in new tab

FIG. S.1: Our predictions for the total number of confirmed cases for August 15.

Dots: The total number of confirmed cases in the presence of decline in the VE. Squares:The total number of confirmed cases without a decline in the VE. The errorbars are due to the uncertainty about the effectiveness of the vaccine and the numerical limitations of the calculation.

• Download figure
• Open in new tab

FIG. S.2:

. Bars: Percentage of verified vaccinated persons aged 60+ among all verified persons as of Merch 2021;Solid line: The daily total confirmed cases.(all the data was taken from Refs. [5, 23]).