51
|
Burg D, Ausubel JH. Trajectories of COVID-19: A longitudinal analysis of many nations and subnational regions. PLoS One 2023; 18:e0281224. [PMID: 37352253 PMCID: PMC10289358 DOI: 10.1371/journal.pone.0281224] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 06/07/2023] [Indexed: 06/25/2023] Open
Abstract
The COVID-19 pandemic is the first to be rapidly and sequentially measured by nation-wide PCR community testing for the presence of the viral RNA at a global scale. We take advantage of the novel "natural experiment" where diverse nations and major subnational regions implemented various policies including social distancing and vaccination at different times with different levels of stringency and adherence. Initially, case numbers expand exponentially with doubling times of ~1-2 weeks. In the nations where interventions were not implemented or perhaps lees effectual, case numbers increased exponentially but then stabilized around 102-to-103 new infections (per km2 built-up area per day). Dynamics under effective interventions were perturbed and infections decayed to low levels. They rebounded concomitantly with the lifting of social distancing policies or pharmaceutical efficacy decline, converging on a stable equilibrium setpoint. Here we deploy a mathematical model which captures this V-shape behavior, incorporating a direct measure of intervention efficacy. Importantly, it allows the derivation of a maximal estimate for the basic reproductive number Ro (mean 1.6-1.8). We were able to test this approach by comparing the approximated "herd immunity" to the vaccination coverage observed that corresponded to rapid declines in community infections during 2021. The estimates reported here agree with the observed phenomena. Moreover, the decay (0.4-0.5) and rebound rates (0.2-0.3) were similar throughout the pandemic and among all the nations and regions studied. Finally, a longitudinal analysis comparing multiple national and regional results provides insights on the underlying epidemiology of SARS-CoV-2 and intervention efficacy, as well as evidence for the existence of an endemic steady state of COVID-19.
Collapse
Affiliation(s)
- David Burg
- Tel Hai Academic College, Qiryhat Shemona, Israel
- Hemdat Academic College, Netivot, Israel
- Ahskelon Academic College, Ashkelon, Israel
- Program for the Human Environment, The Rockefeller University, New York, NY, United States of America
| | - Jesse H. Ausubel
- Program for the Human Environment, The Rockefeller University, New York, NY, United States of America
| |
Collapse
|
52
|
Koyama S. Multivariate epidemic count time series model. PLoS One 2023; 18:e0287389. [PMID: 37327242 PMCID: PMC10275427 DOI: 10.1371/journal.pone.0287389] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2022] [Accepted: 06/05/2023] [Indexed: 06/18/2023] Open
Abstract
An infectious disease spreads not only over a single population or community but also across multiple and heterogeneous communities. Moreover, its transmissibility varies over time because of various factors such as seasonality and epidemic control, which results in strongly nonstationary behavior. In conventional methods for assessing transmissibility trends or changes, univariate time-varying reproduction numbers are calculated without taking into account transmission across multiple communities. In this paper, we propose a multivariate-count time series model for epidemics. We also propose a statistical method for estimating the transmission of infections across multiple communities and the time-varying reproduction numbers of each community simultaneously from a multivariate time series of case counts. We apply our method to incidence data for the novel coronavirus disease 2019 (COVID-19) pandemic to reveal the spatiotemporal heterogeneity of the epidemic process.
Collapse
Affiliation(s)
- Shinsuke Koyama
- Department of Statistical Modeling, The Institute of Statistical Mathematics, Tachikawa, Tokyo, Japan
- Department of Statistical Science, Graduate University for Advanced Studies (SOKENDAI), Tachikawa, Tokyo, Japan
| |
Collapse
|
53
|
Belvis F, Aleta A, Padilla-Pozo Á, Pericàs JM, Fernández-Gracia J, Rodríguez JP, Eguíluz VM, De Santana CN, Julià M, Benach J. Key epidemiological indicators and spatial autocorrelation patterns across five waves of COVID-19 in Catalonia. Sci Rep 2023; 13:9709. [PMID: 37322048 PMCID: PMC10272129 DOI: 10.1038/s41598-023-36169-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Accepted: 05/26/2023] [Indexed: 06/17/2023] Open
Abstract
This research studies the evolution of COVID-19 crude incident rates, effective reproduction number R(t) and their relationship with incidence spatial autocorrelation patterns in the 19 months following the disease outbreak in Catalonia (Spain). A cross-sectional ecological panel design based on n = 371 health-care geographical units is used. Five general outbreaks are described, systematically preceded by generalized values of R(t) > 1 in the two previous weeks. No clear regularities concerning possible initial focus appear when comparing waves. As for autocorrelation, we identify a wave's baseline pattern in which global Moran's I increases rapidly in the first weeks of the outbreak to descend later. However, some waves significantly depart from the baseline. In the simulations, both baseline pattern and departures can be reproduced when measures aimed at reducing mobility and virus transmissibility are introduced. Spatial autocorrelation is inherently contingent on the outbreak phase and is also substantially modified by external interventions affecting human behavior.
Collapse
Affiliation(s)
- Francesc Belvis
- Research Group on Health Inequalities, Environment, and Employment Conditions (GREDS-EMCONET), Department of Political and Social Sciences, Universitat Pompeu Fabra, 08005, Barcelona, Spain.
- Johns Hopkins University-Universitat Pompeu Fabra Public Policy Center (JHU-UPF PPC), 08005, Barcelona, Spain.
| | - Alberto Aleta
- Institute for Biocomputation and Physics of Complex Systems (BIFI), University of Zaragoza, 50018, Zaragoza, Spain
| | - Álvaro Padilla-Pozo
- Research Group on Health Inequalities, Environment, and Employment Conditions (GREDS-EMCONET), Department of Political and Social Sciences, Universitat Pompeu Fabra, 08005, Barcelona, Spain
- Johns Hopkins University-Universitat Pompeu Fabra Public Policy Center (JHU-UPF PPC), 08005, Barcelona, Spain
- Department of Sociology, Cornell University, Ithaca, New York, USA
| | - Juan-M Pericàs
- Research Group on Health Inequalities, Environment, and Employment Conditions (GREDS-EMCONET), Department of Political and Social Sciences, Universitat Pompeu Fabra, 08005, Barcelona, Spain
- Johns Hopkins University-Universitat Pompeu Fabra Public Policy Center (JHU-UPF PPC), 08005, Barcelona, Spain
- Liver Unit, Internal Medicine Department, Vall d'Hebron University Hospital, Vall d'Hebron Institute for Research, CIBERehd, 08035, Barcelona, Spain
- Infectious Disease Department, Hospital Clínic, 08036, Barcelona, Spain
| | - Juan Fernández-Gracia
- Instituto de Física Interdisciplinar Y Sistemas Complejos IFISC (CSIC-UIB), 07122, Palma de Mallorca, Spain
| | - Jorge P Rodríguez
- Instituto de Física Interdisciplinar Y Sistemas Complejos IFISC (CSIC-UIB), 07122, Palma de Mallorca, Spain
- Instituto Mediterráneo de Estudios Avanzados IMEDEA (CSIC-UIB), 07190, Esporles, Spain
| | - Víctor M Eguíluz
- Instituto de Física Interdisciplinar Y Sistemas Complejos IFISC (CSIC-UIB), 07122, Palma de Mallorca, Spain
| | - Charles Novaes De Santana
- Instituto de Física Interdisciplinar Y Sistemas Complejos IFISC (CSIC-UIB), 07122, Palma de Mallorca, Spain
| | - Mireia Julià
- Research Group on Health Inequalities, Environment, and Employment Conditions (GREDS-EMCONET), Department of Political and Social Sciences, Universitat Pompeu Fabra, 08005, Barcelona, Spain
- Johns Hopkins University-Universitat Pompeu Fabra Public Policy Center (JHU-UPF PPC), 08005, Barcelona, Spain
- ESIMar (Mar Nursing School), Parc de Salut Mar, Universitat Pompeu Fabra-Affiliated, 08003, Barcelona, Spain
- SDHEd (Social Determinants and Health Education Research Group), IMIM (Hospital del Mar Medical Research Institute), 08005, Barcelona, Spain
| | - Joan Benach
- Research Group on Health Inequalities, Environment, and Employment Conditions (GREDS-EMCONET), Department of Political and Social Sciences, Universitat Pompeu Fabra, 08005, Barcelona, Spain
- Johns Hopkins University-Universitat Pompeu Fabra Public Policy Center (JHU-UPF PPC), 08005, Barcelona, Spain
- Ecological Humanities Research Group (GHECO), Universidad Autónoma de Madrid, 28049, Madrid, Spain
| |
Collapse
|
54
|
Fornace KM, Topazian HM, Routledge I, Asyraf S, Jelip J, Lindblade KA, Jeffree MS, Ruiz Cuenca P, Bhatt S, Ahmed K, Ghani AC, Drakeley C. No evidence of sustained nonzoonotic Plasmodium knowlesi transmission in Malaysia from modelling malaria case data. Nat Commun 2023; 14:2945. [PMID: 37263994 PMCID: PMC10235043 DOI: 10.1038/s41467-023-38476-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Accepted: 05/02/2023] [Indexed: 06/03/2023] Open
Abstract
Reported incidence of the zoonotic malaria Plasmodium knowlesi has markedly increased across Southeast Asia and threatens malaria elimination. Nonzoonotic transmission of P. knowlesi has been experimentally demonstrated, but it remains unknown whether nonzoonotic transmission is contributing to increases in P. knowlesi cases. Here, we adapt model-based inference methods to estimate RC, individual case reproductive numbers, for P. knowlesi, P. falciparum and P. vivax human cases in Malaysia from 2012-2020 (n = 32,635). Best fitting models for P. knowlesi showed subcritical transmission (RC < 1) consistent with a large reservoir of unobserved infection sources, indicating P. knowlesi remains a primarily zoonotic infection. In contrast, sustained transmission (RC > 1) was estimated historically for P. falciparum and P. vivax, with declines in RC estimates observed over time consistent with local elimination. Together, this suggests sustained nonzoonotic P. knowlesi transmission is highly unlikely and that new approaches are urgently needed to control spillover risks.
Collapse
Affiliation(s)
- Kimberly M Fornace
- School of Biodiversity, One Health and Veterinary Medicine, University of Glasgow, Glasgow, UK.
- Saw Swee Hock School of Public Health, National University of, Singapore, Singapore.
- Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK.
| | - Hillary M Topazian
- MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, UK
| | - Isobel Routledge
- MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, UK
- University of California, San Francisco, San Francisco, USA
| | - Syafie Asyraf
- Faculty of Medicine and Health Sciences, Universiti Malaysia Sabah, Kota Kinabalu, Malaysia
| | - Jenarun Jelip
- Vector-borne Disease Control Division, Ministry of Health Malaysia, Putrajaya, Malaysia
| | - Kim A Lindblade
- Global Malaria Programme, World Health Organization, Geneva, Switzerland
| | | | - Pablo Ruiz Cuenca
- Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
| | - Samir Bhatt
- MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, UK
- Section of Epidemiology, University of Copenhagen, Copenhagen, Denmark
| | - Kamruddin Ahmed
- Faculty of Medicine and Health Sciences, Universiti Malaysia Sabah, Kota Kinabalu, Malaysia
| | - Azra C Ghani
- MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, UK
| | - Chris Drakeley
- Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
| |
Collapse
|
55
|
Romanyukha AA, Karkach AS, Borisov SE, Belilovsky EM, Sannikova TE. Identification of growing tuberculosis incidence clusters in a region with a decrease in tuberculosis prevalence in Moscow (2000-2019). J Glob Health 2023; 13:04052. [PMID: 37224511 DOI: 10.7189/jogh.13.04052] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023] Open
Abstract
Background The control of tuberculosis (TB) may benefit from a prospective identification of areas where the incidence may increase in addition to the traditionally identified foci of high incidence. We aimed to identify residential areas with growing tuberculosis incidence rates and assess their significance and stability. Methods We analysed the changes in TB incidence rates using case data georeferenced with spatial granularity to apartment buildings in the territory of Moscow from 2000 to 2019. We identified sparsely distributed areas with significant increases in the incidence rate inside residential areas. We tested the stability of found growth areas to case underreporting via stochastic modelling. Results For 21 350 cases with smear- or culture-positive pulmonary TB among residents from 2000 to 2019, we identified 52 small-scale clusters of growing incidence rate responsible for 1% of all registered cases. We tested clusters of disease growth for underreporting and found them to be relatively unstable to resampling with case drop-out, but their spatial displacement was small. Territories with a stable increase in TB incidence rate were identified and compared to the rest of the city, which is characterised by a significant decrease in incidence. Conclusions Identified areas with a tendency for an increase in the TB incidence rate may be important targets for disease control services.
Collapse
Affiliation(s)
- Alexei A Romanyukha
- Marchuk Institute of Numerical Mathematics, Russian Academy of Sciences, Moscow, Russia
- Moscow State University, Moscow, Russia
| | - Arseny S Karkach
- Marchuk Institute of Numerical Mathematics, Russian Academy of Sciences, Moscow, Russia
| | - Sergey E Borisov
- Moscow Research and Clinical Center for Tuberculosis Control, Moscow Department of Public Health, Moscow, Russia
| | - Evgeny M Belilovsky
- Moscow Research and Clinical Center for Tuberculosis Control, Moscow Department of Public Health, Moscow, Russia
| | - Tatiana E Sannikova
- Marchuk Institute of Numerical Mathematics, Russian Academy of Sciences, Moscow, Russia
| |
Collapse
|
56
|
Toh KB, Runge M, Richardson RA, Hladish TJ, Gerardin J. Design of effective outpatient sentinel surveillance for COVID-19 decision-making: a modeling study. BMC Infect Dis 2023; 23:287. [PMID: 37142984 PMCID: PMC10158704 DOI: 10.1186/s12879-023-08261-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Accepted: 04/17/2023] [Indexed: 05/06/2023] Open
Abstract
BACKGROUND Decision-makers impose COVID-19 mitigations based on public health indicators such as reported cases, which are sensitive to fluctuations in supply and demand for diagnostic testing, and hospital admissions, which lag infections by up to two weeks. Imposing mitigations too early has unnecessary economic costs while imposing too late leads to uncontrolled epidemics with unnecessary cases and deaths. Sentinel surveillance of recently-symptomatic individuals in outpatient testing sites may overcome biases and lags in conventional indicators, but the minimal outpatient sentinel surveillance system needed for reliable trend estimation remains unknown. METHODS We used a stochastic, compartmental transmission model to evaluate the performance of various surveillance indicators at reliably triggering an alarm in response to, but not before, a step increase in transmission of SARS-CoV-2. The surveillance indicators included hospital admissions, hospital occupancy, and sentinel cases with varying levels of sampling effort capturing 5, 10, 20, 50, or 100% of incident mild cases. We tested 3 levels of transmission increase, 3 population sizes, and conditions of either simultaneous transmission increase or lagged increase in the older population. We compared the indicators' performance at triggering alarm soon after, but not prior, to the transmission increase. RESULTS Compared to surveillance based on hospital admissions, outpatient sentinel surveillance that captured at least 20% of incident mild cases could trigger an alarm 2 to 5 days earlier for a mild increase in transmission and 6 days earlier for a moderate or strong increase. Sentinel surveillance triggered fewer false alarms and averted more deaths per day spent in mitigation. When transmission increase in older populations lagged the increase in younger populations by 14 days, sentinel surveillance extended its lead time over hospital admissions by an additional 2 days. CONCLUSIONS Sentinel surveillance of mild symptomatic cases can provide more timely and reliable information on changes in transmission to inform decision-makers in an epidemic like COVID-19.
Collapse
Affiliation(s)
- Kok Ben Toh
- Department of Preventive Medicine, Institute for Global Health, Northwestern University Feinberg School of Medicine, Chicago, IL, USA.
| | - Manuela Runge
- Department of Preventive Medicine, Institute for Global Health, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
| | - Reese Ak Richardson
- Department of Chemical and Biological Engineering, Northwestern University, Chicago, IL, USA
| | - Thomas J Hladish
- Department of Biology, University of Florida, Gainesville, FL, USA
- Emerging Pathogen Institute, University of Florida, Gainesville, FL, USA
| | - Jaline Gerardin
- Department of Preventive Medicine, Institute for Global Health, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
| |
Collapse
|
57
|
Lison A, Banholzer N, Sharma M, Mindermann S, Unwin HJT, Mishra S, Stadler T, Bhatt S, Ferguson NM, Brauner J, Vach W. Effectiveness assessment of non-pharmaceutical interventions: lessons learned from the COVID-19 pandemic. Lancet Public Health 2023; 8:e311-e317. [PMID: 36965985 PMCID: PMC10036127 DOI: 10.1016/s2468-2667(23)00046-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Revised: 02/07/2023] [Accepted: 02/08/2023] [Indexed: 03/27/2023]
Abstract
Effectiveness of non-pharmaceutical interventions (NPIs), such as school closures and stay-at-home orders, during the COVID-19 pandemic has been assessed in many studies. Such assessments can inform public health policies and contribute to evidence-based choices of NPIs during subsequent waves or future epidemics. However, methodological issues and no standardised assessment practices have restricted the practical value of the existing evidence. Here, we present and discuss lessons learned from the COVID-19 pandemic and make recommendations for standardising and improving assessment, data collection, and modelling. These recommendations could contribute to reliable and policy-relevant assessments of the effectiveness of NPIs during future epidemics.
Collapse
Affiliation(s)
- Adrian Lison
- Department of Biosystems Science and Engineering, ETH Zurich, Zurich, Switzerland
| | - Nicolas Banholzer
- Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland
| | - Mrinank Sharma
- Department of Statistics, University of Oxford, Oxford, UK; Future of Humanity Institute, University of Oxford, Oxford, UK
| | - Sören Mindermann
- Department of Computer Science, University of Oxford, Oxford, UK
| | - H Juliette T Unwin
- Medical Research Council Centre for Global Infectious Disease Analysis, Jameel Institute, Imperial College London, London, UK
| | - Swapnil Mishra
- Department of Public Health, University of Copenhagen, Copenhagen, Denmark; Saw Swee Hock School of Public Health, National University of Singapore, Singapore
| | - Tanja Stadler
- Department of Biosystems Science and Engineering, ETH Zurich, Zurich, Switzerland
| | - Samir Bhatt
- Medical Research Council Centre for Global Infectious Disease Analysis, Jameel Institute, Imperial College London, London, UK; Department of Public Health, University of Copenhagen, Copenhagen, Denmark.
| | - Neil M Ferguson
- Medical Research Council Centre for Global Infectious Disease Analysis, Jameel Institute, Imperial College London, London, UK
| | - Jan Brauner
- Department of Computer Science, University of Oxford, Oxford, UK; Future of Humanity Institute, University of Oxford, Oxford, UK
| | - Werner Vach
- Basel Academy for Quality and Research in Medicine, Basel, Switzerland; Department of Environmental Sciences, University of Basel, Basel, Switzerland
| |
Collapse
|
58
|
Dai C, Zhou D, Gao B, Wang K. A new method for the joint estimation of instantaneous reproductive number and serial interval during epidemics. PLoS Comput Biol 2023; 19:e1011021. [PMID: 37000844 PMCID: PMC10096265 DOI: 10.1371/journal.pcbi.1011021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Revised: 04/12/2023] [Accepted: 03/09/2023] [Indexed: 04/03/2023] Open
Abstract
Although some methods for estimating the instantaneous reproductive number during epidemics have been developed, the existing frameworks usually require information on the distribution of the serial interval and/or additional contact tracing data. However, in the case of outbreaks of emerging infectious diseases with an unknown natural history or undetermined characteristics, the serial interval and/or contact tracing data are often not available, resulting in inaccurate estimates for this quantity. In the present study, a new framework was specifically designed for joint estimates of the instantaneous reproductive number and serial interval. Concretely, a likelihood function for the two quantities was first introduced. Then, the instantaneous reproductive number and the serial interval were modeled parametrically as a function of time using the interpolation method and a known traditional distribution, respectively. Using the Bayesian information criterion and the Markov Chain Monte Carlo method, we ultimately obtained their estimates and distribution. The simulation study revealed that our estimates of the two quantities were consistent with the ground truth. Seven data sets of historical epidemics were considered and further verified the robust performance of our method. Therefore, to some extent, even if we know only the daily incidence, our method can accurately estimate the instantaneous reproductive number and serial interval to provide crucial information for policymakers to design appropriate prevention and control interventions during epidemics.
Collapse
|
59
|
Paireau J, Charpignon ML, Larrieu S, Calba C, Hozé N, Boëlle PY, Thiebaut R, Prague M, Cauchemez S. Impact of non-pharmaceutical interventions, weather, vaccination, and variants on COVID-19 transmission across departments in France. BMC Infect Dis 2023; 23:190. [PMID: 36997873 PMCID: PMC10061408 DOI: 10.1186/s12879-023-08106-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2022] [Accepted: 02/20/2023] [Indexed: 04/01/2023] Open
Abstract
BACKGROUND Multiple factors shape the temporal dynamics of the COVID-19 pandemic. Quantifying their relative contributions is key to guide future control strategies. Our objective was to disentangle the individual effects of non-pharmaceutical interventions (NPIs), weather, vaccination, and variants of concern (VOC) on local SARS-CoV-2 transmission. METHODS We developed a log-linear model for the weekly reproduction number (R) of hospital admissions in 92 French metropolitan departments. We leveraged (i) the homogeneity in data collection and NPI definitions across departments, (ii) the spatial heterogeneity in the timing of NPIs, and (iii) an extensive observation period (14 months) covering different weather conditions, VOC proportions, and vaccine coverage levels. FINDINGS Three lockdowns reduced R by 72.7% (95% CI 71.3-74.1), 70.4% (69.2-71.6) and 60.7% (56.4-64.5), respectively. Curfews implemented at 6/7 pm and 8/9 pm reduced R by 34.3% (27.9-40.2) and 18.9% (12.04-25.3), respectively. School closures reduced R by only 4.9% (2.0-7.8). We estimated that vaccination of the entire population would have reduced R by 71.7% (56.4-81.6), whereas the emergence of VOC (mainly Alpha during the study period) increased transmission by 44.6% (36.1-53.6) compared with the historical variant. Winter weather conditions (lower temperature and absolute humidity) increased R by 42.2% (37.3-47.3) compared to summer weather conditions. Additionally, we explored counterfactual scenarios (absence of VOC or vaccination) to assess their impact on hospital admissions. INTERPRETATION Our study demonstrates the strong effectiveness of NPIs and vaccination and quantifies the role of weather while adjusting for other confounders. It highlights the importance of retrospective evaluation of interventions to inform future decision-making.
Collapse
Affiliation(s)
- Juliette Paireau
- Mathematical Modelling of Infectious Diseases Unit, Institut Pasteur, Université Paris Cité, CNRS UMR 2000, Paris, France.
- Infectious Diseases Department, Santé Publique France, Saint Maurice, France.
| | - Marie-Laure Charpignon
- Institute for Data, Systems, and Society (IDSS), Cambridge, MA, USA
- Computational Health Informatics Program, Boston Children's Hospital, Boston, MA, USA
- University of Bordeaux, Inria Bordeaux Sud-Ouest, Inserm, Bordeaux Population Health Research Center, SISTM Team, UMR1219, Bordeaux, France
| | - Sophie Larrieu
- Regions Department, Regional Office Nouvelle-Aquitaine, Santé publique France, Bordeaux, France
| | - Clémentine Calba
- Regions Department, Regional Office Provence-Alps-French Riviera and Corsica, Santé Publique France, Marseille, France
| | - Nathanaël Hozé
- Mathematical Modelling of Infectious Diseases Unit, Institut Pasteur, Université Paris Cité, CNRS UMR 2000, Paris, France
| | - Pierre-Yves Boëlle
- INSERM, Sorbonne Université, Institut Pierre Louis d'Epidémiologie et de Santé Publique, Paris, France
| | - Rodolphe Thiebaut
- University of Bordeaux, Inria Bordeaux Sud-Ouest, Inserm, Bordeaux Population Health Research Center, SISTM Team, UMR1219, Bordeaux, France
| | - Mélanie Prague
- University of Bordeaux, Inria Bordeaux Sud-Ouest, Inserm, Bordeaux Population Health Research Center, SISTM Team, UMR1219, Bordeaux, France
| | - Simon Cauchemez
- Mathematical Modelling of Infectious Diseases Unit, Institut Pasteur, Université Paris Cité, CNRS UMR 2000, Paris, France
| |
Collapse
|
60
|
Rotejanaprasert C, Lawson AB, Maude RJ. Spatiotemporal reproduction number with Bayesian model selection for evaluation of emerging infectious disease transmissibility: an application to COVID-19 national surveillance data. BMC Med Res Methodol 2023; 23:62. [PMID: 36915077 PMCID: PMC10010957 DOI: 10.1186/s12874-023-01870-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Accepted: 02/20/2023] [Indexed: 03/16/2023] Open
Abstract
BACKGROUND To control emerging diseases, governments often have to make decisions based on limited evidence. The effective or temporal reproductive number is used to estimate the expected number of new cases caused by an infectious person in a partially susceptible population. While the temporal dynamic is captured in the temporal reproduction number, the dominant approach is currently based on modeling that implicitly treats people within a population as geographically well mixed. METHODS In this study we aimed to develop a generic and robust methodology for estimating spatiotemporal dynamic measures that can be instantaneously computed for each location and time within a Bayesian model selection and averaging framework. A simulation study was conducted to demonstrate robustness of the method. A case study was provided of a real-world application to COVID-19 national surveillance data in Thailand. RESULTS Overall, the proposed method allowed for estimation of different scenarios of reproduction numbers in the simulation study. The model selection chose the true serial interval when included in our study whereas model averaging yielded the weighted outcome which could be less accurate than model selection. In the case study of COVID-19 in Thailand, the best model based on model selection and averaging criteria had a similar trend to real data and was consistent with previously published findings in the country. CONCLUSIONS The method yielded robust estimation in several simulated scenarios of force of transmission with computing flexibility and practical benefits. Thus, this development can be suitable and practically useful for surveillance applications especially for newly emerging diseases. As new outbreak waves continue to develop and the risk changes on both local and global scales, our work can facilitate policymaking for timely disease control.
Collapse
Affiliation(s)
- Chawarat Rotejanaprasert
- Department of Tropical Hygiene, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand.
- Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand.
| | - Andrew B Lawson
- Department of Public Health Sciences, Medical University of South Carolina, Charleston, SC, USA
- Usher Institute, University of Edinburgh, Edinburgh, UK
| | - Richard J Maude
- Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand
- Harvard T.H. Chan School of Public Health, Harvard University, Cambridge, MA, USA
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK
- The Open University, Milton Keynes, UK
| |
Collapse
|
61
|
Kerr GH, Badr HS, Barbieri AF, Colston JM, Gardner LM, Kosek MN, Zaitchik BF. Evolving Drivers of Brazilian SARS-CoV-2 Transmission: A Spatiotemporally Disaggregated Time Series Analysis of Meteorology, Policy, and Human Mobility. GEOHEALTH 2023; 7:e2022GH000727. [PMID: 36960326 PMCID: PMC10030230 DOI: 10.1029/2022gh000727] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Revised: 02/17/2023] [Accepted: 03/06/2023] [Indexed: 06/09/2023]
Abstract
Brazil has been severely affected by the COVID-19 pandemic. Temperature and humidity have been purported as drivers of SARS-CoV-2 transmission, but no consensus has been reached in the literature regarding the relative roles of meteorology, governmental policy, and mobility on transmission in Brazil. We compiled data on meteorology, governmental policy, and mobility in Brazil's 26 states and one federal district from June 2020 to August 2021. Associations between these variables and the time-varying reproductive number (R t ) of SARS-CoV-2 were examined using generalized additive models fit to data from the entire 15-month period and several shorter, 3-month periods. Accumulated local effects and variable importance metrics were calculated to analyze the relationship between input variables and R t . We found that transmission is strongly influenced by unmeasured sources of between-state heterogeneity and the near-recent trajectory of the pandemic. Increased temperature generally was associated with decreased transmission and increased specific humidity with increased transmission. However, the impacts of meteorology, policy, and mobility on R t varied in direction, magnitude, and significance across our study period. This time variance could explain inconsistencies in the published literature to date. While meteorology weakly modulates SARS-CoV-2 transmission, daily or seasonal weather variations alone will not stave off future surges in COVID-19 cases in Brazil. Investigating how the roles of environmental factors and disease control interventions may vary with time should be a deliberate consideration of future research on the drivers of SARS-CoV-2 transmission.
Collapse
Affiliation(s)
- Gaige Hunter Kerr
- Department of Environmental and Occupational HealthGeorge Washington UniversityWashingtonDCUSA
| | - Hamada S. Badr
- Department of Civil and Systems EngineeringJohns Hopkins UniversityBaltimoreMDUSA
- Now at Sales, Market, and Global ServicesAmazon Web ServicesSeattleWAUSA
| | - Alisson F. Barbieri
- Demography DepartmentUniversidade Federal de Minas GeraisBelo HorizonteBrazil
| | - Josh M. Colston
- Division of Infectious Diseases and International HealthUniversity of Virginia School of MedicineCharlottesvilleVAUSA
| | - Lauren M. Gardner
- Department of Civil and Systems EngineeringJohns Hopkins UniversityBaltimoreMDUSA
| | - Margaret N. Kosek
- Division of Infectious Diseases and International HealthUniversity of Virginia School of MedicineCharlottesvilleVAUSA
| | - Benjamin F. Zaitchik
- Department of Earth and Planetary SciencesJohns Hopkins UniversityBaltimoreMDUSA
| |
Collapse
|
62
|
Aravamuthan S, Mandujano Reyes JF, Yandell BS, Döpfer D. Real-time estimation and forecasting of COVID-19 cases and hospitalizations in Wisconsin HERC regions for public health decision making processes. BMC Public Health 2023; 23:359. [PMID: 36803324 PMCID: PMC9937741 DOI: 10.1186/s12889-023-15160-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2022] [Accepted: 01/30/2023] [Indexed: 02/19/2023] Open
Abstract
BACKGROUND The spread of the COVID-19 (SARS-CoV-2) and the surging number of cases across the United States have resulted in full hospitals and exhausted health care workers. Limited availability and questionable reliability of the data make outbreak prediction and resource planning difficult. Any estimates or forecasts are subject to high uncertainty and low accuracy to measure such components. The aim of this study is to apply, automate, and assess a Bayesian time series model for the real-time estimation and forecasting of COVID-19 cases and number of hospitalizations in Wisconsin healthcare emergency readiness coalition (HERC) regions. METHODS This study makes use of the publicly available Wisconsin COVID-19 historical data by county. Cases and effective time-varying reproduction number [Formula: see text] by the HERC region over time are estimated using Bayesian latent variable models. Hospitalizations are estimated by the HERC region over time using a Bayesian regression model. Cases, effective Rt, and hospitalizations are forecasted over a 1-day, 3-day, and 7-day time horizon using the last 28 days of data, and the 20%, 50%, and 90% Bayesian credible intervals of the forecasts are calculated. The frequentist coverage probability is compared to the Bayesian credible level to evaluate performance. RESULTS For cases and effective [Formula: see text], all three time horizons outperform the three credible levels of the forecast. For hospitalizations, all three time horizons outperform the 20% and 50% credible intervals of the forecast. On the contrary, the 1-day and 3-day periods underperform the 90% credible intervals. Questions about uncertainty quantification should be re-calculated using the frequentist coverage probability of the Bayesian credible interval based on observed data for all three metrics. CONCLUSIONS We present an approach to automate the real-time estimation and forecasting of cases and hospitalizations and corresponding uncertainty using publicly available data. The models were able to infer short-term trends consistent with reported values at the HERC region level. Additionally, the models were able to accurately forecast and estimate the uncertainty of the measurements. This study can help identify the most affected regions and major outbreaks in the near future. The workflow can be adapted to other geographic regions, states, and even countries where decision-making processes are supported in real-time by the proposed modeling system.
Collapse
Affiliation(s)
- Srikanth Aravamuthan
- Department of Medical Sciences, University of Wisconsin, Madison, WI, USA. .,Department of Statistics, University of Wisconsin, Madison, WI, USA.
| | - Juan Francisco Mandujano Reyes
- grid.28803.310000 0001 0701 8607Department of Medical Sciences, University of Wisconsin, Madison, WI USA ,grid.28803.310000 0001 0701 8607Department of Statistics, University of Wisconsin, Madison, WI USA
| | - Brian S. Yandell
- grid.28803.310000 0001 0701 8607Department of Statistics, University of Wisconsin, Madison, WI USA
| | - Dörte Döpfer
- grid.28803.310000 0001 0701 8607Department of Medical Sciences, University of Wisconsin, Madison, WI USA
| |
Collapse
|
63
|
Dei Bardi L, Acampora A, Cacciani L, Di Martino M, Agabiti N, Davoli M, Cesaroni G. SARS-CoV-2 spread and area economic disadvantage in the italian three-tier restrictions: a multilevel approach. BMC Public Health 2023; 23:329. [PMID: 36788600 PMCID: PMC9926448 DOI: 10.1186/s12889-023-15246-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Accepted: 02/09/2023] [Indexed: 02/16/2023] Open
Abstract
BACKGROUND To face the second wave of COVID-19, Italy implemented a tiered restriction system with different limitation levels (yellow = medium; orange = medium-high, red = high) at the beginning of November 2020. The restrictions systematically reduced the transmission of SARS-CoV-2 with increasing strength for increasing tier. However, it is unknown whether the effect of limitations was equal between provinces with different socioeconomic levels. Therefore, we investigated the association between the province's socioeconomic level and SARS-CoV-2 infection daily reproduction number in each restriction level. METHODS We measured the province's socioeconomic level as the percentage of individuals whose 2019 total yearly income was lower than 10,000€, using the measure as a proxy of economic disadvantage. We estimated the daily reproduction number (Rt) at the province level using the SARS-CoV-2 daily incidence data from November 2020 to May 2021. We then used multilevel linear regression models with random intercepts stratified by restriction level to estimate the association between economic disadvantage and Rt. We also adjusted the analyses for potential confounders of the association between the province's economic disadvantage and the Rt: the percentage of people with 0-5 years, the quartiles of population density, and the geographical repartition. RESULTS Overall, we found increasing Rt in yellow (+ 0.004 p < 0.01, from Rt = 0.99 to 1.08 in three weeks) and containing effects for the orange (-0.005 p < 0.01, from Rt = 1.03 to 0.93) and the red tier (-0.014 p < 0.01, from Rt = 1.05 to 0.76). More economically disadvantaged provinces had higher Rt levels in every tier, although non-significantly in the yellow level (yellow = 0.001 p = 0.19; orange = 0.002 p = 0.02; red = 0.004 p < 0.01). The results showed that the association between economic disadvantage and Rt differed by level of restriction. The number of days into the restriction and the economic disadvantage had statistically significant interactions in every adjusted model. Compared to better off, more economically disadvantaged provinces had slower increasing trends in yellow and steeper Rt reductions in orange, but they showed slower Rt reductions in the highest tier. CONCLUSION Lower tiers were more effective in more economically disadvantaged provinces, while the highest restriction level had milder effects. These results underline the importance of accounting for socioeconomic level when implementing public health measures.
Collapse
Affiliation(s)
- Luca Dei Bardi
- Department of Epidemiology - Regional Health Service, ASL Roma 1, Rome, Italy. .,Sapienza University of Rome, Rome, Italy.
| | - Anna Acampora
- Department of Epidemiology - Regional Health Service, ASL Roma 1, Rome, Italy
| | - Laura Cacciani
- Department of Epidemiology - Regional Health Service, ASL Roma 1, Rome, Italy
| | - Mirko Di Martino
- Department of Epidemiology - Regional Health Service, ASL Roma 1, Rome, Italy
| | - Nera Agabiti
- Department of Epidemiology - Regional Health Service, ASL Roma 1, Rome, Italy
| | - Marina Davoli
- Department of Epidemiology - Regional Health Service, ASL Roma 1, Rome, Italy
| | - Giulia Cesaroni
- Department of Epidemiology - Regional Health Service, ASL Roma 1, Rome, Italy
| |
Collapse
|
64
|
Bilancia M, Vitale D, Manca F, Perchinunno P, Santacroce L. A dynamic causal modeling of the second outbreak of COVID-19 in Italy. ADVANCES IN STATISTICAL ANALYSIS : ASTA : A JOURNAL OF THE GERMAN STATISTICAL SOCIETY 2023:1-30. [PMID: 36776481 PMCID: PMC9904269 DOI: 10.1007/s10182-023-00469-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Accepted: 01/04/2023] [Indexed: 02/10/2023]
Abstract
While the vaccination campaign against COVID-19 is having its positive impact, we retrospectively analyze the causal impact of some decisions made by the Italian government on the second outbreak of the SARS-CoV-2 pandemic in Italy, when no vaccine was available. First, we analyze the causal impact of reopenings after the first lockdown in 2020. In addition, we also analyze the impact of reopening schools in September 2020. Our results provide an unprecedented opportunity to evaluate the causal relationship between the relaxation of restrictions and the transmission in the community of a highly contagious respiratory virus that causes severe illness in the absence of prophylactic vaccination programs. We present a purely data-analytic approach based on a Bayesian methodology and discuss possible interpretations of the results obtained and implications for policy makers.
Collapse
Affiliation(s)
- Massimo Bilancia
- Department of Precision and Regenerative Medicine and Ionian Area (DiMePRe-J), University of Bari Aldo Moro, Policlinic University Hospital – Piazza G. Cesare 11, 70124 Bari, Italy
| | - Domenico Vitale
- MEMOTEF Department, University of Roma La Sapienza, Via del Castro Laurenziano 9, 00161 Rome, Italy
| | - Fabio Manca
- Department of Education, Psychology, Communication (ForPsiCom), University of Bari Aldo Moro, Palazzo Chiaia Napolitano – Via S. Crisanzio 42, 70122 Bari, Italy
| | - Paola Perchinunno
- Department of Business and Law Studies (DEMDI), University of Bari Aldo Moro, Largo Abbazia di Santa Scolastica 53, 70124 Bari, Italy
| | - Luigi Santacroce
- Department of Interdisciplinary Medicine (DIM) and Microbiology and Virology Unit, University of Bari Aldo Moro, Policlinic University Hospital – Piazza G. Cesare 11, 70124 Bari, Italy
| |
Collapse
|
65
|
Creswell R, Robinson M, Gavaghan D, Parag KV, Lei CL, Lambert B. A Bayesian nonparametric method for detecting rapid changes in disease transmission. J Theor Biol 2023; 558:111351. [PMID: 36379231 DOI: 10.1016/j.jtbi.2022.111351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Revised: 10/11/2022] [Accepted: 11/01/2022] [Indexed: 11/15/2022]
Abstract
Whether an outbreak of infectious disease is likely to grow or dissipate is determined through the time-varying reproduction number, Rt. Real-time or retrospective identification of changes in Rt following the imposition or relaxation of interventions can thus contribute important evidence about disease transmission dynamics which can inform policymaking. Here, we present a method for estimating shifts in Rt within a renewal model framework. Our method, which we call EpiCluster, is a Bayesian nonparametric model based on the Pitman-Yor process. We assume that Rt is piecewise-constant, and the incidence data and priors determine when or whether Rt should change and how many times it should do so throughout the series. We also introduce a prior which induces sparsity over the number of changepoints. Being Bayesian, our approach yields a measure of uncertainty in Rt and its changepoints. EpiCluster is fast, straightforward to use, and we demonstrate that it provides automated detection of rapid changes in transmission, either in real-time or retrospectively, for synthetic data series where the Rt profile is known. We illustrate the practical utility of our method by fitting it to case data of outbreaks of COVID-19 in Australia and Hong Kong, where it finds changepoints coinciding with the imposition of non-pharmaceutical interventions. Bayesian nonparametric methods, such as ours, allow the volume and complexity of the data to dictate the number of parameters required to approximate the process and should find wide application in epidemiology. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".
Collapse
Affiliation(s)
- Richard Creswell
- Department of Computer Science, University of Oxford, Oxford, United Kingdom.
| | - Martin Robinson
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - David Gavaghan
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Kris V Parag
- MRC Centre of Global Infectious Disease Analysis, Jameel Institute for Disease and Emergency Analytics, Imperial College London, London, United Kingdom; NIHR Health Protection Research Unit in Behavioural Science and Evaluation, University of Bristol, Bristol, United Kingdom
| | - Chon Lok Lei
- Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Macau, Macao Special Administrative Region of China
| | - Ben Lambert
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter, United Kingdom.
| |
Collapse
|
66
|
Chen K, Klompmaker JO, Roscoe CJ, Nguyen LH, Drew DA, James P, Laden F, Fecht D, Wang W, Gulliver J, Wolf J, Steves CJ, Spector TD, Chan AT, Hart JE. Associations between greenness and predicted COVID-19-like illness incidence in the United States and the United Kingdom. Environ Epidemiol 2023; 7:e244. [PMID: 36788976 PMCID: PMC9916094 DOI: 10.1097/ee9.0000000000000244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2022] [Accepted: 01/20/2023] [Indexed: 02/10/2023] Open
Abstract
Green spaces may be protective against COVID-19 incidence. They may provide outdoor, ventilated, settings for physical and social activities and therefore decrease transmission risk. We examined the association between neighborhood greenness and COVID-19-like illness incidence using individual-level data. Methods The study population includes participants enrolled in the COVID Symptom Study smartphone application in the United Kingdom and the United States (March-November 2020). All participants were encouraged to report their current health condition and suspected risk factors for COVID-19. We used a validated symptom-based classifier that predicts COVID-19-like illness. We estimated the Normalized Difference Vegetation Index (NDVI), for each participant's reported neighborhood of residence for each month, using images from Landsat 8 (30 m2). We used time-varying Cox proportional hazards models stratified by age, country, and calendar month at study entry and adjusted for the individual- and neighborhood-level risk factors. Results We observed 143,340 cases of predicted COVID-19-like illness among 2,794,029 participants. Neighborhood NDVI was associated with a decreased risk of predicted COVID-19-like illness incidence in the fully adjusted model (hazard ratio = 0.965, 95% confidence interval = 0.960, 0.970, per 0.1 NDVI increase). Stratified analyses showed protective associations among U.K. participants but not among U.S. participants. Associations were slightly stronger for White individuals, for individuals living in rural neighborhoods, and for individuals living in high-income neighborhoods compared to individuals living in low-income neighborhoods. Conclusions Higher levels of greenness may reduce the risk of predicted COVID-19-like illness incidence, but these associations were not observed in all populations.
Collapse
Affiliation(s)
- Kelly Chen
- Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - Jochem O. Klompmaker
- Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts
| | - Charlotte J. Roscoe
- Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - Long H. Nguyen
- Clinical & Translational Epidemiology Unit, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts
| | - David A. Drew
- Clinical & Translational Epidemiology Unit, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts
- Division of Gastroenterology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts
| | - Peter James
- Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
- Department of Population Medicine, Harvard Medical School and Harvard Pilgrim Health Care Institute, Boston, Massachusetts
| | - Francine Laden
- Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts
| | - Daniela Fecht
- MRC Centre for Environment and Health, School of Public Health, Imperial College London, London, United Kingdom
| | - Weiyi Wang
- MRC Centre for Environment and Health, School of Public Health, Imperial College London, London, United Kingdom
| | - John Gulliver
- Centre for Environmental Health and Sustainability, George Davies Centre, University of Leicester, Leicester, United Kingdom
| | | | - Claire J. Steves
- Kings College Department of Twin Research and Genetic Epidemiology, King’s College London, London, United Kingdom
| | - Tim D. Spector
- Kings College Department of Twin Research and Genetic Epidemiology, King’s College London, London, United Kingdom
| | - Andy T. Chan
- Clinical & Translational Epidemiology Unit, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts
- Division of Gastroenterology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts
- Department of Immunology and Infectious Disease, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
| | - Jaime E. Hart
- Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts
| |
Collapse
|
67
|
Sun Q, Miyoshi T, Richard S. Analysis of COVID-19 in Japan with extended SEIR model and ensemble Kalman filter. JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS 2023; 419:114772. [PMID: 36061090 PMCID: PMC9420319 DOI: 10.1016/j.cam.2022.114772] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Revised: 08/08/2022] [Indexed: 06/15/2023]
Abstract
We introduce an extended SEIR infectious disease model with data assimilation for the study of the spread of COVID-19. In this framework, undetected asymptomatic and pre-symptomatic cases are taken into account, and the impact of their uncertain proportion is fully investigated. The standard SEIR model does not consider these populations, while their role in the propagation of the disease is acknowledged. An ensemble Kalman filter is implemented to assimilate reliable observations of three compartments in the model. The system tracks the evolution of the effective reproduction number and estimates the unobservable subpopulations. The analysis is carried out for three main prefectures of Japan and for the entire country of Japan. For these four communities, our estimated effective reproduction numbers are more stable than the corresponding ones estimated by a different method (Toyokeizai). We also perform sensitivity tests for different values of some uncertain medical parameters, like the relative infectivity of symptomatic/asymptomatic cases. The regional analysis results suggest the decreasing efficiency of the states of emergency.
Collapse
Affiliation(s)
- Q Sun
- Data Assimilation Research Team, RIKEN Center for Computational Science (R-CCS), Kobe 650-0047, Japan
- Graduate School of Mathematics, Nagoya University, Nagoya 464-8602, Japan
| | - T Miyoshi
- Data Assimilation Research Team, RIKEN Center for Computational Science (R-CCS), Kobe 650-0047, Japan
- Prediction Science Laboratory, RIKEN Cluster for Pioneering Research (CPR), Kobe 650-0047, Japan
- RIKEN Interdisciplinary Theoretical and Mathematical Sciences Program (iTHEMS), Wako 351-0198, Japan
| | - S Richard
- Data Assimilation Research Team, RIKEN Center for Computational Science (R-CCS), Kobe 650-0047, Japan
- Graduate School of Mathematics, Nagoya University, Nagoya 464-8602, Japan
| |
Collapse
|
68
|
Lustig A, Vattiato G, Maclaren O, Watson LM, Datta S, Plank MJ. Modelling the impact of the Omicron BA.5 subvariant in New Zealand. J R Soc Interface 2023; 20:20220698. [PMID: 36722072 PMCID: PMC9890098 DOI: 10.1098/rsif.2022.0698] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 01/06/2023] [Indexed: 02/02/2023] Open
Abstract
New Zealand experienced a wave of the Omicron variant of SARS-CoV-2 in early 2022, which occurred against a backdrop of high two-dose vaccination rates, ongoing roll-out of boosters and paediatric doses, and negligible levels of prior infection. New Omicron subvariants have subsequently emerged with a significant growth advantage over the previously dominant BA.2. We investigated a mathematical model that included waning of vaccine-derived and infection-derived immunity, as well as the impact of the BA.5 subvariant which began spreading in New Zealand in May 2022. The model was used to provide scenarios to the New Zealand Government with differing levels of BA.5 growth advantage, helping to inform policy response and healthcare system preparedness during the winter period. In all scenarios investigated, the projected peak in new infections during the BA.5 wave was smaller than in the first Omicron wave in March 2022. However, results indicated that the peak hospital occupancy was likely to be higher than in March 2022, primarily due to a shift in the age distribution of infections to older groups. We compare model results with subsequent epidemiological data and show that the model provided a good projection of cases, hospitalizations and deaths during the BA.5 wave.
Collapse
Affiliation(s)
| | - Giorgia Vattiato
- Mathematics and Statistics, University of Canterbury, Christchurch, New Zealand
- Department of Physics, The University of Auckland, Auckland, New Zealand
| | - Oliver Maclaren
- Department of Engineering Science, University of Auckland, Auckland, New Zealand
| | - Leighton M. Watson
- School of Earth and Environment, University of Canterbury, Christchurch, New Zealand
| | - Samik Datta
- National Institute of Water and Atmospheric Research, Wellington, New Zealand
| | - Michael J. Plank
- Mathematics and Statistics, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
69
|
Golding N, Price DJ, Ryan G, McVernon J, McCaw JM, Shearer FM. A modelling approach to estimate the transmissibility of SARS-CoV-2 during periods of high, low, and zero case incidence. eLife 2023; 12:78089. [PMID: 36661303 PMCID: PMC9995112 DOI: 10.7554/elife.78089] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Accepted: 01/16/2023] [Indexed: 01/21/2023] Open
Abstract
Against a backdrop of widespread global transmission, a number of countries have successfully brought large outbreaks of COVID-19 under control and maintained near-elimination status. A key element of epidemic response is the tracking of disease transmissibility in near real-time. During major outbreaks, the effective reproduction number can be estimated from a time-series of case, hospitalisation or death counts. In low or zero incidence settings, knowing the potential for the virus to spread is a response priority. Absence of case data means that this potential cannot be estimated directly. We present a semi-mechanistic modelling framework that draws on time-series of both behavioural data and case data (when disease activity is present) to estimate the transmissibility of SARS-CoV-2 from periods of high to low - or zero - case incidence, with a coherent transition in interpretation across the changing epidemiological situations. Of note, during periods of epidemic activity, our analysis recovers the effective reproduction number, while during periods of low - or zero - case incidence, it provides an estimate of transmission risk. This enables tracking and planning of progress towards the control of large outbreaks, maintenance of virus suppression, and monitoring the risk posed by re-introduction of the virus. We demonstrate the value of our methods by reporting on their use throughout 2020 in Australia, where they have become a central component of the national COVID-19 response.
Collapse
Affiliation(s)
- Nick Golding
- Telethon Kids InstituteNedlandsAustralia
- Curtin UniversityPerthAustralia
| | - David J Price
- Peter Doherty Institute for Infection and Immunity, The Royal Melbourne Hospital and The University of MelbourneVictoriaAustralia
- Melbourne School of Population and Global Health, The University of MelbourneVictoriaAustralia
| | - Gerard Ryan
- Telethon Kids InstituteNedlandsAustralia
- Melbourne School of Population and Global Health, The University of MelbourneVictoriaAustralia
- School of Ecosystem and Forest Sciences, The University of MelbourneVictoriaAustralia
| | - Jodie McVernon
- Peter Doherty Institute for Infection and Immunity, The Royal Melbourne Hospital and The University of MelbourneVictoriaAustralia
- Melbourne School of Population and Global Health, The University of MelbourneVictoriaAustralia
- Murdoch Childrens Research Institute, The Royal Children’s HospitalVictoriaAustralia
| | - James M McCaw
- Peter Doherty Institute for Infection and Immunity, The Royal Melbourne Hospital and The University of MelbourneVictoriaAustralia
- Melbourne School of Population and Global Health, The University of MelbourneVictoriaAustralia
- School of Mathematics and Statistics, The University of MelbourneVictoriaAustralia
| | - Freya M Shearer
- Telethon Kids InstituteNedlandsAustralia
- Melbourne School of Population and Global Health, The University of MelbourneVictoriaAustralia
| |
Collapse
|
70
|
Stolerman LM, Clemente L, Poirier C, Parag KV, Majumder A, Masyn S, Resch B, Santillana M. Using digital traces to build prospective and real-time county-level early warning systems to anticipate COVID-19 outbreaks in the United States. SCIENCE ADVANCES 2023; 9:eabq0199. [PMID: 36652520 PMCID: PMC9848273 DOI: 10.1126/sciadv.abq0199] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Accepted: 12/19/2022] [Indexed: 06/17/2023]
Abstract
Coronavirus disease 2019 (COVID-19) continues to affect the world, and the design of strategies to curb disease outbreaks requires close monitoring of their trajectories. We present machine learning methods that leverage internet-based digital traces to anticipate sharp increases in COVID-19 activity in U.S. counties. In a complementary direction to the efforts led by the Centers for Disease Control and Prevention (CDC), our models are designed to detect the time when an uptrend in COVID-19 activity will occur. Motivated by the need for finer spatial resolution epidemiological insights, we build upon previous efforts conceived at the state level. Our methods-tested in an out-of-sample manner, as events were unfolding, in 97 counties representative of multiple population sizes across the United States-frequently anticipated increases in COVID-19 activity 1 to 6 weeks before local outbreaks, defined when the effective reproduction number Rt becomes larger than 1 for a period of 2 weeks.
Collapse
Affiliation(s)
- Lucas M. Stolerman
- Computational Health Informatics Program, Boston Children’s Hospital, Boston, MA, USA
- Department of Pediatrics, Harvard Medical School, Boston, MA, USA
- Department of Mathematics, Oklahoma State University, Stillwater, OK, USA
| | - Leonardo Clemente
- Computational Health Informatics Program, Boston Children’s Hospital, Boston, MA, USA
- Machine Intelligence Group for the Betterment of Health and the Environment, Network Science Institute, Northeastern University, Boston, MA, USA
| | - Canelle Poirier
- Computational Health Informatics Program, Boston Children’s Hospital, Boston, MA, USA
- Department of Pediatrics, Harvard Medical School, Boston, MA, USA
| | - Kris V. Parag
- NIHR Health Protection Research Unit, Behavioural Science and Evaluation, University of Bristol, Bristol, UK
| | | | - Serge Masyn
- Global Public Health, Janssen R&D, Beerse, Belgium
| | - Bernd Resch
- Department of Geoinformatics - Z-GIS, University of Salzburg, Salzburg, Austria
- Center for Geographic Analysis, Harvard University, Cambridge, MA, USA
| | - Mauricio Santillana
- Department of Pediatrics, Harvard Medical School, Boston, MA, USA
- Machine Intelligence Group for the Betterment of Health and the Environment, Network Science Institute, Northeastern University, Boston, MA, USA
- Harvard University, T.H. Chan School of Public Health, Boston, MA, USA
| |
Collapse
|
71
|
Gressani O, Faes C, Hens N. An approximate Bayesian approach for estimation of the instantaneous reproduction number under misreported epidemic data. Biom J 2023:e2200024. [PMID: 36639234 DOI: 10.1002/bimj.202200024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2022] [Revised: 11/04/2022] [Accepted: 11/18/2022] [Indexed: 01/15/2023]
Abstract
In epidemic models, the effective reproduction number is of central importance to assess the transmission dynamics of an infectious disease and to orient health intervention strategies. Publicly shared data during an outbreak often suffers from two sources of misreporting (underreporting and delay in reporting) that should not be overlooked when estimating epidemiological parameters. The main statistical challenge in models that intrinsically account for a misreporting process lies in the joint estimation of the time-varying reproduction number and the delay/underreporting parameters. Existing Bayesian approaches typically rely on Markov chain Monte Carlo algorithms that are extremely costly from a computational perspective. We propose a much faster alternative based on Laplacian-P-splines (LPS) that combines Bayesian penalized B-splines for flexible and smooth estimation of the instantaneous reproduction number and Laplace approximations to selected posterior distributions for fast computation. Assuming a known generation interval distribution, the incidence at a given calendar time is governed by the epidemic renewal equation and the delay structure is specified through a composite link framework. Laplace approximations to the conditional posterior of the spline vector are obtained from analytical versions of the gradient and Hessian of the log-likelihood, implying a drastic speed-up in the computation of posterior estimates. Furthermore, the proposed LPS approach can be used to obtain point estimates and approximate credible intervals for the delay and reporting probabilities. Simulation of epidemics with different combinations for the underreporting rate and delay structure (one-day, two-day, and weekend delays) show that the proposed LPS methodology delivers fast and accurate estimates outperforming existing methods that do not take into account underreporting and delay patterns. Finally, LPS is illustrated in two real case studies of epidemic outbreaks.
Collapse
Affiliation(s)
- Oswaldo Gressani
- Interuniversity Institute for Biostatistics and Statistical Bioinformatics (I-BioStat), Data Science Institute, Hasselt University, Hasselt, Belgium
| | - Christel Faes
- Interuniversity Institute for Biostatistics and Statistical Bioinformatics (I-BioStat), Data Science Institute, Hasselt University, Hasselt, Belgium
| | - Niel Hens
- Interuniversity Institute for Biostatistics and Statistical Bioinformatics (I-BioStat), Data Science Institute, Hasselt University, Hasselt, Belgium.,Centre for Health Economics Research and Modelling Infectious Diseases, Vaxinfectio, University of Antwerp, Antwerp, Belgium
| |
Collapse
|
72
|
Community Mitigation Strategies, Mobility, and COVID-19 Incidence Across Three Waves in the United States in 2020. Epidemiology 2023; 34:131-139. [PMID: 36137192 PMCID: PMC9811991 DOI: 10.1097/ede.0000000000001553] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
BACKGROUND Summarizing the impact of community-based mitigation strategies and mobility on COVID-19 infections throughout the pandemic is critical for informing responses and future infectious disease outbreaks. Here, we employed time-series analyses to empirically investigate the relationships between mitigation strategies and mobility on COVID-19 incident cases across US states during the first three waves of infections. METHODS We linked data on daily COVID-19 incidence by US state from March to December 2020 with the stringency index, a well-known index capturing the strictness of mitigation strategies, and the trip ratio, which measures the ratio of the number of trips taken per day compared with the same day in 2019. We utilized multilevel models to determine the relative impacts of policy stringency and the trip ratio on COVID-19 cumulative incidence and the effective reproduction number. We stratified analyses by three waves of infections. RESULTS Every five-point increase in the stringency index was associated with 2.89% (95% confidence interval = 1.52, 4.26%) and 5.01% (3.02, 6.95%) reductions in COVID-19 incidence for the first and third waves, respectively. Reducing the number of trips taken by 50% compared with the same time in 2019 was associated with a 16.2% (-0.07, 35.2%) decline in COVID-19 incidence at the state level during the second wave and 19.3% (2.30, 39.0%) during the third wave. CONCLUSIONS Mitigation strategies and reductions in mobility are associated with marked health gains through the reduction of COVID-19 infections, but we estimate variable impacts depending on policy stringency and levels of adherence.
Collapse
|
73
|
Chang JT, Kaplan EH. Modeling local coronavirus outbreaks. EUROPEAN JOURNAL OF OPERATIONAL RESEARCH 2023; 304:57-68. [PMID: 34413569 PMCID: PMC8364218 DOI: 10.1016/j.ejor.2021.07.049] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2021] [Accepted: 07/27/2021] [Indexed: 05/25/2023]
Abstract
This article presents an overview of methods developed for the modeling and control of local coronavirus outbreaks. The article reviews early transmission dynamics featuring exponential growth in infections, and links this to a renewal epidemic model where the current incidence of infection depends upon the expected value of incidence randomly lagged into the past. This leads directly to simple formulas for the fraction of the population infected in an unmitigated outbreak, and reveals herd immunity as the solution to an optimization problem. The model also leads to direct and easy-to-understand formulas for aligning observable epidemic indicators such as cases, hospitalizations and deaths with the unobservable incidence of infection, and as a byproduct leads to a simple first-order approach for estimating the effective reproduction number R t . The model also leads naturally to direct assessments of the effectiveness of isolation in preventing the spread of infection. This is illustrated with application to repeat asymptomatic screening programs of the sort utilized by universities, sports teams and businesses to prevent the spread of infection.
Collapse
Affiliation(s)
- Joseph T Chang
- Department of Statistics and Data Science, Yale University, 24 Hillhouse Avenue, New Haven, CT 06511, USA
| | - Edward H Kaplan
- Yale School of Management, 165 Whitney Avenue, New Haven, CT 06511, USA
| |
Collapse
|
74
|
Tegally H, Khan K, Huber C, de Oliveira T, Kraemer MUG. Shifts in global mobility dictate the synchrony of SARS-CoV-2 epidemic waves. J Travel Med 2022; 29:taac134. [PMID: 36367200 PMCID: PMC9793401 DOI: 10.1093/jtm/taac134] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/11/2022] [Revised: 09/20/2022] [Accepted: 09/23/2022] [Indexed: 11/13/2022]
Abstract
BACKGROUND Human mobility changed in unprecedented ways during the SARS-CoV-2 pandemic. In March and April 2020, when lockdowns and large travel restrictions began in most countries, global air-travel almost entirely halted (92% decrease in commercial global air travel in the months between February and April 2020). Initial recovery in global air travel started around July 2020 and subsequently nearly tripled between May and July 2021. Here, we aim to establish a preliminary link between global mobility patterns and the synchrony of SARS-CoV-2 epidemic waves across the world. METHODS We compare epidemic peaks and human global mobility in two time periods: November 2020 to February 2021 (when just over 70 million passengers travelled) and November 2021 to February 2022 (when more than 200 million passengers travelled). We calculate the time interval during which continental epidemic peaks occurred for both of these time periods, and we calculate the pairwise correlations of epidemic waves between all pairs of countries for the same time periods. RESULTS We find that as air travel increases at the end of 2021, epidemic peaks around the world are more synchronous with one another, both globally and regionally. Continental epidemic peaks occur globally within a 20 day interval at the end of 2021 compared with 73 days at the end of 2020, and epidemic waves globally are more correlated with one another at the end of 2021. CONCLUSIONS This suggests that the rebound in human mobility dictates the synchrony of global and regional epidemic waves. In line with theoretical work, we show that in a more connected world, epidemic dynamics are more synchronized.
Collapse
Affiliation(s)
- Houriiyah Tegally
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
- Centre for Epidemic Response and Innovation (CERI), School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa
| | - Kamran Khan
- BlueDot, Toronto, Canada
- Department of Medicine, Division of Infectious Diseases, University of Toronto, Toronto, Canada
| | | | - Tulio de Oliveira
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
- Centre for Epidemic Response and Innovation (CERI), School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa
- Centre for the AIDS Programme of Research in South Africa (CAPRISA), Durban, South Africa
- Department of Global Health, University of Washington, Seattle, WA, USA
| | - Moritz U G Kraemer
- Department of Biology, University of Oxford, Oxford, United Kingdom
- Pandemic Sciences Institute, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
75
|
Salimović-Bešić I, Dedeić-Ljubović A, Zahirović E, Hasanović M, Šehić M, Vukovikj M, Boshevska G, Vegar-Zubović S, Mehmedika-Suljić E, Izetbegović S. The SARS-CoV-2 Delta (B.1.617.2) variant with spike N501Y mutation in the shadow of Omicron emergence. Heliyon 2022; 8:e12650. [PMID: 36590492 PMCID: PMC9789543 DOI: 10.1016/j.heliyon.2022.e12650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Revised: 10/05/2022] [Accepted: 12/20/2022] [Indexed: 12/25/2022] Open
Abstract
Circulation of the Omicron variant with the reemergence of the N501Y mutation along with many others in the spike protein has once again stirred the academic community. Interestingly, tracing the genetic diversity of SARS-CoV-2 shed light on a less frequent N501Y + Delta variant which has been in the global circulation for some time before the Omicron appearance. This paper aims to present the molecular characteristics of the SARS-CoV-2 Spike_N501Y + Delta variant detected in Bosnia and Herzegovina. The study was conducted during November and December 2021. All patients were tested using real-time RT-PCR for detection of SARS-CoV-2. A representative number of SARS-CoV-2 positive samples was pre-screened using VirSNiP SARS-CoV-2 Spike N501Y kit. The characterization of the viruses was carried out with Illumina RNA Prep with enrichment and the Respiratory Virus Oligo Panel kit. Among the analyzed sequences, we found two isolates of the Delta variant that differ from their most related clade- GK AY.4.3 in additional mutations N501Y and L54F. In this study, we described the presence of a rare form of Delta variant with Spike_N501Y mutation in the shadow of the Omicron emergence. Despite the set of mutations in the Spike protein, this form of Delta variant does not indicate the large-scale consequences for the general population. Further functional studies of this form could provide more information about its antigenicity and infectivity.
Collapse
Affiliation(s)
- Irma Salimović-Bešić
- Clinical Center of the University of Sarajevo, Bolnička 25, 71000, Sarajevo, Bosnia and Herzegovina,Corresponding author.
| | - Amela Dedeić-Ljubović
- Clinical Center of the University of Sarajevo, Bolnička 25, 71000, Sarajevo, Bosnia and Herzegovina
| | - Edina Zahirović
- Clinical Center of the University of Sarajevo, Bolnička 25, 71000, Sarajevo, Bosnia and Herzegovina
| | - Medina Hasanović
- Clinical Center of the University of Sarajevo, Bolnička 25, 71000, Sarajevo, Bosnia and Herzegovina
| | - Merima Šehić
- Clinical Center of the University of Sarajevo, Bolnička 25, 71000, Sarajevo, Bosnia and Herzegovina
| | - Maja Vukovikj
- Institute of Public Health of Republic of North Macedonia, 50-ta Divizija 6, 1000, Skopje, Macedonia
| | - Golubinka Boshevska
- Faculty for Medical Sciences, University Goce Delchev, Krste Misirkov No.10-A P.O. Box 201, 2000, Shtip, Macedonia
| | - Sandra Vegar-Zubović
- Clinical Center of the University of Sarajevo, Bolnička 25, 71000, Sarajevo, Bosnia and Herzegovina
| | - Enra Mehmedika-Suljić
- Clinical Center of the University of Sarajevo, Bolnička 25, 71000, Sarajevo, Bosnia and Herzegovina
| | - Sebija Izetbegović
- Clinical Center of the University of Sarajevo, Bolnička 25, 71000, Sarajevo, Bosnia and Herzegovina
| |
Collapse
|
76
|
Inferring time-varying generation time, serial interval, and incubation period distributions for COVID-19. Nat Commun 2022; 13:7727. [PMID: 36513688 PMCID: PMC9747081 DOI: 10.1038/s41467-022-35496-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Accepted: 12/07/2022] [Indexed: 12/15/2022] Open
Abstract
The generation time distribution, reflecting the time between successive infections in transmission chains, is a key epidemiological parameter for describing COVID-19 transmission dynamics. However, because exact infection times are rarely known, it is often approximated by the serial interval distribution. This approximation holds under the assumption that infectors and infectees share the same incubation period distribution, which may not always be true. We estimated incubation period and serial interval distributions using 629 transmission pairs reconstructed by investigating 2989 confirmed cases in China in January-February 2020, and developed an inferential framework to estimate the generation time distribution that accounts for variation over time due to changes in epidemiology, sampling biases and public health and social measures. We identified substantial reductions over time in the serial interval and generation time distributions. Our proposed method provides more reliable estimation of the temporal variation in the generation time distribution, improving assessment of transmission dynamics.
Collapse
|
77
|
Inward RP, Jackson F, Dasgupta A, Lee G, Battle AL, Parag KV, Kraemer MU. Impact of spatiotemporal heterogeneity in COVID-19 disease surveillance on epidemiological parameters and case growth rates. Epidemics 2022; 41:100627. [PMID: 36099708 PMCID: PMC9443927 DOI: 10.1016/j.epidem.2022.100627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Revised: 08/04/2022] [Accepted: 09/03/2022] [Indexed: 02/08/2023] Open
Abstract
SARS-CoV-2 case data are primary sources for estimating epidemiological parameters and for modelling the dynamics of outbreaks. Understanding biases within case-based data sources used in epidemiological analyses is important as they can detract from the value of these rich datasets. This raises questions of how variations in surveillance can affect the estimation of epidemiological parameters such as the case growth rates. We use standardised line list data of COVID-19 from Argentina, Brazil, Mexico and Colombia to estimate delay distributions of symptom-onset-to-confirmation, -hospitalisation and -death as well as hospitalisation-to-death at high spatial resolutions and throughout time. Using these estimates, we model the biases introduced by the delay from symptom-onset-to-confirmation on national and state level case growth rates (rt) using an adaptation of the Richardson-Lucy deconvolution algorithm. We find significant heterogeneities in the estimation of delay distributions through time and space with delay difference of up to 19 days between epochs at the state level. Further, we find that by changing the spatial scale, estimates of case growth rate can vary by up to 0.13 d-1. Lastly, we find that states with a high variance and/or mean delay in symptom-onset-to-diagnosis also have the largest difference between the rt estimated from raw and deconvolved case counts at the state level. We highlight the importance of high-resolution case-based data in understanding biases in disease reporting and how these biases can be avoided by adjusting case numbers based on empirical delay distributions. Code and openly accessible data to reproduce analyses presented here are available.
Collapse
Affiliation(s)
- Rhys P.D. Inward
- Department of Biology, University of Oxford, United Kingdom,Corresponding author
| | - Felix Jackson
- Department of Biology, University of Oxford, United Kingdom,Department of Computer Science, University of Oxford, United Kingdom
| | - Abhishek Dasgupta
- Department of Biology, University of Oxford, United Kingdom,Department of Computer Science, University of Oxford, United Kingdom
| | - Graham Lee
- Department of Biology, University of Oxford, United Kingdom,Department of Computer Science, University of Oxford, United Kingdom
| | | | - Kris V. Parag
- MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, United Kingdom,NIHR Health Protection Research Unit in Behavioural Science and Evaluation, University of Bristol, Bristol, United Kingdom
| | - Moritz U.G. Kraemer
- Department of Biology, University of Oxford, United Kingdom,Reuben College, University of Oxford, United Kingdom,Corresponding author at: Department of Biology, University of Oxford, United Kingdom
| | | |
Collapse
|
78
|
Hotton AL, Ozik J, Kaligotla C, Collier N, Stevens A, Khanna AS, MacDonell MM, Wang C, LePoire DJ, Chang YS, Martinez-Moyano IJ, Mucenic B, Pollack HA, Schneider JA, Macal C. Impact of changes in protective behaviors and out-of-household activities by age on COVID-19 transmission and hospitalization in Chicago, Illinois. Ann Epidemiol 2022; 76:165-173. [PMID: 35728733 PMCID: PMC9212859 DOI: 10.1016/j.annepidem.2022.06.005] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2021] [Revised: 06/02/2022] [Accepted: 06/10/2022] [Indexed: 12/14/2022]
Abstract
PURPOSE Even with an efficacious vaccine, protective behaviors (social distancing, masking) are essential for preventing COVID-19 transmission and could become even more important if current or future variants evade immunity from vaccines or prior infection. METHODS We created an agent-based model representing the Chicago population and conducted experiments to determine the effects of varying adult out-of-household activities (OOHA), school reopening, and protective behaviors across age groups on COVID-19 transmission and hospitalizations. RESULTS From September-November 2020, decreasing adult protective behaviors and increasing adult OOHA both substantially impacted COVID-19 outcomes; school reopening had relatively little impact when adult protective behaviors and OOHA were maintained. As of November 1, 2020, a 50% reduction in young adult (age 18-40) protective behaviors resulted in increased latent infection prevalence per 100,000 from 15.93 (IQR 6.18, 36.23) to 40.06 (IQR 14.65, 85.21) and 19.87 (IQR 6.83, 46.83) to 47.74 (IQR 18.89, 118.77) with 15% and 45% school reopening. Increasing adult (age ≥18) OOHA from 65% to 80% of prepandemic levels resulted in increased latent infection prevalence per 100,000 from 35.18 (IQR 13.59, 75.00) to 69.84 (IQR 33.27, 145.89) and 38.17 (IQR 15.84, 91.16) to 80.02 (IQR 30.91, 186.63) with 15% and 45% school reopening. Similar patterns were observed for hospitalizations. CONCLUSIONS In areas without widespread vaccination coverage, interventions to maintain adherence to protective behaviors, particularly among younger adults and in out-of-household settings, remain a priority for preventing COVID-19 transmission.
Collapse
Affiliation(s)
- Anna L. Hotton
- Department of Medicine, University of Chicago, Chicago, IL,Corresponding author: Department of Medicine, University of Chicago, 5837 N Maryland Ave, Chicago, IL, 60637
| | - Jonathan Ozik
- Argonne National Laboratory, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL, Northwestern Argonne Institute for Science and Engineering, Evanston, IL
| | | | - Nick Collier
- Argonne National Laboratory, Lemont, IL, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL
| | - Abby Stevens
- Department of Statistics, University of Chicago, Chicago, IL
| | - Aditya S. Khanna
- Department of Behavioral and Social Sciences, Brown University School of Public Health, Providence, RI
| | - Margaret M. MacDonell
- Radiological, Chemical and Environmental Risk Analysis (RACER), Environmental Science Division (EVS), Argonne National Laboratory, Lemont, IL
| | - Cheng Wang
- RACER EVS, Argonne National Laboratory, Lemont, IL, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL
| | | | - Young-Soo Chang
- Department of Climate and Earth System Science (CESS), EVS, Argonne National Laboratory, Lemont, IL
| | - Ignacio J. Martinez-Moyano
- Argonne National Laboratory, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL, Northwestern Argonne Institute for Science and Engineering, Evanston, IL
| | | | - Harold A. Pollack
- Crown School of Social Work Policy and Practice, University of Chicago, Chicago, IL
| | - John A. Schneider
- Departments of Medicine and Public Health Sciences, University of Chicago, Chicago, IL
| | - Charles Macal
- Argonne National Laboratory, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL, Northwestern Argonne Institute for Science and Engineering, Evanston, IL
| |
Collapse
|
79
|
Reichmuth ML, Hodcroft EB, Riou J, Neher RA, Hens N, Althaus CL. Impact of cross-border-associated cases on the SARS-CoV-2 epidemic in Switzerland during summer 2020 and 2021. Epidemics 2022; 41:100654. [PMID: 36444785 PMCID: PMC9671612 DOI: 10.1016/j.epidem.2022.100654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Revised: 10/01/2022] [Accepted: 11/12/2022] [Indexed: 11/18/2022] Open
Abstract
During the summers of 2020 and 2021, the number of confirmed cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in Switzerland remained at relatively low levels, but grew steadily over time. It remains unclear to what extent epidemic growth during these periods was a result of the relaxation of local control measures or increased traveling and subsequent importation of cases. A better understanding of the role of cross-border-associated cases (imports) on the local epidemic dynamics will help to inform future surveillance strategies. We analyzed routine surveillance data of confirmed cases of SARS-CoV-2 in Switzerland from 1 June to 30 September 2020 and 2021. We used a stochastic branching process model that accounts for superspreading of SARS-CoV-2 to simulate epidemic trajectories in absence and in presence of imports during summer 2020 and 2021. The Swiss Federal Office of Public Health reported 22,919 and 145,840 confirmed cases of SARS-CoV-2 from 1 June to 30 September 2020 and 2021, respectively. Among cases with known place of exposure, 27% (3,276 of 12,088) and 25% (1,110 of 4,368) reported an exposure abroad in 2020 and 2021, respectively. Without considering the impact of imported cases, the steady growth of confirmed cases during summer periods would be consistent with a value of Re that is significantly above the critical threshold of 1. In contrast, we estimated Re at 0.84 (95% credible interval, CrI: 0.78-0.90) in 2020 and 0.82 (95% CrI: 0.74-0.90) in 2021 when imported cases were taken into account, indicating that the local Re was below the critical threshold of 1 during summer. In Switzerland, cross-border-associated SARS-CoV-2 cases had a considerable impact on the local transmission dynamics and can explain the steady growth of the epidemic during the summers of 2020 and 2021.
Collapse
Affiliation(s)
- Martina L. Reichmuth
- Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland,Correspondence to: Institute of Social and Preventive Medicine, University of Bern, Mittelstrasse 43, CH-3012 Bern, Switzerland
| | - Emma B. Hodcroft
- Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Julien Riou
- Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland,Federal Office of Public Health, Liebefeld, Switzerland
| | - Richard A. Neher
- Swiss Institute of Bioinformatics, Lausanne, Switzerland,Biozentrum, University of Basel, Basel, Switzerland
| | - Niel Hens
- Interuniversity Institute for Biostatistics and statistical Bioinformatics, Data Science Institute, Hasselt University, Hasselt, Belgium,Centre for Health Economics Research and Modelling Infectious Diseases, Vaccine and Infectious Disease Institute, University of Antwerp, Antwerp, Belgium
| | - Christian L. Althaus
- Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland
| |
Collapse
|
80
|
Zhang Y, Britton T, Zhou X. Monitoring real-time transmission heterogeneity from incidence data. PLoS Comput Biol 2022; 18:e1010078. [PMID: 36455043 DOI: 10.1371/journal.pcbi.1010078] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Revised: 12/13/2022] [Accepted: 11/16/2022] [Indexed: 12/03/2022] Open
Abstract
The transmission heterogeneity of an epidemic is associated with a complex mixture of host, pathogen and environmental factors. And it may indicate superspreading events to reduce the efficiency of population-level control measures and to sustain the epidemic over a larger scale and a longer duration. Methods have been proposed to identify significant transmission heterogeneity in historic epidemics based on several data sources, such as contact history, viral genomes and spatial information, which may not be available, and more importantly ignore the temporal trend of transmission heterogeneity. Here we attempted to establish a convenient method to estimate real-time heterogeneity over an epidemic. Within the branching process framework, we introduced an instant-individualheterogenous infectiousness model to jointly characterize the variation in infectiousness both between individuals and among different times. With this model, we could simultaneously estimate the transmission heterogeneity and the reproduction number from incidence time series. We validated the model with data of both simulated and real outbreaks. Our estimates of the overall and real-time heterogeneities of the six epidemics were consistent with those presented in the literature. Additionally, our model is robust to the ubiquitous bias of under-reporting and misspecification of serial interval. By analyzing recent data from South Africa, we found evidence that the Omicron might be of more significant transmission heterogeneity than Delta. Our model based on incidence data was proved to be reliable in estimating the real-time transmission heterogeneity.
Collapse
Affiliation(s)
- Yunjun Zhang
- Department of Biostatistics, School of Public Health, Peking University, Beijing, China.,Center for Statistical Science, Peking University, Beijing, China
| | - Tom Britton
- Department of Mathematics, Stockholm University, Stockholm, Sweden
| | - Xiaohua Zhou
- Department of Biostatistics, School of Public Health, Peking University, Beijing, China.,Center for Statistical Science, Peking University, Beijing, China.,Beijing International Center for Mathematical Research, Peking University, Beijing, China.,School of Mathematical Sciences, Peking University, Beijing, China
| |
Collapse
|
81
|
Linton NM, Akhmetzhanov AR, Nishiura H. Correlation between times to SARS-CoV-2 symptom onset and secondary transmission undermines epidemic control efforts. Epidemics 2022; 41:100655. [PMID: 36413921 PMCID: PMC9661582 DOI: 10.1016/j.epidem.2022.100655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Revised: 10/12/2022] [Accepted: 11/12/2022] [Indexed: 11/16/2022] Open
Abstract
Severe acute respiratory coronavirus 2 (SARS-CoV-2) infections have been associated with substantial presymptomatic transmission, which occurs when the generation interval-the time between infection of an individual with a pathogen and transmission of the pathogen to another individual-is shorter than the incubation period-the time between infection and symptom onset. We collected a dataset of 257 SARS-CoV-2 transmission pairs in Japan during 2020 and jointly estimated the mean incubation period of infectors (4.8 days, 95 % CrI: 4.4-5.1 days), mean generation interval to when they infect others (4.3 days, 95 % credible interval [CrI]: 4.0-4.7 days), and the correlation (Kendall's tau: 0.5, 95 % CrI: 0.4-0.6) between these two epidemiological parameters. Our finding of a positive correlation and mean generation interval shorter than the mean infector incubation period indicates ample infectiousness before symptom onset and suggests that reliance on isolation of symptomatic COVID-19 cases as a focal point of control efforts is insufficient to address the challenges posed by SARS-CoV-2 transmission dynamics.
Collapse
Affiliation(s)
- Natalie M. Linton
- Kyoto University School of Public Health, Yoshidakonoe-cho, Sakyo-ku, Kyoto city, 606-8501, Japan,Graduate School of Medicine, Hokkaido University, Kita 15 Jo Nishi 7 Chome, Kita-ku, Sapporo-shi, Hokkaido 060-8638, Japan
| | - Andrei R. Akhmetzhanov
- College of Public Health, National Taiwan University, 17 Xu-Zhou Road, Taipei 10055, Taiwan
| | - Hiroshi Nishiura
- Kyoto University School of Public Health, Yoshidakonoe-cho, Sakyo-ku, Kyoto city, 606-8501, Japan,Core Research for Evolutional Science and Technology (CREST), Japan Science and Technology Agency, Saitama, Japan,Corresponding author at: Kyoto University School of Public Health, Yoshidakonoe-cho, Sakyo-ku, Kyoto city, 606-8501, Japan
| |
Collapse
|
82
|
Using mechanistic model-based inference to understand and project epidemic dynamics with time-varying contact and vaccination rates. Sci Rep 2022; 12:20451. [PMID: 36443439 PMCID: PMC9702885 DOI: 10.1038/s41598-022-25018-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2022] [Accepted: 11/23/2022] [Indexed: 11/29/2022] Open
Abstract
Epidemiological models range in complexity from relatively simple statistical models that make minimal assumptions about the variables driving epidemic dynamics to more mechanistic models that include effects such as vaccine-derived and infection-derived immunity, population structure and heterogeneity. The former are often fitted to data in real-time and used for short-term forecasting, while the latter are more suitable for comparing longer-term scenarios under differing assumptions about control measures or other factors. Here, we present a mechanistic model of intermediate complexity that can be fitted to data in real-time but is also suitable for investigating longer-term dynamics. Our approach provides a bridge between primarily empirical approaches to forecasting and assumption-driven scenario models. The model was developed as a policy advice tool for New Zealand's 2021 outbreak of the Delta variant of SARS-CoV-2 and includes the effects of age structure, non-pharmaceutical interventions, and the ongoing vaccine rollout occurring during the time period studied. We use an approximate Bayesian computation approach to infer the time-varying transmission coefficient from real-time data on reported cases. We then compare projections of the model with future, out-of-sample data. We find that this approach produces a good fit with in-sample data and reasonable forward projections given the inherent limitations of predicting epidemic dynamics during periods of rapidly changing policy and behaviour. Results from the model helped inform the New Zealand Government's policy response throughout the outbreak.
Collapse
|
83
|
Epidemiology and Transmission Dynamics of Infectious Diseases and Control Measures. Viruses 2022; 14:v14112510. [PMID: 36423119 PMCID: PMC9695084 DOI: 10.3390/v14112510] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Revised: 11/10/2022] [Accepted: 11/10/2022] [Indexed: 11/16/2022] Open
Abstract
The epidemiology and transmission dynamics of infectious diseases must be understood at the individual and community levels to improve public health decision-making for real-time and integrated community-based control strategies. Herein, we explore the epidemiological characteristics for assessing the impact of public health interventions in the community setting and their applications. Computational statistical methods could advance research on infectious disease epidemiology and accumulate scientific evidence of the potential impacts of pharmaceutical/nonpharmaceutical measures to mitigate or control infectious diseases in the community. Novel public health threats from emerging zoonotic infectious diseases are urgent issues. Given these direct and indirect mitigating impacts at various levels to different infectious diseases and their burdens, we must consider an integrated assessment approach, 'One Health', to understand the dynamics and control of infectious diseases.
Collapse
|
84
|
Late Surges in COVID-19 Cases and Varying Transmission Potential Partially Due to Public Health Policy Changes in 5 Western States, March 10, 2020, to January 10, 2021. Disaster Med Public Health Prep 2022; 17:e277. [PMID: 36325878 PMCID: PMC9794457 DOI: 10.1017/dmp.2022.248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
OBJECTIVE This study investigates the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission potential in North Dakota, South Dakota, Montana, Wyoming, and Idaho from March 2020 through January 2021. METHODS Time-varying reproduction numbers, R t , of a 7-d-sliding-window and of non-overlapping-windows between policy changes were estimated using the instantaneous reproduction number method. Linear regression was performed to evaluate if per-capita cumulative case-count varied across counties with different population size or density. RESULTS The median 7-d-sliding-window R t estimates across the studied region varied between 1 and 1.25 during September through November 2020. Between November 13 and 18, R t was reduced by 14.71% (95% credible interval, CrI, [14.41%, 14.99%]) in North Dakota following a mask mandate; Idaho saw a 1.93% (95% CrI [1.87%, 1.99%]) reduction and Montana saw a 9.63% (95% CrI [9.26%, 9.98%]) reduction following the tightening of restrictions. High-population and high-density counties had higher per-capita cumulative case-count in North Dakota on June 30, August 31, October 31, and December 31, 2020. In Idaho, North Dakota, South Dakota, and Wyoming, there were positive correlations between population size and per-capita weekly incident case-count, adjusted for calendar time and social vulnerability index variables. CONCLUSIONS R t decreased after mask mandate during the region's case-count spike suggested reduction in SARS-CoV-2 transmission.
Collapse
|
85
|
Forecasting COVID-19 cases by assessing control-intervention effects in Republic of Korea: A statistical modeling approach. ALEXANDRIA ENGINEERING JOURNAL 2022; 61:9203-9217. [PMCID: PMC8872739 DOI: 10.1016/j.aej.2022.02.037] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 02/02/2022] [Accepted: 02/12/2022] [Indexed: 05/25/2023]
Abstract
The Coronavirus disease of 2019 (COVID-19) is an ongoing public health concern worldwide. COVID-19 infections continue to occur and thus, it is important to assess the effects of various public health measures. This study aims to forecast COVID-19 cases by geographical area in Korea, based on the effects of different control-intervention intensities (CII). Methods involved estimating the effective reproduction number (Rt) by Korean geographical area using the SEIHR model, and the instantaneous reproduction number using statistical model, comparing the epidemic curves and high-, intermediate-, and low-intensity control interventions. Here, short-term four-week forecasts by geographical area were conducted. The mean of delayed instantaneous reproduction number was estimated at 1.36, 1.03, and 0.93 for the low-, intermediate-, and high-intensity control interventions, respectively, in the capital area of Korea from July 16, 2020, to March 4, 2021. The COVID-19 cases were forecasted with an accuracy rate of 11.28%, 13.62%, and 20.19% MAPE in Korea, including both the capital and non-capital areas. High-intensity control measures significantly reduced the reproduction number to be less than one. The proposed model forecasted COVID-19 transmission dynamics with good accuracy and interpretability. High-intensity control intervention, active case detection, and isolation efforts should be maintained to control the pandemic.
Collapse
|
86
|
da Costa Avelar PH, Del Coco N, Lamb LC, Tsoka S, Cardoso-Silva J. A Bayesian predictive analytics model for improving long range epidemic forecasting during an infection wave. HEALTHCARE ANALYTICS (NEW YORK, N.Y.) 2022; 2:100115. [PMID: 37520620 PMCID: PMC9533637 DOI: 10.1016/j.health.2022.100115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Revised: 08/17/2022] [Accepted: 09/26/2022] [Indexed: 11/04/2022]
Abstract
Following the outbreak of the coronavirus epidemic in early 2020, municipalities, regional governments and policymakers worldwide had to plan their Non-Pharmaceutical Interventions (NPIs) amidst a scenario of great uncertainty. At this early stage of an epidemic, where no vaccine or medical treatment is in sight, algorithmic prediction can become a powerful tool to inform local policymaking. However, when we replicated one prominent epidemiological model to inform health authorities in a region in the south of Brazil, we found that this model relied too heavily on manually predetermined covariates and was too reactive to changes in data trends. Our four proposed models access data of both daily reported deaths and infections as well as take into account missing data (e.g., the under-reporting of cases) more explicitly, with two of the proposed versions also attempting to model the delay in test reporting. We simulated weekly forecasting of deaths from the period from 31/05/2020 until 31/01/2021, with first week data being used as a cold-start to the algorithm, after which we use a lighter variant of the model for faster forecasting. Because our models are significantly more proactive in identifying trend changes, this has improved forecasting, especially in long-range predictions and after the peak of an infection wave, as they were quicker to adapt to scenarios after these peaks in reported deaths. Assuming reported cases were under-reported greatly benefited the model in its stability, and modelling retroactively-added data (due to the "hot" nature of the data used) had a negligible impact on performance.
Collapse
Affiliation(s)
- Pedro Henrique da Costa Avelar
- Data Science Brigade, Porto Alegre, Rio Grande do Sul, Brazil
- Institute of Informatics, Federal University of Rio Grande do Sul, Porto Alegre, Rio Grande do Sul, Brazil
- Department of Informatics, King's College London, London, United Kingdom
- Machine Intellection Department, Institute for Infocomm Research, A*STAR, Singapore
| | | | - Luis C Lamb
- Institute of Informatics, Federal University of Rio Grande do Sul, Porto Alegre, Rio Grande do Sul, Brazil
| | - Sophia Tsoka
- Department of Informatics, King's College London, London, United Kingdom
| | - Jonathan Cardoso-Silva
- Data Science Brigade, Porto Alegre, Rio Grande do Sul, Brazil
- Department of Informatics, King's College London, London, United Kingdom
- Data Science Institute, London School of Economics and Political Science, London, United Kingdom
| |
Collapse
|
87
|
Pellis L, Birrell PJ, Blake J, Overton CE, Scarabel F, Stage HB, Brooks‐Pollock E, Danon L, Hall I, House TA, Keeling MJ, Read JM, De Angelis D. Estimation of reproduction numbers in real time: Conceptual and statistical challenges. JOURNAL OF THE ROYAL STATISTICAL SOCIETY. SERIES A, (STATISTICS IN SOCIETY) 2022; 185:S112-S130. [PMID: 37063605 PMCID: PMC10100071 DOI: 10.1111/rssa.12955] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Accepted: 07/06/2022] [Indexed: 06/19/2023]
Abstract
The reproduction numberR has been a central metric of the COVID-19 pandemic response, published weekly by the UK government and regularly reported in the media. Here, we provide a formal definition and discuss the advantages and most common misconceptions around this quantity. We consider the intuition behind different formulations ofR , the complexities in its estimation (including the unavoidable lags involved), and its value compared to other indicators (e.g. the growth rate) that can be directly observed from aggregate surveillance data and react more promptly to changes in epidemic trend. As models become more sophisticated, with age and/or spatial structure, formulatingR becomes increasingly complicated and inevitably model-dependent. We present some models currently used in the UK pandemic response as examples. Ultimately, limitations in the available data streams, data quality and time constraints force pragmatic choices to be made on a quantity that is an average across time, space, social structure and settings. Effectively communicating these challenges is important but often difficult in an emergency.
Collapse
Affiliation(s)
- Lorenzo Pellis
- Department of MathematicsThe University of ManchesterManchesterUK
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- The Alan Turing InstituteLondonUK
| | - Paul J. Birrell
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- MRC Biostatistics Unit, School of Clinical MedicineUniversity of CambridgeCambridgeUK
- Statistics Modelling and Economics DepartmentPublic Health EnglandLondonUK
- Joint Modelling TeamPublic Health EnglandLondonUK
| | - Joshua Blake
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- MRC Biostatistics Unit, School of Clinical MedicineUniversity of CambridgeCambridgeUK
| | - Christopher E. Overton
- Department of MathematicsThe University of ManchesterManchesterUK
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- Manchester University NHS Foundation TrustManchesterUK
| | - Francesca Scarabel
- Department of MathematicsThe University of ManchesterManchesterUK
- Joint UNIversities Pandemic and Epidemiological ResearchUK
| | - Helena B. Stage
- Department of MathematicsThe University of ManchesterManchesterUK
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- Department of PhysicsHumboldt University of BerlinBerlinGermany
- Department of Physics and AstronomyUniversity of PotsdamPotsdamGermany
| | - Ellen Brooks‐Pollock
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- NIHR Health Protection Research Unit (HPRU) in Behavioural Science and Evaluation, Population Health SciencesUniversity of BristolBristolUK
| | - Leon Danon
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- The Alan Turing InstituteLondonUK
- Department of Engineering MathematicsUniversity of BristolBristolUK
| | - Ian Hall
- Department of MathematicsThe University of ManchesterManchesterUK
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- The Alan Turing InstituteLondonUK
- Joint Modelling TeamPublic Health EnglandLondonUK
- School of Health SciencesThe University of ManchesterManchesterUK
| | - Thomas A. House
- Department of MathematicsThe University of ManchesterManchesterUK
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- The Alan Turing InstituteLondonUK
| | - Matt J. Keeling
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- Mathematics Institute and School of Life SciencesUniversity of WarwickCoventryUK
| | - Jonathan M. Read
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- Centre for Health Informatics, Computing and Statistics, Lancaster Medical SchoolLancaster UniversityLancasterUK
| | | | - Daniela De Angelis
- Joint UNIversities Pandemic and Epidemiological ResearchUK
- MRC Biostatistics Unit, School of Clinical MedicineUniversity of CambridgeCambridgeUK
- Statistics Modelling and Economics DepartmentPublic Health EnglandLondonUK
- Joint Modelling TeamPublic Health EnglandLondonUK
| |
Collapse
|
88
|
Bekker‐Nielsen Dunbar M, Hofmann F, Held L. Assessing the effect of school closures on the spread of COVID-19 in Zurich. JOURNAL OF THE ROYAL STATISTICAL SOCIETY. SERIES A, (STATISTICS IN SOCIETY) 2022; 185:S131-S142. [PMID: 38607867 PMCID: PMC9878126 DOI: 10.1111/rssa.12910] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Accepted: 04/08/2022] [Indexed: 04/14/2024]
Abstract
The effect of school closure on the spread of COVID-19 has been discussed intensively in the literature and the news. To capture the interdependencies between children and adults, we consider daily age-stratified incidence data and contact patterns between age groups which change over time to reflect social distancing policy indicators. We fit a multivariate time-series endemic-epidemic model to such data from the Canton of Zurich, Switzerland and use the model to predict the age-specific incidence in a counterfactual approach (with and without school closures). The results indicate a 17% median increase of incidence in the youngest age group (0-14 year olds), whereas the relative increase in the other age groups drops to values between 2% and 3%. We argue that our approach is more informative to policy makers than summarising the effect of school closures with time-dependent effective reproduction numbers, which are difficult to estimate due to the sparsity of incidence counts within the relevant age groups.
Collapse
Affiliation(s)
| | - Felix Hofmann
- Epidemiology, Biostatistics and Prevention Institute (EBPI)University of Zurich (UZH)ZurichSwitzerland
| | - Leonhard Held
- Epidemiology, Biostatistics and Prevention Institute (EBPI)University of Zurich (UZH)ZurichSwitzerland
| | | |
Collapse
|
89
|
Binny RN, Lustig A, Hendy SC, Maclaren OJ, Ridings KM, Vattiato G, Plank MJ. Real-time estimation of the effective reproduction number of SARS-CoV-2 in Aotearoa New Zealand. PeerJ 2022; 10:e14119. [PMID: 36275456 PMCID: PMC9583856 DOI: 10.7717/peerj.14119] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2022] [Accepted: 09/04/2022] [Indexed: 01/21/2023] Open
Abstract
During an epidemic, real-time estimation of the effective reproduction number supports decision makers to introduce timely and effective public health measures. We estimate the time-varying effective reproduction number, Rt , during Aotearoa New Zealand's August 2021 outbreak of the Delta variant of SARS-CoV-2, by fitting the publicly available EpiNow2 model to New Zealand case data. While we do not explicitly model non-pharmaceutical interventions or vaccination coverage, these two factors were the leading drivers of variation in transmission in this period and we describe how changes in these factors coincided with changes in Rt . Alert Level 4, New Zealand's most stringent restriction setting which includes stay-at-home measures, was initially effective at reducing the median Rt to 0.6 (90% CrI 0.4, 0.8) on 29 August 2021. As New Zealand eased certain restrictions and switched from an elimination strategy to a suppression strategy, Rt subsequently increased to a median 1.3 (1.2, 1.4). Increasing vaccination coverage along with regional restrictions were eventually sufficient to reduce Rt below 1. The outbreak peaked at an estimated 198 (172, 229) new infected cases on 10 November, after which cases declined until January 2022. We continue to update Rt estimates in real time as new case data become available to inform New Zealand's ongoing pandemic response.
Collapse
Affiliation(s)
- Rachelle N. Binny
- Manaaki Whenua-Landcare Research, Lincoln, New Zealand,Te Pūnaha Matatini, Auckland, New Zealand
| | - Audrey Lustig
- Manaaki Whenua-Landcare Research, Lincoln, New Zealand,Te Pūnaha Matatini, Auckland, New Zealand
| | - Shaun C. Hendy
- Te Pūnaha Matatini, Auckland, New Zealand,Department of Physics, University of Auckland, Auckland, New Zealand
| | - Oliver J. Maclaren
- Department of Engineering Science, University of Auckland, Auckland, New Zealand
| | - Kannan M. Ridings
- Te Pūnaha Matatini, Auckland, New Zealand,Department of Physics, University of Auckland, Auckland, New Zealand
| | - Giorgia Vattiato
- Te Pūnaha Matatini, Auckland, New Zealand,Department of Physics, University of Auckland, Auckland, New Zealand,School of Mathematics and Statistics, University of Canterbury, Christchurch, New Zealand
| | - Michael J. Plank
- Te Pūnaha Matatini, Auckland, New Zealand,School of Mathematics and Statistics, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
90
|
Gressani O, Wallinga J, Althaus CL, Hens N, Faes C. EpiLPS: A fast and flexible Bayesian tool for estimation of the time-varying reproduction number. PLoS Comput Biol 2022; 18:e1010618. [PMID: 36215319 PMCID: PMC9584461 DOI: 10.1371/journal.pcbi.1010618] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2022] [Revised: 10/20/2022] [Accepted: 09/30/2022] [Indexed: 11/17/2022] Open
Abstract
In infectious disease epidemiology, the instantaneous reproduction number [Formula: see text] is a time-varying parameter defined as the average number of secondary infections generated by an infected individual at time t. It is therefore a crucial epidemiological statistic that assists public health decision makers in the management of an epidemic. We present a new Bayesian tool (EpiLPS) for robust estimation of the time-varying reproduction number. The proposed methodology smooths the epidemic curve and allows to obtain (approximate) point estimates and credible intervals of [Formula: see text] by employing the renewal equation, using Bayesian P-splines coupled with Laplace approximations of the conditional posterior of the spline vector. Two alternative approaches for inference are presented: (1) an approach based on a maximum a posteriori argument for the model hyperparameters, delivering estimates of [Formula: see text] in only a few seconds; and (2) an approach based on a Markov chain Monte Carlo (MCMC) scheme with underlying Langevin dynamics for efficient sampling of the posterior target distribution. Case counts per unit of time are assumed to follow a negative binomial distribution to account for potential overdispersion in the data that would not be captured by a classic Poisson model. Furthermore, after smoothing the epidemic curve, a "plug-in'' estimate of the reproduction number can be obtained from the renewal equation yielding a closed form expression of [Formula: see text] as a function of the spline parameters. The approach is extremely fast and free of arbitrary smoothing assumptions. EpiLPS is applied on data of SARS-CoV-1 in Hong-Kong (2003), influenza A H1N1 (2009) in the USA and on the SARS-CoV-2 pandemic (2020-2021) for Belgium, Portugal, Denmark and France.
Collapse
Affiliation(s)
- Oswaldo Gressani
- Interuniversity Institute for Biostatistics and statistical Bioinformatics (I-BioStat), Data Science Institute, Hasselt University, Hasselt, Belgium,* E-mail:
| | - Jacco Wallinga
- Centre for Infectious Disease Control, National Institute for Public Health and the Environment, Bilthoven, The Netherlands,Department of Biomedical Data Sciences, Leiden University Medical Centre, Leiden, The Netherlands
| | - Christian L. Althaus
- Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland
| | - Niel Hens
- Interuniversity Institute for Biostatistics and statistical Bioinformatics (I-BioStat), Data Science Institute, Hasselt University, Hasselt, Belgium,Centre for Health Economics Research and Modelling Infectious Diseases, Vaxinfectio, University of Antwerp, Antwerp, Belgium
| | - Christel Faes
- Interuniversity Institute for Biostatistics and statistical Bioinformatics (I-BioStat), Data Science Institute, Hasselt University, Hasselt, Belgium
| |
Collapse
|
91
|
Li W, Bulekova K, Gregor B, White LF, Kolaczyk ED. Estimation of local time-varying reproduction numbers in noisy surveillance data. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20210303. [PMID: 35965456 PMCID: PMC9376722 DOI: 10.1098/rsta.2021.0303] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 04/11/2022] [Indexed: 05/04/2023]
Abstract
A valuable metric in understanding local infectious disease dynamics is the local time-varying reproduction number, i.e. the expected number of secondary local cases caused by each infected individual. Accurate estimation of this quantity requires distinguishing cases arising from local transmission from those imported from elsewhere. Realistically, we can expect identification of cases as local or imported to be imperfect. We study the propagation of such errors in estimation of the local time-varying reproduction number. In addition, we propose a Bayesian framework for estimation of the true local time-varying reproduction number when identification errors exist. And we illustrate the practical performance of our estimator through simulation studies and with outbreaks of COVID-19 in Hong Kong and Victoria, Australia. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.
Collapse
Affiliation(s)
- Wenrui Li
- Department of Mathematics and Statistics, Boston University, Boston, MA 02215, USA
| | - Katia Bulekova
- Research Computing Services, Information Services and Technology Boston University, Boston, MA 02215, USA
| | - Brian Gregor
- Research Computing Services, Information Services and Technology Boston University, Boston, MA 02215, USA
| | - Laura F. White
- Department of Biostatistics, Boston University, Boston, MA 02215, USA
| | - Eric D. Kolaczyk
- Department of Mathematics and Statistics, Boston University, Boston, MA 02215, USA
- Hariri Institute for Computing, Boston University, Boston,MA 02215, USA
| |
Collapse
|
92
|
Hinch R, Panovska-Griffiths J, Probert WJM, Ferretti L, Wymant C, Di Lauro F, Baya N, Ghafari M, Abeler-Dörner L, Fraser C. Estimating SARS-CoV-2 variant fitness and the impact of interventions in England using statistical and geo-spatial agent-based models. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022. [PMID: 35965459 DOI: 10.6084/m9.figshare.c.6067650] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
The SARS-CoV-2 epidemic has been extended by the evolution of more transmissible viral variants. In autumn 2020, the B.1.177 lineage became the dominant variant in England, before being replaced by the B.1.1.7 (Alpha) lineage in late 2020, with the sweep occurring at different times in each region. This period coincided with a large number of non-pharmaceutical interventions (e.g. lockdowns) to control the epidemic, making it difficult to estimate the relative transmissibility of variants. In this paper, we model the spatial spread of these variants in England using a meta-population agent-based model which correctly characterizes the regional variation in cases and distribution of variants. As a test of robustness, we additionally estimated the relative transmissibility of multiple variants using a statistical model based on the renewal equation, which simultaneously estimates the effective reproduction number R. Relative to earlier variants, the transmissibility of B.1.177 is estimated to have increased by 1.14 (1.12-1.16) and that of Alpha by 1.71 (1.65-1.77). The vaccination programme starting in December 2020 is also modelled. Counterfactual simulations demonstrate that the vaccination programme was essential for reopening in March 2021, and that if the January lockdown had started one month earlier, up to 30 k (24 k-38 k) deaths could have been prevented. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.
Collapse
Affiliation(s)
- Robert Hinch
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Jasmina Panovska-Griffiths
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
- The Queen's College, and, University of Oxford, Oxford, UK
| | - William J M Probert
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Luca Ferretti
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Chris Wymant
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Francesco Di Lauro
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Nikolas Baya
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Mahan Ghafari
- Department of Zoology, University of Oxford, Oxford, UK
| | - Lucie Abeler-Dörner
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Christophe Fraser
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| |
Collapse
|
93
|
Creswell R, Augustin D, Bouros I, Farm HJ, Miao S, Ahern A, Robinson M, Lemenuel-Diot A, Gavaghan DJ, Lambert BC, Thompson RN. Heterogeneity in the onwards transmission risk between local and imported cases affects practical estimates of the time-dependent reproduction number. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20210308. [PMID: 35965464 PMCID: PMC9376709 DOI: 10.1098/rsta.2021.0308] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/07/2021] [Accepted: 05/04/2022] [Indexed: 05/02/2023]
Abstract
During infectious disease outbreaks, inference of summary statistics characterizing transmission is essential for planning interventions. An important metric is the time-dependent reproduction number (Rt), which represents the expected number of secondary cases generated by each infected individual over the course of their infectious period. The value of Rt varies during an outbreak due to factors such as varying population immunity and changes to interventions, including those that affect individuals' contact networks. While it is possible to estimate a single population-wide Rt, this may belie differences in transmission between subgroups within the population. Here, we explore the effects of this heterogeneity on Rt estimates. Specifically, we consider two groups of infected hosts: those infected outside the local population (imported cases), and those infected locally (local cases). We use a Bayesian approach to estimate Rt, made available for others to use via an online tool, that accounts for differences in the onwards transmission risk from individuals in these groups. Using COVID-19 data from different regions worldwide, we show that different assumptions about the relative transmission risk between imported and local cases affect Rt estimates significantly, with implications for interventions. This highlights the need to collect data during outbreaks describing heterogeneities in transmission between different infected hosts, and to account for these heterogeneities in methods used to estimate Rt. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.
Collapse
Affiliation(s)
- R. Creswell
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - D. Augustin
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - I. Bouros
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - H. J. Farm
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - S. Miao
- Mathematical Institute, University of Oxford, Oxford OX2 6GG, UK
| | - A. Ahern
- Mathematical Institute, University of Oxford, Oxford OX2 6GG, UK
| | - M. Robinson
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - A. Lemenuel-Diot
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Basel CH-4070, Switzerland
| | - D. J. Gavaghan
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - B. C. Lambert
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - R. N. Thompson
- Mathematics Institute, University of Warwick, Coventry CV4 7AL, UK
- Zeeman Institute for Systems Biology and Infectious Disease Epidemiology Research, University of Warwick, Coventry CV4 7AL, UK
| |
Collapse
|
94
|
Hinch R, Panovska-Griffiths J, Probert WJM, Ferretti L, Wymant C, Di Lauro F, Baya N, Ghafari M, Abeler-Dörner L, Fraser C. Estimating SARS-CoV-2 variant fitness and the impact of interventions in England using statistical and geo-spatial agent-based models. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20210304. [PMID: 35965459 PMCID: PMC9376717 DOI: 10.1098/rsta.2021.0304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2021] [Accepted: 02/22/2022] [Indexed: 05/04/2023]
Abstract
The SARS-CoV-2 epidemic has been extended by the evolution of more transmissible viral variants. In autumn 2020, the B.1.177 lineage became the dominant variant in England, before being replaced by the B.1.1.7 (Alpha) lineage in late 2020, with the sweep occurring at different times in each region. This period coincided with a large number of non-pharmaceutical interventions (e.g. lockdowns) to control the epidemic, making it difficult to estimate the relative transmissibility of variants. In this paper, we model the spatial spread of these variants in England using a meta-population agent-based model which correctly characterizes the regional variation in cases and distribution of variants. As a test of robustness, we additionally estimated the relative transmissibility of multiple variants using a statistical model based on the renewal equation, which simultaneously estimates the effective reproduction number R. Relative to earlier variants, the transmissibility of B.1.177 is estimated to have increased by 1.14 (1.12-1.16) and that of Alpha by 1.71 (1.65-1.77). The vaccination programme starting in December 2020 is also modelled. Counterfactual simulations demonstrate that the vaccination programme was essential for reopening in March 2021, and that if the January lockdown had started one month earlier, up to 30 k (24 k-38 k) deaths could have been prevented. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.
Collapse
Affiliation(s)
- Robert Hinch
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Jasmina Panovska-Griffiths
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
- The Queen's College, University of Oxford, Oxford, UK
| | - William J. M. Probert
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Luca Ferretti
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Chris Wymant
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Francesco Di Lauro
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Nikolas Baya
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Mahan Ghafari
- Department of Zoology, University of Oxford, Oxford, UK
| | - Lucie Abeler-Dörner
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | | | - Christophe Fraser
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| |
Collapse
|
95
|
Costa GS, Cota W, Ferreira SC. Data-driven approach in a compartmental epidemic model to assess undocumented infections. CHAOS, SOLITONS, AND FRACTALS 2022; 163:112520. [PMID: 35996714 PMCID: PMC9385215 DOI: 10.1016/j.chaos.2022.112520] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Revised: 07/27/2022] [Accepted: 07/29/2022] [Indexed: 06/15/2023]
Abstract
Nowcasting and forecasting of epidemic spreading rely on incidence series of reported cases to derive the fundamental epidemiological parameters for a given pathogen. Two relevant drawbacks for predictions are the unknown fractions of undocumented cases and levels of nonpharmacological interventions, which span highly heterogeneously across different places and times. We describe a simple data-driven approach using a compartmental model including asymptomatic and pre-symptomatic contagions that allows to estimate both the level of undocumented infections and the value of effective reproductive number R t from time series of reported cases, deaths, and epidemiological parameters. The method was applied to epidemic series for COVID-19 across different municipalities in Brazil allowing to estimate the heterogeneity level of under-reporting across different places. The reproductive number derived within the current framework is little sensitive to both diagnosis and infection rates during the asymptomatic states. The methods described here can be extended to more general cases if data is available and adapted to other epidemiological approaches and surveillance data.
Collapse
Affiliation(s)
- Guilherme S Costa
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, MG, Brazil
| | - Wesley Cota
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, MG, Brazil
| | - Silvio C Ferreira
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, MG, Brazil
- National Institute of Science and technology for Complex Systems, Centro Brasileiro de Pesquisas Físicas, Rua Xavier Sigaud 150, 22290-180, Rio de Janeiro, Brazil
| |
Collapse
|
96
|
Meagher J, Friel N. Assessing epidemic curves for evidence of superspreading. JOURNAL OF THE ROYAL STATISTICAL SOCIETY. SERIES A, (STATISTICS IN SOCIETY) 2022; 185:2179-2202. [PMID: 37066104 PMCID: PMC10092342 DOI: 10.1111/rssa.12919] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/22/2021] [Accepted: 07/10/2022] [Indexed: 06/19/2023]
Abstract
The expected number of secondary infections arising from each index case, referred to as the reproduction or R number, is a vital summary statistic for understanding and managing epidemic diseases. There are many methods for estimating R ; however, few explicitly model heterogeneous disease reproduction, which gives rise to superspreading within the population. We propose a parsimonious discrete-time branching process model for epidemic curves that incorporates heterogeneous individual reproduction numbers. Our Bayesian approach to inference illustrates that this heterogeneity results in less certainty on estimates of the time-varying cohort reproduction number R t . We apply these methods to a COVID-19 epidemic curve for the Republic of Ireland and find support for heterogeneous disease reproduction. Our analysis allows us to estimate the expected proportion of secondary infections attributable to the most infectious proportion of the population. For example, we estimate that the 20% most infectious index cases account for approximately 75%-98% of the expected secondary infections with 95% posterior probability. In addition, we highlight that heterogeneity is a vital consideration when estimating R t .
Collapse
Affiliation(s)
- Joe Meagher
- Insight Centre for Data Analytics, School of Mathematics and StatisticsUniversity College DublinDublinIreland
| | - Nial Friel
- Insight Centre for Data Analytics, School of Mathematics and StatisticsUniversity College DublinDublinIreland
| |
Collapse
|
97
|
Bushaj S, Yin X, Beqiri A, Andrews D, Büyüktahtakın İE. A simulation-deep reinforcement learning (SiRL) approach for epidemic control optimization. ANNALS OF OPERATIONS RESEARCH 2022; 328:1-33. [PMID: 36187178 PMCID: PMC9512996 DOI: 10.1007/s10479-022-04926-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 07/28/2022] [Indexed: 05/12/2023]
Abstract
In this paper, we address the controversies of epidemic control planning by developing a novel Simulation-Deep Reinforcement Learning (SiRL) model. COVID-19 reminded constituents over the world that government decision-making could change their lives. During the COVID-19 pandemic, governments were concerned with reducing fatalities as the virus spread but at the same time also maintaining a flowing economy. In this paper, we address epidemic decision-making regarding the interventions necessary given of the epidemic based on the purpose of the decision-maker. Further, we intend to compare different vaccination strategies, such as age-based and random vaccination, to shine a light on who should get priority in the vaccination process. To address these issues, we propose a simulation-deep reinforcement learning (DRL) framework. This framework is composed of an agent-based simulation model and a governor DRL agent that can enforce interventions in the agent-based simulation environment. Computational results show that our DRL agent can learn effective strategies and suggest optimal actions given a specific epidemic situation based on a multi-objective reward structure. We compare our DRL agent's decisions to government interventions at different periods of time during the COVID-19 pandemic. Our results suggest that more could have been done to control the epidemic. In addition, if a random vaccination strategy that allows super-spreaders to get vaccinated early were used, infections would have been reduced by 32% at the expense of 4% more deaths. We also show that a behavioral change of fully quarantining 10% of the risky individuals and using a random vaccination strategy leads to a reduction of the death toll by 14% and 27% compared to the age-based vaccination strategy that was implemented and the New Jersey reported data, respectively. We have also demonstrated the flexibility of our approach to be applied to other locations by validating and applying our model to the COVID-19 case in the state of Kansas.
Collapse
Affiliation(s)
- Sabah Bushaj
- Department of Management Information Systems and Analytics, School of Business and Economics, SUNY Plattsburgh, Plattsburgh, NY USA
| | | | - Arjeta Beqiri
- Department of Management Information Systems and Analytics, School of Business and Economics, SUNY Plattsburgh, Plattsburgh, NY USA
| | - Donald Andrews
- Trinity College Dublin, School of Natural Sciences, Dublin, Ireland
| | - İ. Esra Büyüktahtakın
- Grado Department of Industrial and Systems Engineering, Virginia Tech, Blacksburg, VA USA
| |
Collapse
|
98
|
Inward RPD, Parag KV, Faria NR. Using multiple sampling strategies to estimate SARS-CoV-2 epidemiological parameters from genomic sequencing data. Nat Commun 2022; 13:5587. [PMID: 36151084 PMCID: PMC9508174 DOI: 10.1038/s41467-022-32812-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Accepted: 08/16/2022] [Indexed: 11/09/2022] Open
Abstract
The choice of viral sequences used in genetic and epidemiological analysis is important as it can induce biases that detract from the value of these rich datasets. This raises questions about how a set of sequences should be chosen for analysis. We provide insights on these largely understudied problems using SARS-CoV-2 genomic sequences from Hong Kong, China, and the Amazonas State, Brazil. We consider multiple sampling schemes which were used to estimate Rt and rt as well as related R0 and date of origin parameters. We find that both Rt and rt are sensitive to changes in sampling whilst R0 and the date of origin are relatively robust. Moreover, we find that analysis using unsampled datasets result in the most biased Rt and rt estimates for both our Hong Kong and Amazonas case studies. We highlight that sampling strategy choices may be an influential yet neglected component of sequencing analysis pipelines.
Collapse
Affiliation(s)
| | - Kris V Parag
- MRC Centre of Global Infectious Disease Analysis, Jameel Institute for Disease and Emergency Analytics, Imperial College London, London, UK.
- NIHR Health Protection Research Unit in Behavioural Science and Evaluation, University of Bristol, Bristol, UK.
| | - Nuno R Faria
- Department of Zoology, University of Oxford, Oxford, UK.
- MRC Centre of Global Infectious Disease Analysis, Jameel Institute for Disease and Emergency Analytics, Imperial College London, London, UK.
- Instituto de Medicina Tropical, Faculdade de Medicina da Universidade de Sao Paulo, Sao Paulo, Brazil.
| |
Collapse
|
99
|
Tracking changes in SARS-CoV-2 transmission with a novel outpatient sentinel surveillance system in Chicago, USA. Nat Commun 2022; 13:5547. [PMID: 36138039 PMCID: PMC9499975 DOI: 10.1038/s41467-022-33317-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Accepted: 09/12/2022] [Indexed: 01/08/2023] Open
Abstract
Public health indicators typically used for COVID-19 surveillance can be biased or lag changing community transmission patterns. In this study, we investigate whether sentinel surveillance of recently symptomatic individuals receiving outpatient diagnostic testing for SARS-CoV-2 could accurately assess the instantaneous reproductive number R(t) and provide early warning of changes in transmission. We use data from community-based diagnostic testing sites in the United States city of Chicago. Patients tested at community-based diagnostic testing sites between September 2020 and June 2021, and reporting symptom onset within four days preceding their test, formed the sentinel population. R(t) calculated from sentinel cases agreed well with R(t) from other indicators. Retrospectively, trends in sentinel cases did not precede trends in COVID-19 hospital admissions by any identifiable lead time. In deployment, sentinel surveillance held an operational recency advantage of nine days over hospital admissions. The promising performance of opportunistic sentinel surveillance suggests that deliberately designed outpatient sentinel surveillance would provide robust early warning of increasing transmission.
Collapse
|
100
|
Epidemiological model based periodic intervention policies for COVID-19 mitigation in the United Kingdom. Sci Rep 2022; 12:15660. [PMID: 36123382 PMCID: PMC9483909 DOI: 10.1038/s41598-022-19630-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Accepted: 08/31/2022] [Indexed: 11/08/2022] Open
Abstract
As the UK, together with numerous countries in the world, moves towards a new phase of the COVID-19 pandemic, there is a need to be able to predict trends in sufficient time to limit the pressure faced by the National Health Service (NHS) and maintain low hospitalisation levels. In this study, we explore the use of an epidemiological compartmental model to devise a periodic adaptive suppression/intervention policy to alleviate the pressure on the NHS. The proposed model facilitates the understanding of the progression of the specific stages of COVID-19 in communities in the UK including: the susceptible population, the infected population, the hospitalised population, the recovered population, the deceased population, and the vaccinated population. We identify the parameters of the model by relying on past data within the period from 1 October 2020 to 1 June 2021. We use the total number of hospitalised patients and the fraction of those infected who are being admitted to hospital to develop adaptive policies: these modulate the recommended level of social restriction measures and realisable vaccination target adjustments. The analysis over the period 1 October 2020 to 1 June 2021 demonstrates our periodic adaptive policies have the potential to reduce the hospitalisation by 58% on average per month. In a further prospective analysis over the period August 2021 to May 2022, we analyse several future scenarios, characterised by the relaxation of restrictions, the vaccination ineffectiveness and the gradual decay of the vaccination-induced immunity within the population. In addition, we simulate the surge of plausible variants characterised by an higher transmission rate. In such scenarios, we show that our periodic intervention is effective and able to maintain the hospitalisation rate to a manageable level.
Collapse
|