Modelling COVID-19 Vaccine Breakthrough Infections in Highly Vaccinated Israel – the effects of waning immunity and third vaccination dose

Abstract

In August 2021, a major wave of the SARS-CoV-2 Delta variant erupted in the highly vaccinated population of Israel. The Delta variant has a transmission advantage over the Alpha variant, and thus replaced it in approximately two months. The outbreak led to an unexpectedly large proportion of breakthrough infections (BTI)-- a phenomenon that received worldwide attention. The BTI proportion amongst cases in the age group of 60+ years reached levels as high as ∼85% in August 2021. Most of the Israeli population, especially those 60+ age, received their second dose of the vaccination, four months before the invasion of the Delta variant. Hence, either the vaccine induced immunity dropped significantly or the Delta variant possesses immunity escaping abilities. In this work, we analyzed and model age-structured cases, vaccination coverage, and vaccine BTI data obtained from the Israeli Ministry of Health, to help understand the epidemiological factors involved in the outbreak. We propose a mathematical model which captures a multitude of factors, including age structure, the time varying vaccine efficacy, time varying transmission rate, BTIs, reduced susceptibility and infectivity of vaccinated individuals, protection duration of the vaccine induced immunity, and the vaccine distribution. We fitted our model to the cases among vaccinated and unvaccinated, for <60 and 60+ age groups, to address the aforementioned factors. We found that the transmission rate was driven by multiple factors including the invasion of Delta variant and the mitigation measures. Through a model reconstruction of the reproductive number R₀(t), it was found that the peak transmission rate of the Delta variant was 1.96 times larger than the previous Alpha variant. The model estimated that the vaccine efficacy dropped significantly from >90% to ∼40% over 6 months, and that the immunity protection duration has a peaked Gamma distribution (rather than exponential). We further performed model simulations quantifying the important role of the third vaccination booster dose in reducing the levels of breakthrough infections. This allowed us to explore “what if” scenarios should the booster not have been rolled out. Application of this framework upon invasion of new pathogens, or variants of concern, can help elucidate important factors in the outbreak dynamics and highlight potential routes of action to mitigate their spread.