Dynamics of SARS-CoV-2 with Waning Immunity in the UK Population ================================================================ * Thomas Crellen * Li Pi * Emma L. Davis * Timothy M. Pollington * Tim C. D. Lucas * Diepreye Ayabina * Anna Borlase * Jaspreet Toor * Kiesha Prem * Graham F. Medley * Petra Klepac * T. Déirdre Hollingsworth ## Abstract The dynamics of immunity are crucial to understanding the long-term patterns of the SARS-CoV-2 pandemic. Several cases of reinfection with SARS-CoV-2 have been documented 48–142 days after the initial infection and immunity to seasonal circulating coronaviruses is estimated to be shorter than one year. Using an age-structured, deterministic model, we explore potential immunity dynamics using contact data from the UK population. In the scenario where immunity to SARS-CoV-2 lasts an average of three months for non-hospitalised individuals, a year for hospitalised individuals, and the effective reproduction number after lockdown ends is 1.2 (our worst case scenario), we find that the secondary peak occurs in winter 2020 with a daily maximum of 387,000 infectious individuals and 125,000 daily new cases; three-fold greater than in a scenario with permanent immunity. Our models suggests that longitudinal serological surveys to determine if immunity in the population is waning will be most informative when sampling takes place from the end of the lockdown in June until autumn 2020. After this period, the proportion of the population with antibodies to SARS-CoV-2 is expected to increase due to the secondary wave. Overall, our analysis presents considerations for policy makers on the longer term dynamics of SARS-CoV-2 in the UK and suggests that strategies designed to achieve herd immunity may lead to repeated waves of infection as immunity to reinfection is not permanent. ## Introduction ### Public health context As of the 1st of July 2020, SARS-CoV-2 has infected at least 10 million people worldwide and resulted in over 500,000 deaths [1, 2]. Following the initial outbreak from a live animal market in Wuhan, China [3], the United Kingdom (UK) has been among the countries most severely affected; reporting over 310,000 cases and 44,000 deaths, which is among the highest per-capita rates [2, 4]. Since the 23rd of March, nationwide non-pharmaceutical interventions (lockdown) have been in place to reduce social contacts by closing schools and shops; encouraging home working; and social distancing in public places. Similar measures have been in place in other European countries since late February 2020 with restrictions easing in France, Germany and Italy from May 2020. Within the European picture of disease control strategies, Sweden has been an outlier by placing fewer restrictions on social mixing while aiming to build up immunity in the population [5]. ### Immune response to SARS-CoV-2 Following infection with the virus, hospitalised patients have an acute immune response where virus-specific IgM and IgG antibodies titres reach a maximum 15–21 and 22–27 days respectively after symptom onset [6, 7]. Antibodies raised in hospitalised patients and animal models against SARS-CoV-2 provide protection for at least several weeks following infection [8, 9], suggesting that immediate reinfection with the virus is unlikely. There is limited evidence that hospitalised patients with more severe symptoms show a greater antibody response [6, 9]. Asymptomatic individuals have a weaker IgG and specific antibody response to SARS-CoV-2 and are more likely to become seronegative following convalescence [10]. Antibody titres raised against related coronaviruses SARS-CoV and MERS-CoV have been shown to decay over time [11, 12]. Furthermore, immunity to seasonal circulating coronaviruses has been estimated to last for less than one year [13] and recovered individuals from coronavirus NL63 can become reinfected [14]. Concerns that immunity to SARS-CoV-2 may also wane therefore motivated the present study [15]. ### Epidemiological modelling Dynamic epidemiological models play a major role in shaping the timing and intensity of interventions against SARS-CoV-2 in the UK and elsewhere [16]. Many models or simulations have assumed that infected individuals recover with permanent immunity [16, 17, 18]. In such models the epidemic reaches extinction after running out of infected individuals, although they do not preclude a second wave of infections after lockdown [19]. If immunity wanes over a period of time, or recovered individuals have only partial immunity to reinfection, this substantially alters the dynamics of the system [20]. In the absence of stochastic extinction and demography (births and deaths) in a population with equal mixing where; *R* is the basic reproduction number; *γ* is the average duration of infection; and *ω* is the reciprocal of the average duration of immunity; the endemic equilibrium proportion of infected in the population *I****, is given by (*R* − 1) *ω/γR* and thus, in the absence of interventions, the infection persists indefinitely when *R* *>* 1 [21]. In dynamic models which make the assumption of homogeneous mixing, the ‘classic’ threshold at which the proportion of the population with immunity prevents invasion by a pathogen is given by 1 1*/R* [21]. As *R* for SARS-CoV-2 is generally estimated between 2.4–5.6 [22, 23, 24], this equates to 58–82% of the population requiring immunity to eventually halt the epidemic. Serological studies conducted in affected countries to-date have reported the proportion of the population with antibodies against SARS-CoV-2 to be much lower than this figure [22, 25]. However, when more realistic non-homogeneous mixing is considered, the observed immunity threshold is lower than the classical threshold [26]. Recent studies have considered this question for SARS-CoV-2 [27, 28], with Britton *et al*. noting that the disease-induced population immunity threshold could be closer to 40% in an age-structured population when *R* is 2.5, rather than the 60% ‘classic’ population immunity threshold [28]. This phenomenon is driven by individuals that have more contacts, or greater susceptibility to the virus, getting infected earlier on and leaving the susceptible population; thus decelerating the growth of the epidemic. ### Models with waning immunity Kissler *et al*. considered the dynamics of SARS-CoV-2 in the United States with seasonal forcing, homogeneous mixing and waning immunity that could be boosted by exposure to seasonal circulating betacoronaviruses [29]. Under these assumptions, the incidence of SARS-CoV-2 was predicted to rebound in winter months. Here we do not consider seasonality, but rather the dynamics of transmission in an age-structured population with different periods of waning immunity in the context of the UK emerging from lockdown. We developed a discrete time gamma delay-distributed (susceptible-exposed-infectious-recovered-susceptible; SEIRS) model, which incorporates current knowledge about the natural history of the virus and the UK population. Our model accounts for symptomatic and asymptomatic transmission, and heterogeneity in both daily contacts and infection susceptibility by age group. We consider different durations of immunity for hospitalised patients (or those with more severe symptoms) compared to non-hospitalised patients (those with less severe symptoms). We use this model to explore a range of scenarios in the UK population in the context of stringent non-pharmaceutical interventions (lockdown) followed by more limited interventions over a two year period from March 2020, and the impact of immunity duration on the longer term disease equilibrium. ### Reinfection with SARS-CoV-2 After this analysis was presented to policymakers (see Accompanying Letter), the first case of reinfection with SARS-CoV-2 was reported in a 33 year old male in Hong Kong in August 2020; 142 days after the initial diagnosis [30]. Following the first infection in March, the patient was discharged from hospital after two negative PCR results. Genomic sequencing confirmed that the viruses in March and August were from distinct clades. Reinfections have since been reported in single cases in the USA [31], Ecuador [32] and Belgium [33] at 48, 63 and 93 days respectively from the first reported infections. Two asymptomatic reinfections were also reported in Indian healthcare workers 108 and 111 days from the initial infection [34]. The case of reinfection in the USA, in a 25 year old male from Nevada, was notable for the second infection causing more severe symptoms than the first, and led to the patient requiring supplemental oxygen and hospitalisation. While these reported cases confirm the biological plausibility of reinfection with SARS-CoV-2, the data are too sparse to infer an expected duration of immunity at the population level. ## Methods ### Model structure We use current knowledge of the natural history of the virus to construct a plausible epidemiological model (Figure 1). We extend a previously published deterministic compartmental model which has provided general insights into the dynamics of the epidemic at a regional level for a range of scenarios [18]. The general framework of the model is given in Figure 1 and parameter values are shown in Table 1. The disease states for each age group are susceptible (*S*), exposed (*E*), symptomatic infectious (*I**S*), asymptomatic infectious (*I**A*), hospitalised recovered (*R**H*), and non-hospitalised recovered (*R**N*). Each state gives the number of individuals in the compartment at a given time *t*, with the sum over all the states at time *t* equal to the size of the UK population in 2018 (66.4 million). View this table: [Table 1:](http://medrxiv.org/content/early/2021/02/04/2020.07.24.20157982/T1) Table 1: Summary of parameter values used in the modelled scenarios of SARS-CoV-2 transmission in the United Kingdom. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/04/2020.07.24.20157982/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/02/04/2020.07.24.20157982/F1) Figure 1: Flow diagram showing SARS-CoV-2 transmission model outline. Within age class *i*, the disease states are susceptible (*S*), exposed (*E*), symptomatic infectious (*I**S*), asymptomatic infectious (*I**A*), hospitalised recovered (*R**H*), and non-hospitalised recovered (*R**N*). The rate parameters *λ* represent the force of infection acting on an individual in age class *i* (Equation 8), *σ* represents the mean latency period, *γ* the mean infectious period, *ω**N* the mean duration of immunity for non-hospitalised recovered individuals, and *ω**H* the mean duration of immunity for hospitalised recovered individuals. The probabilities *φ**i* give the proportion of individuals in age group *i* that are asymptomatic following infection and *p**i* the proportion in age group *i* that require hospitalisation (or have severe symptoms) following infection. All model parameters are given in Table 1 and state transitions are shown in Equations 1 to 6. ### Distributed natural history of infection The mean latent and infectious periods for SARS-CoV-2 have been estimated as 4.5 days and 3.1 days respectively, using viral load data and the timing of known index and secondary case contacts (Figure 2) [35]. The probability mass of the latent and infectious period distributions are centred around the mean, therefore we consider that gamma distributions with an integer shape parameter, also known as Erlang distributions, give more realistic waiting times than exponential distributions which have a mode of zero [36, 37, 38]. We decided to use a discrete time implementation (difference equations), as this permits greater flexibility in implementing interventions compared to differential equations. Using difference equations results in a discretised approximation of the Erlang distribution compared with continuous time models, with the output converging as the time step approaches zero [39, 40]. We used a time step of 0.01 days to ensure that the infection dynamics closely approximate a continuous time model. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/04/2020.07.24.20157982/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/02/04/2020.07.24.20157982/F2) Figure 2: Probabilities for time spent in each state given Erlang (gamma) distributed waiting times. A. Proportion of individuals in exposed and infectious classes since time from infection. Time exposed and time infectious have mean durations of 4.5 and 3.1 days respectively. B. Proportion of individuals immune since recovery, where time immune has mean durations of 90, 180 or 365 days depending on the scenario. ### Transmissibility and infectivity Estimates of the transmissibility of the virus in the UK at the beginning of the epidemic have ranged from 2.4–5.6 [23, 41, 42], therefore we assume that *R* at the beginning of the epidemic in the UK population is 4.0 [23]. Non-pharmaceutical interventions have been shown to bring the effective reproduction number (*R**t*) below one, and in some settings have led to local elimination of the virus [22, 23]. Testing performed in closed populations suggests that 40-50% of SARS-CoV-2 infections may be asymptomatic [43, 44, 45], while data from contact tracing shows transmission can occur from asymptomatic individuals [46]. We make the assumption that asymptomatic individuals (*I**A*) have 0.5 the infectiousness of symptomatic individuals (*I**S*) [6, 16]. The UK population shows variable contact rates by age [47, 48] and, while studies show mixed results, evidence is accruing that children have a lower susceptibility to acquiring the infection than adults [49, 50, 51]. We assume that children (*≤*15 years) have 0.4 times the susceptibility of adults [52]. ### Immunity following infection We allow the duration of immunity to differ for recovered individuals with severe symptoms that are hospitalised (*R**H*) versus those with less severe symptoms that are not hospitalised (*R**N*), as there is evidence from SARS-CoV-2 and other coronaviruses that individuals with milder symptoms may have a lower antibody response [53]. Among the reported cases of reinfections, all individuals showed at most ‘mild’ symptoms during the first infection [54], supporting the notion that infections which cause more severe symptoms could generate longer lasting immunity. Epidemic transitions for age group *i* at time *t* + 1 are given by: ![Formula][1] ![Formula][2] ![Formula][3] ![Formula][4] ![Formula][5] ![Formula][6] The function notation *f**t*(*X*; *k, µ*) represents the output of the discretised Erlang delay distribution from states *E, I**S*, *I**A*, *R**H* and *R**N*. Our model code contains *k* sub-compartments for each Erlang distributed disease state, with a rate *k/µ* between sub-compartments, and an overall waiting time of *µ*. The function notation is used in equations 1–6 for shorthand convenience instead of writing out *k* difference equations for each disease state. While we do not calculate probability distributions within the model, as a simplified example; if *n* individuals enter state *X* at a single time point *t*, by time *t* + *τ* there will be remaining *n*(1 − *g*(*τ, k, µ*)); where *g*(*τ, k, µ*) gives the cumulative Erlang distribution with shape parameter *k* and mean duration *µ* [39]. The next generation matrix (*K* = *k**ij*) gives the expected number of secondary infections in age group *i* resulting from contact with an index case in age group *j*: ![Formula][7] where *η**i* is the relative susceptibility of age group *i*; *υ* the relative infectiousness of asymptomatic cases; and *c**i,j* is the average number of daily contacts between a single individual in age group *j* and all individuals in age group *i* [48] (Table 1). The basic reproduction number (*R*) is given by the spectral radius *ρ*(*K*) which is the largest absolute eigenvalue of *K*. As we specify the value of *R*, the transmission parameter *β* is left as a free parameter which is scaled to the correct value. The force of infection acting on a single individual in age group *i* at time *t* + 1 (*λ**i,t*+1) is given by: ![Formula][8] where *M* is the number of discrete age groups (*M* = 15); and *N**i* gives the population size of age group *i*. ### Immunological scenarios Using data and timing of events from the UK epidemic, we explore four scenarios with varying average durations of immunity to SARS-CoV-2 (Figure 2). S1. **Permanent:** Where immunity is lifelong for both hospitalised (*R**H*) and non-hospitalised (*R**N*) cases. S2. **Waning (12 months):** Where immunity is lifelong for hospitalised cases and has an average duration of 365 days for non-hospitalised cases. S3. **Waning (6 months):** Where immunity is lifelong for hospitalised cases and has an average duration of 180 days for non-hospitalised cases. S4. **Short-lived:** Where immunity lasts, on average, 365 days for hospitalised cases and 90 days for non-hospitalised cases. ### UK-specific parameterisation All scenarios are initialised with 2968 infected individuals in early March 2020 [59], which is when the majority of importation events into the UK are likely to have occurred [59]. Intervention measures are initiated on the 23rd of March (date the UK nationwide lockdown started), with an immediate reduction in the effective reproduction number (expected number of secondary cases from an index case at time *t*; *R**t*) to 1.2 for a three week period, followed by a further reduction in *R**t* to 0.8 until lockdown measures are eased on 15th June [58]. After this time, *R**t* is brought to 0.9, 1, 1.1 or 1.2 until March 2022. We considered the majority of our analysis over a, relatively short, two year period to explore the epidemic up to a secondary peak; beyond this point the dynamics are likely to be altered depending on further interventions or changes to *R**t*. As we simulate disease dynamics over a relatively short period of time, we do not consider demography (births and deaths) or transitions between age classes (ageing). To obtain equilibrium values, we simulated epidemic trajectories for up to five years. The UK contact matrix (average daily contacts between a single individual in age group *j* with individuals in age group *i*) comes from a ‘citizen science’ project for the BBC, in which individuals in the UK population provided detailed information on their daily contacts in the home, in the workplace, at school and in other settings [47, 48]. The contact matrix is altered to account for changes to contact patterns during and after the main intervention period [56]. During the lockdown, home; work; school; and other contacts are reduced to 0.8, 0.3, 0.1 and 0.2 respectively of their baseline values. This reflects the school closures for all children, except for those of key workers, and that workers were encouraged to work from home. Reduction in home contacts accounts for the absence of visitors to the home during the lockdown. In the post-lockdown phase, home; work; school; and other contacts are scaled to 1, 0.8, 0.85 and 0.75, respectively, of their baseline values to reflect limited social distancing measures that are likely to be in place until at least the end of 2021. Analysis was performed in R version 4.0.2. We present figures from model output in the text to the nearest thousand. Code to reproduce the analysis with user-specified parameters is available at [https://github.com/tc13/covid-19-immunity](https://github.com/tc13/covid-19-immunity). ## Results ### Heterogeneity in transmission The epidemic is driven by the rate of infectious contacts between individuals in different age groups. This is described by the next generation matrix in which the average number of secondary cases generated by an index case in age group *j* is the summation of row *j* (Equation 8 & Figure 3). At the beginning of the epidemic, when SARS-CoV-2 is spreading rapidly, all age groups are involved in transmission; in particular those aged 20–39 years. An index case in the 20–24 age group, for instance, is expected to generate an average of 4.3 secondary cases at baseline. As lockdown measures come into force this dramatically reduces the expected number of secondary cases due to fewer contacts and a lower probability of infection given contact. The average number of secondary cases from an individual aged 20–24 during lockdown drops to 0.9 and the transmission parameter *β*, which captures the probability of infection given contact, is decreased from 0.18 at baseline to 0.11. In the post-lockdown period daily contacts are increased to a higher proportion of their baseline values (see Methods); in order to keep the reproduction number equal to the dominant eigenvalue of the next generation matrix, *β* is consequently reduced to 0.05 when *R**t* =0.9 and to 0.07 when *R**t* =1.2. This implies that, to maintain *R**t* below one when more contacts are occurring in the population post-lockdown, the probability that contact results in infection will need to be reduced. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/04/2020.07.24.20157982/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/02/04/2020.07.24.20157982/F3) Figure 3: Next generation matrix (*K* = *k**ij*) showing the number of secondary cases generated by an index case from age group *j* (rows) in age group *i* (columns). The matrices are shown for different time points; at baseline before the implementation of interventions; during the lockdown period; and in the post-lockdown period when the effective reproduction number (*R**t*) rises to 0.9–1.2. The average number of secondary cases generated by an index case from age group *j* is the summation of row *j*. ### Infection dynamics For the first 100 days until the end of the lockdown, the infection dynamics are equivalent across the four immunity scenarios S1–S4 (Figure 4, panels A, C, E & G). After this time the dynamics depend on both the rate at which recovered individuals lose immunity and become susceptible again, and the post-lockdown *R**t*. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/04/2020.07.24.20157982/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/02/04/2020.07.24.20157982/F4) Figure 4: Projections from immunity scenarios S1–4 with post-lockdown *R**t* ranging from 0.9–1.2. Left panels show the number of infected, both asymptomatic and symptomatic (*I**A* + *I**S*), with SARS-CoV-2 in the UK population over time. Right panels show the proportion of the UK population with immunity (compartments *R**H* + *R**N*). Dashed vertical lines indicate the lockdown period; 23rd of March–15th of June 2020. Given our model and parameters, on the first day the intervention is imposed (23rd March 2020) there are 108,000 new SARS-CoV-2 cases, which is within the 95% credible interval (CrI) of new cases estimated for the UK on that day (95% CrI 54,000–155,000 [60]), and 75,000 people are infectious (infected compartments *I**A* + *I**S*) on this date. From the 3rd of March until the 23rd of March there are 559,000 cumulative cases across all age groups and 524,000 in adults ≥19 years, which is within the credible interval for an estimate of cumulative cases in this period (95% CrI 266,000–628,000 [60]). When most of the lockdown measures were eased in June, 5.7% of the total population and 6.8% of adults aged ≥19 years have immunity to SARS-CoV-2 (in recovered classes *R**H* and *R**N*), which is comparable to estimates of antibody levels in the UK population, estimated as 6.8% of blood donors on 24th May 2020 (95% confidence interval 5.2–8.6%; individuals ≥18 years [25]). ### Secondary peak in infections A secondary peak in infections is expected in spring 2021 where *R**t* = 1.1 or winter 2020 where *R**t* = 1.2 (Figure 4, panels E & G). The height of the secondary peak is determined by the rate at which immunity is lost. In our worst case scenario (S4: short-lived immunity) where immunity lasts an average of three months for non-hospitalised patients, a year for hospitalised patients and *R**t* following lockdown is 1.2, then the secondary peak will exceed the initial peak with a maximum of 387,000 infectious individuals and 125,000 daily new cases in December 2020. This is three-fold greater than the number of new cases in scenario S1 where immunity is permanent; the maximum number of infectious individuals in the secondary peak is 126,000 and there are 41,000 daily new cases (Figure 4G). We note that the timing of the secondary peak in infection curves across immunological scenarios are closely synchronised and in winter 2020–21. This synchrony and timing is also observed during the epidemic when values of *R**t* post-lockdown are greater than 1.2 (explored for values of *R**t* from 1.3 to 2.0). When *R**t* following lockdown is 1.1 the differences between the immunity scenarios is even more pronounced with a ten-fold difference in the height of the secondary peak between a scenario of permanent immunity and one of short-lived immunity, however the secondary peak does not exceed the height (maximum number of infectious individuals) of the initial wave in any scenario. When immunity wanes rapidly, the second peak is observed in April 2021 with a maximum of 154,000 infectious individuals and 50,000 daily new cases. By contrast when immunity is permanent, the number of new infections decays slowly rather than accelerates and there are projected to be a maximum of only 15,000 infectious individuals and 5,000 daily new cases in April 2021 (Figure 4E). ### Population immunity The dynamics of population immunity (recovered compartments *R**H* +*R**N*) are similarly shaped by the expected duration of immunity against SARS-CoV-2 and the post-lockdown *R**t*. Immunity decays from midway through the lockdown period in scenarios S2–S4 of waning (12 or 6 months) and short-lived immunity and resurges following a secondary wave of infection if *R**t* *>* 1 (Figure 4, panels F & H). After lockdown, a fall in the proportion of the population immune to the virus is observed until autumn 2020 for all values of *R**t*, after which point the secondary wave, if *R**t* *>* 1, causes the proportion of the population immune to rise again. This suggests that longitudinal serological surveys to detect waning immunity would be most informative when conducted in the period June–October 2020. ### Consequences of age structure The large differences in the heights of the secondary peaks when *R**t* *>* 1 between immunological scenarios (Figure 4, panels A, C, E & G) can be explained by the heterogeneity in transmission (see the next generation matrix in Figure 3). Infectious and immune cases as a proportion of the total age group are shown in Figure 5 for scenarios S1 & S4 of permanent and short-lived immunity where *R**t* = 1.2 following lockdown. A higher proportion of individuals aged between 20–39 are infected early in the epidemic, and this leads to 12.4–13.9% of individuals in these age groups having antibodies by October 2020 whe’n immunity is life-long (Figure 5B). When immunity wanes, however, by October 2020 this drops to 6.7–7.8% (Figure 5D), thus replenishing the pool of susceptible individuals with the age groups that disproportionately drive transmission. This causes the secondary peak of infectious cases to rise more rapidly and to a greater height when immunity wanes rapidly (Figure 5C), compared with permanent immunity (Figure 5A). Our models suggest that the age distribution of cases in the epidemic will not change greatly over time; as seen in Figure 5 the ordering of the proportion of each age group infected remains constant in both scenarios of permanent and short-lived immunity. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/04/2020.07.24.20157982/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2021/02/04/2020.07.24.20157982/F5) Figure 5: Projections from immunity scenarios S1 & S4 with post-lockdown *R**t* of SARS-CoV-2 at 1.2 in the UK population over time. Left panels show the proportion of each age group infected, for both asymptomatic and symptomatic (*I**A* + *I**S*) individuals. Right panels show the proportion of the each age group with immunity (compartments *R**H* + *R**N*). Dashed vertical lines indicate the intervention (lockdown) period; 23rd of March – 15th of June 2020. ### Longer term dynamics: extinction or endemic equilibrium We explored the impact of waning immunity and *R**t* on the equilibrium values for the different simulations over a longer, five year, period until March 2025 (Table 2). If the post-lockdown *R**t* is suppressed below one following lockdown, then the differences in immunity will have less impact on the longer-term infection dynamics, assuming no imported cases, as transmission of SARS-CoV-2 becomes unsustainable and the virus reaches extinction between April–November 2021 depending on the immunity scenario. In simulations where *R**t* equals one, if immunity is permanent then the epidemic becomes extinct in May 2022. When immunity wanes there is no secondary peak (Figure 4C), however the infections persist at a low level of endemicity equivalent to 107, 236 and 1,122 daily cases in immunity scenarios S2–S4, respectively. For larger values of *R**t*, and where immunity wanes, the system oscillates with subsequent peaks of infection over the next five years until a steady state is reached. We find that, if *R**t* = 1.2 post-lockdown and immunity is short-lived, there is the potential for over 73,000 new cases daily; 5,000 hospitalisations; and 913 intensive care unit (ICU) admissions (calculated as 17% of all hospitalised cases [61]) at endemic equilibrium (January 2025), which would be sufficient to overwhelm contact tracing services and ICU capacity [62, 63]. View this table: [Table 2:](http://medrxiv.org/content/early/2021/02/04/2020.07.24.20157982/T2) Table 2: Values at equilibrium from the modelled scenarios for SARS-CoV-2 in the United Kingdom, explored over a five year horizon (March 2020 to March 2025). 1Effective reproduction number of SARS-CoV-2 after lockdown. 2Assumed duration of immunity for hospitalised and non-hospitalised individuals, see Methods for details of immunity scenarios S1–S4. 3Number of individuals newly infected with SARS-CoV-2 that enter the exposed *E* state. 4Number of symptomatic individuals with SARS-CoV-2 that enter the recovered hospitalised (*R**H*) state. 5Number of hospitalised individuals admitted to intensive care units (ICU), under the assumption that 17% of hospitalised cases in the UK require care in high dependency units [61]. 6Either when the number of daily new cases drops below one (extinction), or when the daily new cases are the same integer value over a sustained period (endemic equilibrium). If models take longer than five years to reach a steady state, the values on the 1st of March 2025 are reported. ## Discussion Despite only 6.8% of the adult UK population having immunity against SARS-CoV-2 in our simulation at the end of the lockdown, the modelled scenarios suggest that, if this acquired immunity wanes over time, there are substantive differences to the subsequent infection dynamics. Waning immunity impacts on the height of the secondary peak and, in the absence of future interventions, establishes the virus at levels of endemic equilibrium that could overwhelm contact tracing services and ICU capacity [62, 63]. We predict that surveys to detect waning immunity at the population level would be most effective when carried out in the period between the end of lockdown and autumn 2020, as after this point an upsurge in cases is expected that will increase the proportion of the population with immunity to SARS-CoV-2. In particular, this will allow evaluation of whether specific antibodies generated against the virus are short-lived if reductions in antibody prevalence are observed at the population level. We find that transmission is driven disproportionately by individuals of working age, and subsequently a higher proportion of individuals aged 20–39 years become infected early in the pandemic and subsequently develop immunity (Figures 3 & 5). This prediction is borne out by serological data from Switzerland, which showed that individuals aged 20–49 years were significantly more likely to be seropositive in May 2020 compared with younger and older age groups [64]. We postulate that ‘key workers’ in the UK population who have continued to work during the lockdown are more likely to have antibodies against SARS-CoV-2. Higher immunity among individuals of working age has the effect of slowing the subsequent epidemic when immunity is permanent. Conversely, when immunity wanes, previously infected individuals of working age rejoin the susceptible pool and so contribute again to transmission; leading to a high growth rate and a larger secondary peak of infected cases. In these circumstances, efforts to suppress transmission will be challenging in the absence of a transmission-blocking vaccine [15]. We note that the model structure developed here is capable of simulating the impact of vaccination with a vaccine that provides temporary transmission-blocking immunity, and could be used to predict the optimal timing for booster shots. The projected trajectory of the epidemic after lockdown is highly sensitive to the effective reproduction number, with model behaviour for values of *R**t* slightly above or below one displaying qualitatively different dynamics (Figure 4). Although we do not describe model output for higher values of *R**t*, the dynamics over a wider range of effective reproduction numbers post-lockdown (1.3–2.0) are qualitatively the same as when *R**t* =1.2. This shows the importance of timely and accurate estimates of *R**t* to inform control strategies, and ensuring widespread community testing and contact tracing is in place. Our calculations show that to suppress *R**t* below one when contact rates rise to a higher fraction of baseline (pre-lockdown) values, the probability of infection given contact (represented by the *β* parameter), must drop by around half. Interventions that have the potential to reduce the probability of infection include social distancing; regular hand washing; and the wearing of face masks outside the home [65]. Our study reinforces the importance of better understanding SARS-CoV-2 immunity among recovered individuals of different ages and disease severity. In scenarios where immunity wanes and *R**t* following lockdown is greater than one, the SARS-CoV-2 epidemic never reaches extinction due to herd immunity, but rather the number of infected cases oscillates with subsequent waves of infection before reaching endemic equilibrium (Table 2). Even in simulations where the post-lockdown reproduction number only narrowly exceeds one (*R**t* =1.1), if immunity wanes over an average of one year for severe cases and three months for non-severe cases, this is projected to lead to an equilibrium state of nearly 40,000 daily new cases and 500 daily admissions to intensive care. Policy strategies aiming to achieve herd immunity are therefore risky [5]. As individuals can lose their immunity and become reinfected with SARS-CoV-2, a stable population immunity threshold can never be reached in the absence of a vaccine with lifelong efficacy [21]. The establishment of an endemic equilibrium state is dependent on no future interventions or changes to *R**t*, which we consider unlikely as policy makers and public health agencies are likely to react to future outbreaks with localised control measures. One of the strengths of our study is that the model is calibrated to key features of the UK epidemic. While we did not explicitly fit to data, new cases at the start of the lockdown; cumulative cases in March; and the proportion of the adult population with antibodies to SARS-CoV-2 are highly comparable between our output and current estimates [25, 60]. We used contact matrices from a comprehensive study of contact patterns in the UK population [47] in addition to demographic data from the Office for National Statistics, to give our simulations the best chance of capturing realistic age-specific transmission patterns in the UK population. Plausible estimates on which to base expectations for the duration of immunity are sparse in the current literature. Rosado *et al*. estimated that antibodies could wane in 50% of recovered individuals after one year [57], which is similar to the estimated duration of immunity against seasonal circulating coronaviruses [13, 29]. Even with this consideration, there are many probability distributions that can be used to capture a median duration of immunity, and our selection of an Erlang distribution with a shape parameter of two is somewhat arbitrary. Our assumptions on the duration of the latent and infectious periods are more closely informed by estimates from data [35, 55, 56]. We made the decision to capture the expected duration of these states as a discretised Erlang distribution rather than the, more conventional, exponential distribution. This has the benefit of closely replicating fitted gamma or log-normal distributions within a compartmental model [36], and has important implications for the dynamics of the epidemic [66, 67]. We make a number of assumptions regarding the natural history of the virus, such as the relative susceptibility of children compared with adults and the relative infectiousness of symptomatic versus asymptomatic cases based on the current literature [52, 49]. Future empirical studies are likely to add to and further refine these epidemiological parameters. After we completed the analysis, a study of 37 asymptomatic individuals in China were found to have a longer period of viral shedding when compared with symptomatic individuals [10]. While viral shedding is not necessarily indicative of transmission potential [7], if these findings are replicated in larger studies this may suggest a need to use different durations of infectiousness for asymptomatic and symptomatic infections in subsequent models. We have aimed to capture future infection dynamics at a national level in the UK under a range of scenarios. Our analysis is limited by not considering regional differences in transmission rates, for instance through a patch (metapopulation) model [48], or a stochastic approach that allows for local extinction events [21]. There are no deaths in our model, either from demography or infection. Accounting for mortality would mainly affect dynamics in the oldest age group (over 70 years) [16, 61], and the higher probability of disease-induced mortality would reduce the build up of immunity in this age group (Figure 5D). We also do not explicitly consider transmission in settings such as schools, hospitals or care homes, although such dynamics may be captured indirectly through the contact matrix. Given the simplicity of the model structure, we advise against treating the output as an exact prediction of the future. In addition to the limitations listed above, the epidemic trajectory will be substantially altered by any future interventions such as a return to full lockdown conditions, or intensive contact tracing and isolation [29, 68]. We emphasise that the average duration of immunity following recovery impacts substantially on epidemic growth following the easing of non-pharmaceutical interventions, and that slight differences in *R**t* above or below one can lead to qualitatively different infection dynamics. ## Data Availability This is a modelling study and does not use primary data. Open access code to reproduce the analysis is provided. [https://github.com/tc13/covid-19-immunity](https://github.com/tc13/covid-19-immunity) ## Disclosure statement The authors declare no conflicts of interest. No other authors have been paid to write this article. ## Open access Anyone can share and adapt this article provided they give credit and link to this article’s DOI and CC BY 4.0 licence (creativecommons.org/licences/by/4.0). ## Funding TC is funded by a Sir Henry Wellcome Postdoctoral Fellowship from the Wellcome Trust (reference 215919/Z/19/Z). ELD, TCDL, DA, AB, TMP and TDH gratefully acknowledge funding of the NTD Modelling Consortium by the Bill & Melinda Gates Foundation (BMGF) (grant number OPP1184344). Views, opinions, assumptions or any other information set out in this article should not be attributed to BMGF or any person connected with them. TMP’s PhD was supported by the Engineering & Physical Sciences Research Council, Medical Research Council and University of Warwick (grant number EP/L015374/1). All funders had no role in the study design, collection, analysis, interpretation of data, writing of the report, or decision to submit the manuscript for publication. ## Acknowledgements This work was undertaken as a contribution to the Rapid Assistance in Modelling the Pandemic (RAMP) initiative, coordinated by The Royal Society. The authors thank Prof. Julia Gog and Dr. Rosalind Eggo for their helpful comments and input. ## Footnotes * Revisions made following peer review at Philosophical Transactions B * Received July 24, 2020. * Revision received February 4, 2021. * Accepted February 4, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. [1].Worldometer. Covid-19 data. [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/), June 2020. 2. [2].Esteban Ortiz-Ospina Max Roser, Hannah Ritchie and Joe Hasell. Our World in Data. Coronavirus (COVID-19) cases. [https://ourworldindata.org/covid-cases](https://ourworldindata.org/covid-cases), June 2020. 3. [3]. Kristian G. Andersen, Andrew Rambaut, W. Ian Lipkin, Edward C. Holmes, and Robert F. Garry. The proximal origin of SARS-CoV-2. Nature Medicine, 26(4):450–452, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0820-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 4. [4].UK Government. Coronavirus (COVID-19) in the uk. [https://coronavirus.data.gov.uk/](https://coronavirus.data.gov.uk/), June 2020. 5. [5]. Heba Habib. Has Sweden’s controversial COVID-19 strategy been successful? BMJ, 369, 2020. 6. [6]. Quan-Xin Long, Bai-Zhong Liu, Hai-Jun Deng, Gui-Cheng Wu, Kun Deng, Yao-Kai Chen, Pu Liao, Jing-Fu Qiu, Yong Lin, Xue-Fei Cai, et al. Antibody responses to SARS-CoV-2 in patients with COVID-19. Nature Medicine, 26:845–848, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 7. [7]. Roman Wölfel, Victor M. Corman, Wolfgang Guggemos, Michael Seilmaier, Sabine Zange, Marcel A. Müller, Daniela Niemeyer, Terry C. Jones, Patrick Vollmar, Camilla Rothe, et al. Virological assessment of hospitalized patients with COVID-2019. Nature, 581:465–469, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2196-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 8. [8]. Linlin Bao, Wei Deng, Hong Gao, Chong Xiao, Jiayi Liu, Jing Xue, Qi Lv, Jiangning Liu, Pin Yu, Yanfeng Xu, et al. Lack of reinfection in rhesus macaques infected with SARS-CoV-2. bioRxiv, 2020. 9. [9]. Juanjuan Zhao, Quan Yuan, Haiyan Wang, Wei Liu, Xuejiao Liao, Yingying Su, Xin Wang, Jing Yuan, Tingdong Li, Jinxiu Li, et al. Antibody responses to SARS-CoV-2 in patients of novel coronavirus disease 2019. Clinical Infectious Diseases, 2020. ciaa344. 10. [10]. Quan-Xin Long, Xiao-Jun Tang, Qiu-Lin Shi, Qin Li, Hai-Jun Deng, Jun Yuan, Jie-Li Hu, Wei Xu, Yong Zhang, Fa-Jin Lv, et al. Clinical and immunological assessment of asymptomatic SARS-CoV-2 infections. Nature Medicine, 26:1200–1204, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591020-0965-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 11. [11]. Wu-Chun Cao, Wei Liu, Pan-He Zhang, Fang Zhang, and Jan H. Richardus. Disappearance of antibodies to SARS–associated coronavirus after recovery. New England Journal of Medicine, 357(11):1162–1163, 2007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc070348&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17855683&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249387200031&link_type=ISI) 12. [12]. Pyoeng Gyun Choe, R. A. P. M. Perera, Wan Beom Park, Kyoung-Ho Song, Ji Hwan Bang, Eu Suk Kim, Hong Bin Kim, Long Wei Ronald Ko, Sang Won Park, Nam-Joong Kim, et al. MERS-CoV antibody responses 1 year after symptom onset, South Korea, 2015. Emerging Infectious Diseases, 23(7):1079, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2307.170310&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28585916&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 13. [13]. Arthur WD Edridge, Joanna Kaczorowska, Alexis CR Hoste, Margreet Bakker, Michelle Klein, Katherine Loens, Maarten F Jebbink, Amy Matser, Cormac M Kinsella, Paloma Rueda, et al. Seasonal coronavirus protective immunity is short-lasting. Nature Medicine, 2020. 14. [14]. Patience K. Kiyuka, Charles N. Agoti, Patrick K. Munywoki, Regina Njeru, Anne Bett, James R. Otieno, Grieven P. Otieno, Everlyn Kamau, Taane G. Clark, Lia van der Hoek, et al. Human coronavirus nl63 molecular epidemiology and evolutionary patterns in rural coastal Kenya. The Journal of Infectious Diseases, 217(11):1728–1739, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiy098&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29741740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 15. [15]. Paul Kellam and Wendy Barclay. The dynamics of humoral immune responses following SARS-CoV-2 infection and the potential for reinfection. Journal of General Virology, 101(8):791–797, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/jgv.0.001439&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32430094&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 16. [16]. Neil Ferguson, Daniel Laydon, Gemma Nedjati Gilani, Natsuko Imai, Kylie Ainslie, Marc Baguelin, Sangeeta Bhatia, Adhiratha Boonyasiri, Zulma Cucunuba Perez, Gina Cuomo-Dannenburg, et al. Report 9: Impact of non-pharmaceutical interventions (NPIs) to reduce COVID19 mortality and healthcare demand. 2020. 17. [17]. Adam J. Kucharski, Timothy W. Russell, Charlie Diamond, Yang Liu, John Edmunds, Sebastian Funk, and Rosalind M. Eggo. Early dynamics of transmission and control of COVID-19: a mathematical modelling study. The Lancet Infectious Diseases, 20(5):553–558, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s1473-3099(20)30144-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 18. [18]. Kiesha Prem, Yang Liu, Timothy W. Russell, Adam J. Kucharski, Rosalind M. Eggo, Nicholas Davies, Stefan Flasche, Samuel Clifford, Carl A. B. Pearson, James D. Munday, et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. The Lancet Public Health, 5(5):e261–e270, 2020. 19. [19]. T Déirdre Hollingsworth, Don Klinkenberg, Hans Heesterbeek, and Roy M Anderson. Mitigation strategies for pandemic influenza a: balancing conflicting policy objectives. PLoS Comput Biol, 7(2):e1001076, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1001076&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21347316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 20. [20]. M. Gabriela M. Gomes, Lisa J. White, and Graham F. Medley. Infection, reinfection, and vaccination under suboptimal immune protection: epidemiological perspectives. Journal of Theoretical Biology, 228(4):539–549, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2004.02.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15178201&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222157400009&link_type=ISI) 21. [21]. Matt J Keeling and Pejman Rohani. Modeling infectious diseases in humans and animals. Princeton University Press, 2011. 22. [22]. Henrik Salje, Cécile Tran Kiem, Noémie Lefrancq, Noémie Courtejoie, Paolo Bosetti, Juliette Paireau, Alessio Andronico, Nathanaël Hozé, Jehanne Richet, Claire-Lise Dubost, et al. Estimating the burden of SARS-CoV-2 in France. Science, 369(6500):208–211, 2020. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjkvNjUwMC8yMDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMi8wNC8yMDIwLjA3LjI0LjIwMTU3OTgyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 23. [23]. Seth Flaxman, Swapnil Mishra, Axel Gandy, H Juliette T Unwin, Thomas A Mellan, Helen Coupland, Charles Whittaker, Harrison Zhu, Tresnia Berah, Jeffrey W Eaton, et al. Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature, pages 1–8, 2020. 24. [24]. Alessia Lai, Annalisa Bergna, Carla Acciarri, Massimo Galli, and Gianguglielmo Zehender. Early phylogenetic estimate of the effective reproduction number of SARS-CoV-2. Journal of Medical Virology, 92(6):675–679, 2020. 25. [25].Office for National Statistics. Coronavirus (COVID-19) Infection Survey pilot: 28 May 2020. Technical report, 2020. 26. [26]. Roy M. Anderson and Robert M. May. Infectious Diseases of Humans: Dynamics and Control. Oxford University Press, 1992. 27. [27]. M. Gabriela M. Gomes, Ricardo Aguas, Rodrigo M. Corder, Jessica G. King, Kate E. Langwig, Caetano Souto-Maior, Jorge Carneiro, Marcelo U. Ferreira, and Carlos Penha-Goncalves. Individual variation in susceptibility or exposure to SARS-CoV-2 lowers the herd immunity threshold. medRxiv, 2020. 28. [28]. Tom Britton, Frank Ball, and Pieter Trapman. A mathematical model reveals the influence of population heterogeneity on herd immunity to sars-cov-2. Science, 369(6505):846–849, 2020. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjkvNjUwNS84NDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMi8wNC8yMDIwLjA3LjI0LjIwMTU3OTgyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 29. [29]. Stephen M. Kissler, Christine Tedijanto, Edward Goldstein, Yonatan H. Grad, and Marc Lipsitch. Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science, 368(6493):860–868, 2020. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5My84NjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMi8wNC8yMDIwLjA3LjI0LjIwMTU3OTgyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 30. [30]. Kelvin Kai-Wang To, Ivan Fan-Ngai Hung, Jonathan Daniel Ip, Allen Wing-Ho Chu, Wan-Mui Chan, Anthony Raymond Tam, Carol Ho-Yan Fong, Shuofeng Yuan, Hoi-Wah Tsoi, Anthony Chin-Ki Ng, Larry Lap-Yip Lee, Polk Wan, Eugene Yuk-Keung Tso, Wing-Kin To, Dominic Ngai-Chong Tsang, Kwok-Hung Chan, Jian-Dong Huang, Kin-Hang Kok, Vincent Chi-Chung Cheng, and Kwok-Yung Yuen. Coronavirus Disease 2019 (COVID-19) Re-infection by a Phylogenetically Distinct Severe Acute Respiratory Syndrome Coronavirus 2 Strain Confirmed by Whole Genome Sequencing. Clinical Infectious Diseases, 2020. 31. [31]. Richard L Tillett, Joel R Sevinsky, Paul D Hartley, Heather Kerwin, Natalie Crawford, Andrew Gorzalski, Chris Laverdure, Subhash C Verma, Cyprian C Rossetto, David Jackson, et al. Genomic evidence for reinfection with SARS-CoV-2: a case study. The Lancet Infectious Diseases, 2020. 32. [32]. Belen Prado-Vivar, Monica Becerra-Wong, Juan Jose Guadalupe, Sully Marquez, Bernardo Gutierrez, Patricio Rojas-Silva, Michelle Grunauer, Gabriel Trueba, Veronica Barragan, and Paul Cardenas. COVID-19 Re-Infection by a Phylogenetically Distinct SARS-CoV-2 Variant, First Confirmed Event in South America. SSRN, 2020. 33. [33]. Jan Van Elslande, Pieter Vermeersch, Kris Vandervoort, Tony Wawina-Bokalanga, Bert Vanmechelen, Elke Wollants, Lies Laenen, Emmanuel André, Marc Van Ranst, Katrien Lagrou, and Piet Maes. Symptomatic SARS-CoV-2 reinfection by a phylogenetically distinct strain. Clinical Infectious Diseases, 2020. ciaa1330. 34. [34]. Vivek Gupta, Rahul C Bhoyar, Abhinav Jain, Saurabh Srivastava, Rashmi Upadhayay, Mohamed Imran, Bani Jolly, Mohit Kumar Divakar, Disha Sharma, Paras Sehgal, Gyan Ranjan, Rakesh Gupta, Vinod Scaria, and Sridhar Sivasubbu. Asymptomatic reinfection in two healthcare workers from India with genetically distinct SARS-CoV-2. Clinical Infectious Diseases, 2020. 35. [35]. Xi He, Eric H. Y. Lau, Peng Wu, Xilong Deng, Jian Wang, Xinxin Hao, Yiu Chung Lau, Jessica Y. Wong, Yujuan Guan, Xinghua Tan, et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nature Medicine, 26:672–675, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-3012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 36. [36]. Helen J. Wearing, Pejman Rohani, and Matt J. Keeling. Appropriate models for the management of infectious diseases. PLoS Medicine, 2(7), 2005. 37. [37]. James O. Lloyd-Smith, Alison P. Galvani, and Wayne M. Getz. Curtailing transmission of severe acute respiratory syndrome within a community and its hospital. Proceedings of the Royal Society of London. Series B: Biological Sciences, 270(1528):1979–1989, 2003. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2003.2481&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14561285&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185739600001&link_type=ISI) 38. [38]. Steven Riley, Christophe Fraser, Christl A Donnelly, Azra C Ghani, Laith J Abu-Raddad, Anthony J Hedley, Gabriel M Leung, Lai-Ming Ho, Tai-Hing Lam, Thuan Q Thach, et al. Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions. Science, 300(5627):1961–1966, 2003. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMDAvNTYyNy8xOTYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDIvMDQvMjAyMC4wNy4yNC4yMDE1Nzk4Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 39. [39]. Wayne M Getz and Eric R Dougherty. Discrete stochastic analogs of erlang epidemic models. Journal of biological dynamics, 12(1):16–38, 2018. 40. [40]. Paul J Hurtado and Adam S Kirosingh. Generalizations of the ‘linear chain trick’: incorporating more flexible dwell time distributions into mean field ode models. Journal of mathematical biology, 79(5):1831–1883, 2019. 41. [41]. Christopher I Jarvis, Kevin Van Zandvoort, Amy Gimma, Kiesha Prem, Petra Klepac, G James Rubin, and W John Edmunds. Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK. BMC Medicine, 18:1–10, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12916-020-01597-8&link_type=DOI) 42. [42]. J Kucharski Adam, Klepac Petra, JK Andrew, M Kissler Stephen, L Tang Maria, Fry Hannah, R Julia, CMMID COVID-19 working group, et al. Effectiveness of isolation, testing, contact tracing, and physical distancing on reducing transmission of SARS-CoV-2 in different settings: A mathematical modelling study. The Lancet Infectious Diseases, pages S1473–3099, 2020. 43. [43]. Enrico Lavezzo, Elisa Franchin, Constanze Ciavarella, Gina Cuomo-Dannenburg, Luisa Barzon, Claudia Del Vecchio, Lucia Rossi, Riccardo Manganelli, Arianna Loregian, Nicolò Navarin, et al. Suppression of COVID-19 outbreak in the municipality of Vo, Italy. medRxiv, 2020. 44. [44]. Robert Verity, Lucy C Okell, Ilaria Dorigatti, Peter Winskill, Charles Whittaker, Natsuko Imai, Gina Cuomo-Dannenburg, Hayley Thompson, Patrick GT Walker, Han Fu, et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. The Lancet Infectious Diseases, 20(6):669–677, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30243-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32240634&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 45. [45]. Kenji Mizumoto, Katsushi Kagaya, Alexander Zarebski, and Gerardo Chowell. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Eurosurveillance, 25(10):2000180, 2020. 46. [46]. C. Li, F. Ji, L. Wang, J. Hao, M. Dai, Y. Liu, X. Pan, J. Fu, L. Li, G. Yang, et al. Asymptomatic and human-to-human transmission of SARS-CoV-2 in a 2-family cluster, Xuzhou, China. Emerging Infectious Diseases, 26(7), 2020. 47. [47]. Petra Klepac, Stephen Kissler, and Julia Gog. Contagion! The BBC Four Pandemic—the model behind the documentary. Epidemics, 24:49–59, 2018. 48. [48]. Petra Klepac, Adam J. Kucharski, Andrew J. K. Conlan, Stephen Kissler, Maria Tang, Hannah Fry, and Julia R. Gog. Contacts in context: large-scale setting-specific social mixing matrices from the BBC Pandemic project. medRxiv, 2020. 49. [49]. Nicholas G. Davies, Petra Klepac, Yang Liu, Kiesha Prem, Mark Jit, CMMID COVID-19 working group Eggo, Rosalind M., et al. Age-dependent effects in the transmission and control of COVID-19 epidemics. Nature Medicine, 26:1205–1211, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0962-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32546824&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 50. [50]. Juanjuan Zhang, Maria Litvinova, Yuxia Liang, Yan Wang, Wei Wang, Shanlu Zhao, Qianhui Wu, Stefano Merler, Cecile Viboud, Alessandro Vespignani, et al. Age profile of susceptibility, mixing, and social distancing shape the dynamics of the novel coronavirus disease 2019 outbreak in China. medRxiv, 2020. 51. [51]. Daniel F. Gudbjartsson, Agnar Helgason, Hakon Jonsson, Olafur T. Magnusson, Pall Melsted, Gudmundur L. Norddahl, Jona Saemundsdottir, Asgeir Sigurdsson, Patrick Sulem, Arna B. Agustsdottir, et al. Spread of SARS-CoV-2 in the Icelandic population. New England Journal of Medicine, 382:2302–2315, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2006100&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 52. [52]. Qin-Long Jing, Ming-Jin Liu, Zhou-Bin Zhang, Li-Qun Fang, Jun Yuan, An-Ran Zhang, Natalie E Dean, Lei Luo, Meng-Meng Ma, Ira Longini, et al. Household secondary attack rate of COVID-19 and associated determinants in Guangzhou, China: a retrospective cohort study. The Lancet Infectious Diseases, 20(10):1141–1150, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 53. [53]. Jiuxin Qu, Chi Wu, Xiaoyong Li, Guobin Zhang, Zhaofang Jiang, Xiaohe Li, Lei Liu, et al. Profile of IgG and IgM antibodies against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Clinical Infectious Diseases, 2020. ciaa489. 54. [54]. Akiko Iwasaki. What reinfections mean for COVID-19. The Lancet Infectious Diseases, 2020. 55. [55]. Qun Li, Xuhua Guan, Peng Wu, Xiaoye Wang, Lei Zhou, Yeqing Tong, Ruiqi Ren, Kathy SM Leung, Eric HY Lau, Jessica Y Wong, et al. Early transmission dynamics in Wuhan, China, of novel coronavirus– infected pneumonia. New England Journal of Medicine, 382:1199–1207, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) 56. [56]. Nicholas G Davies, Adam J Kucharski, Rosalind M Eggo, Amy Gimma, W John Edmunds, Thibaut Jombart, Kathleen O’Reilly, Akira Endo, Joel Hellewell, Emily S Nightingale, et al. Effects of non-pharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: a modelling study. The Lancet Public Health, 5(7):e375–e385, 2020. 57. [57]. Jason Rosado, Charlotte Cockram, Sarah Helene Merkling, Caroline Demeret, Annalisa Meola, Solen Kerneis, Benjamin Terrier, Samira Fafi-Kremer, Jerome de Seze, Marija Backovic, et al. Serological signatures of SARS-CoV-2 infection: Implications for antibody-based diagnostics. medRxiv, 2020. 58. [58].Centre for the Mathematical Modelling of Infectious Disases, London School of Hygiene and Tropical Medicine. Temporal variation in transmission during the COVID-19 outbreak. [https://epiforecasts.io/covid/](https://epiforecasts.io/covid/), June 2020. 59. [59]. Louis du Plessis, John T McCrone, Alexander E Zarebski, Verity Hill, Christopher Ruis, Bernardo Gutierrez, Jayna Raghwani, Jordan Ashworth, Rachel Colquhoun, Thomas R Connor, et al. Establishment & lineage dynamics of the sars-cov-2 epidemic in the uk. medRxiv, 2020. 60. [60]. Mark Jit, Thibaut Jombart, Emily S. Nightingale, Akira Endo, Sam Abbott, L.S.H.T.M. Centre for Mathematical Modelling of Infectious Diseases COVID-19 Working Group, and W. John Edmunds. Estimating number of cases and spread of coronavirus disease (COVID-19) using critical care admissions, United Kingdom, February to March 2020. Eurosurveillance, 25(18), 2020. 61. [61]. Annemarie B. Docherty, Ewen M. Harrison, Christopher A. Green, Hayley E. Hardwick, Riinu Pius, Lisa Norman, Karl A. Holden, Jonathan M. Read, Frank Dondelinger, Gail Carson, et al. Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: prospective observational cohort study. BMJ, 369, 2020. 62. [62].NHS England. Monthly critical care beds and cancelled urgent operations (MSITREP) data, England. [https://www.england.nhs.uk/statistics/wp-content/uploads/sites/2/2020/04/MSitRep-SPN-FEBRUARY-2020-auY71.pdf](https://www.england.nhs.uk/statistics/wp-content/uploads/sites/2/2020/04/MSitRep-SPN-FEBRUARY-2020-auY71.pdf), April 2020. 63. [63]. Emma L Davis, Tim CD Lucas, Anna Borlase, Timothy M Pollington, Sam Abbott, Diepreye Ayabina, Thomas Crellen, Joel Hellewell, Li Pi, Graham F Medley, et al. An imperfect tool: COVID-19 ‘test & trace’ success relies on minimising the impact of false negatives and continuation of physical distancing. medRxiv, 2020. 64. [64]. Silvia Stringhini, Ania Wisniak, Giovanni Piumatti, Andrew S Azman, Stephen A Lauer, Héléne Baysson, David De Ridder, Dusan Petrovic, Stephanie Schrempft, Kailing Marcus, et al. Seroprevalence of anti-SARS-CoV-2 IgG antibodies in Geneva, Switzerland (SEROCoV-POP): a population-based study. The Lancet, 396(10247):313–319, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)31304-0&link_type=DOI) 65. [65]. Benjamin J Cowling and Allison E Aiello. Public health measures to slow community spread of coronavirus disease 2019. The Journal of Infectious Diseases, 221(11):1749–1751, 2020. 66. [66]. Alun L Lloyd. Realistic distributions of infectious periods in epidemic models: changing patterns of persistence and dynamics. Theoretical population biology, 60(1):59–71, 2001. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/tpbi.2001.1525&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11589638&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F04%2F2020.07.24.20157982.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171555800004&link_type=ISI) 67. [67]. AC Fowler and T Déirdre Hollingsworth. Simple approximations for epidemics with exponential and fixed infectious periods. Bulletin of Mathematical Biology, 77(8):1539–1555, 2015. 68. [68]. Joel Hellewell, Sam Abbott, Amy Gimma, Nikos I Bosse, Christopher I Jarvis, Timothy W Russell, James D Munday, Adam J Kucharski, W John Edmunds, Fiona Sun, et al. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. The Lancet Global Health, 8(4):e488–e496, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2214-109X(20)30074-7&link_type=DOI) [1]: /embed/graphic-4.gif [2]: /embed/graphic-5.gif [3]: /embed/graphic-6.gif [4]: /embed/graphic-7.gif [5]: /embed/graphic-8.gif [6]: /embed/graphic-9.gif [7]: /embed/graphic-10.gif [8]: /embed/graphic-11.gif