Super-Spreaders Out, Super-Spreading In: The Effects of Infectiousness Heterogeneity and Lockdowns on Herd Immunity =================================================================================================================== * Jhonatan Tavori * Hanoch Levy ## Abstract Recently, [8] has proposed that heterogeneity of infectiousness (and susceptibility) across individuals in infectious diseases, plays a major role in affecting the Herd Immunity Thresh-old (HIT). Such heterogeneity has been observed in COVID-19 and is recognized as *overdis-persion* (or “super-spreading”). The model of [8] suggests that super-spreaders contribute significantly to the effective reproduction factor, *R*, and that they are likely to get infected and immune early in the process. Consequently, under *R* ≈ 3 (attributed to COVID-19), the Herd Immunity Threshold (HIT) is as low as 5%, in contrast to 67% according to the traditional models [1, 2, 4, 10]. This work follows up on [8] and proposes that heterogeneity of infectiousness (susceptibility) has two “faces” whose mix affects dramatically the HIT: (1) *Personal-Trait-*, and (2) *Event-Based-*Infectiousness (Susceptibility). The former is a personal trait of specific individuals (*super-spreaders*) and is nullified once those individuals are immune (as in [8]). The latter is event-based (e.g cultural *super-spreading* events) and remains effective throughout the process, even after the super-spreaders immune. We extend [8]’s model to account for these two factors, analyze it and conclude that the HIT is very sensitive to the mix between (1) and (2), and under *R* ≈ 3 it can vary between 5% and 67%. Preliminary data from COVID-19 suggests that herd immunity is not reached at 5%. We address operational aspects and analyze the effects of lockdown strategies on the spread of a disease. We find that herd immunity (and HIT) is very sensitive to the lock-down type. While some lockdowns affect positively the disease blocking and increase herd immunity, others have *adverse* effects and *reduce* the herd immunity. ## 1 Introduction In 1923 Topley and Wilson described experimental epidemics in which the rising prevalence of immune individuals would end an epidemic. They named this phenomenon as *“Herd-Immunity”* [14]. Once the Herd Immunity Threshold (HIT, measured in fractions of the population that got immune) is surpassed, then the effective reproduction number, *R*, reduces below 1 and the number of infection cases decreases. The exact value of this threshold is an important measure used in infectious disease control and immunization, and its estimation for the COVID-19 disease are used by governments worldwide in determining policies to fight against the current pandemic. A recent study of Oz, Rubinstein and Safra [8] proposed a new model for the spreading of infectious diseases, such as COVID-19. Under that model and the assumption that the basic reproduction number, *R*, of COVID-19 is approximately 3, the HIT is approximately 5%, namely, when 5% of the population is infected herd immunity is reached. This estimation was in contrast to the allegedly “axiomatic” cutpoint of *HIT* ≈ 67% assumed for COVID-19 [10]. Preliminary data from COVID-19 suggests that herd immunity is not reached at the approximate 5% fraction expected by [8]. For Example, the US states of North Dakota, South Dakota, Iowa, Utah and Tennessee with at least 12%, 11.5%, 9.5%, 9.5% and 9% infected, respectively. Four of them have current *R >* 1 [12, 15]. The model of [8] is based on the observation [9, 11, 13] that the epidemic spread network is not homogeneous, where distinct individuals are infectious (likely to infect others) and susceptible (likely to become infected themselves) in various degrees. The heterogeneity of these values among individuals is recognized as *overdispersion* or *super-spreading* (*super-spreaders* are a class of individuals whose secondary infection rate is very high [5]). The estimates for the COVID-19 pandemic fits this property and asserts that between 5% to 10% of the infected individuals cause 80% of the secondary infections [3, 7]. Furthermore, a correlation between the infectiousness and susceptibility of each individual have a drastic effect on the over-time reduction of the effective reproduction number, *R*, under the spreading model of [8]. The heterogeneity and correlation of these parameters yields that the “super-spreaders” are extremely likely to get infected and develop immunity in an early stage of the pandemic process. [8]’s 5% estimation of the percentage of the population that contract the disease before herd immunity is reached was based on these properties. In this work, we follow-up on [8] regarding the effect of infectiousness heterogeneity on the HIT. However, we propose that infectiousness (and susceptibility) should be classified into two inherently different types: (1) *Personal-Trait Infectiousness (Susceptibility)*, and (2) *Event-Based Infectiousness (Susceptibility)*. We will refer to the combination of the infectiousness and the susceptibility as *Spreading*. The first type stems from traits of an i ndividual. The second type relates to social events in which every individual may participate, regardless of its personal traits. To demonstrate these two types, consider, for example, a package-delivery person and compare it with an academic researcher. During a single day, the delivery person has interactions with tens or hundreds of people, and therefore has high personal-trait infectiousness. On the other hand, the researcher may work most of the time in his/her office or interact with a small research group and as a result has a lower personal-trait infectiousness. Yet, both of them may participate in a social-gathering event (such as a concert, a wedding, or “just” a family birthday party). During such an event, both have approximately the same amount of interaction (which may be quite large), and therefore have the same event-based spreading degree, regardless of their personal traits. It is important to note that personal-trait infectiousness relate not only to the social behaviour of the individual; it may also relate to his/hers biological properties (e.g., his/her body reproduces a virus faster and therefore he/she is more infectious). We assume that the likelihood parameters of each individual consist each of the sum of two parameters. (1) *S**p*(*a*) and *I**p*(*a*) which are the personal-trait susceptibility and infectiousness parameters of *a*, respectively. Those values reflect p ersonal t raits a nd a re d rawn o nce (pandemic beginning) and remain constant throughout, exactly as in [8]. (2) *S**e* and *I**e* which are the event-based (cultural) susceptibility and infectiousness parameters. Those values reflect occasional event-based spreading which is probabilistically redrawn for each individual at every step of the pandemic. We assume that the likelihood of *a* to be infected is: ![Formula][1] and the likelihood of *a* to infect others is ![Formula][2] where *p* determines the weight of each infection type. The symmetry between S(a) and I(a) is similar to that of [8] and stems from the assumption that infectiousness level and susceptibility level are proportional to the level of interaction the individual has with others or to its biological properties. We use this model to analyze the progression of an infectious disease and the value of *R*(*n*), the effective reproduction number, as a function of the fraction of population infected. We show that the contribution of the personal-trait spreading drops sharply at early stages of the disease, as the super-spreaders are likely to contract the disease and develop immune at early stages of the process. On the other hand, the contribution of the event-based spreading drops much more slowly, and is affected very little at early stages, as its reduction is proportional to the decrease of the susceptible population size. In other words – when super-spreaders are out (event-based) super-spreading is in. Hence, *R*(*n*) may remain at high values even after the super-spreaders population is fully immune, as opposed to [8]. This results in a slower decay of *R*(*n*) and leads to a higher value of the Herd Immunity Threshold. We show that for COVID-19 the Herd Immunity Threshold depends on the *mix* (*p*) between the weights of the personal-trait and event-based spreading. In particular, we prove that the *HIT* ≈ 5% estimate of [8] holds when assuming only personal-trait spreading, and that the traditional prediction of *HIT* ≈ 67% holds when assuming only event-based spreading. Having established a formula expressing *R*(*n*), we address operational aspects and analyze the effects of lockdowns on the Herd Immunity Threshold. Lockdowns, of a variety of variants, have been enforced worldwide in order to fight the COVID-19 pandemic. While lockdowns might have immediate impact such as collapse of the effective reproduction and suppression of infections and mortality, they have long-term impact as well. We discuss two different lockdown policies: (1) An Event-based spreading targeted lockdown (e.g shut down of cultural events). (2) A Personal-trait spreading targeted lockdown (e.g. restricting daily/professional activities). We analyze the effect of these policies on the composition of the infected population, on the effective reproduction number, and on the Herd Immunity Threshold (HIT). We show that while a lockdown (or a sequence of lockdowns) which is targeted at event-based spreading reduces the disease spread by decreasing the HIT, a lockdown (or a sequence of lockdowns) which is targeted at personal-trait spreading will act adversely and will *increase* the disease spread by increasing the HIT. The rest of the paper in organized as follows: In Section 2 we formally describe our model, and present the effective reproduction number which plays a major role in the analysis. In Section 3 we develop an expression for *R*(*n*), and calculate the Herd Immunity Threshold for a general-case disease. We examine the result in a numerical discussion based on COVID-19 spreading distributions. Then, in Section 4, we study various lockdown policies, and analyze their effect on the Herd Immunity Threshold. Finally, concluding remarks are given in Section 5. ## 2 The Disease Spread Model In this section, we present our model for the spreading of infectious disease accounting for heterogeneity of infectiousness/susceptibility. We extend the model of [8] and propose that there exists an event-based infection factor in addition to the factors described in their model. Our analysis begins with a certain number of infected individuals. We measure the spread of the disease as a function of the number of individuals who got infected. Specifically, we index the individuals by the order they are infected and have *R*(*n*) denote the effective reproduction number associated with the *n*th infected individual. Namely, the event whereby the *n*th individual gets infected is the *n*th event (or step *n*). We use *n* also to denote the *step of the disease*. Measuring the spread of the disease as a function of the infected population size will be useful in deriving the Herd Immunity Threshold (HIT) of the disease, namely the fraction of the population that gets infected prior to reaching *R*(*n*) ≤ 1. ### 2.1 Susceptibility and Infectiousness We follow the model of [8] and assign to each individual *a* a personal-trait-susceptibility parameter *S**p*(*a*) and a personal-trait-infectiousness parameter *I**p*(*a*) drawn from some probability distributions. Those values quantify how likely *a* is to be infected and infect others, respectively, according to its *personal* traits. The values of *S**p*(*a*) and *I**p*(*a*) accompany *a* throughout the entire progress of the disease, and remain at the same values. We follow [8] and define the average conditional infectiousness *φ* (*s*). In our case, it is logical to parametrize *φ* (*s*) only by the personal-trait susceptibility and infectiousness: ![Formula][3] As was discussed, the heterogeneity of the spreading values of the population will play a major role in our analysis. Hence, we will measure: 1 ![Formula][4] where *H**n* is the healthy population at step *n*. In addition, and beyond the model of [8], we assign an event-based infectiousness parameter and event-based susceptibility parameter to each individual. Those values are subject to change through the progress of the disease. At step *i* we assign to *a* ![Graphic][5], the event-based-susceptibility parameter and ![Graphic][6], the event-based-infectiousness parameter, both are random variables whose realizations hold only for iteration *i*. Since the values of *S**e*, *I**e* measure the event-based-spreading of the society (assigned to its individuals at a given time), they are drawn from probability distributions that are common for the entire population, denote them by Λ*S* and Λ*I*. The *susceptibility* of *a* at step *i*, which is the likelihood of *a* to be infected, is: ![Formula][7] and the *infectiousness* of *a*, which is the likelihood of *a* to infect others is ![Formula][8] where 0 *< p <* 1 and *q* = 1−*p*. The value of *p* determines the *mix* between the spreading types. We call *p* (and respectively, *q*) the weights of the personal-trait (and respectively, event-based) spreading of the disease. As will be seen later, the value of *p* will have a drastic effect on the Herd Immunity Threshold. Note that the special case where *p* = 1 gives exactly the model of [8]. Under this model, the probability that *a* will be infected at step *n*, assuming that *a* was healthy at step *n* − 1 is: ![Formula][9] ### 2.2 Basic and Effective Reproduction Number The *basic reproduction number, R*, is a measure of how transferable a disease is. It is defined as the expected number of secondary cases produced by a single (typical) infection in a completely susceptible population (whose size is *N*). In reality, varying proportions of the population are immune to any given disease at any given time. Hence, as in [8], we will measure the *effective reproduction number, R*(*n*), which is defined as the expected number of infections directly generated by the *n*th infected individual. ![Formula][10] where the expectation is taken over the *n*th individual to be infected. As in [8], Eq. (6) equals to2 ![Formula][11] this can be approximated by: ![Formula][12] Using Eq. (3), (4) and (7) we have: ![Formula][13] where *N* (*n*) is the size of the healthy population at step *n* and *λ**I*, *λ**S* are the means of Λ*I*, Λ*S*, respectively. ## 3 Reaching Herd Immunity In this section, we analyze the changes in the composition of the population through the spread of the disease, and the decrease of *R*(*n*) as the fraction of the population that contracted with the disease increases. We prove the following theorem: Theorem 3.1 (General Case Herd Immunity Threshold). *For any δ when* ![Formula][14] *fraction of the population is infected, the effective reproduction number, R*(), *will be reduced by a factor of* ![Formula][15] *relatively to the basic reproduction number, R*. *The threshold for herd immunity is when the value of the effective reproduction number is* 1. Having this expression for the change in the effective reproduction number for a general distribution, we will later use the special case of the Gamma distribution with estimated parameters for COVID-19 [3, 5, 6] and inspect the HIT values for different *p* values. ### 3.1 Proof of Theorem 3.1 for A General Spreading Distribution In order to prove Theorem 3.1, we establish the following claim: **Claim 3.2** (The likelihood of an individual to be infected). *For a person a*, ![Formula][16] *Where* ![Formula][17] *Proof of Claim 3.2*. The proof follows the proof of Claim I provided in [8] with modifications required for our extended model. The proof is based on using Eq. (5) and obtaining: ![Formula][18] The rest of the proof consists of algebraic manipulations of (13) and the full details are given in Appendix A. □ We next establish two supporting lemmas, corresponding to equations (3.6) and (3.7) in [8], and conclude with the proof of Theorem 3.1. Lemma 3.3 (Heterogeneity of the population during the process). *For any s* ∈ Supp(*S**p*), ![Formula][19] *Proof of Lemma 3.3*. By definition ![Formula][20] Using Eq. (11), ![Formula][21] . □ Lemma 3.4 (The size of the susceptible population). *For any n* ∈ [*N*], ![Formula][22] *where N* (*n*) *is the size of the susceptible population at step n, and N* *is the total size of the population*. *Proof of Lemma 3.4*. By definition, ![Formula][23] Using Eq. 11, we have Eq. (15). □ *Proof of Theorem 3.1*. Denote ![Formula][24] Using Eq. (8) we know that: ![Formula][25] We develop the ratio: ![Formula][26] According to Lemma 3.3 and Lemma 3.4 we can replace the values of *ρ*(*s, n*) and the ration *N* (*n*)*/N* and have: ![Formula][27] Since the value of expression (9) is given by: ![Formula][28] we have that Eq. (10) holds.. □ ### 3.2 COVID-19: Analysis and Discussion We move to demonstrate the results of Theorem 3.1 on a Gamma distribution with shape and scale parameters *k* and *θ*, respectively. The Gamma distribution was previously attributed to the infectiousness of COVID-2 [5]. We substitute the estimates for COVID-19: *R* ≈ 3 and *k* ≈ 0.1. [3, 6, 8]. We assume that the personal-trait-susceptibility and personal-trait-infectiousness of the population are highly correlated. I.e., we set *φ* (*s*) = *s*. This stems from assuming that both personal-trait infectiousness and susceptibility levels are correlated to social interaction levels of the individual or to its biological properties. In Figure 1 we demonstrate the decay of the effective reproduction number, *R*(*n*), and its contributing factors, classified by their spreading types. We plot the following values as a function of the fraction of the infected population. ![Formula][29] ![Formula][30] ![Formula][31] ![Formula][32] ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/26/2021.01.23.21250242/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/01/26/2021.01.23.21250242/F1) Figure 1: The over-time reduction in the effective reproduction number, *R*(*n*), and its contributing factors as a function of *n*, assuming *p* = 0.5, *k* = 0.1 and *R* = 3. Note that *n* (horizontal-axis) is normalized to percentage. In red -*R*(*n*) (Eq. (17)); In blue - *R**p*(*n*) (Eq. (18)); In Green - *R**e*(*n*) (Eq. (19)); In yellow - *R**mix*(*n*) (Eq. (20)); The blue curve in Figure 1 depicts the contribution of the personal-trait spreading to *R*(*n*), while the green curve depicts the contribution of the event-based spreading. As can be seen, the contribution of the personal-trait spreading drops at early stages of the disease. This results from the assumption that infectiousness is positively correlated with susceptibility, since if more infectious people are also more susceptible, then they have higher probability to be infected and develop natural immunity much sooner (than the less infectious individuals). On the other hand, the contribution of the event-based spreading is affected very little at early stages. Its reduction is proportional to the decrease of the susceptible population, which is linear in *n*. The weight of each spreading type, which is determined by *p*, determines the combined behaviour of *R*(*n*). In Figure 2 we plot the Herd Immunity Threshold as a function of *p*. This demonstrates the effect of the weight of the personal-trait spreading on the decay of *R*(*n*). ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/26/2021.01.23.21250242/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/01/26/2021.01.23.21250242/F2) Figure 2: The Herd-Immunity Threshold (HIT) as a function of *p* assuming Gamma distribution with shape parameter *k* = 0.1. In green - *R* = 3. In red - *R* = 9. When assuming only event-based-spreading (i.e., *p* = 0) the HIT is approximately ![Formula][33] i.e., 67% when *R* ≈ 3 and 90% when *R* ≈ 9. This follows the classical models [1, 2, 4, 10]. On the other hand, assuming only personal-trait-spreading (i.e., *p* = 1) the HIT is approximately 5% when *R* ≈ 3 and 10% when *R* ≈ 9. The 5% matches the expected threshold given by [8]. In reality, the two types of spreading contribute to infections and hence 0 *< p <* 1; In order to predict the Herd Immunity Threshold, one has to estimate the value of *p*, given a population. In Figure 3 we plot the effective reproduction number, *R*(*n*), throughout the spread of the disease for a number of *p* values. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/26/2021.01.23.21250242/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/01/26/2021.01.23.21250242/F3) Figure 3: The effective reproduction number *R*(*n*) as a function of *n* throughout the disease spread for different *p* values where *R* = 3 and *k* = 0.1. Note that *n* (horizontal-axis) is normalized to percentage. Figure 3 and the shapes of the various curves can be used to assist in predicting the value of *p* at *early stages* of a disease. As can be seen, even when a small percentage of the population is infected (*<* 5%), there is a significant difference between the *R*(*n*) values for different *p*’s. Such prediction can assist in predicting the Herd Immunity Threshold. The current data available to us for the current-spreading COVID-19 is biased due to lockdowns of many types, which makes it challenging to estimate *p* and *q*. This estimation problem is left open for forthcoming research. The effect of lockdown policies on the value of *R*(*n*) and the Herd Immunity Threshold is the subject of our next section. ## 4 The Effect of Lockdowns Lockdown strategies, of a variety of variants, have been used worldwide as a major means to fight an epidemy, specifically COVID-19. The question addressed in this section is how such strategies affect the scope of the disease, namely how they affect the Herd Immunity Threshold. We show that, depending on its type, a lockdown may either increase or decrease the size of the population infected prior to reaching herd immunity. In particular, we show that a personal-trait spreading targeted lockdown acts *adversely* on the efforts to reduce the spread of a disease since it increases the threshold. In contrast, an event-based spreading targeted lockdown affect positively the disease blocking and reduces the Herd Immunity Threshold. ### 4.1 Lockdown Policies Large scale physical distancing measures and moving-around restrictions, often referred to as *lockdowns*, can slow disease transmission by limiting contacts between people. These days, many restrictions are used worldwide in order to slow down the spread of the COVID-19 pandemic. We classify the restrictions into two inherently different types: (1) a Personal-trait spreading targeted lockdown, for example by closing or restricting workplaces. (2) an Event-based spreading targeted lockdown, for example by prohibition on cultural events. Formally, these lockdowns will be defined as follows: Definition 4.1 (Personal-Trait Spreading Targeted Lockdown). *During a personal-trait spreading targeted lockdown, the infectiousness and susceptibility of an individual a is:* ![Formula][34] *where q remains the same and p**L* = 0. Definition 4.2 (Event-Based Spreading Targeted Lockdown). *During an event-based spreading targeted lockdown, the infectiousness and susceptibility of an individual a is:* ![Formula][35] *where p remains the same and q**L* = 0. We will use the following notation. Assume that a lockdown starts at the *n**b*th step of the disease and ends at step *n**e*. For any *n* ∈ [*n**b*, *n**e*], let *R**L*(*n*) be the expected value of the effective reproduction number during the lockdown; For any *n > n**e* let *R**L*(*n*) be its expected value after the lockdown is released and *p* (or *q*) returns to its original value. Similarly, let *HIT* be the Herd Immunity Threshold assuming a “natural” spread of the disease (i.e., no lockdown), and *HIT**L* be the threshold assuming a lockdown was performed. Note that the comparison between the natural evolution and the lockdown evolution is based on the number of individuals that contract the disease during those evolutions. I.e., the coupling between the evolutions is done user by user by the order of infection rather than by the time period in which the policies are compared. This coupling method is useful in deriving *HIT* and *HIT**L*. At *n**b* we set *p**L* = 0 or *q**L* = 0. Hence, the value of the effective reproduction number is likely to drop during the lockdown. The question we answer in this section is *if* – and *what* – will be the long-term impact of the lockdown, after it is released (at *n > n**e*). Since at *n**e* the value of *p* (or *q*) will increase, *R**L*(*n*) is most likely to increase as well. Yet, will it pass *R*(*n*)? Will it stay below it? Consequently, what will be the impact on the Herd Immunity Threshed? The next two theorems establish that *any* personal-trait spreading targeted lockdown will increase the HIT while *any* event-based spreading targeted lockdown will decrease the HIT. Theorem 4.3 (HIT of Personal-Trait Spreading Targeted Lockdown). *Assume that a personal-trait spreading targeted lockdown was performed for n* ∈ [*n**b*, *n**e*]. *Then for any n > n**e*, ![Formula][36] *Consequently*, ![Formula][37] Theorem 4.4 (HIT of Event-Based Spreading Targeted Lockdown). *Assume that an event-based spreading targeted lockdown was performed for n* ∈ [*n**b*, *n**e*]. *Then for any n > n**e*, ![Formula][38] *Consequently*, ![Formula][39] Due to the generality of Theorems 4.3 and 4.4, we state the following corollary: Corollary 4.5. *Any sequence of personal-trait spreading targeted lockdowns will result with HIT < HIT**L*. *Similarly, any sequence of event-Based spreading targeted lockdowns will result with HIT > HIT**L*. The following claim will be useful in proving Theorems 4.3 and 4.4. Claim 4.6. *Let S*1, *S*2 *be continuous random variables. Let ρ*1(), *ρ*2() *be their pdfs, and* P1(), P2() *their CDFs, respectively. Let r*() *be an injective monotone function. If for any s*1 *> s*2: ![Formula][40] *then* ![Formula][41] *Proof of Claim 4.6*. The proof is given in Appendix A. □ ### 4.2 Personal-Trait Lockdown: Proving Theorem 4.3 Recall that *ρ*(*s, n*) is the normalized susceptibility distribution (density) at step *n* assuming a “natural” spread of the disease (i.e., no lockdown). Let *ρ**L*(*s, n*) be the normalized susceptibility distribution (density) at step *n* assuming that a lockdown was performed. *Proof of Theorem 4.3*. Resulting from Lemma 3.3 and by the assumption that *p**L* = 0 for *n**b* ≤ *n* ≤ *n**e*, at the end of the (personal-trait) lockdown, i.e. at *n**e*, it holds that ![Formula][42] However, ![Formula][43] Therefore, ![Formula][44] where *c* = ∫ *ρ*(*σ, n**b*) · exp (−*β*(*n**b*, *n**e*) · (*p* · *σ*)) *dσ*. Therefore, for any *s*1 *> s*2: ![Formula][45] Note that: *r*(*s*) = (*p* · *s* + *q* · *λ**S*)2 is an injective monotone function. Hence, by the definition of *R*(*n*) (Eq. (8)) and from Claim 4.6: ![Formula][46] This is demonstrated in Figure 4. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/26/2021.01.23.21250242/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/01/26/2021.01.23.21250242/F4) Figure 4: The expected value of *R**L*(*n*) under a personal-trait spreading targeted lockdown (vs. *R*(*n*) in a natural evolution) assuming *p* = 0.5, *R* = 3 and *k* = 0.1. The lockdown begins at *n**b* = 5% and ends at *n**e* = 30%. More specifically, Since *ρ**L*(*s, n**e*) = *ρ*(*s, n**b*) (Eq. (24)), and by Eq. (8), it holds that ![Formula][47] We move to measure the value of the effective reproduction number after the lockdown ends. Let *n > n**e*, and denote ![Formula][48] Let ![Formula][49] (note that *n* = *n**e* + *x* · (*N* − *n**e*)). *n* and *n*′ are depicted in Figure 5. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/26/2021.01.23.21250242/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2021/01/26/2021.01.23.21250242/F5) Figure 5: The expected value of *R**L*(*n*) under a personal-trait spreading targeted lockdown (vs. *R*(*n*) in a natural evolution) assuming *p* = 0.5, *R* = 3 and *k* = 0.1. The lockdown begins at *n**b* = 5% and ends at *n**e* = 30%. We use *x* = 10%. The dashed purple line demonstrates a linear reduction in the value of *R*(*n*′). As proven, the purple and red curves intersect at *n*. At the end of the lockdown *p* returns to its original value and the disease spreads naturally, with its original parameters. At that point (*n**e*), according to Eq. (24), the distribution of the population is the same as it was before the lockdown began, where only the size of the susceptible population was changed. In other words, *R**L* (red) behave to the right of *n**e* exactly as *R* (blue) behaves to the right of *n**b* (see Figure 5). Therefore, ![Formula][50] This is demonstrated in Figure 5. Combining Eq. (30) with Eq. (29), ![Formula][51] Let us look at ![Formula][52] Hence, ![Formula][53] Note that: ![Formula][54] ![Formula][55] Using Eq. (31) ![Formula][56] In Eq. (32) ![Graphic][57] was calculated using *ρ*(*σ, n*′). Hence, as in Eq. (27) and (28), ![Formula][58] and we have that *R*(*n*) *< R**L*(*n*). This holds for any *n > n**e*, and we conclude that *HIT < HIT**L*. □ Hence, while performing a personal-trait spreading targeted lockdown will reduce the value of the effective reproduction number during the lockdown, it might increase the number of individuals that contract with the disease, prior to reaching herd immunity. In addition, note that even if “herd immunity” is reached during the lockdown (i.e., having *R**L*(*n*) *<* 1 for *n* ∈ [*n**b*, *n**e*]), it might not hold once the lockdown releases as the effective reproduction number might become greater than 1. ### 4.3 Event-Based Lockdown: Proving Theorem 4.4 In the proof of Theorem 4.4 we use the following notation. Let ![Formula][59] I.e., *β*(*n*1, *n*2) = *β*(*n*2) − *β*(*n*1) according to Eq. (12). Let *β**L*(*n*1, *n*2) be defined assuming that a lockdown was performed. We prove the following claim: Claim 4.7. *Under an event-base lockdown, for any n* ≥ *n**e* *and for any s*1 *> s*2, ![Formula][60] *Proof of Claim 4.7*. When the lockdown begins, we have: ![Formula][61] For any *n**b* ≤ *n* ≤ *n**e*, ![Formula][62] and ![Formula][63] The ratio between the density functions is: ![Formula][64] Hence: ![Formula][65] Where ![Graphic][66]. We establish the following claim. Claim 4.8. *For any n**b* *< n* ≤ *n**e*, ![Formula][67] *Proof of Claim 4.8*. We will prove the claim using induction on *n*. * (Base case). Since *q**L* = 0 (*< q*) we have that: ![Graphic][68] for any *a*. Hence, *β**L*(*n**b*, *n**b* + 1) *> β*(*n**b*, *n**b* + 1). * (Inductive step). Let *n* ∈ (*n**b*, *n**e*). Assume that ![Formula][69] By Eq. (37), for any *s*1 *> s*2, ![Formula][70] and using Claim 4.6, ![Formula][71] Using the definition of *β* and by our assumption, ![Formula][72] we conclude that for any *n**b* *< n* ≤ *n**e*, *β**L*(*n**b*, *n*) *< β*(*n**b*, *n*).. □ Next, we apply Claim 4.8 for *n**e* and have that *β**L*(*n**b*, *n**e*) − *β*(*n**b*, *n**e*) *>* 0. Therefore, by Eq. (37) for any *s*1 *> s*2 ![Formula][73] and according to Claim 4.6, *R**L*(*n**e*) *< R*(*n**e*), where *R* and *R**L* are calculated using the original values of *p* and *q*. This is demonstrated in Figure 6. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/26/2021.01.23.21250242/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2021/01/26/2021.01.23.21250242/F6) Figure 6: The expected *R**L*(*n*) under an event-based spreading targeted lockdown (vs. *R*(*n*) in a natural evolution) assuming *p* = 0.5, *R* = 3 and *k* = 0.1. The lockdown begins at *n**b* = 5% and ends at *n**e* = 30%. In a similar develop as in Eq. (37), we have that for any *n > n**e*, ![Formula][74] where ![Graphic][75]. We establish the following claim, similar to Claim 4.8: Claim 4.9. *For any n**e* *< n*, ![Formula][76] *Proof of Claim 4.9*. The proof follows the same idea of the proof of Claim 4.8.The full proof is given in Appendix A. □ We conclude the proof of the claim using Claim 4.9 and Eq. (39). It holds that ∀*s*1 *> s*2, ∀*n* ≥ *n**e*, ![Formula][77] □ *Proof of Theorem 4.4*. Combining Claim 4.7 with Claim 4.6 we have that for any *n* ≥ *n**e*: ![Formula][78] Consequently, *HIT**L* *< HIT*, and Theorem 4.4 follows. □ Note that the proof of Theorem 4.4 was not based on the assumption that *q**L* = 0, but on the fact that *q**L* *< q* (see the proof of Claim 4.8). Hence, Theorem 4.4 can be generalized to hold for any 0 ≤ *q**L* *< q*, namely for *Partial Lockdown*. Corollary 4.10 (Generalization of Theorem 4.4). *Any event-based spreading targeted partial lockdown, namely where p remains at the same level and q**L* *< q, results in HIT**L* *< HIT*. Remark 4.11 (Generalization of Theorem 4.3). *We conjecture that the generalization of Theorem 4.3 for personal-trait spreading targeted partial lockdown holds as well*. ## 5 Conclusions and Discussion In this work we studied the effects of infectiousness heterogeneity (overdispersion) and lock-downs on herd immunity. Recent literature suggests that COVID-19 is characterized by such heterogeneity which affects dramatically herd immunity. We proposed that infectiousness (and susceptibility) should be classified into two inherently different types which we called *personal-trait* and *event-based* spreading. We followed up [8] and proposed a model that accounts for both of these types. Under this new model we showed that herd immunity and the HIT strongly depend on the mix between the two types of spreading. We analyzed the decay of the effective reproduction number, *R*(), and showed that the contribution of the personal-trait spreading drops sharply at early stages of the disease while the contribution of the event-based spreading drops much more slowly. That is -the *super-spreaders* “leave the game” at early stages, while the *super-spreading* events remain active. We demonstrated the results on a (Gamma) distribution which previously attributed to the infectiousness of COVID-2 [5] and COVID-19 [3, 8]. We showed that in order to predict the mix between the spreading types, only a small fraction of the population is needed to get infected. Such prediction can help in calculating the HIT in an early stage of the disease. This estimation problem for COVID-19 is left open for forthcoming research. We addressed operational aspects of disease blocking and analyzed the effect of lockdowns on the HIT. We showed that different lockdown strategies, targeting different spreading types, result in opposite effect on the HIT. In particular, a lock-down which targets personal-trait spreading would act adversely and *reduces* herd immunity. This seems to fit a lockdown that focuses on daily/professional activities. In contrast, a lock-down which targets event-based spreading will *increase* herd immunity. This may fit a lockdown that focuses on sports/social events. Of course, a lockdown may have other objectives such as achieving temporary slow down of the disease spreading (to allow handling the patients masses) which may justify the lockdown strategy. Yet – the effect on herd immunity requires consideration, especially due to the opposite effects of various lockdowns. ## Data Availability Online data sources and cited papers/books can be found at the list of references. ## A Proofs *Proof of Claim 3.2 (continued)*. The proof follows the proof of Claim I provided in [8]. Taking natural log of Eq. (13), ![Formula][79] It holds that: ![Formula][80] We attempt to find the herd-immunity threshold, and hence can bound number of steps, *n*, by (1 − 1*/R*) · *N*. Hence, ![Formula][81] Counting over the steps 1, …, *n* we have: ![Formula][82] And as long as ![Graphic][83] we have that: ![Formula][84] Note that for any *i*, ![Formula][85] Since our calculation is done by taking an expectation over all possible scenarios of infections, we have: ![Formula][86] Hence, ![Formula][87] And the proof is complete. □ *Proof of Claim 4.6*. Let *s*′ *>* 0. By our assumption, for any *s < s*′: ![Formula][88] Hence, integrating *s* over (0, *s*′) we have: ![Formula][89] Therefore: ![Formula][90] By the same ideas we have that: ![Formula][91] Using Eq. (45) and Eq. (46) we have: ![Formula][92] and hence: ![Formula][93] Since *r* is an injective monotone function, it holds that: ![Formula][94] In the same way: ![Formula][95] Using Eq. (47) combined with Eq. (48) and 49): ![Formula][96] For any non-negative random variable *X* it holds that: ![Formula][97] Therefore, by Eq. (50) combined with 51 we know that: ![Formula][98] and the proof is complete. □ *Proof of Claim 4.9*. We will prove the claim using induction on *n*. * (Base case). Using Eq. (38), we have that: ![Formula][99] Hence, ![Formula][100] * (Inductive step). Let *n**e* *< n*. Assume that ![Formula][101] By Claim 4.8 and Eq. (39), for any *s*1 *> s*2, ![Formula][102] and using Claim 4.6, ![Formula][103] Using the definition of *β* and by our assumption, ![Formula][104] we conclude that for any *n**e* *< n, β**L*(*n**e*, *n*) *< β*(*n**e*, *n*). □ ## Acknowledgment This research was supported in part by Len Blavatnik and the Blavatnik Family foundation. ## Footnotes * * E-mail: jhonatan.tavori{at}cs.tau.ac.il * † E-mail: hanoch{at}tauex.tau.ac.il * 1 For continuous *s*, Eq. (2) should be considered as a density function. * 2 The expectation is taken over all possible scenarios of infection. * Received January 23, 2021. * Revision received January 23, 2021. * Accepted January 26, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. [1].The second (1926) milroy lecture on experimental epidemiology. The Lancet, 207(5350): 531 – 537, 1926. ISSN 0140-6736. doi:>[https://doi.org/10.1016/S0140-6736(00)92941-6](https://doi.org/10.1016/S0140-6736(00)92941-6).Originally published as Volume 1, Issue 5350. 2. [2]. R. M. Anderson and R. M. May. Infectious diseases of humans: dynamics and control. Oxford university press, 1992. 3. [3]. A. Endo, S. Abbott, A. J. Kucharski, S. Funk, et al. Estimating the overdispersion in covid-19 transmission using outbreak sizes outside china. Wellcome Open Research, 5(67): 67, 2020. 4. [4]. P. Fine, K. Eames, and D. L. Heymann. “herd immunity”: a rough guide. Clinical infectious diseases, 52(7):911–916, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cir007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21427399&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F26%2F2021.01.23.21250242.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000288802600016&link_type=ISI) 5. [5]. J. O. Lloyd-Smith, S. J. Schreiber, P. E. Kopp, and W. M. Getz. Superspreading and the effect of individual variation on disease emergence. Nature, 438(7066):355–359, 2005. doi:10.1038/nature04153. URL [https://doi.org/10.1038/nature04153](https://doi.org/10.1038/nature04153). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature04153&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16292310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F26%2F2021.01.23.21250242.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233300200048&link_type=ISI) 6. [6]. J. Lourenco, F. Pinotti, C. Thompson, and S. Gupta. The impact of host resistance on cumulative mortality and the threshold of herd immunity for sars-cov-2. MedRxiv, 2020. 7. [7]. D. Miller, M. A. Martin, N. Harel, O. Tirosh, T. Kustin, M. Meir, N. Sorek, S. Gefen-Halevi, S. Amit, O. Vorontsov, et al. Full genome viral sequences inform patterns of sars-cov-2 spread into and within israel. Nature communications, 11(1):1–10, 2020. 8. [8]. Y. Oz, I. Rubinstein, and M. Safra. Heterogeneity and superspreading effect on herd immunity. arXiv preprint arXiv:2009.01923v2, 2020. 9. [9]. R. Pastor-Satorras, C. Castellano, P. Van Mieghem, and A. Vespignani. Epidemic processes in complex networks. Reviews of modern physics, 87(3):925, 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/RevModPhys.87.925&link_type=DOI) 10. [10]. H. E. Randolph and L. B. Barreiro. Herd immunity: Understanding covid-19. Immunity, 52(5):737–741, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.immuni.2020.04.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32433946&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F26%2F2021.01.23.21250242.atom) 11. [11]. K. Rock, S. Brand, J. Moir, and M. J. Keeling. Dynamics of infectious diseases. Reports on Progress in Physics, 77(2):026602, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/0034-4885/77/2/026602&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24444713&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F26%2F2021.01.23.21250242.atom) 12. [12].Rt.live. Current r value in the us, 2020. URL [https://rt.live](https://rt.live). 13. [13]. A. V. Tkachenko, S. Maslov, A. Elbanna, G. N. Wong, Z. J. Weiner, and N. Goldenfeld. Persistent heterogeneity not short-term overdispersion determines herd immunity to covid-19. arXiv preprint arXiv:2008.08142, 2020. 14. [14]. W. Topley and G. Wilson. The spread of bacterial infection. the problem of herd-immunity. Epidemiology & Infection, 21(3):243–249, 1923. 15. [15].Worldometers.info. Covid-19 coronavirus pandemic, 2020. URL [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/). [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/inline-graphic-1.gif [6]: /embed/inline-graphic-2.gif [7]: /embed/graphic-5.gif [8]: /embed/graphic-6.gif [9]: /embed/graphic-7.gif [10]: /embed/graphic-8.gif [11]: /embed/graphic-9.gif [12]: /embed/graphic-10.gif [13]: /embed/graphic-11.gif [14]: /embed/graphic-12.gif [15]: /embed/graphic-13.gif [16]: /embed/graphic-14.gif [17]: /embed/graphic-15.gif [18]: /embed/graphic-16.gif [19]: /embed/graphic-17.gif [20]: /embed/graphic-18.gif [21]: /embed/graphic-19.gif [22]: /embed/graphic-20.gif [23]: /embed/graphic-21.gif [24]: /embed/graphic-22.gif [25]: /embed/graphic-23.gif [26]: /embed/graphic-24.gif [27]: /embed/graphic-25.gif [28]: /embed/graphic-26.gif [29]: /embed/graphic-27.gif [30]: /embed/graphic-28.gif [31]: /embed/graphic-29.gif [32]: /embed/graphic-30.gif [33]: /embed/graphic-33.gif [34]: /embed/graphic-35.gif [35]: /embed/graphic-36.gif [36]: /embed/graphic-37.gif [37]: /embed/graphic-38.gif [38]: /embed/graphic-39.gif [39]: /embed/graphic-40.gif [40]: /embed/graphic-41.gif [41]: /embed/graphic-42.gif [42]: /embed/graphic-43.gif [43]: /embed/graphic-44.gif [44]: /embed/graphic-45.gif [45]: /embed/graphic-46.gif [46]: /embed/graphic-47.gif [47]: /embed/graphic-49.gif [48]: /embed/graphic-50.gif [49]: /embed/graphic-51.gif [50]: /embed/graphic-53.gif [51]: /embed/graphic-54.gif [52]: /embed/graphic-55.gif [53]: /embed/graphic-56.gif [54]: /embed/graphic-57.gif [55]: /embed/graphic-58.gif [56]: /embed/graphic-59.gif [57]: /embed/inline-graphic-3.gif [58]: /embed/graphic-60.gif [59]: /embed/graphic-61.gif [60]: /embed/graphic-62.gif [61]: /embed/graphic-63.gif [62]: /embed/graphic-64.gif [63]: /embed/graphic-65.gif [64]: /embed/graphic-66.gif [65]: /embed/graphic-67.gif [66]: /embed/inline-graphic-4.gif [67]: /embed/graphic-68.gif [68]: /embed/inline-graphic-5.gif [69]: /embed/graphic-69.gif [70]: /embed/graphic-70.gif [71]: /embed/graphic-71.gif [72]: /embed/graphic-72.gif [73]: /embed/graphic-73.gif [74]: /embed/graphic-75.gif [75]: /embed/inline-graphic-6.gif [76]: /embed/graphic-76.gif [77]: /embed/graphic-77.gif [78]: /embed/graphic-78.gif [79]: /embed/graphic-79.gif [80]: /embed/graphic-80.gif [81]: /embed/graphic-81.gif [82]: /embed/graphic-82.gif [83]: /embed/inline-graphic-7.gif [84]: /embed/graphic-83.gif [85]: /embed/graphic-84.gif [86]: /embed/graphic-85.gif [87]: /embed/graphic-86.gif [88]: /embed/graphic-87.gif [89]: /embed/graphic-88.gif [90]: /embed/graphic-89.gif [91]: /embed/graphic-90.gif [92]: /embed/graphic-91.gif [93]: /embed/graphic-92.gif [94]: /embed/graphic-93.gif [95]: /embed/graphic-94.gif [96]: /embed/graphic-95.gif [97]: /embed/graphic-96.gif [98]: /embed/graphic-97.gif [99]: /embed/graphic-98.gif [100]: /embed/graphic-99.gif [101]: /embed/graphic-100.gif [102]: /embed/graphic-101.gif [103]: /embed/graphic-102.gif [104]: /embed/graphic-103.gif