Immunological characteristics govern the changing severity of COVID-19 during the transition to endemicity ========================================================================================================== * Jennie S Lavine * Ottar N Bjornstad * Rustom Antia ## 1 Abstract As prospects for eradicating CoV-2 dwindle, we are faced with the question of how the severity of CoV-2 disease may change in the years ahead. Will CoV-2 continue to be a pathogenic scourge that, like smallpox or measles, can be tamed only by ongoing vaccination, or will it join the ranks of mild endemic human coronaviruses (HCoVs)? Our analysis of immunological and epidemiological data on HCoVs shows that infection-blocking immunity wanes rapidly, but disease-reducing immunity is long-lived. We estimate the relevant parameters and incorporate them into a new epidemiological model framework which separates these different components of immunity. Our model recapitulates both the current severity of CoV-2 and the relatively benign nature of HCoVs; suggesting that once the endemic phase is reached, CoV-2 may be no more virulent than the common cold. The benign outcome at the endemic phase is contingent on the virus causing primary infections in children. We predict a very different outcome were a CoV like MERS (that causes severe disease in children) to become endemic. These results force us to re-evaluate control measures that rely on identifying and isolating symptomatic infections, and reconsider ideas regarding herd immunity and the use of immune individuals as shields to protect vulnerable groups. ## 2 Introduction Emerging pathogens have been invading the human host niche through history (such as smallpox emerging from an ancestral animal pox virus and rinderpest evolving to measles) to kill a substantial fraction of all people born. Recent decades have seen multiple challenges that cause severe viral disease, such as SARS, MERS, Hendra, Nipha and Ebola. Fortunately, all these are examples of emerging diseases that were locally contained. When containment is not successful, as is likely for the novel betacoronavirus SARS CoV-2 (CoV-2) [1, 2], we need to understand and plan for the transition to endemicity and continued circulation, with possible changes in disease severity due to virus evolution and build-up of host immunity and resistance. COVID-19 is an emerging infectious disease caused by CoV-2 that has a high basic re-productive number (*R*) and significant transmission during the asymptomatic phase, both of which make it hard to control [3]. However, CoV-2 is not the first coronavirus identified with human-to-human transmission: there are six other coronaviruses with known human chains of transmission, and these can provide clues to future scenarios for the current pandemic. There are four human coronaviruses (HCoVs) that circulate endemically across the globe; they cause only mild symptoms and are not a significant public health burden [4]. Another two HCoV strains, SARS CoV-1 and MERS, emerged in the past few decades and had higher case fatality ratios (CFRs) than COVID-19 but were contained and never spread widely [5, 6]. We propose a model to explore the potential changes in both transmission and disease severity of emerging HCoVs through the transition to endemicity. We focus on CoV-2, and discuss how the conclusions would differ for emerging coronaviruses more akin to SARS and MERS. Our hypothesis is that all HCoVs elicit immunity with similar characteristics, and the current acute public health problem is a consequence of epidemic emergence into an immunologically naive population in which frail age-groups have no previous exposure. We use our estimates of immunological and epidemiological parameters for endemic HCoVs to develop a quantitative transmission model and incorporate the age-dependence of severity from primary infections of emerging HCoVs. Our model explicitly considers three separate measures of the efficacy of immunity which can wane at different rates. These novel features allow us new insights into potential changes in both transmission and disease severity in the coming years (Fig 1). ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F1) Figure 1: Schematic of a compartmental model for the epidemiology of CoV infections. The model incorporates the following features of natural HCoV infections: recovered individuals (*R*1) have sterilizing immunity (high IE*S*) immediately after infection, and this immunity wanes over time at rate *ω*; individuals whose sterilizing immunity has waned (*R*2) can become reinfected (*I*2), albeit with mild symptoms (IE*P* is high) and reduced transmission to others (compared with transmission from primary infections, i.e., *ρ ≤* 1). We use an age-structured version of this model with parameters for demography of the US population in our simulations of transient dynamics. See SI Sec 2 for equations and details. Building on ideas from the vaccine modeling literature, immunity may be thought to provide protection in three ways [7]. In its most robust form, “sterilizing” immunity can prevent a pathogen from replicating, thereby rendering the host refractory to reinfection. We term this property immune efficacy with respect to susceptibility, IE*S*. If immunity does not prevent reinfection, it may still attenuate the pathology due to reinfection (IE*P*) and/or reduce transmissibility or infectiousness (IE*I*). Experimental re-exposure studies on endemic HCoVs provide evidence that the three IE’s do not wane at the same rate [8, 9]. Callow’s experimental study [8] shows that reinfection is possible within one year (relatively short IE*S*); however, upon reinfection symptoms are mild (high IE*P*) and the virus is cleared more quickly (moderate IE*I*). Details the derivation of the model are in SI section 2. ## 3 Results ### 3.1 Endemic dynamics and age structure We use available data on age-specific seroprevalence of IgM (acute response) and IgG (long-term memory) against the four circulating HCoVs to estimate parameter ranges for transmission and waning of immunity [10]. The data are plotted and further analyzed in Fig 2A. The rapid rise in seroprevelance indicates that primary infection with all four endemic HCoV strains happens early in life; and our analysis of the data gives us an estimate for the mean age of primary infection (MAPI) between 3.4 and 5.1 years, with almost everyone infected by age 15 (see SI Sec 1 for details). For most people to be infected so early in life – younger even than measles in the pre-vaccine era – there must be a very high attack rate, too high to be due only to transmission from primary infections. Steady state analysis of the model shows how this attack rate arises from a combination of high transmissibility from primary infections (i.e., high *R*) and waning of sterilizing immunity and significant transmission from reinfections in older individuals. The rapid waning of sterilizing immunity is also documented in experimental HCoV infections of humans which showed that reinfection is possible one year after an earlier infection, albeit with milder symptoms (IE*P*) and a shorter duration (IE*I*) [11]. Figure 2B shows the plausible combinations of waning immunity and transmission from reinfected individuals that is required to produce the MAPI observed in Figure 2A. Table 1 shows the ranges of the parameters used in our simulations. ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F2) Figure 2: **Panel A.** Using data published by Zhou et al (2013) [10], we plotted the mean proportion seropositive for IgG (green, top lines) and IgM (purple, bottom lines) against the four endemic HCoV strains (dots connected by dashed lines) along with the 95% CIs for the proportions (vertical lines going through each dot). Additionally, we calculate the mean age of primary infection (MAPI) based on the IgM data, and show the mean and its 95% CI in text inside each panel (see SI for details). **Panel B.** We plot how the MAPI depends on the rate of waning of sterilizing immunity (*ω*, on the y-axis) and the relative transmissibility of reinfections (compared with primary infections) (*ρ*, on the x-axis). The MAPI was calculated from the equilibrium dynamics of the model shown in Fig 1 and SI Eqns 3–9 with a plausible basic reproductive number (*R* = 5) and 0 *< ω* 2 (i.e., sterilizing immunity is as short as *≥* 0.5 yr) and 0 *< ρ* 1 (i.e., reinfections transmit less than primary infections). See SI Sec 2.1 for details. The white band in Panel B indicates the plausible combinations of values of *ρ* and *ω* consistent with the MAPI for HCoVs estimated in Panel A. If additionally we restrict the duration of sterilizing immunity to at most two years (*ω >* 0.5, plot area above the black line in Panel B), then the plausible range for the relative transmissibility of reinfections is 0.14 *< ρ* 0.23. Reinfected individuals cause between 54 and 68% of primary infections. (See SI Fig 1 for parallel figures calculated at the extremes of plausible values for *R* as defined in Table 1 (i.e., *R* = 2 and *R* = 10).) View this table: [Table 1:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/T1) Table 1: Characteristics of coronavirus-immune interactions and relevant parameter ranges ## 3.2 Transition from pandemic to endemic dynamics At the beginning of an outbreak, the age distribution of cases mirrors that of the population (Fig 3A). Once the demographics of infection reaches a steady state, however primary cases occur almost entirely in babies and young children, who in the case of COVID-19, experience a low CFR. Reinfections in older individuals are predicted to be common and contribute to transmission, but in this steady-state population, older individuals, who would be at risk for severe disease from a *primary* infection, have acquired disease-reducing immunity following infection during childhood. ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F3) Figure 3: Simulations for the transition from epidemic to endemic dynamics for emerging CoVs. Simulations from an extension of the model presented in Fig 1 that includes age structure. Demographic characteristics such as the population’s age distribution, birth, and age-specific death rates are taken from the USA, and seasonality is incorporated via a sinusoidal forcing function (see SI Sec 2.2). **Panel A**. The black line representing the number of infections per day shows an initial peak and subsequent transition to endemic dynamics with much lower numbers of infections (with the endemic dynamics from years 5 to 10 shown in the inset). The three plots in A show that a higher *R* results in a larger and faster initial epidemic and more rapid transition to endemic dynamics. We also show (on the same plot) how the proportion of primary cases in different age groups changes over time (plotted in different colors), and how the transition from epidemic to endemic dynamics results in primary cases being restricted to the lower age groups. Parameters for simulations: *ω* = 1 and *ρ* = 0.7. In **Panel B**, we plot how long it takes for the average CFR (calculated as a 6 month moving average) to fall to 0.001, the CFR associated with seasonal influenza. Grey areas represent simulations where the CFR did not reach 0.001 within ten years. We see that the time taken for the CFR to decline to that of influenza decreases as the transmissibility (*R*) increases and the duration of sterilizing immunity becomes shorter. Results are shown for *ρ* = 0.7. See SI Sec 2.3 and SI Figs 3–5 for sensitivity analyses and model specifications. The overall CFR for CoV-2 is predicted to drop dramatically, eventually falling below that of seasonal influenza (approximately 0.001, see Fig 4 B top panel) due to the age-structure of primary infection. The time it takes to complete the shift in CFR en route to endemicity depends on both transmission (*R*) and loss of immunity (*ω* and *ρ*) (Fig 3B and SI Fig 1). ![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F4) Figure 4: The age dependence of the CFR determines how the overall CFR changes during the transition from epidemic to endemic dynamics for emerging CoVs. **Panel A** shows the age dependence of the CFRs for the three emerging CoVs. CoV-1 and CoV-2 have a J shaped profile, with a monotonic increase in CFR with age. In contrast, the age dependence of the CFR for MERS is U shaped, with high mortality in the younger as well as older age groups. Details of the statistical smoothing are described in SI Sec 5. **Panel B** shows how the overall CFR (computed from model-predicted infections in each age group multiplied by the age-specific CFRs for the different CoVs) changes with the transition from epidemic to endemic dynamics. We see the overall CFR declines as CoV-2 transitions to endemicity, and we predict a similar trend should a CoV-1 (SARS)-like viurs spread widely. In contrast, the model predicts the CFR for a MERS-like virus would increase with time; this increase is a consequence of the shape of the age-dependence of the CFR for MERS. In the absence of vaccination or treatment, high transmission is predicted to lead to a high case load and death rate in the initial epidemic phase during the first few years following emergence (Fig 3 and SI Figs 3 & 4). The total number of deaths is largely independent of *R* (see SI Fig 6). Slowing down the epidemic through social distancing measures that lower *R* flattens the curve, thus delaying infections and preventing the majority of deaths from happening early on, affording critical time for the development of an effective vaccine or treatment. However, as the time scales of disease spread and waning of immunity meet, transmission from reinfected individuals (likely asymptomatic or mild) to the vulnerable is predicted to increasingly drive excess deaths (SI Fig 5). Our prediction of mild CoV-2-associated pathology in the long run stems from the low CFR for primary infections in young children. Were a more pathogenic emerging coronavirus to become endemic, for example one with an age-dependent CFR similar to CoV-1, we still expect a low disease burden in the long run (Fig 4). However, data suggest not all emerging HCoVs follow this optimistic pattern; the overall CFR of an endemic MERS-like virus would increase during the transition to endemicity because disease severity is high in children, the age group expected to experience the bulk of primary cases during the endemic phase. ## 4 Discussion The key result from our new model framework that explicitly recognizes that functional immunity to reinfection, disease and shedding are different is that, in contrast with infections that are severe in childhood, CoV-2 could join the ranks of mild, cold-causing endemic human coronaviruses in the long run. A critical prediction is that the severity of emergent CoVs once they reach endemicity depends only on the severity of primary infection in children (Fig 4) because all available evidence suggests immunity to HCoVs has short IE*S*, moderate IE*I* but strong IE*P* such that childhood infection provides protection from pathology upon reinfection in adulthood. Strain-specific virulence factors, such as the shared cellular receptor, ACE-2, to which CoV-1, CoV-2 and the endemic strain NL63 all bind [18, 19, 20, 21], may affect the CFR during the emergence phase but have little impact on the severity of disease in the endemic phase. Because the four endemic HCoVs have been globally circulating for a long time and almost everyone is infected at a young age, we cannot ascertain how much pathology would result from a primary case of any of these in an elderly person. The key insights come from how our model explicitly incorporates different components of immunological protection with respect to susceptibility, pathology and infectivity (IE*S*, IE*P* and IE*I*) and their different rates of waning. In our analysis we hypothesized that the rates for CoV-2 are comparable to endemic HCoVs. Testing this assumption for CoV-2 will follow from ongoing and future longitudinal studies that combine clinical observations with measures of virus titers, virus shedding and antibody and T cell immunity following primary, secondary and subsequent infections of individuals. In our projections, we make the simplifying assumption that immune memory following a single infection generates protection from pathology (high IE*P*). However, we expect our results to be robust to a more gradual build-up of immunity in childhood when disease is mild. The details of the change in overall CFR through the transient period will be impacted by a wide array of factors, such as age-specific human contact rates [22] and susceptibility to infection [23], as well as improvement in treatment protocols, hospital capacity, and virus evolution. The qualitative result of mild disease in the endemic phase is robust to these complexities, but quantitative predictions for the transient phase will depend on a careful consideration of these realities and how they interact with the dynamics of infection and components of immunity. The changes in the CFR over time predicted by the model have implications for vaccination strategy against current and future emerging HCoVs. The widespread (pandemic) circulation of any novel pathogenic HCoV will invariably require rapid generation and deployment of a vaccine. In the longer term, however, the necessity for continual vaccination will depend on the age-dependence of the CFR. If primary infections of children are mild (CoV-1 and CoV-2), continued vaccination may not be needed as primary cases recede to mild childhood sniffles. If, on the other hand, primary infection is severe in children (as for MERS), then vaccination of children will need to be continued. Our model provides strategic insights into plausible transitions from pandemic to endemic dynamics for CoV-2, which have implications for control measures during a likely high-CFR transient period. These insights are outlined in Table 2. The findings presented here suggest that using symptoms as a surveillance tool to curb the virus’s spread will become more difficult, as milder reinfections increasingly contribute to chains of transmission and population level attack rates. In addition, infection or vaccination may protect against disease but not provide the type of transmission-blocking immunity that allows for shielding [24] or the generation of long-term herd immunity [2]. View this table: [Table 2:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/T2) Table 2: Effects of separated IEs model on control strategies for SARS-CoV-2 From an ecological and evolutionary perspective, our study opens the door to questions regarding the within-host and between-host dynamics of human immunity and pathogen populations in the face of IE’s with different kinetics. It also opens the question of how these IE’s interplay with strain cross-immunity, which is likely relevant within the alpha- and beta-coronaviruses. Considering data and model predictions from emergence through endemicity of HCoVs revealed a framework for understanding immunity and vaccination that may apply to a variety of infections, such as RSV and seasonal influenza, which share similar age distributions and immune responses. ## Data Availability Analysis involves data from literature. It is available in the supplementary information ## Supplement ### 1 Serostudy analysis details We used data from a 2013 study [1] on the seroprevalence of IgG and IgM antibodies against the Spike (S) protein of the four endemic HCoV strains in a cross-section of the population in Beijing, none of whom exhibited symptoms. A striking feature of these data is that IgM titers are undetectable in all of the more than 140 subjects ages 15 yr and older. This observation combined with general knowledge of the the IgM response [2] suggests that S-specific IgM is only elicited during primary infection. Additionally, IgM titers decay quickly, which makes IgM seropositivity a useful marker of recent primary infection. We calculate error bars on the seroprevalence for both IgM and IgG and estimate the mean age of primary infection (MAPI) for each of the four strains from the IgM data. We assume the data for both IgM and IgG stem from a binomial process where the probability of seropositivity in age group *j* is *pj*, estimated by ^*pj* and the sample size for each age group is *Nj*. The 95% confidence interval around the mean proportion seropositive for each age group is then ![Formula][1] We further estimate the MAPI using only the IgM data (we assume that the cases are uniformly distributed within each age group – a more accurate estimate could be obtained from the raw data with smaller age bins; unfortunately we were unable to gain access to it). For each HCoV strain, *s*, we create a vector, *As* with length *L* containing the ages ages of IgM seropositivity (using the midpoints of each age range). We assume this is a reasonable reflection of ages of first infection, as IgM titers increase only during primary infection and decay in a matter of weeks [3]. The 95% CI for the MAPI for each strain, *s*, is therefore estimated by ![Formula][2] ### 2 Model derivation We derived the model presented in Fig 1 by combining three key sources of information: (1) classic SIRS disease transmission models, (2) separated functional immune efficacies (IE’s, [4]), and (3) a human reinfection experiment with HCoV 229e in fifteen healthy adult volunteers [5]. ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F5.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F5) Figure 1: Mean age of primary infection In aforementioned reinfection experiment, all subjects had serum specific antibodies at the start, suggesting that participants had already experienced a primary infection. The first experimental exposure resulted in viral replication and a boosting of IgG titers in ten of the fifteen participants; the group that supported viral replication had lower serum specific IgG, IgA and nasal IgA levels prior to exposure and shed virus for on average 5.6 days; eight out of ten had cold symptoms. Antibody titers increased significantly approximately a week after infection and then slowly decayed over the course of the next year. Among the five participants who did not get infected following exposure at the beginning of the trial, antibody titers remained relatively constant following exposure; however their IgG titers did drop significantly by the end of the year. One year after the initial exposure, fourteen of the fifteen participants were re-exposed. All five who had not sustained an infection the first time became infected (their *IE*S waned) and one developed cold symptoms (*IE*P was still strong for four of the five). Six of the nine who had supported viral replication a year earlier sustained an infection following exposure the second time (their *IE*S waned within the span of year) but none developed cold symptoms (they retained *IE*P); the other three did not become infected following exposure. Among all who got infected at the second time point, the mean duration of viral shedding was only 2.0 days, suggesting that *IE*I had not waned completely. Based on these observations, the following are the basic equations that correspond to Fig 1. ![Formula][3] ![Formula][4] ![Formula][5] ![Formula][6] ![Formula][7] #### 2.1 Steady-state analysis The above equations are used for steady state analysis of the predicted mean ages of primary infection (MAPIs). The model-predicted MAPI for a given set of parameters is calculated as the waiting time from birth to first infection according to the following equation: ![Formula][8] where *Îi* is the equilibrium proportion of the population in class *Ii*. The equilibrium values are calculated by first running a short (100 iterations = 1/10 yr) numerical simulation using a wrapper for the R function lsoda and then using the final values as estimates to start the Newton-Raphson method to find equilibria as implemented in the R package rootSolve, function stode. We additionally calculate the proportion of cases caused by reinfections as follows: ![Formula][9] In addition to the results shown in the main text for *R* = 5, we here show figures parallel to Fig 3b for *R* = 2 and *R* = 10 (SI Fig 1). ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F6.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F6) Figure 2: Age specific death rate in the US ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F7.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F7) Figure 3: Time to CFR = 0.001, parallel to main text Fig 3b #### 2.2 Transient to endemic dynamic simulations To incorporate seasonality with a peak in early January (modeling on influenza and seasonal coronaviruses) and the introduction of the virus in early March (as was approximately observed with CoV-2 in the US), cases are introduced at t = 0 and we allow *β* to fluctuate annually according to ![Formula][10] *β* is the mean value of *β*, and *β*1 is the amplitude of the sin wave. All results shown here use *β*1 = 0.2 (SI Figs 3–6 and Fig 3 in the main text). Additionally, to incorporate the age-specific case fatality rates, death rates (*δ*), and the current age distribution of the US population, we separate each immune state, *X*, into nine age classes *Xj* with 10-yr widths for the first eight. The aging rates correspond to the width of the age classes; the death rates and age distribution are taken from US data. This yields the following equations: ![Formula][11] ![Formula][12] ![Formula][13] ![Formula][14] ![Formula][15] where ![Graphic][16] is the total population size at time *t*. The birth rate, *µ*, is a vector of length nine, containing the overall population birthrate (based on demographic data) followed by zeros for all subsequent age classes (i.e., people are only born into the youngest age class). The age-specific death rates, *δj*, are fixed at values estimated from demographic data. The aging rates, *λj*, are contained in a vector of length ten: (0, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.04, 0). We fix *γ* at ![Graphic][17] corresponding to an infectious period that lasts on average nine days. *β* is calculated according to ![Graphic][18]. We consider a range of values of *R* (2-10), *ω* (0-2), and *ρ* (0-1). The age-specific death rates were inferred from CDC data Health Statistics [6], and were calculated as follows: ![](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F8/graphic-25.medium.gif) [](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F8/graphic-25) ![](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F8/graphic-26.medium.gif) [](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F8/graphic-26) The ages used in the model prediction for death rates were the midpoints of the age classes (e.g. 5, 15, 25yr, etc). We used the death rate for 85 yr olds as the highest death rate. The initial conditions were set so the population was distributed into the nine susceptible classes according to the age distribution of the US population [7]. One infected individuals was seeded into each age group’s *I*1 class. #### 2.3 Calculating infections and the case fatality rate from simulations The following steps were used to calculate the number of daily infections and 6-month moving average case fatality rate (CFR): 1. Numerically integrate equations with chosen parameters and initial conditions as described above using the R function lsoda with a time step of one day (1/365 yr). 2. Calculate the probability of staying in *I*1 for a timestep of one day given that you’re already there. In the simplest version of the model we consider here, this can be calculated using the cumulative distribution function, *F*: ![Formula][19] ![Figure9](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F9.medium.gif) [Figure9](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F9) ## [1] 0.8948393 3. For each age group, calculate the number of new primary infections in each time step, *Xt* ![Formula][20] ![Formula][21] ![Figure10](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F10.medium.gif) [Figure10](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F10) 4. Calculate the projected number of HCoV-induced deaths in each age group based on the age-specific CFR’s. ![Formula][22] 5. For every 6-month window, calculate the overall CFR beginning six months into the pandemic ![Formula][23] The key result that the overall CFR drops to something akin to seasonal influenza (0.001), is robust across a wide range of values for *R*, *ω* and *ρ* for COVID-like CFRs. We additionally show results for simulations with different initial conditions, in which 2000 infections were seeded into each age groups approximating the situation when the pandemic was under control in the U.S. (i.e., there had been already been many infections and *R* was close to 1). The CFR evolves very similarly to the case where only nine cases are introduced. Additionally, we see that keeping *R* below a threshold value (in these simulations approximately 2, e.g., by social distancing and the use of Personal Protective Equipment) allows us to stop the majority of deaths from happening early on, buying time for the development of an effective vaccine and/or treatment (SI Fig 5). ![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F11.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F11) Figure 4: Effect of initial outbreak on time to CFR = 0.001 ![Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F12.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F12) Figure 5: The time it takes (in yr) for 75% of the initial population size to become infected. For *R* *>* 2 the vast majority of infections occur in the first year. By slowing the epidemic down (i.e., decreasing *R*, infections are delayed buying time for the development of a vaccine and/or treatment. However, as we slow down the epidemic by decreasing *R*, the time scales of disease spread and immune waning meet, and reinfections increasingly drive excess deaths, and the total number of infection-induced deaths within the first decade after emergence is very similar regardless of how fast the infection is spreading (SI Fig 6). ![Figure 6:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F13.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F13) Figure 6: Number of infection-induced deaths in first ten years after emergence. When reinfection occurs quickly and is transmissible (high *ω* and *ρ*), the total number of infection-induced deaths is close to independent of *R*. ### 3 Loss of immunity kernel We estimate the duration of sterilizing immunity based on Callow et al. [5], in which six out of nine people who got infected at the start of the experiment were susceptible to reinfection one year later. We can calculate the mean of *ω* given this for exponentially distributed waning times setting the CDF equal to 2/3. ![Formula][24] ![Figure14](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F14.medium.gif) [Figure14](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F14) Given this, the estimate for *ω* is 1.099 and the mean waiting time is ![Graphic][25]. We can also find the average duration of immunity for normally distributed waning times. We assume that the variance is the same as for the exponential model, ![Graphic][26] and the standard deviation is ![Graphic][27]. ![Figure15](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F15.medium.gif) [Figure15](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F15) Here, we find the estimate for the mean of *ω* to be 0.61/yr, so the mean waiting time is ![Graphic][28]. ### 4 Quarantine analysis To better understand how waning of immunity impacts symptom-based quarantining in the first few years after disease emergence, we consider the following three phases. pexp(1, rate = lambdas) ![Figure 7:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F16.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F16) Figure 7: Estimated duration of immunity given exponentially distributed waning times ### 4.1 Phase 1 The initial epidemic, before control measures are introduced, burns through a proportion of the population, *p*1. At some point, let’s call it *t*, we introduce highly effective quarantining, which breaks chains of transmission. Because the initial epidemic burns through the population quickly, we approximate the time of entry into *R*1 as being uniform. Additionally, since chains of transmission are broken upon the introduction of effective quarantining and contact tracing, we assume no new cases arise after *t* during Phase1. Therefore at *t* (or shortly after): ![Formula][29] ![Formula][30] ![Formula][31] ### 4.2 Phase 2 There is a time period during which the occasional immigrant case enters, but no chains of transmission start. Even if the immigrant with a primary or secondary infection infects a Susceptible person before being quarantined, that *I*1 is immediately quarantined and does not infect others. If an *R*2 person is infected, they are not quarantined but neither can they start a chain of transmission because any infection to S results in quarantine, and a transmission chain all among secondary (or later) cases will not happen because ![Graphic][32]. During this time period, people who have recovered from a primary infection are losing their immunity. During this phase: ![Formula][33] ![Formula][34] ![Formula][35] ![Formula][36] ### 4.3 Phase 3 There is a threshold value of *R*2, let’s call it ![Graphic][37], above which an immigrant case can lead to an outbreak of asymptomatic secondary infections, which acts as a reservoir for symptomatic primary infections. This occurs even in the presence of a strong quarantining program for symptomatic primary cases. If ![Graphic][38], quarantining can continue working indefinitely. However, if ![Graphic][39], immune waning will lead to a sufficient build up of secondary susceptibles (in *R*2) to sustain an outbreak. In this scenario, the ‘safe time’ during which strong quarantining of primary cases can prevent outbreaks is influenced by: how fast immunity wanes (*ω*), how transmissible secondary cases are (*ρβ*) and the proportion of the population that was infected in the initial outbreak (*p*1). ![Figure 8:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F17.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F17) Figure 8: Estimated duration of immunity given normally distributed waning times ![Figure 9:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F18.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F18) Figure 9: Phase 1 transmission During this phase: ![Formula][40] ![Formula][41] ![Formula][42] ![Formula][43] ![Formula][44] ### 4.4 Window of safety Here, we show results for the window of safety during which strict quarantining of primary cases will prevent an immigrant infection from leading to an outbreak for *R* = 6 and *R* = 1.5 to mimic a high but reasonable transmissionrate and a rate that may be more commensurate with the current scenario given social distancing. ![Figure 10:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F19.medium.gif) [Figure 10:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F19) Figure 10: Phase 2 transmission For *R* = 2 we see ![Figure20](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F20.medium.gif) [Figure20](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F20) ![Figure 11:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F21.medium.gif) [Figure 11:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F21) Figure 11: Phase 3 transmission ![Figure 12:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F22.medium.gif) [Figure 12:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F22) Figure 12: Safe time for *R* = 2 For *R* = 6, a smaller *p*1 can support an outbreak (as low as p1 = 0.3). The amount of time it takes for people’s immunity to wane sufficiently has a broader range. ![](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F23/graphic-59.medium.gif) [](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F23/graphic-59) ![](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F23/graphic-60.medium.gif) [](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F23/graphic-60) ![Figure 13:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F24.medium.gif) [Figure 13:](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F24) Figure 13: Safe time for *R* = 6 In reality, many primary cases have also been found to be asymptomatic, making symptom-based quarantine substantially more challenging and highlighting the need for asymptomatic surveillance to protect against transmission to vulnerable individuals [8]. ## 5 Age-severity curves We estimate the age-severity curves for the three HCoVs to have emerged in the past few decades using published data (CoV-2 [9], SARS CoV-1 [10], and MERS [11]). We then fit a generalized linear models to each data set to estimate a smoothed CFR as a function of age. We use a binomial model in which total cases in an age group is considered the number of trials, and the number of deaths is considered the number of ‘successes’. For MERS, we allow the function to be a third degree polynomial to account for the non-monotonicity of the data. We consider a second degree polynomial for SARS CoV-1 and −2, since the relationship between age and CFR appears monotonic in the data for these. Our code is included below. ![Figure25](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F25.medium.gif) [Figure25](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F25) ![Figure26](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/05/2020.09.03.20187856/F26.medium.gif) [Figure26](http://medrxiv.org/content/early/2020/09/05/2020.09.03.20187856/F26) ## 6 Code availability and reproducibility All code will be available as an RMarkdown file. * Received September 3, 2020. * Revision received September 3, 2020. * Accepted September 5, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].Cobey, S. Modeling infectious disease dynamics. Science 368, 713–714 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5Mi83MTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8wNS8yMDIwLjA5LjAzLjIwMTg3ODU2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 2. [2].Randolph, H. E. & Barreiro L. B. Herd immunity: Understanding COVID-19. Immunity 52, 737–741 (2020). URL [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.immuni.2020.04.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32433946&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 3. [3].Fraser, C., Riley, S., Anderson, R. M. & Ferguson, N. M. Factors that make an infectious disease outbreak controllable. Proceedings of the National Academy of Sciences 101, 6146–6151 (2004). URL [http://dx.doi.org/10.1073/pnas.0307506101](http://dx.doi.org/10.1073/pnas.0307506101). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAxLzE2LzYxNDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8wNS8yMDIwLjA5LjAzLjIwMTg3ODU2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 4. [4].Su, S. et al. Epidemiology, genetic recombination, and pathogenesis of coronaviruses. Trends in Microbiology 24, 490–502 (2016). URL [http://dx.doi.org/10.1016/j.tim.2016.03.003](http://dx.doi.org/10.1016/j.tim.2016.03.003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tim.2016.03.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27012512&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 5. [5].Salamatbakhsh, M., Mobaraki, K., Sadeghimohammadi, S. & Ahmadzadeh, J. The global burden of premature mortality due to the middle east respiratory syndrome (mers) using standard expected years of life lost, 2012 to 2019. BMC Public Health 19 (2019). URL [http://dx.doi.org/10.1186/s12889-019-7899-2](http://dx.doi.org/10.1186/s12889-019-7899-2). 6. [6].Ruan, S. Likelihood of survival of coronavirus disease 2019. The Lancet Infectious Diseases 20, 630–631 (2020). URL [http://dx.doi.org/10.1016/S1473-3099(20)30257-7](http://dx.doi.org/10.1016/S1473-3099(20)30257-7). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30257-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32240633&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 7. [7].Halloran, M. E., Longini, I. M. & Struchiner, C. J. Design and analysis of vaccine studies. Statistics for biology and health (Springer, New York, 2010). 8. [8].[8] Callow, K. A., Parry, H. F., Sergeant, M. & Tyrrell, D. A. J The time course of the immune response to experimental coronavirus infection of man. Epidemiology and Infection 105, 435–446 (1990). URL [http://dx.doi.org/10.1017/s0950268800048019](http://dx.doi.org/10.1017/s0950268800048019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268800048019&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2170159&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 9. [9].Bradburne, A. F., Bynoe, M. L. & Tyrrell, D. A. Effects of a “new” human respiratory virus in volunteers. BMJ 3, 767–769 (1967). URL [http://dx.doi.org/10.1136/bmj.3.5568.767](http://dx.doi.org/10.1136/bmj.3.5568.767). [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czozOiJibWoiO3M6NToicmVzaWQiO3M6MTA6IjMvNTU2OC83NjciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8wNS8yMDIwLjA5LjAzLjIwMTg3ODU2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 10. [10].Zhou, W., Wang, W., Wang, H., Lu, R. & Tan, W. First infection by all four non-severe acute respiratory syndrome human coronaviruses takes place during childhood. BMC Infectious Diseases 13 (2013). URL [http://dx.doi.org/10.1186/1471-2334-13-433](http://dx.doi.org/10.1186/1471-2334-13-433). 11. [11].Galanti, M. & Shaman, J. Direct observation of repeated infections with endemic coronaviruses. The Journal of Infectious Diseases (2020). URL [http://dx.doi.org/10.1093/infdis/jiaa392](http://dx.doi.org/10.1093/infdis/jiaa392). 12. [12].He, X. et al. Temporal dynamics in viral shedding and transmissibility of covid-19. Nature Medicine 26, 672–675 (2020). URL [http://dx.doi.org/10.1038/s41591-020-0869-5](http://dx.doi.org/10.1038/s41591-020-0869-5). [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%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 13. [13].Sanche, S. et al. High contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2. Emerging Infectious Diseases 26, 1470–1477 (2020). URL [http://dx.doi.org/10.3201/eid2607.200282](http://dx.doi.org/10.3201/eid2607.200282). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 14. [14].Kissler, S.M., Tedijanto, C., Goldstein, E., Grad, Y.H. & Lipsitch, M. Projecting the transmission dynamics of sars-cov-2 through the postpandemic period. Science 368, 860–868 (2020). URL [http://dx.doi.org/10.1126/science.abb5793](http://dx.doi.org/10.1126/science.abb5793). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5My84NjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8wNS8yMDIwLjA5LjAzLjIwMTg3ODU2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 15. [15].Chan-Yeung, M. & Xu, R.-H. SARS: epidemiology. Respirology 8, S9–S14 (2003). URL [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1440-1843.2003.00518.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15018127&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000187345100004&link_type=ISI) 16. [16].Salamatbakhsh, M., Mobaraki, K., Sadeghimohammadi, S. & Ahmadzadeh, J. The global burden of premature mortality due to the middle east respiratory syndrome (mers) using standard expected years of life lost, 2012 to 2019. BMC Public Health 19, 1523 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12889-019-7899-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31727042&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 17. [17].Verity, R. et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect Dis 20, 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%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 18. [18].Yan, R. et al. Structural basis for the recognition of sars-cov-2 by full-length human ace2. Science 367, 1444–1448 (2020). URL [http://dx.doi.org/10.1126/science.abb2762](http://dx.doi.org/10.1126/science.abb2762). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjcvNjQ4NS8xNDQ0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDUvMjAyMC4wOS4wMy4yMDE4Nzg1Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 19. [19].Kuba, K. et al. A crucial role of angiotensin converting enzyme 2 (ace2) in sars coronavirus–induced lung injury. Nature Medicine 11, 875–879 (2005). URL [http://dx.doi.org/10.1038/nm1267](http://dx.doi.org/10.1038/nm1267). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm1267&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16007097&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230964700033&link_type=ISI) 20. [20].Imai, Y. et al. Angiotensin-converting enzyme 2 protects from severe acute lung failure. Nature 436, 112–116 (2005). URL [http://dx.doi.org/10.1038/nature03712](http://dx.doi.org/10.1038/nature03712). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature03712&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16001071&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230296600049&link_type=ISI) 21. [21].Fung, T.S. & Liu, D.X. Human coronavirus: Host-pathogen interaction. Annual Review of Microbiology 73, 529–557 (2019). URL [http://dx.doi.org/10.1146/annurev-micro-020518-115759](http://dx.doi.org/10.1146/annurev-micro-020518-115759). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-micro-020518-115759&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31226023&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 22. [22].Wallinga, J., Teunis, P. & Kretzschmar, M. Using data on social contacts to estimate age-specific transmission parameters for respiratory-spread infectious agents. American Journal of Epidemiology 164, 936–944 (2006). URL [http://dx.doi.org/10.1093/aje/kwj317](http://dx.doi.org/10.1093/aje/kwj317). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwj317&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16968863&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241958900002&link_type=ISI) 23. [23].Zhang, J. et al. Changes in contact patterns shape the dynamics of the covid-19 outbreak in china. Science 368, 1481–1486 (2020). URL [http://dx.doi.org/10.1126/science.abb8001](http://dx.doi.org/10.1126/science.abb8001). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjgvNjQ5OC8xNDgxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDUvMjAyMC4wOS4wMy4yMDE4Nzg1Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 24. [24].Weitz, J. S. et al. Modeling shield immunity to reduce covid-19 epidemic spread. Nature Medicine (2020). URL [https://doi.org/10.1038/s41591-020-0895-3](https://doi.org/10.1038/s41591-020-0895-3). ## References 1. [1].Zhou, W., Wang, W., Wang, H., Lu, R. & Tan, W. First infection by all four non-severe acute respiratory syndrome human coronaviruses takes place during childhood. BMC Infectious Diseases 13 (2013). URL [http://dx.doi.org/10.1186/1471-2334-13-433](http://dx.doi.org/10.1186/1471-2334-13-433). 2. [2].Schroeder, H. W. & Cavacini, L. Structure and function of immunoglobulins. Journal of Allergy and Clinical Immunology 125, S41–S52 (2010). URL [http://dx.doi.org/10.1016/j.jaci.2009.09.046](http://dx.doi.org/10.1016/j.jaci.2009.09.046). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2009.09.046&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20176268&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280170600005&link_type=ISI) 3. [3].Wu, L.-P. et al. Duration of antibody responses after severe acute respiratory syndrome. Emerging Infectious Diseases 13, 1562–1564 (2007). URL [https://pubmed.ncbi.nlm.nih.gov/18258008/](https://pubmed.ncbi.nlm.nih.gov/18258008/). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1310.070576&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18258008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 4. [4].Halloran, M. E., Longini, I. M. & Struchiner, C. J. Design and analysis of vaccine studies. Statistics for biology and health (Springer, New York, 2010). 5. [5].Callow, K. A., Parry, H. F., Sergeant, M. & Tyrrell, D. A. J. The time course of the immune response to experimental coronavirus infection of man. Epidemiology and Infection 105, 435–446 (1990). URL [http://dx.doi.org/10.1017/s0950268800048019](http://dx.doi.org/10.1017/s0950268800048019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268800048019&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2170159&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 6. [6].for Health Statistics, N. C. Mortality in the united states, 2018. 7. [7].Howden L, M. J. Age and sex composition: 2010: 2010 census briefs. 8. [8].Lee, S. et al. Clinical Course and Molecular Viral Shedding Among Asymptomatic and Symptomatic Patients With SARS-CoV-2 Infection in a Community Treatment Center in the Republic of Korea. JAMA Internal Medicine (2020). URL [https://doi.org/10.1001/jamainternmed.2020.3862](https://doi.org/10.1001/jamainternmed.2020.3862). [https://jamanetwork.com/journals/jamainternalmedicine/articlepdf/2769235/jamainternal\_lee\_2020\_oi\_200057](https://jamanetwork.com/journals/jamainternalmedicine/articlepdf/2769235/jamainternal\_lee_2020_oi_200057). 9. [9].Verity, R. et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect Dis 20, 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%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) 10. [10].Chan-Yeung, M. & Xu, R.-H. SARS: epidemiology. Respirology 8, S9–S14 (2003). URL [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1440-1843.2003.00518.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15018127&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000187345100004&link_type=ISI) 11. [11].Salamatbakhsh, M., Mobaraki, K., Sadeghimohammadi, S. & Ahmadzadeh, J. The global burden of premature mortality due to the middle east respiratory syndrome (mers) using standard expected years of life lost, 2012 to 2019. BMC Public Health 19, 1523 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12889-019-7899-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31727042&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F05%2F2020.09.03.20187856.atom) [1]: /embed/graphic-7.gif [2]: /embed/graphic-8.gif [3]: /embed/graphic-10.gif [4]: /embed/graphic-11.gif [5]: /embed/graphic-12.gif [6]: /embed/graphic-13.gif [7]: /embed/graphic-14.gif [8]: /embed/graphic-15.gif [9]: /embed/graphic-16.gif [10]: /embed/graphic-19.gif [11]: /embed/graphic-20.gif [12]: /embed/graphic-21.gif [13]: /embed/graphic-22.gif [14]: /embed/graphic-23.gif [15]: /embed/graphic-24.gif [16]: /embed/inline-graphic-1.gif [17]: /embed/inline-graphic-2.gif [18]: /embed/inline-graphic-3.gif [19]: /embed/graphic-27.gif [20]: /embed/graphic-29.gif [21]: /embed/graphic-30.gif [22]: /embed/graphic-32.gif [23]: /embed/graphic-33.gif [24]: /embed/graphic-37.gif [25]: /embed/inline-graphic-4.gif [26]: /embed/inline-graphic-5.gif [27]: /embed/inline-graphic-6.gif [28]: /embed/inline-graphic-7.gif [29]: /embed/graphic-41.gif [30]: /embed/graphic-42.gif [31]: /embed/graphic-43.gif [32]: /embed/inline-graphic-8.gif [33]: /embed/graphic-44.gif [34]: /embed/graphic-45.gif [35]: /embed/graphic-46.gif [36]: /embed/graphic-47.gif [37]: /embed/inline-graphic-9.gif [38]: /embed/inline-graphic-10.gif [39]: /embed/inline-graphic-11.gif [40]: /embed/graphic-50.gif [41]: /embed/graphic-51.gif [42]: /embed/graphic-52.gif [43]: /embed/graphic-53.gif [44]: /embed/graphic-54.gif