Abstract
Background Societal separation of unvaccinated people from public spaces has been a novel and controversial COVID-era public health practice in many countries. Models exploring potential consequences of vaccination-status-based separation have not considered how separation influences the contact frequencies in the separated groups; we systematically investigate implementing effects of separation on person-to-person contact frequencies and show this critically determines the predicted epidemiological outcomes, focusing on the attack rates in the vaccinated and unvaccinated populations and the share of infections among vaccinated people that were due to contacts with infectious unvaccinated people.
Methods We describe a susceptible-infectious-recovered (SIR) two-population model for vaccinated and unvaccinated groups of individuals that transmit an infectious disease by person-to-person contact. The degree of separation between the two groups, ranging from zero to complete separation, is implemented using the like-to-like mixing approach developed for sexually-transmitted diseases [1-3], adapted for presumed SARS-CoV-2 transmission. We allow the contact frequencies for individuals in the two groups to be different and depend, with variable strength, on the degree of separation.
Results Separation can either increase or decrease the attack rate among the vaccinated, depending on the type of separation (isolating or compounding), and the contagiousness of the disease. For diseases with low contagiousness, separation can cause an attack rate in the vaccinated, which does not occur without separation.
Interpretation There is no blanket epidemiological advantage to separation, either for the vaccinated or the unvaccinated. Negative epidemiological consequences can occur for both groups.
Introduction
Models can be used to investigate infectious disease dynamics under different hypotheses about the characteristics of a disease and the effects of health policy. In this endeavour, there are advantages to working with simplest-possible but sufficiently realistic models [3,4], where one should exclude simple models that are not sufficiently realistic for the intended application, either because of their structure or because of incorrect assumptions about the underlying mechanisms. Following this approach, researchers have extended the foundational simple susceptible-infectious-recovered (SIR)-type model to explore diseases with birth and death dynamics, maternal- or vaccine-derived immunity, latency of infection, and so on [5-7]. Recently, SIR models of epidemic dynamics have been implemented with two interacting societal groups (vaccinated and unvaccinated) to examine epidemic outcomes for variable degrees of interaction between the two groups, including whether the unvaccinated put the vaccinated unduly or disproportionately at risk, using epidemiological parameters intended to be representative of SARS-CoV-2 [8-10]. These prior implementations take the person-to-person contact frequencies of the majority and socially-excluded groups to be equal and held constant, irrespective of the degree of separation (or exclusion or “like-to-like mixing”), which is not realistic.
Here, we implement person-to-person contact frequencies that can be different for the two groups and that can either increase or decrease with increasing separation. This is necessary because, for example, in many actual regulatory policies the excluded unvaccinated group is barred from public venues or services where people gather and from public transport where people are in close proximity for various durations. In general, the person-to-person contact frequency of the excluded group decreases with increasing separation if isolation is in effect, and increases with increasing separation if the excluded individuals are crowded together. Implementing this essential model feature gives rise to more complex behaviour of the attack rates in the vaccinated and unvaccinated populations (Avand Au, respectively), which can increase or decrease, or rise to a maximum before decreasing, as the two groups are increasingly separated. This is also true for the share of infections among vaccinated people that are due to contacts with infectious unvaccinated people (Bv).
Methods
Model Design
We adopt the standard SIR framework in a structure with two sub-populations. Following the usual SIR model design, a person can be in one of three states: susceptible to infection (S), infectious (I), or recovered and immune (R). If a susceptible person comes into contact with an infectious person, the susceptible person can become infectious, and infectious people eventually recover and become permanently immune.
We divide the population into two groups: vaccinated and unvaccinated. Vaccination is “all or nothing”, such that a proportion VE of the vaccinated population is immune (are in the R state from the outset of the simulation), where the parameter VE represents vaccine efficacy. The model also includes a natural immunity parameter, NI, equal to the proportion of unvaccinated that are immune from the outset due to previous infection [8].
The model parameter η controls the degree of separation between vaccinated and unvaccinated people. When η = 0, there is no separation, and the two groups mix randomly. When η = 1, there is complete separation, such that vaccinated only come into contact with other vaccinated, and unvaccinated only come into contact with other unvaccinated.
The parameter η follows from Garnett and Anderson [3], who modeled sexually-transmitted disease spread in a population divided into groups with different frequencies of sexual contacts. They take the contact frequency to be a constant characteristic of the individuals within a group. However, contact frequency is not generally and solely an intrinsic individual characteristic [11], and separation based on vaccination status may increase or decrease contact frequencies, depending on how the separation is implemented.
In our model, the contact frequencies of the vaccinated and unvaccinated individuals (cvand cu, respectively) can increase, decrease, or remain constant as the two groups are separated. We implement a new approach to achieve this: we keep the first two terms in Taylor expansions of cvand cuversus η and ; see Eqs. A3 in Appendix A1). Thus, mv and mu determine the degree of increase or decrease of the contact frequency in either group, as η is increased.
For example, when mu < 0, as separation is increased, the contact frequency of unvaccinated people decreases. This corresponds to separation policy that excludes unvaccinated people from public spaces such as restaurants, cinemas, workplaces, airplanes, trains, etc. [12-15]. Conversely, mu > 0 corresponds to separation policy that increases contacts between unvaccinated people; for example, by requiring returning unvaccinated travelers to stay in designated facilities [16-18].
The model of Fisman et al. [8] is the special case of our model with mu= mv= 0 and , in which case the equal contact frequencies of both vaccinated and unvaccinated remain constant regardless of the level of separation. Such separation implies two parallel societies, where both populations have their own public spaces and within these spaces maintain the same contact frequencies they would have with no separation. Such an implementation does not represent how separation has been applied during the COVID era in Canada and many countries [12-15,19], since unvaccinated people were excluded from public spaces while vaccinated people were allowed access.
Model Parameterization
The parameters of our model are listed in Table 1; calculated quantities in Table 2. Technical details of the model are in Appendix 1.
Analysis
The attack rate among the vaccinated population is defined as the proportion of initially-susceptible vaccinated people who become infected during the epidemic: where Sv(t0) is the number of susceptible vaccinated people at the beginning of the epidemic and Sv(tf) is the number of susceptible vaccinated people remaining once there are no longer any infectious people in the entire (vaccinated and unvaccinated) population. Au is defined equivalently, for the unvaccinated, replacing the v subscripts with u in Eq. 1. The overall attack rate for the full (vaccinated plus unvaccinated population) is: We also define Bv as the share of infections among vaccinated people that were due to contacts with infectious unvaccinated people (see Eq. A4 of Appendix 1).
We focus on separation types that are targeted at the unvaccinated group. We assume, for simplicity, that separation has no impact on the contact frequency of vaccinated people (mv = 0). We also assume that the contact frequencies in both groups are the same when there is no separation . We use the same values as used by Fisman et al. [8] for the remaining parameters: Pv = 0.8, VE=0.8, NI=0.2, γ = 73 yr-1, and N = 107. These values were intended to be representative for COVID-19 and vaccination; in particular, the recovery rate of 73 yr-1 is equivalent to a recovery time of 5 days [20,21], and assumed to be the same for vaccinated and unvaccinated people.
Appendix 2 contains supplementary figures with results for different parameter combinations, including Pv ≠ 0.8, mv ≠ 0 and . In all results in this paper, simulations were initiated with a seed number of 100 infectious individuals distributed proportionately among the two sub-populations.
Results
Fig. 1 shows simulation results for a range of model parameters for different epidemiological conditions and degrees and types of societal separation. Each row of panels is for a fixed value of . Moving from the top row (Figs. 1a.i to a.iv) to the bottom (Figs. 1e.i to e.iv), . The left column of panels shows how the attack rate among the vaccinated population, Av, changes with the degree of separation, η. The second and third columns show Au and At as functions of η, and the right column shows how Bv, the share of vaccinated infections that were due to contacts with unvaccinated people, varies with η.
Figs. 1c.i to c.iv show results for a moderate value of . For reference, in a single-population (no vaccination) model, c = 300 and γ = 73 corresponds to a basic reproduction number R0= c/γ = 4.1.
In Fig. 1c.i, when mu=-1 and mu= -0.5 (reflecting large and moderate degrees of exclusion and isolation of unvaccinated people) the vaccinated attack rate, Av decreases with increasing separation. However, when mu > 0 (compounding of unvaccinated people) or mu= 0 (separation has no influence on contact frequency of unvaccinated people), there is a maximum in Av for moderate values of η. Therefore, with compounding separation, large values of η are required for Av to be lower than its value for no separation (η = 0). Fig. 1c.ii shows that the unvaccinated attack rate, Au, increases with separation for anything other than strong isolating separation (mu approaching -1). This produces a maximum in the overall attack rate At at moderate degrees of separation, even for values of mu for which Av decreases monotonically (mu = -0.5).
Fig. 1c.iv shows that Bv, the share of vaccinated infections that are due to unvaccinated people, has a shape similar to Av(η, mu). In all panels, 20% of the total population is unvaccinated (Pv= 0.8; Table 1).
Figs. 1c.i to c.iv therefore demonstrate that whether separation increases or decreases the vaccinated-population attack rate depends on both the degree of separation and how separation affects contact frequency.
Figs. 1a.i to a.iv and 1b.i to b.iv show results for larger . Compared to Fig. 1c.i, Av in Figs. 1a.i and 1b.i does not increase much with η when mu > 0, and Av no longer has a maximum when mu = 0. It can also be seen that Av increases with increasing when there is no separation (η = 0).
Reducing (Figs. 1d.i to d.iv and 1e.i to e.iv), decreases Av(η = 0), and larger η can dramatically increase Av. Even with an isolating separation policy (mu = −0.5 in Fig. 1d.i), Av is increased for moderate values of η.
When are small enough ( in Fig. 1e.i, corresponding to R0 = 2.7 in a single population (no vaccination) model), there is no epidemic among the vaccinated in the absence of separation (Av(η = 0) = 0). However, a non-zero vaccinated-population attack rate (Av> 0) occurs if η is sufficiently large, and emerges regardless of whether one isolates or compounds the unvaccinated. Therefore, for small enough values of , any separation could increase infections among the vaccinated.
The main qualitative features of the above results for Pv= 0.8 hold for other values of Pv. Appendix 2 provides a detailed exploration of results for Pv = 0.1 through 0.99; and for two values of VE (0.4 and 0.8). When VE is decreased, Av is not strongly influenced by η, regardless of mu; therefore, any beneficial effect of separation on Av is reduced as VE decreases.
Appendix 2 also explores . For example, when , the unvaccinated contact frequency is reduced even when there is no separation; increasing η can then increase Av substantially compared to the case of , holding all other parameter values constant (see panels a.i and b.i in Figs. A2.28 and A2.31).
Interpretation
Separation can have substantially different and negative impacts on the outcome of an epidemic, depending on the type and degree of separation, and depending on cultural and population-density factors, for example, that co-determine and .
Separation that compounds the unvaccinated (mu > 0 and mv = 0) generally causes an increase in the vaccinated-population attack rate, Av, for small and intermediate degrees of separation, η, while for large η, Av decreases below its value in an unseparated society. Separation that isolates and excludes the unvaccinated (mu < 0 and mv = 0) decreases Av for “more contagious viruses” (i.e. large , large R0); however, for “less contagious viruses” (smaller , smaller R0), both isolating and compounding types of separation can increase Av beyond its value in an unseparated society. For “viruses that are not very contagious” (small , small R0), applying separation can cause a sizeable epidemic among the vaccinated even though virtually no vaccinated people would be infected in an unseparated society. Separation increases the unvaccinated attack rate, Au, for compounding and moderately isolating types of separation, and Au is only decreased for strongly isolating separation (mu approaching -1).
Except for large negative values of mu, applying separation has the effect of increasing the overall probability of a susceptible-infectious interaction, since the unvaccinated population has a higher fraction of susceptibles. This creates a form of core group dynamics [22-24]. At the same time, increasing separation shields the vaccinated population from the increased prevalence of infection in the unvaccinated population. This trade-off causes the non-monotonic relationship between Av and η. The same dynamic causes the emergence of an epidemic for large η when (and thus R0) is small.
We find that Bv, the share of vaccinated infections that are due to contact with unvaccinated people, follows a similar trend to Av as a function of the degree of separation, when separation has no impact on the vaccinated contact frequency (mv = 0). For this type of separation, Av and Bv either increase or decrease simultaneously with increasing η, depending on the value of mu, and Bv is minimized for complete separation.
In contrast, when mv ≠ 0, such that separation affects the contact frequencies of vaccinated people, increasing separation can cause Av to increase while Bv decreases and vice-versa (see Appendix 2, Figs. A2.25 and A2.26).
The impact of separation on contact frequencies has not previously been considered to our knowledge, including in network-based models in which unvaccinated people cluster together in “cliques” or households [25-27].
Limitations
Our model assumes only two risk populations (vaccinated and unvaccinated), considers only the attack rates on epidemic completion (Av and Au), and takes the degree of separation η to be time-independent, without variation due to public holidays and such. It does not consider other outcomes such as death or hospitalization. Our model assumes an all-or-nothing VE, without waning immunity or influence on infectiousness; and no possibility of reinfection. SIR models and their variations are based on the paradigm of transmission due to pairwise contact between a recently infected and a susceptible individual. However, this paradigm is unable to account for important features of viral respiratory disease incidence, in particular its rapid emergence and disappearance occurring at essentially the same time at widely dispersed locations [28]. Air-borne transmission via suspended aerosol particles is not directly compatible with pairwise transmission, since it occurs in built environments where many people may transit or be present [29]. A related and unavoidable limitation is the lack of reliable empirical evaluations of needed infectious contact frequencies, which is important because our calculated outcomes are sensitive to the chosen contact frequency values.
Conclusion
Vaccination-status-based societal separation can lead to substantially different and counter-intuitive epidemic outcomes depending on the type and degree of separation. Negative epidemiological consequences can occur for either separated group, irrespective of the expected deleterious health impacts of the separation policies themselves [30-32].
Data Availability
This is a modeling paper, therefore all data is from simulations. Information for reproducing the simulations is contained in the manuscript and appendices.
Appendix 1 Elaboration of the model
A1.1: Model differential equations and “mixing” rule
The model is a susceptible-infectious-recovered (SIR) model with two populations: vaccinated (subscript v) and unvaccinated (subscript u) people, consisting of the following six differential equations: Su, Iu, and Ru represent the number of susceptible, infectious, and recovered unvaccinated people, at time t. Nu represents the total number of unvaccinated people. cu represents the contact frequency (number of contacts per unit time) of unvaccinated people. βu is the probability that a susceptible unvaccinated person becomes infected upon contact with an infectious person (regardless of whether the infectious person is vaccinated or unvaccinated). γu is the rate at which infected unvaccinated people recover from infection. The quantities Sv, Iv, Rv, Nv, cv, βv, and γv are defined equivalently, for vaccinated people.
There are thus two “β parameters”, two “c parameters” and two “γ parameters” in our model. Since each β parameter always occurs as part of a product with its respective c parameter, the β parameters can freely be set equal to 1: this imposes that the “contacts” considered in the model are, by definition, only those contacts that are of sufficiently close proximity and long duration that an infection is guaranteed to occur when a susceptible and an infectious person meet. We set βv= βu = 1 in the main text, without any loss of generality.
fij is the probability that a person of type i (either u or v) has a contact with a person of type j (either u or v), and is defined as follows:, where 0 ≤ η ≤ 1 and δij is the Kronecker delta, such that δuu = δvv = 1 and δuv = δvu = 0. η is therefore a parameter that controls the degree of separation between the u and v sub-populations. For example, when η = 1, then fuu = fvv = 1 and fuv = fvu = 0, such that u people only ever have contacts with other u people and likewise for v people (complete separation). At the other extreme, when η = 0, then the probability that a u person has a contact with a v person is entirely determined by the relative proportions of u and v people, weighted by their respective contact frequencies (no separation, or “random mixing”).
A1.2: Variation of contact frequency with degree of separation
To allow contact frequency to vary with degree of separation, we define: where and represent the contact frequencies for unvaccinated and vaccinated people when there is no separation (η = 0), and mu ≥ −1 and mv ≥ −1 are two parameters that control how non-zero separation impacts the contact frequencies of unvaccinated and vaccinated people.
Fig. A1.1 shows cu as a function of η, for different values of the parameter mu and for . As can be seen, when mu = 0, , such that the contact frequency is constant regardless of the degree of separation.
When mu < 0, the contact frequency of unvaccinated people decreases with increasing η. This represents a separation policy that excludes unvaccinated people from public spaces while also isolating them from themselves to some degree. For example, a pair of intermediate values of mu < 0 and 0 < η < 1 could represent a separation policy (such as with “vaccination passports”) that excludes unvaccinated people from recreational venues like restaurants and cinemas but not more essential services such as grocery stores and hospitals. On the other hand, in the extreme case of mu = −1 and η = 1, cu = 0 such that unvaccinated people are completely separated and isolated, having no contacts with anyone.
When mu > 0, the contact frequency of unvaccinated people increases with increasing η. This represents a separation policy that compounds unvaccinated people by placing them together in close quarters, for example in designated facilities for returning unvaccinated travelers.
Via the two parameters mu and η, the model therefore spans the full range of contact frequencies, from cu = 0 in the isolating extreme of mu = −1 and η = 1 to arbitrarily high contact frequency in the compounding extreme with η = 1 and a mu > 0.
The impact of separation on the vaccinated sub-population is independently controlled via the parameters and mv.
A1.3: Quantities calculated from simulation results
The attack rate among the vaccinated population, Av, is defined as the proportion of initially-susceptible vaccinated people who become infected during the epidemic (Eq. 1 of the main text, reproduced below): where Sv(t0) is the number of susceptible vaccinated people at the beginning of the epidemic and Sv(tf) is the number of susceptible vaccinated people remaining once there are no longer any infectious people in the entire (vaccinated and unvaccinated) population. Au is defined equivalently, for the unvaccinated.
The quantity Bv is equal to the share of infections among vaccinated people that were due to contacts with infectious unvaccinated people, i.e.: Similarly, Bu represents the share of infections among unvaccinated people that were due to contacts with infectious vaccinated people, and is defined in the same way as Bv, (interchanging the v and u subscripts in Eq. A5).
Appendix 2 Supplementary figures
A2.2: Epidemic outcomes for different values of Pv
A2.2.1: VE = 0.8
A2.3: Epidemic outcomes for different values of mv
A2.4: Epidemic outcomes for different values of
A2.4.1: Fixed , vary
A2.4.2: Fixed weighted sum
Footnotes
Revisions following reviewer comments, to improve clarity and precision, add new references and discussion, and add new figures to Appendix 2.