Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Transmission dynamics reveal the impracticality of COVID-19 herd immunity strategies

View ORCID ProfileTobias S. Brett and View ORCID ProfilePejman Rohani
PNAS October 13, 2020 117 (41) 25897-25903; first published September 22, 2020; https://doi.org/10.1073/pnas.2008087117
Tobias S. Brett
aOdum School of Ecology, University of Georgia, Athens, GA, 30602;
bCenter for the Ecology of Infectious Diseases, University of Georgia, Athens, GA, 30602;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Tobias S. Brett
  • For correspondence: tsbrett@uga.edu
Pejman Rohani
aOdum School of Ecology, University of Georgia, Athens, GA, 30602;
bCenter for the Ecology of Infectious Diseases, University of Georgia, Athens, GA, 30602;
cDepartment of Infectious Diseases, University of Georgia, Athens, GA, 30602
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Pejman Rohani
  1. Edited by Alan Hastings, University of California, Davis, CA, and approved August 27, 2020 (received for review April 26, 2020)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Confronted with escalating COVID-19 outbreaks, countries at the leading edge of the pandemic have had to resort to imposing drastic social distancing measures which have serious societal and economic repercussions. Establishing herd immunity in a population by allowing the epidemic to spread, while mitigating the negative health impacts of COVID-19, presents a tantalizing resolution to the crisis. Our study simulating SARS-CoV-2 spread in the United Kingdom finds that achieving herd immunity without overwhelming hospital capacity leaves little room for error. Intervention levels must be carefully manipulated in an adaptive manner for an extended period, despite acute sensitivity to poorly quantified epidemiological factors. Such fine-tuning of social distancing renders this strategy impractical.

Abstract

The rapid growth rate of COVID-19 continues to threaten to overwhelm healthcare systems in multiple countries. In response, severely affected countries have had to impose a range of public health strategies achieved via nonpharmaceutical interventions. Broadly, these strategies have fallen into two categories: 1) “mitigation,” which aims to achieve herd immunity by allowing the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus to spread through the population while mitigating disease burden, and 2) “suppression,” aiming to drastically reduce SARS-CoV-2 transmission rates and halt endogenous transmission in the target population. Using an age-structured transmission model, parameterized to simulate SARS-CoV-2 transmission in the United Kingdom, we assessed the long-term prospects of success using both of these approaches. We simulated a range of different nonpharmaceutical intervention scenarios incorporating social distancing applied to differing age groups. Our modeling confirmed that suppression of SARS-CoV-2 transmission is possible with plausible levels of social distancing over a period of months, consistent with observed trends. Notably, our modeling did not support achieving herd immunity as a practical objective, requiring an unlikely balancing of multiple poorly defined forces. Specifically, we found that 1) social distancing must initially reduce the transmission rate to within a narrow range, 2) to compensate for susceptible depletion, the extent of social distancing must be adaptive over time in a precise yet unfeasible way, and 3) social distancing must be maintained for an extended period to ensure the healthcare system is not overwhelmed.

  • infectious diseases
  • mathematical modeling
  • dynamical systems

Caused by a novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (1), COVID-19 is an infectious disease capable of severe respiratory illness and death (2). Since its identification in Wuhan, China, COVID-19 has become an on-going and rapidly expanding global pandemic that is causing substantial mortality and healthcare system strain in multiple countries (3). While older individuals and those with underlying conditions are most at risk (4), infection has been seen across age groups (5, 6). Worryingly, detection of viral loads in the upper respiratory tract suggests potential for presymptomatic and ogliosymptomatic transmission (7⇓–9). Due to the absence of a vaccine, current attempts at controlling SARS-CoV-2 spread are focused on social measures that reduce rates of viral transmission: social distancing (a generalized reduction of contact rates between individuals in the population) and self-isolation by symptomatic individuals (10).

Broadly speaking, two distinct approaches to controlling the spread of SARS-CoV-2 have received much attention. The first aims to suppress transmission in the target population (referred to hereafter as “suppression”) (10). Under this objective, control measures reduce viral transmission to such a degree that sustained endogenous transmission is no longer possible. By maintaining control measures in place for a sufficient period of time, the virus will be eliminated in the focal population. The focus will then shift to preventing subsequent reintroduction to prevent resurgence. The second approach aims to manage or mitigate the negative health impacts (referred to hereafter as “mitigation”) (10). While suppression aims to ultimately halt local transmission, mitigation aims to reduce the growth rate of the epidemic to ensure disease burden does not overwhelm healthcare systems (3). By reducing (rather than halting) transmission, this strategy allows the susceptible pool to diminish, with the population potentially able to achieve herd immunity (whereby sustained local transmission is impossible, even without social distancing) (11). In practice, both approaches require the rollout of the same types of control measures (social distancing and self-isolation), although the necessary intensities and durations vary.

At the time of writing, trends in incidence data suggest multiple countries, including China, South Korea, Spain, and Italy, have successfully implemented suppression strategies (12). Sweden is the poster child for mitigation strategies, and appears to be aiming for herd immunity (13). Meanwhile, other countries continue to experience sustained transmission, and strategic intention is less clear [including the United States (14) and United Kingdom (15)]. The severe economic costs and acute societal pressures associated with social distancing measures have led to a push for their relaxation (10). Given the potentially long wait until a vaccine is available, the UK government appears to have considered following Sweden’s example and attempt to achieve herd immunity in the country (16).

The consequences of failure to either adequately mitigate or suppress COVID-19 are potentially catastrophic. Due to the many uncertainties surrounding SARS-CoV-2 transmission, evolution, and immunity, public health decision makers are presented with an unenviable task. To help inform control policies under uncertainty, mathematical modeling can assist in evaluating the viability of mitigation and suppression as objectives (17), by simulating the impacts of control strategies on viral transmission, hospital burden, fatalities, and population-level immunity.

Recent studies have modeled impacts of both mitigation and suppression strategies, including for China (18), low-income countries (19), and the United Kingdom (20, 21). Crucially, we are not aware of any systematic studies that focus on 1) walking the tightrope of achieving herd immunity without overburdening healthcare systems and 2) the control effort (e.g., reduction in contacts) required for successful mitigation relative to suppression. These two knowledge gaps motivate our study. We use an age-stratified disease transmission model, taking the United Kingdom as an example, to simulate SARS-CoV-2 spread controlled by individual self-isolation and widespread social distancing. We simulated various levels of self-isolation effectiveness and three distinct types of social-distancing measures: 1) school (including university) closures, 2) work and social place closures, and 3) effective isolation by older individuals.

Our modeling confirms that suppression is possible with plausible levels of social distancing and self-isolation, consistent with experience in multiple countries. Our research does not, however, support attempting to mitigate COVID-19 with the aim of building herd immunity. Achieving herd immunity while simultaneously maintaining hospital burden at manageable levels requires adaptive fine-tuning of mitigation efforts, in the face of imperfect epidemiological intelligence—something that is impractical.

Results

In the absence of any intervention measures, our modeling suggests SARS-CoV-2 will spread rapidly through the United Kingdom (Fig. 1A), and ultimately infect approximately 77% of the population (Fig. 1B). Using data on the age-specific fatality rate of COVID-19 (22), our results show around 350,000 fatalities among individuals aged over 60 y, and around 60,000 aged below 60 y (Fig. 1C). While we caution that our model makes a number of simplifying assumptions [e.g., no spatial dependence in transmission (11)], the total of fatalities is comparable to predictions made in other UK modeling studies (e.g., within the 95% prediction interval of ref. 20).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Simulated examples of SARS-CoV-2 spread in the United Kingdom using an age-structured Susceptible Exposed Infectious Recovered (SEIR) model under different control scenarios. (A) Daily new cases assuming no intervention measures enacted. (B) Cumulative proportion of the population exposed to SARS-CoV-2 over the course of the epidemic. (C) Cumulative fatalities assuming fixed age-specific case fatality rates (see Materials and Methods). (D–F) Same as A–C, but assuming that specific control measures are introduced when daily cases reached 10,000 (day 57): Older individuals social distance, and symptomatic individuals self-isolate (at 20% effectiveness). (G–I) Same as D–F but, in addition, schools close on day 57. Reopening of schools after 100 d results in a resurgence.

Sustained social distancing by older individuals (assumed to result in a 90% reduction in contacts with individuals under 25 y old, a 70% reduction with 25- to 59-y-olds, and a 50% reduction between one another), and moderately effective self-isolation by symptomatic individuals (at 20% efficacy) results in a shallower epidemic curve (Fig. 1D) and a much smaller outbreak size among individuals aged 60+ y (Fig. 1E). The attendant mortality burden among 60+-y-old individuals is also substantially reduced (to 62,000), with a smaller reduction in fatalities in those aged <60 y (to 43,000; Fig. 1F).

The addition of school (and university) closures, corresponding to a 70% reduction in contacts among school-aged individuals and a 20% reduction with 25- to 59-y-olds, dramatically reduces the rate of epidemic growth (Fig. 1G), although such levels of control are insufficient to suppress the epidemic (i.e., the number of daily cases still rises after implementation). The premature reopening of schools after 100 d (while the virus is still circulating) triggers a second wave of infection, with only a moderately reduced peak in daily new cases, largely eroding any additional gains made (23). The final proportion of the population exposed (Fig. 1H) and the number of fatalities (Fig. 1I) are largely unaltered compared to if schools had not been closed (compare Fig. 1 E and F).

Our modeling indicates that, if sustained, such control measures can lead to the suppression of COVID-19 in the United Kingdom by reducing R0 to <1 (Fig. 2 A and B). The effectiveness of self-isolation by symptomatic individuals (i.e., its impact on reducing overall transmission) is a product of two factors. Firstly, the proportion of infections generated while the primary case is exhibiting symptoms (including mild symptoms) and, secondly, the engagement of symptomatic individuals with self-isolation policies (see Materials and Methods). As both of these parameters decrease, the self-isolation efficacy drops, and greater social distancing measures are necessary to achieve suppression (Fig. 2A). At present, there is a large uncertainty in the relationship between symptoms and viral shedding (9). For the social distancing strengths considered, suppression is possible if over 14% of infections are caused while the primary case is exhibiting symptoms (ps). The associated self-isolation observance necessary to achieve suppression is inversely proportional to ps, decreasing from 100% if ps=0.14 to 50% if ps=0.28. Given the uncertainty surrounding asymptomatic transmission, the likelihood of successful suppression is greatest if all social distancing measures are enacted (Fig. 2B).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Prospects for disease suppression. Simulations of the age-structured SEIR model performed assuming control measures are initiated when there is a total of 10,000 infectious individuals in the population. The different control measures and strengths are listed in Materials and Methods. (A) Whether suppression is possible depends on both the self-isolation observance rate and the proportion of infections due to symptomatic individuals. More extensive social distancing measures increase the ranges of these two parameters for which suppression is possible. (B) Increases in self-isolation effectiveness drive down the reproductive number, which also depends on the social distancing measures employed. (C) The time taken for COVID-19 to be suppressed (modeled as a 100-fold reduction in infectious individuals) depends on the amount the reproductive number is decreased below one.

The time taken for suppression to be achieved (modeled as a 100-fold reduction in infectious individuals) once control measures are implemented is shown in Fig. 2C. If self-isolation effectiveness is high (>70% reduction in transmission), then suppression can be achieved in 2 mo regardless of any additional social distancing measures. There is little additional decrease in the necessary duration of social distancing unless schools and workplaces are both closed, in which case suppression can be achieved within 2 mo at much lower levels of self-isolation effectiveness (≳45%).

If suppression cannot be achieved (due to unfeasibility or lack of political will to reduce transmission sufficiently), then the objective of control measures is mitigation. Social distancing by 60+-y-old individuals results in a marked reduction of the final fraction of this age group that are exposed; however, unless both schools and workplaces are closed, additional social distancing measures do not lead to much further reduction (Fig. 3A).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Outcomes of disease mitigation attempts, simulated using the age-structured SEIR model. As in Fig. 2, control measures are implemented when there are 10,000 cases in the population. (A) The final fraction of 60+-y-old individuals exposed to COVID-19 for each of the control strategies simulated, assuming that social distancing measures can be maintained at the same strength indefinitely (“without fatigue”). (B) Size and (C) timing of the peak in daily new cases among 60+-y-old individuals. (D) Peak hospital burden, assuming age-specific hospitalization rates (see Materials and Methods) and a mean hospital stay of 7 d (24). (E–H) Simulation results using the same control strategies as in A–D, but assuming that, due to fatigue, schools and workplaces closures last 100 d.

These results are also mirrored in the impacts of social distancing on the daily cases in 60+-y-old individuals (Fig. 3B). Unfortunately, the hospital burden remains high for most intervention strategies, unless self-isolation is very effective (≳50%; Fig. 3D). Taking around 100,000 hospitalized cases to be the upper limit of hospital capacity, we find that there is a relatively small range of parameters where mitigation is successful at preventing hospitals being overwhelmed, but the disease is not also successfully suppressed (Fig. 3D; compare to Fig. 2B). If social distancing is applied to all age groups, this range is 0 to 14% self-isolation effectiveness, whereas, if just the 60+-y age group socially distance, the range is 41 to 54%.

As mentioned previously, if schools and workplaces reopen simultaneously after 100 d [e.g., due to social distancing “fatigue” (10)] and the disease has not been successfully suppressed, then much of the benefit of their closure is lost (Fig. 3 E–H) due to a resurgent second wave. In this scenario, the principle effect of school and workplace closures is in delaying the peak, buying more time for preparations (Fig. 3G).

It has been suggested that children might have reduced susceptibility to infection with SARS-CoV-2 (25). We repeated our analysis of suppression and mitigation assuming 50% reduction in susceptibility of individuals under 20 y old, with transmission rates among adults increased to ensure the same reproductive number, R0=2.3. As might be expected, there is a decreased impact of school closures; however, this is compensated by the increased impact of social distancing by over-20-y-olds, resulting in little net change in our findings if both work and school closures occur (SI Appendix, Figs. S1 and S2).

To summarize, aiming to build herd immunity to SARS-CoV-2 in a population while mitigating the burden on hospitals requires initially reducing the reproductive number to ensure available hospital capacity is not exceeded (Fig. 4A). If social distancing is maintained at a fixed level, hospital capacity needs to be much larger than presently available to achieve herd immunity without exceeding capacity; otherwise, the final outbreak size will be insufficient to achieve herd immunity (Fig. 4B). Relaxing social distancing measures linearly is also insufficient (SI Appendix, Fig. S3).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Summary of prospects for achieving herd immunity. (A) The peak hospital burden (shown on a log scale) is highly sensitive to the reproductive number. There is a narrow window of reproductive number values (shaded) where either 1) the number of COVID-19 cases requiring hospitalization does not overwhelm hospital capacity (modeled at the average hospital burden for April, 17,800 beds) or 2) circulation is suppressed. This window depends subtly on the exact age-specific social distancing configuration; however, all strategies studied fall between the two curves shown. (B) None of the simulated control scenarios shown in Figs. 2 and 3 achieved herd immunity while also keeping cases below hospital capacity. For the parameters considered, a hospital capacity in excess of 300,000 is required for this to be possible—almost 3 times the total UK NHS hospital beds (around 125,000 beds; see Materials and Methods), and around 15 times the average hospital burden of April.

Instead, achieving herd immunity without exceeding hospital capacity requires that social distancing is implemented nonlinearly. An idealized optimal strategy for achieving social distancing, calculated using a two-age group (<60 y old and ≥60 y old) SEIR model, is shown in Fig. 5 A and B. We assumed that individuals ≥60 y old were able to perfectly self-isolate, while herd immunity is built up in the <60-y age group. The epidemic in <60-y-old individuals is initially allowed to grow unimpeded until the hospital burden caused by the number of sick individuals reaches the hospital capacity. At this point, social distancing among those younger than 60 y commences (Fig. 5B), at a level that ensures the epidemic growth stalls (i.e., Reff=1). To achieve herd immunity in minimal time, the rate of new infections must then remain at the maximum afforded by the healthcare capacity. This requires reducing social distancing (i.e., increasing contacts) at the exact rate to balance out reductions in transmission stemming from the decreasing susceptible population (Fig. 5B)—too quickly and the epidemic grows above hospital capacity; too slowly and the epidemic with die out without reaching herd immunity. Over time, social distancing in <60-y-olds drops to zero, after which there are no longer enough susceptible <60-y-old individuals for the epidemic to be sustained. New infections continue but decline in frequency (meanwhile, further decreasing the effective reproductive number), before reaching zero.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Proof-of-principle strategy for achieving herd immunity. This strategy 1) achieves herd immunity with minimal social distancing duration 2) without exceeding hospital capacity, and 3) prevents infection among the most vulnerable. We simulated a reduced two-age group (<60 y old and ≥60 y old) model, assuming that all ≥60-y-old individuals completely isolate. (A) To achieve herd immunity in the minimum time requires control measures that fix the rate of new hospitalized cases to ensure hospital beds are at the maximum acceptable capacity until herd immunity is achieved (results showing taking capacity to be 1× and 2× the average hospital burden of April). (B) The reduction in contacts (via the product of self-isolation and social distancing) among <60-y-olds needs to vary nonlinearly. Until the hospital burden hits the hospital capacity, there is no social distancing; subsequently, it is tuned to ensure Reff=1; that is, the epidemic neither grows nor shrinks. Contact rates have to then gradually increase to exactly balance the reduction in Reff due to susceptible depletion. Eventually, there are no longer enough <60-y-old susceptible individuals to sustain the epidemic, and Reff<1. (C) While this strategy ensures that the final reproductive number drops below 1 if individuals ≥60 y old remain in isolation, if they return to pre-COVID-19 contact rates, then whether the reduction in susceptibles is sufficient to achieve herd immunity depends on the hospital capacity. The greater the proportion of the population that needs protecting, the greater the hospital capacity needed to achieve herd immunity. (D) The duration of social distancing also depends on the available hospital capacity. For the parameters considered, if hospital capacity is less than around 20,000, then it is not possible to achieve herd immunity without individuals aged ≥60 y remaining socially distanced.

While, at this point, the effective reproductive number is below 1 without any social distancing in the <60-y age group, full herd immunity has not necessarily been achieved (Fig. 5C). If ≥60-year-olds cease self-isolating, then the increased contact rates cause the effective reproductive number to increase, and may take it above 1. Achieving herd immunity requires that, while ≥60-y-olds are isolated, the additional new infections after Reff<1 sufficiently further reduce Reff (the “overshoot”) such that, when ≥60-y-olds cease self-isolation, it remains below 1. The amount of overshoot increases with the hospital capacity (Fig. 5C). In addition to increasing the prospects of achieving herd immunity (and also its robustness; see also SI Appendix, Fig. S3C), the necessary duration of social distancing is inversely proportional to hospital capacity (Fig. 5D).

Our modeling required estimates for the case hospitalization probability (22) and mean hospitalization stay (24). We explored sensitivity of the duration of social distancing to these two parameters, and also performed a comparison with a rough calculation of the duration based on estimates of seroprevalance at the end of April (26) (SI Appendix, Fig. S4). Our estimated necessary social distancing duration of 12 mo is close to the estimate from serology of 7 mo to 11 mo, especially considering the uncertainty in the case hospitalization probability (22).

Throughout our analysis, we used R0=2.3, broadly in line with most early estimates; however, some estimates were higher. We repeated our analysis using R0=5.7 (27), finding that, as might be expected, the range of intervention parameters that led to successful suppression was reduced (SI Appendix, Fig. S5). However, our main finding was reinforced: Higher values of R0 narrow the window for successful mitigation (SI Appendix, Fig. S6).

Discussion

Various governments have entertained the idea of achieving herd immunity through natural infection as a means of ending the long-term threat of COVID-19. This has provoked alarm in sections of the public health community (16, 28). Our work confirms that this alarm is well founded.

Attempting to achieve herd immunity while simultaneously mitigating the impact of COVID-19 on hospital burden is an extremely challenging task. In order to ensure the hospital burden does not exceed levels comparable with that of the United Kingdom in April 2020, R0 needs to be reduced from its initial value (assumed to be R0=2.3) to about 1.2. Suppression is possible if R0 is reduced below 1. Due to the fine margins (in terms of control effectiveness) between successful disease suppression and overwhelming hospitals, making herd immunity the primary objective (rather than applying maximal social distancing and aiming for suppression) is not supported by our modeling. Put another way, mitigation (via “flattening the curve”) is not a practical objective: If mitigation efforts are sufficient to prevent hospitals from being overwhelmed, only a comparatively small further increase in control measures will drive R0 below 1, and make suppression possible.

If herd immunity is the objective, then, in addition to the narrow range of R0 that must be aimed for, social distancing measures must be subsequently relaxed gradually in a highly controlled manner over a period of months to years. We were able to find a mathematical solution for building herd immunity with the minimal duration of social distancing using a reduced two-age group model. This solution requires knowledge of unobserved epidemiological determinants, namely, the remaining susceptible population, fraction exposed, and hospitalization probability. Additionally, complexities neglected by the model (e.g., failure of the ≥60-y age group to completely isolate, spatial structure) will alter the exact social distancing function. We, instead, view our results as serving as an indicator of the general shape and duration of social distancing necessary. Developing methods for implementable social distancing strategies capable of building toward herd immunity and that rely only on observable epidemiological data (e.g., the incidence curve and the hospital bed occupancy) requires further research.

The estimates of hospitalization probability and fatalities were calculated using results from a study of cases in Wuhan, China (22). For this study, we took the point estimates; however, these had uncertainties associated with them and are unlikely to be the same across regions. Furthermore, as treatment of COVID-19 continues to improve, fatality rates will fall. While there is an obvious feedback between fatality rates and healthcare system burden, the extent to which the Wuhan healthcare system (used in estimation) was overwhelmed is unclear. We therefore assumed the fatality rates fixed regardless of hospital burden. For these reasons, we have avoided attaching confidence intervals to estimates of fatalities, and they should be interpreted as plausible projections and not predictions.

A major unknown remains the nature, duration, and effectiveness of natural immunity. Here, we made the pragmatic assumption that, over the time scales under consideration, infection confers perfect long-lasting immunity (the best-case scenario for mitigation strategies). If immunity is not perfect, and there is a moderate to high chance of reinfection, then prospects for achieving herd immunity via natural infection are slim (29). To shed light on the kinetics of immunity, mass longitudinal antibody testing is necessary. This would both permit the identification of previously infected individuals, and provide information regarding immunity through time (30). We submit that models such as the one explored here, when integrated into statistical inference algorithms (31), provide a powerful means of integrating parallel serological and epidemiological data streams to quantify population-level immunity. Further, such models can be central to the development of efficient age-stratified serological testing schemes.

Finally, we stress that our study only explored the epidemiological impacts of nonpharmaceutical interventions (social distancing and self-isolation). Ultimately, any comprehensive public health policy needs to take into account the concomitant and wide-ranging societal and economic consequences of control measures.

Materials and Methods

Model.

We used a deterministic age-structured SEIR transmission model to simulate COVID-19 transmission in the United Kingdom. Contact rates ci,j, the number of daily contacts an individual of age i y makes with individuals of age j y, were taken from the POLYMOD study of social mixing patterns for the United Kingdom (32) corrected for reciprocity (33). The simulated age groups were matched to those of the POLYMOD study: 14 5-y increments from 0 y to 69 y and then 70+ y. Age-stratified population sizes (Nj) were taken from 2018 UK demographic data.

The mean latent and infectious periods were set to 1/ρ=3 and 1/γ=3 d, respectively, consistent with various estimates of the serial interval (34, 35) and incubation period (4, 36, 37), assuming infectiousness starts 1 d to 2 d before symptoms develop.

Both latent and infectious periods were assumed to be gamma distributed and modeled using the method of stages (11, 38, 39), by dividing the exposed (Ei) and infectious (Ii) compartments for each age group into four subcompartments, Ei=∑k=14Ei(k) and Ii=∑k=14Ii(k), where the superscript labels the subcompartment. The transmission dynamics for the age groups were governed by a system of ordinary differential equations,dSidt=−λi(t)Si,[1]dEi1dt=λi(t)Si−4ρEi(1),[2]dEikdt=4ρEi(k−1)−4ρEi(k)fork=2,3,4,[3]dIi1dt=4ρEi(4)−4γIi(1),[4]dIikdt=4γIi(k−1)−4γIi(k)fork=2,3,4,[5]λi(t)=β∑jci,j(t)Ij(t)Nj.[6]The transmission rate β was tuned using the next-generation matrix (40) to give a value of R0=2.3, consistent with estimates (34, 35). Simulations were initialized with one initial introduction in a fully susceptible population (Si=Ni). The resulting doubling time was observed to be about 3 d, broadly consistent with early observations from the United Kingdom.

Modeling Nonpharmaceutical Interventions.

Two types of nonpharmaceutical intervention were modeled: 1) self-isolation by symptomatic infectious individuals and 2) mass social distancing by differing age groups (Fig. 6). The effectiveness of self-isolation of symptomatic individuals is dependent on the product of two factors: 1) the proportion of infections that occur due to symptomatic individuals (excluding both presymptomatic and asymptomatic transmission), ps, and 2) the observance rate of social isolation among symptomatic individuals, k. The fractional reduction of contacts between age groups i and j due to social distancing is given by qi,j.

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Modeling the impact of nonpharmaceutical interventions on disease transmission. (A) Different social distancing measures (e.g., school closures, work and social place closures, and older individuals distancing) reduce contact rates between individuals in different ways. To reduce the complexity of the model, and to understand their differential impacts, we assume that individuals are only affected by one of these measures, dependent on their age. Individuals aged 0 y to 24 y are affected by school closures (included in this are university closures). School closures were assumed to result in a 70% reduction in contacts among school-aged individuals (qYY) and 20% reduction in their contacts with individuals aged 25 y to 59 y (qYA). Work and social place closures were assumed to reduce contacts among adults (qAA) by 50%. Finally, older individuals distancing reduced contacts by 60+-y-old individuals with 0- to 24-y-olds by 90% (qYO), with 25- to 59-y-olds by 70% (qAO), and among one another by 50% (qOO). The effectiveness of symptomatic individuals self-isolating is dependent on two factors: 1) the observance by symptomatic individuals, k, and 2) the proportion of transmission due to individuals who are symptomatic, ps. (B) We modeled five distinct combinations of social distancing measures, assuming that older individuals’ social distancing will always be prioritized.

Both of these interventions take the form of modifications to the contact matrix between infectious and susceptible individuals,c̃i,j=(1−kps)(1−qi,j)ci,j.[7]This expression for c̃i,j is inserted in place of ci,j in Eq. 6.

Age groups in the model are divided into whether they are young (Y; corresponding to 0- to 24-y-olds and age groups i=1 to 4), adults (A; 25- to 59-y-olds, age groups i=5 to 11) and older (O; 60+-y-olds, age groups i=12 to 15). The reduction in contacts due to social distancing, qi,j, is then determined by which of these three categories the contacter and contactee fall into, given by the block matrixq=qY YqY AqY OqY AqAAqAOqY OqAOqOO.[8]We assume school closures reduce contact rates between young individuals by a factor of qY Y=0.7 and between young people and adults by qY A=0.2. Social distancing among adults (e.g., due to workplace closures and reduction in social events) was modeled as a reduction of qAA=0.5. Social distancing of older individuals was represented by qY O=0.9, qAO=0.7, and qOO=0.5. For simulations with social distancing fatigue, qY Y, qY A, and qAA were modeled as linearly decreasing from these initial values to 0 over the periods indicated in SI Appendix, Fig. S3.

Estimating Hospital Burden and Case Fatalities.

Age-specific hospitalization probabilities, hi, and fatality rates were taken from point estimates calculated in a study of cases in Wuhan, China (22). Due to differences in the final age group of our model (70+ y) and those of the Wuhan study (70- to 79-y-olds and 80+-y-olds), the hospitalization and fatality rates for 70+-y-old individuals were calculated by summing the estimated 70- to 79-y-old and 80+-y-old rates weighted by their relative UK population sizes. Based on data from the United Kingdom, we assumed the average duration of hospitalization with COVID-19, d, was 7 d (24).

Pre–COVID-19 hospital beds and occupancy rates in the UK National Health Systems (NHS) were taken from the most recent (autumn 2019) published numbers for Northern Ireland, Wales, Scotland, and England.

Hospital burden was calculated as H=d∑ihiΛiSi. Although there is a lag between exposure and hospitalization, this has no bearing on our mathematical results.

Optimal Strategy to Achieve Herd Immunity.

We used a two-age group model (<60-y-olds and ≥60-y-olds) to demonstrate an optimal social distancing strategy while shielding ≥60-y-old individuals. To preserve the data-informed contact structure, we made the approximation that the age-specific incidence was constant for age groups <60 y and similarly for ≥60 y. The two-age group contact matrix is then given by a weighted sum of the POLYMOD-derived contact matrix, ca,b(2)=∑i∈a∑j∈b(Ni/Nb)ci,j, where a and b index age groups {<60,≥60}. Aside from the contact matrix, the model is identical in structure to that given in Eqs. 1–6.

If ≥60-y-old individuals are shielded (contact rates reduced to zero), the <60-y-old susceptible population is depleted at the fastest rate if the epidemic is allowed to grow unimpeded until the hospital burden (as defined above) reaches the hospital capacity, Hc. At this point, contact rates are reduced to maintain hospital burden at hospital capacity. Mathematically, this translates into setting the reduction in contacts, q(t), such that Eqs. 2–5 all equal 0. Solving for q(t) givesq(t)=1−1/Reffu(t)ifH≥HcandReffu(t)>1,0otherwise.[9]where Reffu(t)=γN0/(βc0,0(2)S0(t)) is the effective reproductive number without contact reduction by <60-y-olds. Substituting Eq. 9 into Eq. 1 and solving gives an expression for the control duration, td. For the parameters considered, there is limited susceptible depletion prior to hospital capacity being reached, and sotd≈dh0N0Hc1−γβc0,0(2).[10]

Data Availability.

Source code data have been deposited in Zenodo (10.5281/zenodo.400093) (41).

Acknowledgments

This work is funded by the NIH through Models of Infectious Disease Agent Study Grant 5R01GM123007.

Footnotes

  • ↵1To whom correspondence may be addressed. Email: tsbrett{at}uga.edu.
  • Author contributions: T.S.B. and P.R. designed research; T.S.B. performed research; T.S.B. analyzed data; and T.S.B. and P.R. wrote the paper.

  • The authors declare no competing interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.2008087117/-/DCSupplemental.

  • Copyright © 2020 the Author(s). Published by PNAS.

This open access article is distributed under Creative Commons Attribution License 4.0 (CC BY).

References

  1. ↵
    1. N. Zhu et al.
    , A novel coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 328, 727–733 (2020).
    OpenUrl
  2. ↵
    1. C. Huang et al.
    , Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 395, 497–506 (2020).
    OpenUrlCrossRefPubMed
  3. ↵
    1. J. Bedford et al.
    , COVID-19: Towards controlling of a pandemic. Lancet 395, 1015–1018 (2020).
    OpenUrlCrossRefPubMed
  4. ↵
    1. W.-J. Guan et al.
    , Clinical characteristics of coronavirus disease 2019 in China. N. Engl. J. Med. 382 1708–1720 (2020).
    OpenUrlCrossRefPubMed
  5. ↵
    1. W. Liu et al.
    , Detection of COVID-19 in children in early January 2020 in Wuhan, China. N. Engl. J. Med. 382, 1370–1371 (2020).
    OpenUrlCrossRefPubMed
  6. ↵
    1. Q. Bi et al.
    , Epidemiology and transmission of COVID-19 in Shenzhen China: Analysis of 391 cases and 1,286 of their close contacts. medRxiv:10.1101/2020.03.03.20028423 (27 March 2020).
  7. ↵
    1. L. Zou et al.
    , SARS-CoV-2 viral load in upper respiratory specimens of infected patients. N. Engl. J. Med. 382, 1177–1179 (2020).
    OpenUrlCrossRefPubMed
  8. ↵
    1. W. E. Wei et al.
    , Presymptomatic transmission of SARS-CoV-2—Singapore, January 23–March 16, 2020. MMWR (Morb. Mortal. Wkly. Rep.) 69, 411–415 (2020).
    OpenUrlCrossRefPubMed
  9. ↵
    1. R. Wölfel et al.
    , Virological assessment of hospitalized patients with COVID-2019. Nature 581, 465–469 (2020).
    OpenUrlCrossRefPubMed
  10. ↵
    1. R. M. Anderson,
    2. H. Heesterbeek,
    3. D. Klinkenberg,
    4. T. D. Hollingsworth
    , How will country-based mitigation measures influence the course of the COVID-19 epidemic?. Lancet 395, 931–934 (2020).
    OpenUrlCrossRefPubMed
  11. ↵
    1. M. J. Keeling,
    2. P. Rohani
    , Modeling Infectious Diseases in Humans and Animals (Princeton University Press, 2008).
  12. ↵
    1. E. Dong,
    2. H. Du,
    3. L. Gardner
    , An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect. Dis. 20, 533–534 (2020).
    OpenUrlCrossRefPubMed
  13. ↵
    1. H. Habib
    , Has Sweden’s controversial COVID-19 strategy been successful? BMJ 369, m2376 (2020).
    OpenUrlFREE Full Text
  14. ↵
    1. R. L. Haffajee,
    2. M. M. Mello
    , Thinking globally, acting locally—The U.S. response to COVID-19. N. Engl. J. Med. 382, e75 (2020).
    OpenUrl
  15. ↵
    1. G. Scally,
    2. B. Jacobson,
    3. K. Abbasi
    , The UK’s public health response to COVID-19. BMJ 369, m1932 (2020).
    OpenUrlFREE Full Text
  16. ↵
    1. D. J. Hunter
    , COVID-19 and the stiff upper lip—The pandemic response in the United Kingdom. N. Engl. J. Med. 382, e31 (2020).
    OpenUrlPubMed
  17. ↵
    1. R. M. Anderson,
    2. B. Anderson,
    3. R. M. May
    , Infectious Diseases of Humans: Dynamics and Control (Oxford University Press, 1992).
  18. ↵
    1. K. Prem et al.
    , 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 (2020).
    OpenUrl
  19. ↵
    1. P. G. T. Walker et al.
    , The impact of COVID-19 and strategies for mitigation and suppression in low- and middle-income countries. Science 369, 413–422 (2020).
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. N. G. Davies et al.
    , 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 (2020).
    OpenUrl
  21. ↵
    1. N. Ferguson et al.
    , “Impact of non-pharmaceutical interventions (NPIs) to reduce COVID19 mortality and healthcare demand” (Rep. 9, Imperial College London, 2020).
  22. ↵
    1. R. Verity et al.
    , Estimates of the severity of coronavirus disease 2019: A model-based analysis. Lancet Infect. Dis. 20, P669–P677 (2020).
    OpenUrl
  23. ↵
    1. T. D. Hollingsworth,
    2. D. Klinkenberg,
    3. H. Heesterbeek,
    4. R. M. Anderson
    , Mitigation strategies for pandemic influenza A: Balancing conflicting policy objectives. PLoS Comput. Biol. 7, e1001076 (2011).
    OpenUrlCrossRefPubMed
  24. ↵
    1. A. B. Docherty et al.
    , Features of 16,749 hospitalised UK patients with COVID-19 using the ISARIC WHO clinical characterisation protocol. medRxiv:10.1101/2020.04.23.20076042 (28 April 2020).
  25. ↵
    1. S. Stringhini et al.
    , Seroprevalence of anti-SARS-CoV-2 IgG antibodies in Geneva, Switzerland (SEROCoV-POP): A population-based study. Lancet 396, P313–P319 (2020).
    OpenUrl
  26. ↵
    1. Office for National Statistics
    , Coronavirus (COVID-19) infection survey pilot: 28 May 2020. https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/bulletins/coronaviruscovid19infectionsurveypilot/28may2020. Accessed 3 June 2020.
  27. ↵
    1. S. Sanche et al.
    , High contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2. Emerg. Infect. Dis. 26, 1470–1477 (2020).
    OpenUrlPubMed
  28. ↵
    1. R. Horton
    , Offline: COVID-19—A reckoning. Lancet 395, 935 (2020).
    OpenUrlPubMed
  29. ↵
    1. S. M. Kissler,
    2. C. Tedijanto,
    3. E. Goldstein,
    4. Y. H. Grad,
    5. M. Lipsitch
    , Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science 368, 860–868 (2020).
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. A. Antia et al.
    , Heterogeneity and longevity of antibody memory to viruses and vaccines. PLoS Biol. 16, e2006601 (2018).
    OpenUrlCrossRef
  31. ↵
    1. M. D. de Cellès,
    2. F. M. Magpantay,
    3. A. A. King,
    4. P. Rohani
    , The impact of past vaccination coverage and immunity on pertussis resurgence. Sci. Transl. Med. 10, eaaj1748 (2018).
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. J. Mossong et al.
    , Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 5, e74 (2008).
    OpenUrlCrossRefPubMed
  33. ↵
    1. M. A. Riolo,
    2. P. Rohani
    , Combating pertussis resurgence: One booster vaccination schedule does not fit all. Proc. Natl. Acad. Sci. U.S.A. 112, E472–E477 (2015).
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Q. Li et al.
    , Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. N. Engl. J. Med. 382, 1199–1207 (2020).
    OpenUrlCrossRefPubMed
  35. ↵
    1. J. Zhang et al.
    , Evolving epidemiology and transmission dynamics of coronavirus disease 2019 outside Hubei province, China: A descriptive and modelling study. Lancet Infect. Dis. 20, 793–802 (2020).
    OpenUrlPubMed
  36. ↵
    1. N. M. Linton et al.
    , Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: A statistical analysis of publicly available case data. J. Clin. Med. 9, 538 (2020).
    OpenUrl
  37. ↵
    1. S. A. Lauer et al.
    , The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: Estimation and application. Ann. Intern. Med. 172, 577–582 (2020).
    OpenUrlCrossRefPubMed
  38. ↵
    1. A. L. Lloyd
    , Realistic distributions of infectious periods in epidemic models: Changing patterns of persistence and dynamics. Theor. Popul. Biol. 60, 59–71 (2001).
    OpenUrlCrossRefPubMed
  39. ↵
    1. H. Wearing,
    2. P. Rohani,
    3. M. Keeling
    , Appropriate models for the management of infectious diseases. PLoS Med. 2, e174 (2005).
    OpenUrlCrossRefPubMed
  40. ↵
    1. P. Van den Driessche,
    2. J. Watmough
    , Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 180, 29–48 (2002).
    OpenUrlCrossRefPubMed
  41. ↵
    1. T. Brett
    , Data for: Transmission dynamics reveal the impracticality of COVID-19 herd immunity strategies (Version 2.0.0). Zenodo. http://doi.org/10.5281/zenodo.4000983. Deposited 25 August 2020.
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Transmission dynamics reveal the impracticality of COVID-19 herd immunity strategies
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Transmission dynamics reveal the impracticality of COVID-19 herd immunity strategies
Tobias S. Brett, Pejman Rohani
Proceedings of the National Academy of Sciences Oct 2020, 117 (41) 25897-25903; DOI: 10.1073/pnas.2008087117

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Transmission dynamics reveal the impracticality of COVID-19 herd immunity strategies
Tobias S. Brett, Pejman Rohani
Proceedings of the National Academy of Sciences Oct 2020, 117 (41) 25897-25903; DOI: 10.1073/pnas.2008087117
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 117 (41)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Biological Sciences
  • Population Biology

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Data Availability.
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Surgeons hands during surgery
Inner Workings: Advances in infectious disease treatment promise to expand the pool of donor organs
Despite myriad challenges, clinicians see room for progress.
Image credit: Shutterstock/David Tadevosian.
Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Double helix
Journal Club: Noncoding DNA shown to underlie function, cause limb malformations
Using CRISPR, researchers showed that a region some used to label “junk DNA” has a major role in a rare genetic disorder.
Image credit: Nathan Devery.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Multi-color molecular model
Enzymatic breakdown of PET plastic
A study demonstrates how two enzymes—MHETase and PETase—work synergistically to depolymerize the plastic pollutant PET.
Image credit: Aaron McGeehan (artist).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490