Age-dependent effects in the transmission and control of COVID-19 epidemics =========================================================================== * Nicholas G. Davies * Petra Klepac * Yang Liu * Kiesha Prem * Mark Jit * CMMID COVID-19 working group * Rosalind M Eggo ## Abstract The COVID-19 pandemic has shown a markedly low proportion of cases among children. Age disparities in observed cases could be explained by children having lower susceptibility to infection, lower propensity to show clinical symptoms, or both. We evaluate these possibilities by fitting an age-structured mathematical model to epidemic data from six countries. We estimate that clinical symptoms occur in 25% (95% CrI: 19-32%) of infections in 10-19-year-olds, rising to 76% (68-82%) in over-70s, and that susceptibility to infection in under-20s is approximately half that of older adults. Accordingly, we find that interventions aimed at children may have a relatively small impact on total cases, particularly if the transmissibility of subclinical infections is low. The age-specific clinical fraction and susceptibility we have estimated has implications for the expected global burden of COVID-19 because of demographic differences across settings: in younger populations, the expected clinical attack rate would be lower, although it is likely that comorbidities in low-income countries will affect disease severity. Without effective control measures, regions with older populations may see disproportionally more clinical cases, particularly in the later stages of the pandemic. View this table: [Extended Data Table 1.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/T1) Extended Data Table 1. Posterior estimates for the consensus susceptibility and clinical fraction from 6 countries. Note that susceptibility is a relative measure. View this table: [Extended Data Table 2.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/T2) Extended Data Table 2. Model parameters. ![Extended Data Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.03.24.20043018/F1.medium.gif) [Extended Data Figure 1.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/F1) Extended Data Figure 1. Prior distributions (grey dotted lines) and posterior distributions (coloured histograms) for model parameters fitting to the early epidemic in Wuhan (Fig. 1, main text); seed_start is measured in days after November 1st, 2019. **(a)** Model 1 (age-varying contact patterns and susceptibility); **(b)** Model 2 (age-varying contact patterns and clinical fraction); **(c)** Model 3 (age-varying contact patterns only). See also Supplementary Table 3. ![Extended Data Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.03.24.20043018/F2.medium.gif) [Extended Data Figure 2.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/F2) Extended Data Figure 2. Analysis showing how the inferred age-varying susceptibility (first column) and age-varying clinical fraction (second column) depend upon the additional data sources used. ![Extended Data Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.03.24.20043018/F3.medium.gif) [Extended Data Figure 3.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/F3) Extended Data Figure 3. Prior and posterior distributions for the epidemics in **(a)** Beijing and Shanghai, **(b)** South Korea and **(c)** Lombardy using the “consensus” fit for age-specific clinical fraction and assuming subclinical infections are 50% as infectious as clinical infections (see Fig. 2c, main text). For **(a)**, times are in days after December 1st, 2019; for **(b)** and **(c)**, times are in days after January 1st, 2019. Note, seed_d is the inferred duration of the seeding event. See also Supplementary Table 3. ![Extended Data Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.03.24.20043018/F4.medium.gif) [Extended Data Figure 4.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/F4) Extended Data Figure 4. Global projections assuming greater severity in lower-income countries. **(a)** Schematic age-specific clinical fraction for higher-income and lower-income countries. **(b-f)** Illustrative results of the projections for 146 capital cities assuming a higher age-varying clinical fraction in lower-income countries. See Fig. 4 (main text) for details. ![Extended Data Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.03.24.20043018/F5.medium.gif) [Extended Data Figure 5.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/F5) Extended Data Figure 5. Consensus age-specific clinical fraction and susceptibility, assuming subclinical infections are 0%, 25%, 50%, 75%, or 100% as infectious as clinical infections. ![Extended Data Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.03.24.20043018/F6.medium.gif) [Extended Data Figure 6.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/F6) Extended Data Figure 6. **(a)** Projected total and peak clinical case attack rate for 146 capital cities, under different assumptions for the infectiousness of subclinical infections. **(b)** Projected total and peak subclinical infection attack rate for 146 capital cities, under different assumptions for the infectiousness of subclinical infections. **(c)** Projected differences in R0 among 146 capital cities, under different assumptions for the infectiousness of subclinical infections. ![Extended Data Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.03.24.20043018/F7.medium.gif) [Extended Data Figure 7.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/F7) Extended Data Figure 7. Comparison of school closures in three exemplar cities when susceptibility *ui* is fixed across settings instead of R0. ![Extended Data Figure 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.03.24.20043018/F8.medium.gif) [Extended Data Figure 8.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/F8) Extended Data Figure 8. Contact matrices used for Figs. 1-3 of the main text. We have not shown matrices for all 12 regions of Italy modelled, nor for all 13 provinces of China modelled, as these show similar patterns to the matrices for Milan and for Wuhan, Beijing and Shanghai, respectively. ## Supplementary Information View this table: [Supplementary Table 1.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/T3) Supplementary Table 1. Details of model fitting. View this table: [Supplementary Table 2.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/T4) Supplementary Table 2. Details on mixing matrices used in the study. View this table: [Supplementary Table 3.](http://medrxiv.org/content/early/2020/05/03/2020.03.24.20043018/T5) Supplementary Table 3. Posterior means and 95% HDIs from fitting the dynamic transmission model (Figs. 1 and 2, main text). ## Data Availability All data used for the analysis are publicly available. Analysis code will be made available on Github. ## Footnotes * 1 Carl A B Pearson, Billy J Quilty, Adam J Kucharski, Hamish Gibbs, Samuel Clifford, Amy Gimma, Kevin van Zandvoort, James D Munday, Charlie Diamond, W John Edmunds, Rein MGJ Houben, Joel Hellewell, Timothy W Russell, Sam Abbott, Sebastian Funk, Nikos I Bosse, Fiona Sun, Stefan Flasche, Alicia Rosello & Christopher I Jarvis. Order of working group determined at random. * Received March 24, 2020. * Revision received April 29, 2020. * Accepted May 3, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)