Estimating local outbreak risks and the effects of non-pharmaceutical interventions in age-structured populations: SARS-CoV-2 as a case study

Francesca A. Lovell-Read, Silvia Shen, Robin N. Thompson
doi: https://doi.org/10.1101/2021.04.27.21256163

Abstract

During the COVID-19 pandemic, non-pharmaceutical interventions (NPIs) including school closures, workplace closures and social distancing policies have been employed worldwide to reduce transmission and prevent local outbreaks. However, transmission and the effectiveness of NPIs depend strongly on age-related factors including heterogeneities in contact patterns and pathophysiology. Here, using SARS-CoV-2 as a case study, we develop a branching process model for assessing the risk that an infectious case arriving in a new location will initiate a local outbreak, accounting for the age distribution of the host population. We show that the risk of a local outbreak depends on the age of the index case, and we explore the effects of NPIs targeting individuals of different ages. Social distancing policies that reduce contacts outside of schools and workplaces and target individuals of all ages are predicted to reduce local outbreak risks substantially, whereas school closures have a more limited impact. In the scenarios considered here, when different NPIs are used in combination the risk of local outbreaks can be eliminated. We also show that heightened surveillance of infectious individuals reduces the level of NPIs required to prevent local outbreaks, particularly if enhanced surveillance of symptomatic cases is combined with efforts to find and isolate nonsymptomatic infected individuals. Our results reflect real-world experience of the COVID-19 pandemic, during which combinations of intense NPIs have reduced transmission and the risk of local outbreaks. The general modelling framework that we present can be used to estimate local outbreak risks during future epidemics of a range of pathogens, accounting fully for age-related factors.

1. Introduction

Throughout the COVID-19 pandemic, policy makers worldwide have relied on non- pharmaceutical interventions (NPIs) to limit the spread of SARS-CoV-2. Commonly introduced NPIs have included school closures, workplace closures and population-wide social distancing policies, all of which aim to reduce the numbers of contacts between individuals and disrupt potential chains of transmission [1–4]. Similar measures have previously been adopted for countering other infectious diseases such as Ebola and pandemic influenza [5–7], and are likely to remain a key line of defence against emerging pathogens that are directly transmitted between hosts. NPIs are particularly important when no effective treatment or vaccine is available, and they are also beneficial when vaccination programmes are being rolled out [8–10]. If vaccines do not prevent transmission completely, then NPIs may be important even when vaccination is widespread [11]. However, the negative economic, social and non-disease health consequences of NPIs have been widely discussed, with the impact of school closures on the academic progress and wellbeing of school-aged individuals a particular concern [7, 12–16]. Therefore, assessing the effectiveness of different NPIs at reducing transmission is critical for determining whether or not they should be used.

Since NPIs such as school and workplace closures affect distinct age groups within the population, when evaluating their effectiveness it is important to account for age-dependent factors that influence transmission. Multiple studies have documented marked heterogeneities in the patterns of contacts between individuals in different age groups, with school-aged individuals tending to have more contacts each day than older individuals [17–23]. Since close contact between individuals is a key driver of transmission for respiratory pathogens such as influenza viruses and SARS-CoV-2, these contact patterns influence transmission dynamics and consequently the effects of interventions that target different age groups [18, 19, 24–27]. Additionally, many diseases are characterised by significant age-related variations in pathophysiology. For example, for SARS-CoV-2, children may be less susceptible to infection than adults [27–31], and more likely to experience asymptomatic or subclinical courses of infection [28, 31–36]. Since the secondary attack rate (the proportion of close contacts that lead to new infections) from asymptomatic or subclinical hosts is lower than from hosts with clinical symptoms [37–42], children are likely to be less infectious on average than older individuals who are at increased risk of developing symptoms [43–45].

Previous studies have used age-stratified deterministic transmission models to investigate the effects of NPIs on COVID-19 epidemic peak incidence and timing. Prem et al. [46] projected the outbreak in Wuhan, China, over a one year period under different control scenarios, and demonstrated that a period of intense control measures including school closures, a 90% reduction in the workforce and a significant reduction in other social mixing could delay the epidemic peak by several months. Zhang et al. [27] predicted that eliminating all school contacts during the outbreak period would lead to a noticeable decrease in the peak incidence and a later peak; however, they did not take differences between symptomatic and asymptomatic cases into account explicitly. In contrast, Davies et al. [31] used estimates of age-dependent susceptibility and clinical fraction fitted to the observed age distribution of cases in six countries to demonstrate that school closures alone were unlikely to reduce SARS-CoV-2 transmission substantially. Davies et al. [47] subsequently concluded that a combination of several strongly enforced NPIs would be necessary to avoid COVID-19 cases exceeding available healthcare capacity in the UK.

Rather than considering the entire epidemic curve, here we focus on estimating the probability that cases introduced to a new location trigger a local outbreak as opposed to fading out with few cases. Localised clusters of transmission have been a feature of the COVID-19 pandemic [48–50], and assessing the risk that such local outbreaks occur requires a stochastic model in which the pathogen can either invade or fade out. Stochastic branching process models have been applied previously to assess outbreak risks for many pathogens without considering different age groups explicitly [51–55], and extended to consider adults and children as two distinct groups [56]. However, the significant heterogeneities in contact patterns and pathophysiology between individuals across the full range of ages have never previously been considered in estimates of local outbreak risks. Here, we develop an age-structured branching process model that can be used to estimate the probability of a local outbreak occurring for index cases of different ages, and demonstrate how the age-dependent risk profile changes when susceptibility to infection and clinical fraction vary with age.

We use the model to investigate the effects on the local outbreak probability of NPIs that reduce the numbers of contacts between individuals. Specifically, we use location-specific contact data for the UK detailing the average numbers of daily contacts occurring in school, in the workplace and elsewhere [17] to model the impacts of school closures, workplace closures and broader social distancing policies. We demonstrate that, for SARS-CoV-2, contacts occurring outside schools and the workplace are a key driver of sustained transmission. Thus, population-wide social distancing policies that affect individuals of all ages lead to a substantial reduction in the risk of local outbreaks. In contrast, since school-aged individuals only make up around one quarter of the UK population and tend to have large numbers of contacts outside school, school closures are predicted to have only a limited effect when applied as the sole NPI.

We then go on to consider the impacts of mixed strategies made up of multiple NPIs, as well as additional NPIs that do not only reduce numbers of contacts. Specifically, we show that rigorous surveillance and effective isolation of infected hosts can reduce the level of contact-reducing NPIs required to achieve substantial reductions in the risk of local outbreaks. Although we use SARS-CoV-2 as a case study, our approach can be applied more generally to explore the effects of NPIs on the risk of outbreaks of any pathogen for which age-related heterogeneities play a significant role in transmission dynamics.

2. Methods

2.1 Mathematical model

We considered a branching process model in which the population was divided into 16 age groups, denoted G1, G2, … , G1#. The first 15 groups represent individuals aged 0-74, divided into five-year intervals (0-4, 5-9, 10-14 etc.). The final group represents individuals aged 75 and over. The total number of individuals in age group Gk is denoted Nk. Infected individuals in each age group Gk are classified into compartments representing asymptomatic (Ak), presymptomatic (Pk) or symptomatic (Sk) hosts, where an individual in the Ak compartment does not develop symptoms at any time during their course of infection.

An infected individual of any type in group Gk may generate new infections in any age group. In our model, the rate at which a single infected symptomatic individual in group Gk generates infections in group Gj is given by Embedded Image Here, τk represents the infectivity of individuals in group Gk, ωj represents the susceptibility to infection of individuals in group Gj, Ckj represents the daily number of unique contacts a single individual in group Gk has with individuals in group Gj, and B is a scaling factor that can be used to set the reproduction number of the pathogen being considered (see Section 2.2). Since the initial phase of potential local outbreaks are the focus of this study, we did not account for depletion of susceptible hosts explicitly. The relative transmission rates from presymptomatic and asymptomatic individuals compared to symptomatic individuals are given by the scaled quantities ηβkjand θβkj, respectively, where η and θ were chosen so that the proportions of transmissions generated by presymptomatic and asymptomatic hosts were in line with literature estimates [57]. The parameter ξk represents the proportion of asymptomatic infections in group Gk, so that a new infection in group Gk either increases Ak by one (with probability ξk) or increases Pk by one (with probability 1 − ξk).

A presymptomatic individual in group Gk may go on to develop symptoms (transition from Pk to Sk) or be detected and isolated (so that Pk decreases by one). A symptomatic individual in group Gk may be detected and isolated as a result of successful surveillance, or may be removed due to self-isolation, recovery or death (so that Sk decreases by one in either case). Similarly, an asymptomatic individual in group Gk may be detected and isolated or recover (so that Ak decreases by one). A schematic of the different possible events in the model is shown in Fig 1.

Fig 1.
Fig 1. The branching process model used in our analyses.

Schematic showing the different possible events in the branching process model and the rates at which they occur. The parameters of the model are described in the text and in Tables 1, 2 and 3.

View this table:
Table 1. Baseline values of age-dependent parameters.

Values used for the age-dependent relative susceptibility to infection (ωk) and the proportion of infections that are asymptomatic (ξk) for each of the scenarios A, B and C [31].

View this table:
Table 2. Baseline values of scenario-independent parameters.

Values used for the parameters that were assumed not to vary between scenarios A, B and C. We also considered strategies involving enhanced surveillance (ρk, σk > 0) – see Fig 6.

View this table:
Table 3. Baseline values of scenario-dependent scaling parameters.

Values used for the scaling parameters B, η and θ for each of the scenarios A, B and C.

The parameter λ represents the rate at which presymptomatic individuals develop symptoms, so that the expected duration of the presymptomatic infectious period is 1/λ days in the absence of surveillance of nonsymptomatic infected individuals. Similarly, the expected duration of the asymptomatic infectious period in the absence of surveillance is 1/λ days. The parameter μ represents the rate at which symptomatic individuals are removed as a result of self-isolation, recovery or death, so that the duration of time for which they are able to infect others is 1/μ days.

For each group Gk, the rate at which symptomatic individuals are detected and isolated as a result of enhanced surveillance is determined by the parameter ρk. Analogously, the parameter σk governs the rate at which presymptomatic and asymptomatic individuals in Gk are detected and isolated. We assumed that surveillance measures targeting nonsymptomatic hosts are equally effective for those who are presymptomatic and those who are asymptomatic, and therefore used the same rate of isolation due to surveillance for both of these groups.

2.2 Reproduction number

The effective reproduction number, R, represents the expected number of secondary infections generated by a single infected individual during their entire course of infection, accounting for interventions that are in place. Here, we take a heuristic approach to derive the following expression for R: Embedded Image where N = N1 + ⋯ + N1# is the total population size. To obtain this expression, we first consider the expected number of secondary infections an infected individual in age group Gk will generate in age group Gl. If an individual in age group Gk experiences a fully asymptomatic course of infection, which occurs with probability ξk, they will generate new infections in age group Gl at rate θβkl and recover or be isolated at rate λ + σk. Therefore, the total number of infections they are expected to cause in age group Gl is θβkl/(λ + σk). If instead the individual in age group Gk experiences a symptomatic course of infection, which occurs with probability 1 − ξk, whilst presymptomatic they will generate new infections in age group Gl at rate ηβkl and be isolated or develop symptoms at rate λ + σk. Thus, the total number of infections they are expected to cause in age group Gl whilst presymptomatic is ηβkl/(λ + σk). If they go on to develop symptoms before being isolated, which occurs with probability λ/(λ + σk), applying similar reasoning they are expected to cause βkl/(μ + ρk) new infections in age group Gl whilst symptomatic.

Combining these possibilities leads to the term in square brackets in expression (1), which is then summed over all possible age groups Gl of the infectee. Finally, to obtain the full expression (1) we take a weighted average across all possible age groups Gk of the infector, where the weights Nk/N represent the proportions of the population belonging to each age group. This corresponds to the assumption that the initial infected host is more likely to belong to an age group containing more individuals than an age group with fewer individuals.

In the absence of interventions, i.e. when σk = ρk = 0 (representing no enhanced isolation as a result of surveillance) and βkl is calculated using contact patterns that are characteristic of normal behaviour, the effective reproduction number, R, is equal to the basic reproduction number, R(.

2.3 Model parameterisation

The numbers of individuals in each age group (values of Nk) were chosen according to United Nations age demographic data for the UK [58] (Fig 2A). The daily numbers of contacts between individuals in each age group (values of Ckj) were set according to the 16x16 ‘contact matrix’ for the UK, in which the (k, j)th entry represents the expected daily number of unique contacts an individual in age group Gk has with individuals in age group Gj [17]. In addition to matrices representing ‘all’ contacts (Figure 2B), we also considered matrices detailing a breakdown into ‘school’, ‘work’, ‘home’ and ‘other’ contacts (Figures 2C-F), allowing us to investigate the effects of control interventions that reduce contacts in each of these settings.

Fig 2.
Fig 2. Age demographic and age-structured contact patterns for the United Kingdom.

A. United Nations age demographic data for the UK in 2020, split into five-year age groups (where the final age group contains all ages 75+) [58]. B. Heat map of UK ‘all contacts’ matrix, representing the expected daily number of unique contacts that an individual in each age group Gk has with individuals in each other age group G& [17]. C. The analogous figure to B, but showing only the subset of ‘all’ contacts that occur in schools (‘school’ contacts). D. The analogous figure to C, but showing only ‘workplace’ contacts. E. The analogous figure to C, but showing only ‘home’ contacts. F. The analogous figure to C, but showing only ‘other’ contacts (i.e. all contacts outside schools, workplaces or homes).

Since we considered SARS-CoV-2 as a case study, we used studies conducted during the COVID-19 pandemic to inform the epidemiological parameters of our model. Despite previous research assessing the relationships between age and factors such as susceptibility to SARS- CoV-2 infection or the propensity to develop symptoms, there is some variation in estimated parameters between different studies. To test the robustness of our results to this uncertainty, we conducted our analyses under three different scenarios (A, B and C). In scenario A, we assumed that susceptibility to infection (values of ωk) and the proportion of hosts who experience a fully asymptomatic course of infection (values of ξk) are independent of age. In scenario B, susceptibility was assumed to vary with age but the proportion of asymptomatic infections is independent of age. In scenario C, we allowed both susceptibility and the asymptomatic proportion to vary with age. The values used for the parameters ωk and ξk in each of these three scenarios are shown in Table 1 (see also [31]).

In all scenarios considered, the inherent infectivity was not assumed to be age-dependent (i.e. τk = 1 for all values of k). In other words, the expected infectiousness of infected hosts in different age groups was governed solely by the proportion of asymptomatic infections in that age group. We chose the scaling factors η and θ for the relative transmission rates from presymptomatic and asymptomatic individuals compared to symptomatic individuals so that the proportions of infections arising from each of these groups were in line with literature estimates (see Table 3 and [57]).

In the absence of enhanced isolation, we set the expected duration of the presymptomatic infectious period and the time for which symptomatic individuals are able to infect others to be 1/λ = 2 days and 1/μ = 8 days, respectively [59–62]. The asymptomatic infectious period was then chosen so that all infected individuals are expected to be infectious for the same period (i.e. 1/λ = 10 days). In our initial analysis, we set the isolation rates ρk and σk equal to 0 for all k; later, we considered the effects of increasing these rates.

Initially, we fixed R( = 3 (in line with initial estimates of SARS-CoV-2 transmissibility [63–66], before the emergence of more transmissible variants) and used expression (1) to determine the appropriate corresponding value of the scaling factor B. Later, when considering the impact of NPIs on the probability of a local outbreak, we retained this value of B and used expression (1) to determine how the reproduction number changes as a result of the control implemented.

2.4 Probability of a local outbreak

The probability that an infected individual in a particular age group initiates a local outbreak when they are introduced into the population was calculated using the branching process model. One possible approach for approximating the age-dependent local outbreak probability using a branching process model is to run a large number of stochastic simulations of the model starting from a single infected individual in a particular age group, and record the proportion of simulations in which the pathogen does not fade out after only a small number of infections [67]. This would then need to be repeated for index cases of different ages. Here, we instead take an analytic approach, and derive a nonlinear system of simultaneous equations that determine the age-dependent outbreak probabilities, as described below. The local outbreak probabilities are then obtained by solving these equations numerically, and are analogous to the probabilities that would be derived from the simulation approach in the limit of infinitely many simulations. The benefit of our analytic approach is that it does not require a large number of stochastic simulations to be run.

The probability of a local outbreak not occurring (i.e. pathogen fadeout occurs), starting from a single symptomatic (or presymptomatic, asymptomatic respectively) infectious individual in age group Gk, was denoted by xk (yk, zk). Beginning with a single symptomatic individual in Gk, the possibilities for the next event are as follows:

  1. The infected individual in Gk infects a susceptible individual in Gj, so that either Aj increases by one (with probability ξj) or Pj increases by one (with probability (1 − ξj). This occurs with probability Embedded Image

  2. The infected individual in Gk recovers, dies or is isolated before infecting anyone else, so that Sk decreases to zero (and there are no infected individuals left in the population). This occurs with probability Embedded Image

If there are no infectious hosts present in the population, then a local outbreak will not occur. Therefore, assuming that chains of transmission arising from infectious individuals are independent, the probability that no local outbreak occurs beginning from a single symptomatic individual in Gk is Embedded Image Similarly, beginning instead with a single presymptomatic individual in Gk, the possibilities for the next event are:

  1. The presymptomatic infected individual in Gk infects a susceptible individual in Gj, so that (as before) either Aj increases by one (with probability ξj) or Pj increases by one (with probability (1 − ξj)). This occurs with probability Embedded Image

  2. The infected individual in Gk develops symptoms (transitions from Pk to Sk). This occurs with probability Embedded Image

  3. The infected individual in Gk is isolated before infecting anyone else, so that Sk decreases by one. This occurs with probability Embedded Image

Therefore, the probability that no local outbreak occurs beginning from a single presymptomatic individual in Gk is Embedded Image Similarly, the probability zk that a local outbreak does not occur starting from a single asymptomatic individual in Gk satisfies the equation Embedded Image Where Embedded Image

The system of simultaneous equations (2) − (4) can be solved numerically to obtain xk, xk and zk(here, we did this using the MATLAB nonlinear system solver ‘fsolve’). Specifically, we take the minimal non-negative solution, as is standard when calculating extinction probabilities using branching process models [55, 68]. Then, for each k, the probability of a local outbreak occurring beginning from a single symptomatic (or presymptomatic, asymptomatic respectively) individual in group Gk is given by 1 − xk (1 − xk, 1 − zk).

Throughout, we consider the probability pk of a local outbreak occurring beginning from a single nonsymptomatic individual in group Gk arriving in the population at the beginning of their infection: Embedded Image The average local outbreak probability, P, which is defined as the probability of a local outbreak when the index case is chosen randomly from the population, is also considered. The value of P is therefore a weighted average of the pk values, where the weights correspond to the proportion of the population represented by each group: Embedded Image This reflects an assumption that the index case is more likely to come from an age group with more individuals than an age group with fewer individuals.

All computing code used to implement the above methods was written in MATLAB version R2019a, and is available at https://github.com/francescalovellread/age-dependent-outbreak-risks.

3. Results

3.1 Effect of the age of the index case on the risk of a local outbreak

We first considered the probability that a single infected individual in a particular age group Gk initiates a local outbreak when introduced into a new host population. This quantity was calculated for each of the three scenarios A, B and C (Fig 3).

Fig 3.
Fig 3. The probability of a local outbreak depends on the age of the index case.

A. The probability that a single infected individual in any given age group triggers a local outbreak (grey bars) for scenario A, in which susceptibility and clinical fraction are assumed constant across all age groups. The weighted average local outbreak probability P is shown by the black horizontal line. B. The analogous figure to A but for scenario B, in which clinical fraction is assumed constant across all age groups but susceptibility varies with age (Table 1). C. The analogous figure to A but for scenario C, in which both susceptibility and clinical fraction vary with age (Table 1).

In scenario A, the variation in the local outbreak risk for introduced cases of different ages is driven solely by the numbers of contacts between individuals. As a result, due to their higher numbers of daily contacts, school- and working-age individuals are more likely to trigger a local outbreak than children under five or adults over 60, with index cases aged 15-19 posing the highest risk (Fig 3A). These findings do not change significantly when susceptibility is allowed to vary with age in scenario B (Fig 3B). However, in scenario C, assuming that the clinical fraction also varies between age groups alters the age-dependent risk profile substantially. This is because asymptomatic individuals are assumed to be less infectious than symptomatic individuals, and therefore an index case in an age group with a high proportion of asymptomatic infections is less likely to initiate a local outbreak. In this scenario, index cases aged 40 or over had a disproportionately high probability of generating a local outbreak, with individuals aged 70-74 presenting the highest risk (Fig 3C). These individuals are more likely to develop symptoms than younger individuals (Table 1), leading to a higher expected infectiousness. In contrast, individuals under the age of 40 had a below average probability of generating a local outbreak, with individuals aged 10-14 presenting the lowest risk. Noticeably, individuals aged 5- 19 presented relatively low risks, despite the high numbers of contacts occurring among these age groups (Fig 2B). In this scenario, the large number of contacts was offset by the fact that individuals in these age groups are more likely to be asymptomatic and consequently less infectious than older individuals (Table 1). Therefore, an index case in one of these age groups is likely to lead to fewer secondary transmissions. Furthermore, the contact patterns between individuals in these age groups are highly assortative with respect to age (Fig 2B). Therefore, in addition to the index case being less infectious, a high proportion of the contacts they make are with individuals who are also likely to be less infectious, as well as being less susceptible to infection in the first place.

We performed our subsequent analyses for each of the three scenarios A, B and C, with qualitatively similar results. The figures shown in the main text are for scenario C, since we deem this scenario to be the most realistic for SARS-CoV-2 transmission, but the analogous results for scenarios A and B are presented in Supplementary Figs S1-6.

3.2 Effect of the target age group on NPI effectiveness

We next considered the effects of NPIs that reduce the number of contacts between individuals on the probability that an introduced case will lead to a local outbreak. To approximate the relative effects of school closures, workplace closures and population-wide social distancing policies, we calculated the age-dependent risk profiles when each of these types of contact were excluded from the overall contact matrix.

First, we removed all ‘school’ contacts from the total contact matrix (Fig 4A). For scenario C, removing ‘school’ contacts led to a 4.2% reduction in the average probability of a local outbreak (from 0.449 to 0.430). This small reduction is unsurprising for scenario C, since in that scenario school-aged infected individuals are assumed to be more likely to be asymptomatic than other infected individuals, and therefore their expected infectiousness is lower. However, even for scenarios A and B, in which school-aged individuals present the greatest risk of triggering a local outbreak, the effectiveness of removing ‘school’ contacts alone at reducing the local outbreak probability was limited (reductions of 7.2% and 4.75% respectively; see Supplementary Figs S1A, S4A). In each scenario, the reduction in risk was predominantly for school-aged index cases, with the risk from index cases of other ages only slightly reduced. Second, we considered the effects of removing ‘work’ contacts from the total contact matrix (Fig 4B). This led to a more substantial 25.4% reduction in the average probability of a local outbreak for scenario C (with corresponding reductions of 19.0% and 24.0% for scenarios A and B respectively; see Supplementary Figs S1B, S4B). As well as reducing the risk of a local outbreak from an index case of working age, removing ‘work’ contacts also reduced the probability of a local outbreak occurring starting from a school-aged individual. This is because closing workplaces helps to block chains of transmission that begin with an infected child. For example, a transmission chain involving a child transmitting to an adult at home, followed by subsequent spread around the adult’s workplace, will be less likely to occur. Third, we investigated the effect of removing all ‘other’ contacts, reflecting perfect social distancing being observed outside of the home, school or workplace (Fig 4C). This had the most significant effect of the three types of contact-reducing intervention considered, reducing the probability of a local outbreak by 41.7% for scenario C (and 30.7% or 33.2% for scenarios A and B, respectively).

Fig 4.
Fig 4. The effects of interventions that reduce contacts between individuals on the probability of a local outbreak.

A. The effect of removing all ‘school’ contacts on the probability of a local outbreak. Pale grey bars and black dash-dotted line represent the local outbreak probabilities without any contacts removed (as in Fig 3C). Red bars and the solid red line represent the local outbreak probabilities and their weighted average when ‘school’ contacts are removed. B. The analogous figure to A, but with all ‘work’ contacts removed. C. The analogous figure to A, but with all ‘other’ contacts removed. D. Partial reductions in ‘school’, ‘work’ and ‘other’ contacts, and the resulting reductions in the average local outbreak probability (solid red, dashed blue and dotted green lines respectively).

In the three cases described above, we considered complete reductions in ‘school’, ‘work’ and ‘other’ contacts, respectively. In practice, such complete elimination of contacts is unlikely. We therefore also considered partial reductions in ‘school’, ‘work’ and ‘other’ contacts, and compared the resulting reductions in the local outbreak probability (Fig 4D). For any given percentage reduction in contacts, reducing ‘other’ contacts always led to the largest reduction in the local outbreak probability (see also Supplementary Figs S1D, S4D). This suggests that reducing social contacts outside schools and workplaces can be an important component of strategies to reduce the risk of local outbreaks of SARS-CoV-2. However, this alone is not enough to eliminate the risk of local outbreaks entirely. For greater risk reductions using contact- reducing NPIs, a mixed approach involving combinations of reductions in ‘school’, ‘workplace’ and ‘other’ contacts is needed.

3.3 Mixed strategies for reducing the local outbreak risk

Next, we considered the effects of combining reductions in ‘school’, ‘work’ and ‘other’ contacts on the local outbreak probability (Fig 5; analogous results for scenarios A and B are shown in Supplementary Figs S2 and S5). We allowed reductions in ‘school’ and ‘work’ contacts to vary between 0% and 100% whilst ‘other’ contacts were reduced by 25%, 50% or 75% (Fig 5A,B,C, respectively).

Fig 5.
Fig 5. The effects of intervention strategies that combine reductions in ‘school’, ‘work’ and ‘other’ contacts.

A. The effect of reducing ‘school’ and ‘work’ contacts on the weighted average probability of a local outbreak (P), when ‘other’ contacts are reduced by 25% across all age groups. Red dotted lines indicate contours along which the local outbreak probability is constant. B. The analogous figure to A, but with a 50% reduction in ‘other’ contacts. C. The analogous figure to A, but with a 75% reduction in ‘other’ contacts. D. The effect of reducing ‘other’ contacts on the average local outbreak probability when ‘school’ and ‘work’ contacts are not reduced at all (dotted line) and when ‘school’ and ‘work’ contacts are reduced by 100% (solid line).

Since NPIs have negative economic, social and non-disease health consequences, policy makers may choose to implement public health measures in which the risk of local outbreaks is not eliminated completely. These results provide contact reduction targets for mixed strategies in which the local outbreak probability is reduced to a pre-specified ‘acceptable’ level. For example, to reduce the local outbreak probability to 0.25, ‘other’ contacts could be reduced by 25% from the baseline level, and ‘school’ and ‘home’ contacts reduced as indicated by the red dotted contour marked ‘0.25’ in Fig 5A. Alternatively, to achieve the same local outbreak risk, ‘other’ contacts can instead be reduced by 50% or by 75% with the degree of ‘school’ and ‘work’ reductions chosen according to the contours marked ‘0.25’ in Figs 5B,C respectively.

If a policy maker wishes to eliminate the local outbreak risk entirely using contact-reducing NPIs, for the model parameterisation considered very significant reductions in multiple types of contacts are needed in combination. For example, even if all ‘school’ and ‘work’ contacts are removed, ‘other’ contacts must be reduced by 66% for the overall average local outbreak probability to fall below 0.01 (Fig 5D). Since such substantial reductions in multiple types of contacts are unlikely to be possible, this suggests that contact-reducing NPIs must be combined with other interventions, such as effective surveillance and isolation strategies, to eliminate local outbreaks.

3.4 Effect of surveillance on the contact-reducing NPIs required for local outbreak control

We considered whether or not low local outbreak probabilities can be achieved using limited contact-reducing NPIs in combination with other interventions. Specifically, the effects of surveillance and isolation of infected individuals (through e.g. contact tracing) as well as reducing contacts in schools, workplaces and other locations, were assessed. While results are shown for scenario C in Fig 6, analogous results for scenarios A and B are presented in Supplementary Figs S3 and S6.

Fig 6.
Fig 6. Surveillance as part of a mixed strategy to reduce the local outbreak probability.

A. The effect of increasing the isolation rate of symptomatic (red line) or nonsymptomatic infected hosts (blue line) on the average probability of a local outbreak (P), in the absence of contact-reducing NPIs. The isolation rates ρk and σk are varied in turn between 0 days’1 and 1 days’1. B. The effect of simultaneously varying the isolation rate of symptomatic and nonsymptomatic hosts on the average probability of a local outbreak (P), again without contact-reducing NPIs. C. The age-dependent probability of a local outbreak when the isolation rate of symptomatic individuals is ρk = 1/2 days’1, without contact-reducing NPIs or surveillance of nonsymptomatic infected individuals (purple bars and solid line). Pale grey bars and black dash-dotted line represent the original local outbreak probabilities without any contact-reducing NPIs or enhanced surveillance (as in Fig 3C). D. The effect of reducing ‘work’ and ‘other’ contacts when the isolation rate of symptomatic infected individuals is ρk = 1/2 days’1, as in C, without surveillance of nonsymptomatic infected individuals. E,F. The analogous figures to C,D, with enhanced surveillance of both symptomatic and nonsymptomatic infected hosts (ρk = 1/2 days’1 and σk = 1/7 days’1).

Initially, we considered the effect of increasing the rate at which symptomatic and/or nonsymptomatic infected individuals are detected and isolated as a result of surveillance, in the absence of contact-reducing NPIs (i.e. with no reduction in the number of contacts between individuals compared to the baseline case in Fig 2B) (Fig 6A,B). For symptomatic hosts, this represents an enhanced rate of isolation compared to the baseline rate of self-isolation already present in the model. Isolation of nonsymptomatic hosts was more effective at reducing the local outbreak probability than isolation of symptomatic hosts (Figs 6A,B), although of course this is more challenging to achieve [55]. However, if fast isolation of nonsymptomatic hosts could be achieved through efficient large-scale testing (potentially in combination with contact tracing [69]), the probability of local outbreaks could be reduced substantially through this measure alone.

We then demonstrated the effects of combining contact-reducing NPIs with enhanced isolation of infected hosts due to infection surveillance. First, we increased the enhanced isolation rate of symptomatic individuals to ρk = 1/2 days)1. In the absence of other interventions, this reduced the local outbreak probability by 22.0% (Figure 6C). With this level of surveillance, the local outbreak risk could be reduced below 0.01 with a reduction in ‘work’ and ‘other’ contacts of around 73% each, for example.

Finally, keeping the enhanced isolation rate of symptomatic individuals equal to ρk = 1/2 days)1, we increased the isolation rate of nonsymptomatic individuals to σk = 1/7 days)1. In this case, the local outbreak probability without contact-reducing NPIs fell by 59.4% compared to a situation without enhanced surveillance (Fig 6E), and the reductions in ‘work’ and ‘other’ contacts needed to bring the local outbreak probability below 0.01 were significantly smaller (Figure 6F). For example, if ‘work’ contacts can be reduced by 50%, then ‘other’ contacts only need to be reduced by 43%. This indicates that effective surveillance of both symptomatic and nonsymptomatic individuals can substantially lower the extent of contact- reducing NPIs that are required to achieve substantial reductions in local outbreak risks.

4. Discussion

During the COVID-19 pandemic, public health measures that reduce the numbers of contacts between individuals have been implemented in countries globally. These measures include school closures, workplace closures and population-wide social distancing policies. Contact- reducing NPIs have been shown to be effective at reducing SARS-CoV-2 transmission, and have also been used previously during influenza pandemics [5, 70–72]. However, long-term implementation of these measures has negative social, psychological and economic consequences [7, 12–16]. It is therefore important to assess the effectiveness of different contact- reducing NPIs at lowering transmission and preventing local outbreaks, in order to design effective targeted control strategies that avoid unnecessarily strict measures.

Here, we constructed a branching process model to estimate the risk of local outbreaks under different contact-reducing NPIs and different levels of surveillance for symptomatic and nonsymptomatic infected individuals. Unlike previous approaches for estimating outbreak risks using branching processes [51–56], we considered the effects of age-related heterogeneities affecting transmission for infected individuals of a wide range of ages, including age-dependent variations in social mixing patterns, susceptibility to infection and clinical fraction. Using SARS- CoV-2 as a case study, we demonstrated that the risk that an introduced case initiates a local outbreak depends on these age-related factors and on the age of the introduced case (Fig 3), as well as the age-structure of the local population.

We used our model to assess the effects of reducing the numbers of contacts that occur in school, in the workplace and elsewhere. Of the three contact-reducing NPIs considered, removing ‘school’ contacts had the smallest effect on the probability of observing a local outbreak, even when age-dependent variations in susceptibility and clinical fraction were ignored (Figs 4A,D, Supplementary Figs S1A,D, and S4,A,D). This can be attributed to the fact that school closures predominantly reduce contacts between individuals aged 5-19, who only account for approximately 23% of the total population [58]. Additionally, these individuals tend to have large numbers of contacts outside of the school environment (Figs 2E,F). Therefore, compared to other measures, interrupting within-school transmission may have only a limited effect on transmission in the wider population, particularly when school-aged individuals are less susceptible to infection and more likely to experience subclinical courses of infection. In contrast, reducing contacts that occur outside schools or workplaces was the most effective intervention, significantly lowering the local outbreak risk across all age groups, and for those aged over 60 in particular (Figs 4C,D, and Supplementary Figs S1C,D and S4C,D). This could explain the success of social distancing strategies worldwide for reducing observed COVID-19 cases and deaths.

Mixed strategies combining reductions in ‘school’, ‘work’ and ‘other’ contacts led to greater reductions in the local outbreak probability than individual interventions (Figs 5A-C), but very large reductions in all three types of contact were required to eliminate the risk of local outbreaks entirely (Fig 5D). However, implementing effective surveillance to identify infected hosts (followed by isolation) led to substantial reductions in the risk of local outbreaks even in the absence of other control measures (Figs 6A,B). In the scenarios considered here, with an efficient surveillance strategy in place, significantly smaller reductions in ‘work’ and ‘other’ contacts were needed to render the local outbreak probability negligible, even when ‘school’ contacts were not reduced at all (Figs 6C-F). This supports the use of surveillance that targets both symptomatic and nonsymptomatic individuals, such as contact tracing and isolation strategies or population-wide diagnostic testing, to help prevent local outbreaks [55].

Although here we used SARS-CoV-2 as a case study, our model provides a framework for estimating the risk of local outbreaks in age-structured populations that can be adapted for other pathogens, provided sufficient data are available to parametrise the model appropriately. The effects of age-structure on local outbreak risks may vary for pathogens with different epidemiological characteristics. For influenza-A viruses, for example, susceptibility to infection tends to decrease with age, whilst the risk of an infection leading to severe symptoms is greater both for the elderly and for the very young [31, 73, 74]. This is in contrast to SARS-CoV-2, for which children are more likely to experience subclinical courses of infection. In this study, we used age demographic and contact data for the UK, but equivalent data for other countries are available and can easily be substituted into our model to estimate outbreak risks elsewhere [17, 58].

One caveat of the results for SARS-CoV-2 presented here is that, although the epidemiological parameters of our model were chosen to be consistent with reported literature estimates, there is considerable variation between studies. In particular, the precise age-dependent variation in susceptibility and clinical fraction remains unclear, and the relative infectiousness of asymptomatic, presymptomatic and symptomatic hosts has not been determined exactly.

Furthermore, the inherent transmissibility of SARS-CoV-2 is now higher than in the initial stage of the pandemic, due to the appearance of novel variants. To explore ongoing local outbreak risks due to SARS-CoV-2, it would be necessary to update the model to reflect the increased transmissibility of the Delta variant [75]. Due to the uncertainty in model parameter values, we conducted sensitivity analyses to explore the effects of varying the parameters of the model on our results (Supplementary Figs S1-12). In each case that we considered, our main conclusions were unchanged: the probability that an introduced case initiates a local outbreak depends on age-dependent factors affecting pathogen transmission and control, with widespread interventions and combinations of NPIs reducing the risk of local outbreaks most significantly.

An important limitation of our approach to modelling contact-reducing NPIs is that we made a standard assumption in our main analyses that ‘school’, ‘work’ and ‘other’ contacts are independent [27, 31, 47]. In other words, reducing the numbers of contacts in one location did not affect the numbers of contacts occurring in another. In reality, this is unlikely to be the case. For example, closing schools also affects workplace contacts, as adults may then work from home in order to fulfil childcare requirements. Additionally, the contact data that we used represent the number of unique contacts per day and do not reflect the numbers of repeated contacts with the same person, which affect the risk of transmission between individuals. These assumptions could in principle be removed, if relevant data become available – for example, data describing the effects of school closures on numbers of contacts in other locations. To demonstrate how changes in multiple types of contact due to NPIs could be implemented in the model, we conducted a supplementary analysis in which we investigated the effects of removing all ‘school’ contacts and allowing for concurrent changes in ‘work’, ‘home’ and ‘other’ contacts (Supplementary Figs S13-14). These results support our conclusion that school closures are unlikely to have a substantial impact on SARS-CoV-2 transmission when applied as the sole NPI. The benefit of school closures could even potentially be outweighed by the possible secondary effects on other types of contacts. An improved understanding of how NPIs affect different types of contact is important for more accurate assessments of interventions in future.

Despite the simplifications made, our model provides a useful framework for estimating the risk of local outbreaks and the effects of NPIs. Different measures can be considered in combination in the model to develop strategies for lowering local outbreak risks. Our results emphasise the importance of quantifying age-dependent factors that affect transmission dynamics, such as susceptibility to infection and the proportion of hosts who develop clinical symptoms, for individuals of different ages. As we have shown, it is crucial to take age-dependent factors into account when assessing local outbreak risks and designing public health measures.

Declaration of interests

We declare no competing interests.

Authors’ contributions

Conceptualisation: All authors. Methodology: FALR, RNT. Investigation: FALR, SS. Writing - original draft: FALR, RNT. Writing - review and editing: All authors. Supervision: RNT.

Data availability

All computer code used in this paper was written in MATLAB version R2019a, and is available at the following GitHub repository: https://github.com/francescalovellread/age-dependent-outbreak-risks.

Funding

FALR acknowledges funding from the Biotechnology and Biological Sciences Research Council (BBSRC), grant number BB/M011224/1. SS was funded by a BBSRC Research Experience Placement and a Rokos award for undergraduate research from Pembroke College Oxford. RNT was supported by UK Research and Innovation (UKRI) via grant EP/V053507/1. The funders had no role in the study design; the collection, analysis and interpretation of data; the writing of the report; or the decision to submit the article for publication.

Data Availability

All computer code used in this paper was written in MATLAB version R2019a, and is available at the following GitHub repository: https://github.com/francescalovellread/age-dependent-outbreak-risks.

https://github.com/francescalovellread/age-dependent-outbreak-risks

References

  1. 1.↵
    Perra, N. 2021 Non-pharmaceutical interventions during the COVID-19 pandemic: A review. Phys. Rep. 913, 1–52. (doi:10.1016/j.physrep.2021.02.001).
    OpenUrlCrossRef
  2. 2.
    Regmi, K. & Lwin, C.M. 2020 Impact of non-pharmaceutical interventions for reducing transmission of COVID-19: a systematic review and meta-analysis protocol. BMJ Open 10, e041383. (doi:10.1136/bmjopen-2020-041383).
    OpenUrlAbstract/FREE Full Text
  3. 3.
    Imai, N., Gaythorpe, K.A.M., Abbott, S., Bhatia, S., van Elsland, S., Prem, K., Liu, Y. & Ferguson, N.M. 2020 Adoption and impact of non-pharmaceutical interventions for COVID-19. Wellcome Open Res 5, 59–59. (doi:10.12688/wellcomeopenres.15808.1).
    OpenUrlCrossRef
  4. 4.↵
    Thompson, R.N., Hollingsworth, T.D., Isham, V., Arribas-Bel, D., Ashby, B., Britton, T., Challenor, P., Chappell, L.H.K., Clapham, H., Cunniffe, N.J., et al. 2020 Key questions for modelling COVID-19 exit strategies. Proc. R. Soc. Lond. B. Biol. Sci. 287, 20201405. (doi:10.1098/rspb.2020.1405).
    OpenUrlCrossRefPubMed
  5. 5.↵
    Hatchett, R.J., Mecher, C.E. & Lipsitch, M. 2007 Public health interventions and epidemic intensity during the 1918 influenza pandemic. Proc. Natl. Acad. Sci. U. S. A. 104, 7582. (doi:10.1073/pnas.0610941104).
    OpenUrlAbstract/FREE Full Text
  6. 6.
    Kirsch, T.D., Moseson, H., Massaquoi, M., Nyenswah, T.G., Goodermote, R., Rodriguez- Barraquer, I., Lessler, J., Cumings, D.A.T. & Peters, D.H. 2017 Impact of interventions and the incidence of ebola virus disease in Liberia—implications for future epidemics. Health Policy Plan. 32, 205–214. (doi:10.1093/heapol/czw113).
    OpenUrlCrossRefPubMed
  7. 7.↵
    Berkman, B.E. 2008 Mitigating Pandemic Influenza: The Ethics of Implementing a School Closure Policy. J. Public Health Manag. Pract. 14, 372–378. (doi:10.1097/01.PHH.0000324566.72533.0b).
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    Moore, S., Hill, E.M., Tildesley, M.J., Dyson, L. & Keeling, M.J. 2021 Vaccination and non- pharmaceutical interventions for COVID-19: a mathematical modelling study. Lancet Infect Dis 21, 793–802. (doi:10.1016/S1473-3099(21)00143-2).
    OpenUrlCrossRefPubMed
  9. 9.
    Spinelli, M.A., Glidden, D.V., Gennatas, E.D., Bielecki, M., Beyrer, C., Rutherford, G., Chambers, H., Goosby, E. & Gandhi, M. 2021 Importance of non-pharmaceutical interventions in lowering the viral inoculum to reduce susceptibility to infection by SARS-CoV-2 and potentially disease severity. Lancet Infect Dis 21, E296–E301. (doi:10.1016/S1473-3099(20)30982-8).
    OpenUrlCrossRef
  10. 10.↵
    Thompson, R.N., Hill, E.M. & Gog, J.R. 2021 SARS-CoV-2 incidence and vaccine escape. Lancet Infect Dis 21, 913–914. (doi:10.1016/S1473-3099(21)00202-4).
    OpenUrlCrossRefPubMed
  11. 11.↵
    Sachak-Patwa, R., Byrne, H.M., Dyson, L. & Thompson, R.N. 2021 The risk of SARS- CoV-2 outbreaks in low prevalence settings following the removal of travel restrictions. Communications Medicine 1, 39. (doi:10.1038/s43856-021-00038-8).
    OpenUrlCrossRef
  12. 12.↵
    Donohue, J.M. & Miller, E. 2020 COVID-19 and School Closures. JAMA 324, 845–847. (doi:10.1001/jama.2020.13092).
    OpenUrlCrossRef
  13. 13.
    Silverman, M., Sibbald, R. & Stranges, S. 2020 Ethics of COVID-19-related school closures. Can. J. Public Health 111, 462–465. (doi:10.17269/s41997-020-00396-1).
    OpenUrlCrossRef
  14. 14.
    Sadique, M.Z., Adams, E.J. & Edmunds, W.J. 2008 Estimating the costs of school closure for mitigating an influenza pandemic. BMC Public Health 8, 135. (doi:10.1186/1471-2458-8-135).
    OpenUrlCrossRefPubMed
  15. 15.
    Brooks, S.K., Webster, R.K., Smith, L.E., Woodland, L., Wessely, S., Greenberg, N. & Rubin, G.J. 2020 The psychological impact of quarantine and how to reduce it: rapid review of the evidence. The Lancet 395, 912–920. (doi:10.1016/S0140-6736(20)30460-8).
    OpenUrlCrossRefPubMed
  16. 16.↵
    Loades, M.E., Chatburn, E., Higson-Sweeney, N., Reynolds, S., Shafran, R., Brigden, A., Linney, C., McManus, M.N., Borwick, C. & Crawley, E. 2020 Rapid Systematic Review: The Impact of Social Isolation and Loneliness on the Mental Health of Children and Adolescents in the Context of COVID-19. J. Am. Acad. Child Adolesc. Psychiatry 59, 1218–1239.e1213. (doi:10.1016/j.jaac.2020.05.009).
    OpenUrlCrossRefPubMed
  17. 17.↵
    Prem, K., Cook, A.R. & Jit, M. 2017 Projecting social contact matrices in 152 countries using contact surveys and demographic data. PLoS Comp. Biol. 13, e1005697. (doi:10.1371/journal.pcbi.1005697).
    OpenUrlCrossRefPubMed
  18. 18.↵
    Mossong, J., Hens, N., Jit, M., Beutels, P., Auranen, K., Mikolajczyk, R., Massari, M., Salmaso, S., Tomba, G.S., Wallinga, J., et al. 2008 Social Contacts and Mixing Patterns Relevant to the Spread of Infectious Diseases. PLoS Med. 5, e74. (doi:10.1371/journal.pmed.0050074).
    OpenUrlCrossRefPubMed
  19. 19.↵
    Danon, L., Read, J.M., House, T.A., Vernon, M.C. & Keeling, M.J. 2013 Social encounter networks: characterizing Great Britain. Proc. R. Soc. Lond. B. Biol. Sci. 280, 20131037. (doi:10.1098/rspb.2013.1037).
    OpenUrlCrossRefPubMed
  20. 20.
    Leung, K., Jit, M., Lau, E.H.Y. & Wu, J.T. 2017 Social contact patterns relevant to the spread of respiratory infectious diseases in Hong Kong. Sci. Rep. 7. (doi:10.1038/s41598-017-08241-1).
    OpenUrlCrossRefPubMed
  21. 21.
    Béraud, G., Kazmercziak, S., Beutels, P., Levy-Bruhl, D., Lenne, X., Mielcarek, N., Yazdanpanah, Y., Boëlle, P.-Y., Hens, N. & Dervaux, B. 2015 The French Connection: The First Large Population-Based Contact Survey in France Relevant for the Spread of Infectious Diseases. PLoS One 10, e0133203. (doi:10.1371/journal.pone.0133203).
    OpenUrlCrossRefPubMed
  22. 22.
    Ibuka, Y., Ohkusa, Y., Sugawara, T., Chapman, G.B., Yamin, D., Atkins, K.E., Taniguchi, K., Okabe, N. & Galvani, A.P. 2016 Social contacts, vaccination decisions and influenza in Japan. J. Epidemiol. Community Health 70, 162–167. (doi:10.1136/jech-2015-205777).
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Read, J.M., Lessler, J., Riley, S., Wang, S., Tan, L.J., Kwok, K.O., Guan, Y., Jiang, C.Q. & Cummings, D.A.T. 2014 Social mixing patterns in rural and urban areas of southern China. Proc. R. Soc. Lond. B. Biol. Sci. 281, 20140268. (doi:10.1098/rspb.2014.0268).
    OpenUrlCrossRefPubMed
  24. 24.↵
    Read, J.M., Edmunds, W.J., Riley, S., Lessler, J. & Cummings, D.A.T. 2012 Close encounters of the infectious kind: methods to measure social mixing behaviour. Epidemiol. Infect. 140, 2117–2130. (doi:10.1017/S0950268812000842).
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.
    Wallinga, J., Teunis, P. & Kretzschmar, M. 2006 Using Data on Social Contacts to Estimate Age-specific Transmission Parameters for Respiratory-spread Infectious Agents. Am. J. Epidemiol. 164, 936–944. (doi:10.1093/aje/kwj317).
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.
    Kucharski, A.J., Kwok, K.O., Wei, V.W.I., Cowling, B.J., Read, J.M., Lessler, J., Cummings, D.A. & Riley, S. 2014 The Contribution of Social Behaviour to the Transmission of Influenza A in a Human Population. PLoS Pathog. 10, e1004206. (doi:10.1371/journal.ppat.1004206).
    OpenUrlCrossRefPubMed
  27. 27.↵
    Zhang, J., Litvinova, M., Liang, Y., Wang, Y., Wang, W., Zhao, S., Wu, Q., Merler, S., Viboud, C., Vespignani, A., et al. 2020 Changes in contact patterns shape the dynamics of the COVID-19 outbreak in China. Science 368, 1481–1486. (doi:10.1126/science.abb8001).
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    Munro, A.P. & Faust, S.N. 2020 COVID-19 in children: current evidence and key questions. Curr. Opin. Infect. Dis. 33, 540–547. (doi:10.1097/QCO.0000000000000690).
    OpenUrlCrossRef
  29. 29.
    Viner, R.M., Mytton, O.T., Bonell, C., Melendez-Torres, G.J., Ward, J., Hudson, L., Waddington, C., Thomas, J., Russell, S., van der Klis, F., et al. 2021 Susceptibility to SARS- CoV-2 Infection Among Children and Adolescents Compared With Adults: A Systematic Review and Meta-analysis. JAMA Pediatrics 175, 143–156. (doi:10.1001/jamapediatrics.2020.4573).
    OpenUrlCrossRefPubMed
  30. 30.
    Goldstein, E., Lipsitch, M. & Cevik, M. 2020 On the effect of age on the transmission of SARS-CoV-2 in households, schools and the community. J. Infect. Dis. 223, 362–369. (doi:10.1093/infdis/jiaa691).
    OpenUrlCrossRefPubMed
  31. 31.↵
    Davies, N.G., Klepac, P., Liu, Y., Prem, K., Jit, M., Pearson, C.A.B., Quilty, B.J., Kucharski, A.J., Gibbs, H., Clifford, S., et al. 2020 Age-dependent effects in the transmission and control of COVID-19 epidemics. Nat. Med. 26, 1205–1211. (doi:10.1038/s41591-020-0962-9).
    OpenUrlCrossRefPubMed
  32. 32.
    Ludvigsson, J.F. 2020 Systematic review of COVID-19 in children shows milder cases and a better prognosis than adults. Acta paediatrica (Oslo, Norway: 1992) 109, 1088-1095. (doi:10.1111/apa.15270).
    OpenUrlCrossRefPubMed
  33. 33.
    Dong, Y., Mo, X., Hu, Y., Qi, X., Jiang, F., Jiang, Z. & Tong, S. 2020 Epidemiology of COVID-19 Among Children in China. Pediatrics 145, e20200702. (doi:10.1542/peds.2020-0702).
    OpenUrlAbstract/FREE Full Text
  34. 34.
    Siebach, M.K., Piedimonte, G. & Ley, S.H. 2021 COVID-19 in childhood: Transmission, clinical presentation, complications and risk factors. Pediatr. Pulmonol. 56, 1342–1356. (doi:10.1002/ppul.25344).
    OpenUrlCrossRef
  35. 35.
    Mehta, N.S., Mytton, O.T., Mullins, E.W.S., Fowler, T.A., Falconer, C.L., Murphy, O.B., Langenberg, C., Jayatunga, W.J.P., Eddy, D.H. & Nguyen-Van-Tam, J.S. 2020 SARS-CoV-2 (COVID-19): What Do We Know About Children? A Systematic Review. Clin. Infect. Dis. 71, 2469–2479. (doi:10.1093/cid/ciaa556).
    OpenUrlCrossRefPubMed
  36. 36.↵
    Han, M.S., Choi, E.H., Chang, S.H., Jin, B.-L., Lee, E.J., Kim, B.N., Kim, M.K., Doo, K., Seo, J.-H., Kim, Y.-J., et al. 2021 Clinical Characteristics and Viral RNA Detection in Children With Coronavirus Disease 2019 in the Republic of Korea. JAMA Pediatrics 175, 73–80. (doi:10.1001/jamapediatrics.2020.3988).
    OpenUrlCrossRefPubMed
  37. 37.↵
    Pollock, A.M. & Lancaster, J. 2020 Asymptomatic transmission of covid-19. BMJ 371, m4851. (doi:10.1136/bmj.m4851).
    OpenUrlFREE Full Text
  38. 38.
    Buitrago-Garcia, D., Egli-Gany, D., Counotte, M.J., Hossmann, S., Imeri, H., Ipekci, A.M., Salanti, G. & Low, N. 2020 Occurrence and transmission potential of asymptomatic and presymptomatic SARS-CoV-2 infections: A living systematic review and meta-analysis. PLoS Med. 17, e1003346. (doi:10.1371/journal.pmed.1003346).
    OpenUrlCrossRefPubMed
  39. 39.
    Byambasuren, O., Cardona, M., Bell, K., Clark, J., McLaws, M.-L. & Glasziou, P. 2020 Estimating the extent of asymptomatic COVID-19 and its potential for community transmission: Systematic review and meta-analysis. JAMMI 5, 223–234. (doi:10.3138/jammi-2020-0030).
    OpenUrlCrossRef
  40. 40.
    Koh, W.C., Naing, L., Chaw, L., Rosledzana, M.A., Alikhan, M.F., Jamaludin, S.A., Amin, F., Omar, A., Shazli, A., Griffith, M., et al. 2020 What do we know about SARS-CoV-2 transmission? A systematic review and meta-analysis of the secondary attack rate and associated risk factors. PLoS One 15, e0240205. (doi:10.1371/journal.pone.0240205).
    OpenUrlCrossRefPubMed
  41. 41.
    Madewell, Z.J., Yang, Y., Longini, I.M., Jr.., Halloran, M.E. & Dean, N.E. 2020 Household Transmission of SARS-CoV-2: A Systematic Review and Meta-analysis. JAMA Network Open 3, e2031756–e2031756. (doi:10.1001/jamanetworkopen.2020.31756).
    OpenUrlCrossRef
  42. 42.↵
    Qiu, X., Nergiz, A.I., Maraolo, A.E., Bogoch, I.I., Low, N. & Cevik, M. 2021 Defining the role of asymptomatic and pre-symptomatic SARS-CoV-2 transmission - a living systematic review. Clin. Microbiol. Infect. 27, 511–519. (doi:10.1016/j.cmi.2021.01.011).
    OpenUrlCrossRefPubMed
  43. 43.↵
    Pijls, B.G., Jolani, S., Atherley, A., Derckx, R.T., Dijkstra, J.I.R., Franssen, G.H.L., Hendriks, S., Richters, A., Venemans-Jellema, A., Zalpuri, S., et al. 2021 Demographic risk factors for COVID-19 infection, severity, ICU admission and death: a meta-analysis of 59 studies. BMJ Open 11, e044640. (doi:10.1136/bmjopen-2020-044640).
    OpenUrlAbstract/FREE Full Text
  44. 44.
    Verity, R., Okell, L.C., Dorigatti, I., Winskill, P., Whittaker, C., Imai, N., Cuomo- Dannenburg, G., Thompson, H., Walker, P.G.T., Fu, H., et al. 2020 Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect Dis 20, 669–677. (doi:10.1016/S1473-3099(20)30243-7).
    OpenUrlCrossRefPubMed
  45. 45.↵
    Chen, T., Dai, Z., Mo, P., Li, X., Ma, Z., Song, S., Chen, X., Luo, M., Liang, K., Gao, S., et al. 2020 Clinical Characteristics and Outcomes of Older Patients with Coronavirus Disease 2019 (COVID-19) in Wuhan, China: A Single-Centered, Retrospective Study. J. Gerontol. A Biolo. Sci. Med. Sci. 75, 1788–1795. (doi:10.1093/gerona/glaa089).
    OpenUrlCrossRef
  46. 46.↵
    Prem, K., Liu, Y., Russell, T.W., Kucharski, A.J., Eggo, R.M., Davies, N., Flasche, S., Clifford, S., Pearson, C.A.B., Munday, J.D., et al. 2020 The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. Lancet Public Health 5, e261–e270. (doi:10.1016/S2468-2667(20)30073-6).
    OpenUrlCrossRefPubMed
  47. 47.↵
    Davies, N.G., Kucharski, A.J., Eggo, R.M., Gimma, A., Edmunds, W.J. & Centre for the Mathematical Modelling of Infectious Diseases COVID-19 working group. 2020 Effects of non- pharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: a modelling study. Lancet Public Health 5, e375–e385. (doi:10.1016/S2468-2667(20)30133-X).
    OpenUrlCrossRefPubMed
  48. 48.↵
    Liu, T., Gong, D., Xiao, J., Hu, J., He, G., Rong, Z. & Ma, W. 2020 Cluster infections play important roles in the rapid evolution of COVID-19 transmission: A systematic review. Int. J. Infect. Dis. 99, 374–380. (doi:10.1016/j.ijid.2020.07.073).
    OpenUrlCrossRef
  49. 49.
    Leclerc, Q., Fuller, N., Knight, L., CMMID COVID-19 Working Group, Funk, S. & Knight, G. 2020 What settings have been linked to SARS-CoV-2 transmission clusters? [version 2; peer review: 2 approved]. Wellcome Open Res 5. (doi:10.12688/wellcomeopenres.15889.2).
    OpenUrlCrossRefPubMed
  50. 50.↵
    Rossman, H., Keshet, A., Shilo, S., Gavrieli, A., Bauman, T., Cohen, O., Shelly, E., Balicer, R., Geiger, B., Dor, Y., et al. 2020 A framework for identifying regional outbreak and spread of COVID-19 from one-minute population-wide surveys. Nat. Med. 26, 634–638. (doi:10.1038/s41591-020-0857-9).
    OpenUrlCrossRef
  51. 51.↵
    Althaus, C.L., Low, N., Musa, E.O., Shuaib, F. & Gsteiger, S. 2015 Ebola virus disease outbreak in Nigeria: Transmission dynamics and rapid control. Epidemics 11, 80–84. (doi:10.1016/j.epidem.2015.03.001).
    OpenUrlCrossRefPubMed
  52. 52.
    Thompson, R.N., Jalava, K. & Obolski, U. 2019 Sustained transmission of Ebola in new locations: more likely than previously thought. Lancet Infect Dis 19, 1058–1059. (doi:10.1016/S1473-3099(19)30483-9).
    OpenUrlCrossRef
  53. 53.
    Thompson, R.N. 2020 Novel Coronavirus Outbreak in Wuhan, China, 2020: Intense Surveillance Is Vital for Preventing Sustained Transmission in New Locations. J. Clin. Med. 9, 498. (doi:10.3390/jcm9020498).
    OpenUrlCrossRef
  54. 54.
    Hellewell, J., Abbott, S., Gimma, A., Bosse, N.I., Jarvis, C.I., Russell, T.W., Munday, J.D., Kucharski, A.J., Edmunds, W.J., Sun, F., et al. 2020 Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. Lancet Glob. Health 8, e488–e496. (doi:10.1016/S2214-109X(20)30074-7).
    OpenUrlCrossRefPubMed
  55. 55.↵
    Lovell-Read, F.A., Funk, S., Obolski, U., Donnelly, C.A. & Thompson, R.N. 2021 Interventions targeting nonsymptomatic cases can be important to prevent local outbreaks: SARS-CoV-2 as a case-study. J. R. Soc. Interface 18, 20201014. (doi:10.1098/rsif.2020.1014).
    OpenUrlCrossRef
  56. 56.↵
    Nishiura, H., Cook, A.R. & Cowling, B.J. 2011 Assortativity and the Probability of Epidemic Extinction: A Case Study of Pandemic Influenza A (H1N1-2009). Interdiscip. Perspect. Infect. Dis. 2011, 194507. (doi:10.1155/2011/194507).
    OpenUrlCrossRefPubMed
  57. 57.↵
    Ferretti, L., Wymant, C., Kendall, M., Zhao, L., Nurtay, A., Abeler-Dorner, L., Parker, M., Bonsall, D. & Fraser, C. 2020 Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing. Science 368, eabb6936. (doi:10.1126/science.abb6936).
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    United Nations, Department of Economic and Social Affairs, Population Division. 2019 World Population Prospects 2019, Online Edition, Rev. 1. Available at: https://population.un.org/wpp/Download/Standard/Population/.
  59. 59.↵
    Wei, W.E., Li, Z., Chiew, C.J., Yong, S.E., Toh, M.P. & Lee, V.J. 2020 Presymptomatic Transmission of SARS-CoV-2 — Singapore, January 23–March 16, 2020. Morb. Mortal. Weekly Rep. 69, 411–415. (doi:10.15585/mmwr.mm6914e1externalicon).
    OpenUrlCrossRef
  60. 60.
    Arons, M.M., Hatfield, K.M., Reddy, S.C., Kimball, A., James, A., Jacobs, J.R., Taylor, J., Spicer, K., Bardossy, A.C., Oakley, L.P., et al. 2020 Presymptomatic SARS-CoV-2 Infections and Transmission in a Skilled Nursing Facility. N. Engl. J. Med. 382, 2081–2090. (doi:10.1056/NEJMoa2008457).
    OpenUrlCrossRefPubMed
  61. 61.
    Bullard, J., Dust, K., Funk, D., Strong, J.E., Alexander, D., Garnett, L., Boodman, C., Bello, A., Hedley, A., Schiffman, Z., et al. 2020 Predicting Infectious Severe Acute Respiratory Syndrome Coronavirus 2 From Diagnostic Samples. Clin. Infect. Dis. 71, 2663–2666. (doi:10.1093/cid/ciaa638).
    OpenUrlCrossRefPubMed
  62. 62.↵
    Wolfel, R., Corman, V.M., Guggemos, W., Seilmaier, M., Zange, S., Muller, M.A., Niemeyer, D., Jones, T.C., Vollmar, P., Rothe, C., et al. 2020 Virological assessment of hospitalized patients with COVID-2019. Nature 581, 465–469. (doi:10.1038/s41586-020-2196-x).
    OpenUrlCrossRefPubMed
  63. 63.↵
    Lai, C.-C., Shih, T.-P., Ko, W.-C., Tang, H.-J. & Hsueh, P.-R. 2020 Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges. Int. J. Antimicrob. Agents 55, 105924. (doi:10.1016/j.ijantimicag.2020.105924).
    OpenUrlCrossRefPubMed
  64. 64.
    Liu, Y., Gayle, A.A., Wilder-Smith, A. & Rocklov, J. 2020 The reproductive number of COVID-19 is higher compared to SARS coronavirus. J. Travel Med. 27, taaa021. (doi:10.1093/jtm/taaa021).
    OpenUrlCrossRefPubMed
  65. 65.
    Zhai, P., Ding, Y., Wu, X., Long, J., Zhong, Y. & Li, Y. 2020 The epidemiology, diagnosis and treatment of COVID-19. Int. J. Antimicrob. Agents 55, 105955–105955. (doi:10.1016/j.ijantimicag.2020.105955).
    OpenUrlCrossRefPubMed
  66. 66.↵
    Zhao, S., Lin, Q., Ran, J., Musa, S.S., Yang, G., Wang, W., Lou, Y., Gao, D., Yang, L., He, D., et al. 2020 Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak. Int. J. Infect. Dis. 92, 214–217. (doi:10.1016/j.ijid.2020.01.050).
    OpenUrlCrossRefPubMed
  67. 67.↵
    Thompson, R.N., Gilligan, C.A. & Cunniffe, N.J. 2020 Will an outbreak exceed available resources for control? Estimating the risk from invading pathogens using practical definitions of a severe epidemic. J. R. Soc. Interface 17, 20200690. (doi:10.1098/rsif.2020.0690).
    OpenUrlCrossRef
  68. 68.↵
    Norris, J.R. 1997 Markov Chains (Cambridge Series in Statistical and Probabilistic Mathematics). Cambridge, Cambridge University Press.
  69. 69.↵
    Hart, W.S., Maini, P.K. & Thompson, R.N. 2021 High infectiousness immediately before COVID-19 symptom onset highlights the importance of continued contact tracing. Elife 10. (doi:10.7554/eLife.65534).
    OpenUrlCrossRef
  70. 70.↵
    Huang, Q.S., Wood, T., Jelley, L., Jennings, T., Jefferies, S., Daniells, K., Nesdale, A., Dowell, T., Turner, N., Campbell-Stokes, P., et al. 2021 Impact of the COVID-19 nonpharmaceutical interventions on influenza and other respiratory viral infections in New Zealand. Nat. Commun. 12, 1001. (doi:10.1038/s41467-021-21157-9).
    OpenUrlCrossRef
  71. 71.
    Flaxman, S., Mishra, S., Gandy, A., Unwin, H.J.T., Mellan, T.A., Coupland, H., Whittaker, C., Zhu, H., Berah, T., Eaton, J.W., et al. 2020 Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature 584, 257–261. (doi:10.1038/s41586-020-2405-7).
    OpenUrlCrossRefPubMed
  72. 72.↵
    Borse, R.H., Behravesh, C.B., Dumanovsky, T., Zucker, J.R., Swerdlow, D., Edelson, P., Choe-Castillo, J. & Meltzer, M.I. 2011 Closing Schools in Response to the 2009 Pandemic Influenza A H1N1 Virus in New York City: Economic Impact on Households. Clin. Infect. Dis. 52, S168–S172. (doi:10.1093/cid/ciq033).
    OpenUrlCrossRefPubMed
  73. 73.↵
    Greer, A.L., Tuite, A. & Fisman, D.N. 2010 Age, influenza pandemics and disease dynamics. Epidemiol. Infect. 138, 1542–1549. (doi:10.1017/S0950268810000579).
    OpenUrlCrossRefPubMedWeb of Science
  74. 74.↵
    Clohisey, S. & Baillie, J.K. 2019 Host susceptibility to severe influenza A virus infection.Crit. Care 23, 303. (doi:10.1186/s13054-019-2566-7).
    OpenUrlCrossRefPubMed
  75. 75.↵
    Campbell, F., Archer, B., Laurenson-Schafer, H., Jinnai, Y., Konings, F., Batra, N., Pavlin, B., Vandemaele, K., Van Kerkhove, M.D., Jombart, T., et al. 2021 Increased transmissibility and global spread of SARS-CoV-2 variants of concern as at June 2021. Eurosurveillance 26, 2100509. (doi:10.2807/1560-7917.ES.2021.26.24.2100509).
    OpenUrlCrossRef
Back to top
Posted December 02, 2021.
Download PDF

Supplementary Material

Data/Code
Email
Estimating local outbreak risks and the effects of non-pharmaceutical interventions in age-structured populations: SARS-CoV-2 as a case study
Francesca A. Lovell-Read, Silvia Shen, Robin N. Thompson
medRxiv 2021.04.27.21256163; doi: https://doi.org/10.1101/2021.04.27.21256163
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area