Abstract
Background The severity of the Covid-19 pandemic has led to the use of extreme control measures which has halted the spread, but at enormous socioeconomic cost. We set out to explore whether a model that includes ‘superspreader’ could explain aspects of mitigation strategies that explain the dramatic effect.
Methods We developed an age structured agent-based model that includes persons who spread the disease far more widely than others, acting in a society in which transmission occurs in three sectors: home, school/work and a third category representing other social contacts. We employ the model to study the impact of lock-down when superspreaders are present.
Findings As expected, without mitigation imposed, the inclusion of superspreaders in the population does not change the epidemic trajectory. However in a structured society we find that superspreaders made a substantial difference. Our simulations demonstrate that workplaces and schools may open without much effect on the epidemic, as long as the other social contacts are drastically limited.
Interpretation The recent observations of profound pandemic control that cannot be captured by standard SEIR models, require disease transmission models that consider superspreaders. We have found that transmission can be controlled simply by limiting contacts such as public transportation and large events. Indeed, eliminating superspreader opportunities can uniquely explain the observed control of the epidemic with Sweden’s relaxed approach and with the moderate lock-down used in Denmark.
Funding KS received funding from the EU Horizon 2020 research program under ERC grant No. 740704, and LS from the Carlsberg Foundation in Denmark.
Research in context:
Evidence before this study:
Since the emergence of SARS-CoV-2 in China, superspreaders are emerging as a theme, as it was also the case with SARS-CoV in 2003. Analyzing COVID-19 data, substantial individual-level variation in transmissibility has been inferred, and it has been estimated that 10% of infected persons cause 80% of all infections.
Added value of this study:
We developed an agent-based socially structured model to simulate the effect of superspreaders on the epidemic, and in the context of country lockdowns and reopenings. Our simulations demonstrate that COVID-19 is effectively mitigated if one limits social activity where many strangers meet, as superspreaders will infect mainly among the many people they rarely meet.
Implications of all available evidence:
Modelling heterogeneity and super spreading is critical in terms of understanding observed COVID-19 epidemic patterns. Indeed, mitigating superspreading opportunities is an effective way to mitigate Covid-19 and prevent second waves because the remaining population cannot sustain the epidemic. The superspreader phenomenon may also explain the observation of huge COVID-19 epidemic variability between different areas in the same country.
Introduction
The emerging SARS-CoV-2 virus that is causing the COVID-19 pandemic which emerged in late 2019 has lead countries around the world to an unprecedented lock-down strategy. We have seen that these lock-downs work extremely well in terms of halting the spread, at enormous socioeconomic cost. However, we do not know what aspects of the mitigation efforts are causing the effect. As many countries are already in the process of opening up, the relative contribution of different aspects – reducing contacts at home, schools, workplaces and other sectors of society – should be better understood.
For this pandemic, a pattern of “superspreader” events has been documented to occur [1]. For example, as South Korea reopened the country carefully recently, a single infected person visiting a night club led to at least 50 new infections. In the US, a 2,5 hour long choir rehearsal tragically led to 75% of the participants getting infected and several deaths. These, as well as outbreaks in prisons and hospitals and following carnival parties are vivid examples of SARS-CoV-2 superspreader events that illustrate the importance of these events in the spread of this novel virus. It is also reminiscent of documented superspreader events in the dissemination of SARS-CoV in the 2003 outbreak in Asia and Canada. One recent study estimated that 80% of infections may be caused by 10% of the infected population, pointing to transmissibility (R0) being highly heterogeneous [2, 3]. Thus, the practice of relying on an average R0 in dynamic disease models can obscure considerable individual variation in infectiousness when extreme transmitters are included [4, 5, 6, 7]. It is also understood that such heterogeneity will not affect the epidemic spread by much in standard dynamic SEIR models. One might ask why are we often seeing superspreading events as a precursor to outbreaks? Also, why do we see marked heterogeneity in the geographical distribution and intensity of outbreaks across countries and regions?
Such heterogeneity in transmission risk is a well understood phenomenon in infectious diseases. In 2003 Riley et al. pointed out that ‘superspreading events’ could explain 80% of transmission events for a wide variety of diseases, thus being a defining characteristic of many epidemics, including vector-borne diseases due to large variability in parasite density in blood in infected individuals [8]. In the 2014 Ebola outbreak in west Africa, superspreader events are thought to have played a key role in sustaining onward transmission of the epidemy [9]. Overall, this means that certain individuals infect far greater numbers of persons than the average.
For SARS-CoV-2 the observed patient-to-patient variation in virus load appears to be several orders of magnitude [10]. Similarly, observed infectivity varies greatly [2, 3]. Therefore we choose to focus on the variation in personal specific infectivity, while assuming that there is no variation in individual exposure.
We ask here whether the phenomenon of superspreading might play an important role in our ability to control the ongoing COVID-19/SARS-CoV-2 pandemic with non-pharmaceutical mitigation strategies. To investigate the effect of heterogeneity in various spaces we needed to modify our disease transmission model. One can easily do this in a agent-based model in which each person can be assigned individual infection properties. Our findings suggest a rating of the importance of limiting spread in the home, the school-or work place versus in other parts of the social environment where many new contacts are encountered and superspreaders events may occur.
Methods
We developed an age stratified version of the agent-based model previously described in [15], see Fig. 1. The model is simulated in 6 hour time-steps, during which each infected agent has a probability of infecting another agent. In the structured model infections attempts are weighted equally between the three social sectors outlined in Fig. 1. When simulating mitigation scenarios, the infections in one or more sectors are reduced. Infections from superspreaders are simulated by repeating their infectious activity 50 times within the 6 hour time window. When there is 10% of these superspreaders they cause between 75 and 85% of all infections.
The top panel shows various stages of the disease, each simulated with a corresponding rate. People can infect from blue and red stages, and their relative duration corresponds to 45% pre-symptomatic infection [11]. The bottom panel illustrates the schematic society considered in our model, where interactions of each individual are related to three sectors: contacts in“Home”, “Work/school” and “Other” is adjusted to occur equally often (idealized simplification from ref. [12, 13]). Home is the smallest unit (average cluster of size 2.1), while Work/school is larger (each person has an average of 10 connections at work and 26 at school). All ages between 20 and 70 are mixed in the work places. For each agent the Home and Work/school interactions are fixed throughout the epidemic. The “Other” group is treated as consisting of random contacts. are further assigned a social activity as function of age, as shown in the insert, adjusted to fit the data from ref. [14].
In Figs. 2-4 we simulate 200.000 agents, while Fig. 5 considers one million. Simulations are started with 50 infected people. Throughout this work we calibrate the infection rate to fit an exponential growth of the unconstrained epidemics of 23%/day, as reported by [16]. A more complete description, that also includes the age dependent hospitalization and intensive care unit (ICU) probabilities, can be found in the supplement.
We finally explore the real-life situations in Denmark and Sweden (Fig. 5), now fitting the model to daily patients using data from the Danish [17] and Swedish [18] health authorities. The Swedish data were averaged over a 5 days time window.
Panel a: Classic (SEIR) type model where all agents infect at same rate (R0 = 2.7). The 100% curve marks the unmitigated epidemic, the other two curves mark two levels of lockdown in which 50% and 25% of contacts, respectively, remain. Panel b: Same random contacts as in a), but now assuming that 10% of infected agents are superspreaders, that transmit 50 times more effectively than others. Here the lock-down only imposes a limit on the maximum number of contacts, rather than a restriction in the number of contacts for all agents.
a) Unmitigated epidemic scenario (gray curve) using Norwegian age-specific data for ICU utilization (red curve). b) Model simulation without Home contacts, and c) without Work/school contacts, and d) without “Other” contacts, that is, contacts outside Home and Work/school. A three quarters reduction in peak daily ICU usage is achieved when limiting the “Other” contacts, that is, contacts outside home and work/school. e) Final population attack rate as a function of age for each of the four simulations in panel a,b,c,d. We find that mitigating “Other” contacts results in dramatic reductions in ICU. This is explained by the fact that the older generations having contact to the remaining society largely through contacts of the “Other” type.
Panels a-d) recapitulate the simulations done for the social-structured model of Fig.3. a) Unmitigated epidemic scenario (gray curve) using Norwegian age-specific data for ICU utilization (red curve). b) Model simulation without Home contacts, and c) without Work/school contacts, and d) without “Other” contacts. e) Sensitivity analysis where the light blue area marks the fraction of “Work/school” and “Other” sectors that need to be isolated to end the Covid-19 epidemics. Here the superspreaders still infect a factor of 50 times more than other agents but their number is varied between 0 and 20%. f) Sensitivity analysis of model with a distribution of individual infectivity drawn from p(s) ∝ 1/sγ, s < 1000 with γ varying between 1.2 and 2.5. In the sensitivity analysis we used a population size of one million.
Using 1000000 agents and seeding with 100 infected persons, we fitted our model to Danish daily time series of new hospital patients (blue circles). The initial lock down is assumed to reduce contacts in the Work/school sector to 25%, and then the lock down is partially lifted one month later, to 50%. a) Simulation without superspreaders where the “Work/school” and “Other” sectors are reduced as indicated by the width of the black and red regions, respectively, shown at the bottom of the graph. b) Trajectory including superspreaders in the model, where “Other” is a factor 1/6 of pre lockdown value. c) Alternative outlook of simulation from b) if social distancing in the “Other” sector is only partly maintained after the lock down is fully lifted. d) Trajectory of the model with superspreaders compared to Swedish data. The shown fit uses an “Other” sector that is set to 1/6 of the pre-mitigation value. In each scenario we show 3 different stochastic trajectories with same parameters and initial infections. This illustrates that the epidemic when superspreaders are present can be stochastic, even when considering populations sizes of 1 million.
Role of the funding source: The funders of the study had no role whatsoever in the study design, data collection, data interpretation, or writing of this paper. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Results
The agent-based perspective on the spread of an epidemic in a population allows us to address different types of questions and assess different mitigation strategies than what can be done with the more commonly used SEIR models. Naturally an agent based model resembles its corresponding SEIR variant when the agents are simulated without any social structure.
Model without social structure
First we modeled a COVID-19 like epidemic without “superspreaders” (SS) and social structure. Figure 2a) at 100% contact rate corresponds to such an idealization. This closely resembles the corresponding predictions for a SEIR model without social structure. When we reduced the contact rate so that only 50% and 25% of contacts remain, then the Covid-19 epidemic was strongly mitigated or halted, respectively. The largest reduction was enough to bring the effective R0 (Re) below 1.
Next we introduced “superspreaders” into the model to examine the effect of this phenomenon. The simulations in Fig. 2b assume that 10% of the agents are 50 times more infectious than the rest. Using the model we find that 85% of all transmission events are from these highly infectious agents. Although the unmitigated epidemic progresses in a qualitative way similar to the classical model in panel a), one is now able to halt the epidemic by only reducing the maximal number of encounters. In fact, if only 25% of all contacts remain for the superspreader, one can halt the epidemic. Thus, when they are prevented from acting as such, a dramatic level of mitigation can be achieved while everybody else can live their normal lives. This is because, as superspreaders infect many times more than the other agents, a reduction in the maximum number of contacts per agent only affects them. Of course, this scenario is difficult to achieve in real life since we do not know the identify of the superspreaders.
Model with social structure
Next we considered a structured society, by introducing into the model three sectors where interactions can take place (Fig. 1). This will allow us to introduce an indirect cap on the maximal social activity without knowing the identity of the superspreaders.
In Fig. 3 we show simulated epidemics for a socially structured model without superspreaders. The infection rate is again adjusted to fit the overall growth rate of 23%, and we observe an epidemic that resembles the basic model in Fig 2a). We examined, one at a time, the effect of lock-downs on each of the three contact sectors (“Home”, “Work/school”, and “Other”) (see Fig. 3b-d). In panel b), when disallowing any social contact in the “Home” sector, we find that the epidemic peak size and ICU usage are reduced to 60% of their non-mitigated value. Next, when closing down all contacts in the “Work/schools” sector, the peak reduced to 40% of the unmitigated size; ICU usage is lower due the fact that the elderly are not being directly affected by this measure. Mitigating “Work/school” contacts has a larger impact than reducing “Home” contacts because the associated network is greater. Finally, when closing the “Other” sector the epidemic peak again reduces to 40% of its original size but we found an even larger reduction in ICU utilization (Fig 3d). This reects the fact that the elderly are largely connected to the society through \Other” contacts.
We found that the model predicts final epidemic attack rates of between 80% and 100% of the population, with or without mitigation in the: “Home” and Work/school” sectors (Figure 3e). However, when closing down the “Other” sector, the attack rate is further reduced and contagion to the elderly is nearly eliminated.
Finally we introduced superspreaders into the socially structured model. Fig. 4 recapitulates the simulations shown in Fig. 3, but this time with 10% of the population being superspreaders that are 50 times more infective. These 10% superspreaders caused 75% of all infections.
Figure 4 panels b) and c) show that in this superspreader model the elimination of Home or Work/school contacts, respectively, had less influence on the epidemic peak size than in the model without superspreaders (see Figure 3). Even if one eliminated all contacts in both “Home” and “Work/school” sectors, the peak epidemic size would be substantial (not shown), corresponding to an Re value of about 1.7. In contrast, Fig. 4 shows that the epidemic is halted by eliminating “Other” contacts.
We thus conclude that superspreaders matter greatly in terms of the effect of mitigation strategies. Epidemics driven by a small proportion of highly infectious people are less sensitive to reductions in close contacts at Home or Work/school, but are highly affected by changes in their random contacts represented by the “Other” sector.
Sensitivity analyses of infection heterogeneity
Next, we asked ourselves whether this striking prediction was sensitive to assumptions about the frequency and intensity of superspreaders. To simplify our analysis, we considered lockdown in each sector in prioritized order: first “Other”, then “Work/ school” until the epidemic was contained. We do not consider “home” lockdown as a realistic option.
We first allowed the superspreader fraction to vary between 0% and 20% while retaining their 50x greater infection intensity. In Figure 4e), the reduction in contacts needed to halt the epidemic is indicated in light grey. Thus, the boundary of the light grey region marks epidemic control (Re = 1). It can be seen that when the fraction of superspreaders decreases below 8% the epidemic cannot be stopped by closing down the “Other” sector alone. Furthermore, the required lockdown increases substantially when the super-spreader fraction is small. Below ~4%; the superspreaders only contribute little and are no longer important.
We next explored the effect of considering the infection intensity distribution rather than an either-or scenario (Panel f). Each person was at the start of the simulation assigned a particular infectivity drawn from a power law distribution, limited by an upper limit of 1000x of the minimum. When this distribution is wider than about 1/s1.8, corresponding to 25% of superspreaders causing 80% of the infections, the Covid-19 epidemic can be mitigated by a full lockdown of the “Other” sector. If the distribution is wider, for example 1/s1.5, corresponding to 16% of infected causing 80% of infections, one can stop the epidemic while still allowing 14% of the “Other” contacts.
The supplement include sensitivity to variation in infectivity described by a Gamma distribution, and finds that mitigation of Covid-19 remains possible within the other sector. We even find this to be possible if we assume a society that from the beginning had a smaller “other” sector.
Simulations of real life lockdown in Denmark
To further explore the importance of superspreaders, we applied our model to the real life situation that took place in Denmark after a lock-down on March 13. Here, locking down a population of 5,8M was achieved early in the epidemic, and was followed one month later by a limited reopening in the Schools/work sector. After an additional 40 days, a second more extensive reopening took place. Throughout this whole period, the population continued general social distancing and hygiene measures.
First, we considered the model without superspreaders (Fig. 5a) and assume that the first lockdown imposed a reduction in the Work/school sector of about 25% of the pre-pandemic level, and that a month later this is eased to about 50% of the pre-pandemic level. Then in the second reopening we assumed all restrictions on the Work/school sector were dropped. We furthermore assumed that the “Other” sector would be reduced and reinstated in same proportions first down to 25%, then up to 50% and finally back to the original 100%). Even so these reductions are not enough to reproduce the observation that the epidemic virtually came to a halt after the initial lockdown. In the simulations in figure 5a), in order to fit the observed inpatient data we introduced an additional social distancing factor of 1/3 in both the Work/school and “Other” sectors. This additional factor 1/3 can be understood as the effect of the general hygiene measures and social distancing. Now the model traced the observed inpatient data (from early March to mid-May). Importantly, this was not enough to halt the epidemic in a scenario of complete reopening, in which case we predict a large epidemic surge (Figure 5a). These predictions are very sensitive to the fudge factor of 1/3 employed to account for the effects of social distancing and hygienic measures.
Next we introduced superspreaders into the model. We found that the trajectory is more stochastic (see 3 simulation with same parameters in Figure 5b)). With the presence of superspreaders it is sufficient to assume that the general distance/hygiene factor only applies to the “Other” sector with a factor 1/6. Thus one can explain the trajectory of the epidemic without imposing social distancing in the Work/school sector. This emphasizes that the “Other” sector reduction could be enough to control the epidemic i.e. making Re < 1. More realistic simulations should of course include a combination of these factors. Nevertheless, we have demonstrated a clear effect of having super spreaders in the model.
Focusing the reduction on the “Other” sector only allows us to interpret these results as the reduction that could be accomplished by closing down venues where many strangers meet (festivals, public transport, soccer games). Panel c) illustrates this effect of opening the “Other” sector to half the pre-pandemic level after the 2nd reopening: the epidemic takes off again.
The case of Sweden
Finally, we considered the unusual case of Sweden, a country famous for never having locked down their country, but rather depended on encouraging social distancing and hygiene procedures in the general population, as well as a ban on large gatherings (>50 persons), and recommendations to avoid non-essential travel within the country. Interestingly, Sweden has with this strategy been able to mitigate the epidemics with an effective Re that has been below 1 since mid-April [18].
In Fig. 5d) we introduced a reduction of contacts in the “Other” sector by a factor 1/6 early on (see red band in panel d). We found that this assumption in principle allows us to contain the epidemics and bring Re sufficiently down to start seeing effects of herd immunity. We emphasize that this herd-immunity can, according to our model, only function in an environment where the “Other” sector is small. Therefore, an eventual increase in this type of contact will probably lead to a second epidemics wave.
Our model can fit the observed epidemic in Sweden with a factor 1/6 reduction in “Other”, or, alternatively, by a factor 1/4 reduction in both the “Work/school” and “Other” sectors (See Figure S1). This second alternative is rather unappealing, given that schools in Sweden are largely open.
Discussion
SARS-CoV-2 is an intriguing pathogen that has challenged our understanding of disease transmission already from its beginning, in Wuhan. How could the epidemic be so slow to spread in the epicenter in Wuhan, with an estimated growth of 10% a day [19], while it later expanded at a rate of 30% per day. Why did it spread so slowly in Japan, and why did it not yet cause a major epidemic in Africa, while causing disaster with growth in death rates of up to 40%/day in northern Italy and Spain. Are we missing something?
Superspreading is a phenomenon known to occur for coronavira outbreaks of SARS, MERS and now has also been described for COVID-19. Bursts of infections after a choir rehearsal in the United States and South Korean religious meetings are examples of such events. Was it the elimination of superspreading opportunities that ended the emerging SARS-CoV outbreak in 2003? Such questions are not easily addressable with the standard SEIR-type dynamic disease models.
In the present work we have presented an agent-based model to investigate the effect of superspreaders in in the evolution of the epidemics and in various proposed mitigation schemes [1, 3]. Since we do not know how to identify a superspreader individual, we chose to model superspreading in terms of a fixed 10% proportion of individuals transmitting 50-fold more effectively than others. Using this model we could demonstrate that the success of a lockdown and the effects of subsequent re-openings are highly dependent on this heterogeneity of infectivity in the population. The prevalence of 10% superspreaders in our base model causing about 75% of infections is in agreement with recent findings by others [2, 3] where they estimate that ~10% of individuals cause 80% of infections. In sensitivity analyses, we explored other possibilities of infection heterogeneity and found that at least half of all transmissions must be from superspreaders.
Assuming arbitrarily that interactions are evenly distributed in the three sectors – Home, Work/school and “Other” – we investigated the effect of superspreaders following a lockdown of such a socially structured environment. We found that in a mitigated epidemic superspreaders matter tremendously. Our agent-based model allowed us to compare the effect of mitigation in each sector separately. Most importantly we found that mitigating the “Other” sector can have a profound impact, even halt the epidemics. This we interpret to mean that limiting random social contacts at large gatherings and in public transportation is what effectively can control the COVID-19 epidemic.
As a caveat, in our simulations we relegated all non-familiar contacts to an “Other” sector, and assumed all other contacts (with known persons) occurred through fixed social networks at home and work/school. In reality part of the interactions in the “Other” sector is of familiar persons, for example friends and extended family. On the other hand, part of the interactions in the “Work/school” sector would involve unknown persons, such as casual contacts at cafeterias, conferences and workshops. Importantly, all of this is not too concerning, because what really matters is how many different people a superspreader agent meets during the period where it infects excessively.
When we fit our model to actual observed Danish COVID-19 new inpatient data from early March to mid-May we found that our one-month forecast of the epidemic trajectory following the reopening of the country looked far more promising with than without superspreading, providing continued limits on contacts in the “Other” sector continued being applied (see Fig 5, b,c). It is mysterious that the epidemic did not begin again after the rather extensive phase one-reopening in mid-April; we suggest this is due to the continued ban on large gatherings and reduced use of public transportation. If this factor is relaxed in the future, we forecast a return to epidemic growth. In our simulations of the Swedish scenario (Fig. 5d), our super-spreader model fits well with the observed data under realistic assumptions. We predict that if the current strategy is maintained, the epidemic will come to a halt within a few months, owing to few superspreader opportunities. Importantly, we can only explain the control of the Swedish epidemic with superspreaders in the model; without them one has to assume unrealistically large levels of general social distancing.
About 4 months into the COVID-19 pandemic we see a heterogeneous epidemic patterns across the world, in which some cities are devastated while other regions are nearly free of disease. This could be due to the action of superspreaders. An epidemic take-off may well require a superspreader because the R0 is below 1 for the rest of the population. But if a superspreader is infected, the disease may spread to other superspreaders. Although hard to prove, it would not surprise us if the observed major differences across cities, regions and countries could be understood as a stochastic phenomenon modulated by a larger chance for the epidemic to percolate among superspreaders in large cities.
It is in this perspective that we propose that COVID-19 models used to describe the epidemic trajectory and forecast the effect of mitigation strategies and re-openings are inaccurate unless they include superspreaders. Without this element, such models will easily overestimate the epidemic size as well as the mitigation intensity needed in order to regain control.
We conclude that there is an urgent need of including heterogeneity and superspreading phenomena in our understanding of the present COVID-19 epidemic, and a enormous need for experimentally determining the abundance, mechanism and strength of interaction of the superspreader phenomenon.
Data Availability
No data.
Acknowledgements
KS thanks Raul Donangelo and Andreas Eilersen for enlightening discussions and corrections to the manuscript, and LS thanks Viggo Andresen for insights.