First indication of the effect of COVID-19 vaccinations on the course of the COVID-19 outbreak in Israel

Concomitantly with rolling out its rapid COVID-19 vaccine program, Israel is experiencing its third, and so far largest, surge in morbidity. We aimed to estimate whether the high vaccine coverage among individuals aged over 60 years old creates an observable change in disease dynamics. Using observed and simulated data, we suggest that the shape of the outbreak as measured by daily new moderate and severe cases, and in particular of patients aged over 60, has changed because of vaccination, bringing the decline in new moderate and severe cases earlier than expected, by about a week. Our analyses is consistent with the assumption that vaccination lead to higher than 50% protection in preventing clinical disease and with at least some effectiveness in blocking transmission of elderly population, and supports the importance of prioritizing vulnerable population . This is the first indication of the effectivity of COVID-19 vaccine in changing the course of an ongoing pandemic outbreak.

vaccine and prioritizing individuals aged 60+ years, immune deficient patients, and healthcare workers to be the first vaccinated. Since then, and up to date (January 25, 2021), about 2.6 million Israelis (29% of the population) have been vaccinated with at least the first dose, and the second dose is rolling out, yielding vaccine coverage of 80% in the >60 years age group, as published by the Israeli Ministry of Health (4). Thus, monitoring epidemic evolution in Israel may allow assessment of the effectiveness of vaccines in mitigating the outbreak on a national level. A major challenge in this respect is to identify a distinct indicator of the effect of vaccines differentiating it from the effect of lockdown.
We used data analysis tools and simulations to model epidemic dynamics in Israel, comparing different scenarios of lockdown duration and effectiveness. Each scenario is modeled with and without the up-to-date vaccine coverage. We studied the daily trends of new confirmed cases as well as the trends of new moderate or severe COVID-19 cases in Israel. Moderate COVID-19 is defined as a case positive for SARS-COV-2 in RT-PCR test, with bilateral pneumonia. A severe case of COVID-19 is defined in Israel as a SARS-COV-2 confirmed case with oxygen saturation below 93% while breathing ambient air (these include also critical cases). Therefore, as long as hospitals are not overwhelmed, numbers of moderate and severe cases are independent of either testing policy or availability. In addition, these cases are at increased risk for further deterioration, with 15-22% fatality rate within 10 days (See Methods Section B). On average, moderate and severe admissions take place about 5 days following a positive RT-PCR test, and an estimated 10 days from the actual infection. Thus, the number of new moderate and severe cases makes a clear and robust indication for the dynamics of disease, that has the advantage of being evident quite soon after infection. Finally, elderly individuals are at increased risk for having moderate or severe COVID-19. Therefore, we studied the trends of moderate and severe cases as candidate measures of disease dynamics and the effect of vaccine on the outbreak.
We start by analyzing the status of the pandemic in Israel, using the effective reproduction number, , which is approximated using the 4-day growth factor of the weekly average of new COVID-19 cases (5). In Fig-1(a) we show that the dynamics leading to the January lockdown, as described by changes in Re are quantitatively similar to those leading to the second lockdown in Israel, which was implemented on September 19, 2020. The effective , based on the daily admissions of moderate or severe cases shifted by five days (in agreement with clinical course), are highly correlated with the effective based on confirmed cases and show similar correlation with the reproductive numbers shown in September (Fig S-2). The agreement of reproductive numbers based on different measures supports the notion of similar dynamics in both waves, which is not surprising considering the similarities in the non-pharmaceutical intervention restrictions during the lockdowns. It also minimizes the concerns of systematic errors, such as under-detection of cases. At the onset of both nationwide lockdowns, the effective ranged 1.1-1.3.
This similarity is observed also in Fig-1(b), where the absolute number of the 7-day moving average of the daily-confirmed cases, shows the same dynamics starting about three weeks prior to lockdowns, albeit the third wave is about 35% larger, due to late entrance to the January lockdown. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 3, 2021. ;  One sees that while the evolution of the outbreaks is extremely similar, the number of cases in January is bigger by 35%, than the corresponding date in September. Fig-2 shows that similarities between the second and third lockdowns discontinue. While the number of new moderate and severe cases among patients older than 60 has peaked on the 6 th day of the third lockdown and decreased thereafter, in the second lockdown it reached its peak on the 14 th day of the second lockdown, and decreased thereafter. Such a time difference is not found when studying the daily admissions of moderate or severe cases among younger patients. In fact, as of today (January 25, 2021), the current outbreak has not seen a decrease in the latter daily number, as can be seen in Fig. 3 is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 3, 2021. ; We hypothesize that having an early peak and a decline is a result of vaccinations, and suggest the time difference between peaks with and without vaccination as an indicator for vaccination effectiveness. In order to test this hypothesis, we employed homogenous random movement models as introduced by De-Leon and Pederiva (6) to study the expected effects of vaccinations on the dynamics of the moderate and severe hospitalizations during the January lockdown. The models take as input the effective repruduction number as estimated by the detected cases, and use them to model the dynamics of infections, and the resulting moderate and severe cases. Since the estimates of the infection level include some uncertainties, e.g., due to testing policies, we use various scenarios to estimate the lockdown effect, ranging from 2-5 weeks of length, and assuming either prompt or delayed compliance of the population leading to prompt or delayed effect on . We additionally assume that the third lockdown yields similar effects as second lockdown, i.e., reaching of 0.7 (see Fig.1-a). The details of the scenarios are given in Section A of the Methods. We employed models with and without vaccines; the models with vaccines assume coverage according to vaccine rollout in Israel. The effectiveness of vaccines in preventing symptomatic disease was taken as 0% in the first week after the first dose, 30% on days 7-10, 60% on days 11-20 and 95% from day 21 (second dose) onwards. We first assume that vaccinations prevent infections, though we check the sensitivity of the results to this assumption later on, as this was not reported in the Pfizer-BioNTech trial. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 3, 2021. ; https://doi.org/10.1101/2021.02.02.21250630 doi: medRxiv preprint as shown in Fig.3-a, are expected to peak on similar timing in all scenarios, though the rate of the following decline changes according to the scenario, and more importantly, according to whether vaccines were administered, with a steeper decline when the population is vaccinated. This is expected, under the assumption that vaccinations prevent transmission.
An important observation appears when studying the predicted numbers of moderate and severe cases according to these scenarios. In Fig. 3-b, a time difference of 6-11 days is predicted between the peaks of new daily moderate and severe cases of vaccinated and unvaccinated scenarios. Stratifying the models by age (Fig 3-c and 3-d), it is evident that this bringing forward of the decline in moderate and severe cases is limited to patients older than 60, while it is not expected to affect the timing of peak and decline in those younger than 60. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 3, 2021. ; https://doi.org/10.1101/2021.02.02.21250630 doi: medRxiv preprint As aforementioned, there are still open questions regarding the effectiveness of vaccines, in particular, whether the 95% efficacy achieved in clinical trials translates to similar effectiveness in preventing clinical disease in a national vaccination campaign. In addition, the effectiveness in preventing transmission is unknown. We therefore employed simulation models with sensitivity analyses (Fig 4). Trivially, the confirmed cases are not expected to be affected by vaccinations that do not prevent infection ( fig.4-a). However, we see that the indicator we suggested here, i.e., the time difference in the peaks in the new daily moderate and severe cases between vaccinated and unvaccinated scenarios is shortened by about 1-2 days, if the effectiveness of the vaccinations is 95%, even if transmission is not prevented by the vaccines. In addition, we find that the indicator changes by about a day if the 95% effectiveness is reached only on day 28 after the first dose.
The indicator loses its significance, when assuming vaccine effectiveness decrease to 50% in the latter scenario. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted February 3, 2021. ; previous wave when no vaccines were available. Both observed data and models suggest that indeed, prioritizing to vaccinate the elderly was effective in reducing the influx of moderate and severe cases, and bringing the tipping point earlier than expected, by about a week, has there not been any vaccines. Since length of stay is relatively long (10-20 days), the change in admitted patients is not detected immediately in the total numbers of hospitalized patients.
Previous modeling of vaccines looked into a possible impact of vaccines on the outbreak in the US. According to this study, vaccine coverage that will reach 40% in the American population within nearly a year will mitigate morbidity effectively, even though this is far from reaching herd immunity (7). Nevertheless, the modeling ignores current trends in the dynamic of COVID-19 in the US, except for the cumulative incidence of the disease so far. Moreover, this study assumes that vaccination completely blocks transmission, and that asymptomatic disease is not associated with transmission at all. Another modelling study that compared different strategies in prioritizing vaccines (8), suggested that prioritizing individuals aged 60+ years (direct protection) is more effective when transmission is high ( > 1.5), whereas prioritizing younger individuals is more effective when transmission is low ( < 1.2). In all scenarios they assumed rollout rate of 0.1-0.76%. In Israel, those aged 60+ were prioritized in the first month of vaccines and are now gradually being offered to younger individuals. The rollout rate achieved in some days 2% of the population (3).
Our study has some limitations. First, our models assume that the third lockdown would have similar effectiveness in curbing the third wave as the second lockdown if no vaccines were employed. It should be noted that previous lockdowns took place around holidays where people usually work less. Whether because of no holiday, "pandemic fatigue", or loss of trust in the government or healthcare systems (9,10), it is quite possible that the public adherence to the current lockdown is reduced, impairing its effectiveness. Second, our models assume equal chance encounters between individuals. We could not model social interactions, as there is no such data available for Israel. Israel is a country with the highest fertility rate in the OECD countries (11) which translates into big households and crowded classrooms and different social interactions than those observed in European or North American countries. This minimized the possibility to adopt social connections surveys conducted in other countries. Moreover, like any other country, there are sectors in the population that co-interact within the sector than with members of other sectors. This was not taken into account, based on the notion that outbreaks that started within specific population sects (e.g., Israeli Arabs, ultra-orthodox Jews) expanded to the general public and to national outbreaks, leading both to the second and third waves.
To the best of our knowledge, this is the first study to estimate the effect of mass vaccinations in a tome of a pandemic using a new indicator. While our study cannot accurately estimate the effectiveness of the Pfizer-Biontech vaccine, it suggests that the effectiveness is greater than 50% and that there is a considerable level of prevention of transmission following vaccinations.