Abstract:Accurately estimating the effective reproduction number (R t ) of a circulating pathogen is a fundamental challenge in the study of infectious disease. The fields of epidemiology and pathogen phylodynamics both share this goal, but to date, methodologies and data employed by each remain largely distinct. Here we present EpiFusion: a joint approach that can be used to harness the complementary strengths of each field to improve estimation of outbreak dynamics for large and poorly sampled epidemics, such as arboviral or respiratory virus outbreaks, and validate it for retrospective analysis. We propose a model of R t that estimates outbreak trajectories conditional upon both phylodynamic (time-scaled trees estimated from genetic sequences) and epidemiological (case incidence) data. We simulate stochastic outbreak trajectories that are weighted according to epidemiological and phylodynamic observation models and fit using particle Markov Chain Monte Carlo. To assess performance, we test EpiFusion on simulated outbreaks in which transmission and/or surveillance rapidly changes and find that using EpiFusion to combine epidemiological and phylodynamic data maintains accuracy and increases certainty in trajectory and R t estimates, compared to when each data type is used alone. We benchmark EpiFusion's performance against existing methods to estimate R t and demonstrate advances in speed and accuracy. Importantly, our approach scales efficiently with dataset size. Finally, we apply our model to estimate R t during the 2014 Ebola outbreak in Sierra Leone. EpiFusion is designed to accommodate future extensions that will improve its utility, such as explicitly modelling population structure, accommodations for phylogenetic uncertainty, and the ability to weight the contributions of genomic or case incidence to the inference. Understanding infectious disease spread is fundamental to protecting public health, but can be challenging as disease spread is a phenomenon that cannot be directly observed. So, epidemiologists use data in conjunction with mathematical models to estimate disease dynamics. Often, combinations of different models and data can be used to answer the same questions–for example 'traditional' epidemiology commonly uses case incidence data (the number of people who have tested positive for a disease during a certain time period) whereas phylodynamic models use pathogen genomic sequence data and our knowledge of the way their genomes evolve to model disease population dynamics. Each of these approaches have strengths and limitations, and data of each type can be sparse or biased, particularly during rapidly developing outbreaks or in countries with poor pathogen surveillance. An increasing number of approaches attempt to fix this problem by incorporating diverse concepts and data types together in their models. We aim to contribute to this movement by introducing EpiFusion, a modelling framework that improves the efficiency and precision at which we can monitor important changes in pathogen transmission (specifically, in the effective reproduction number). EpiFusion uses particle filtering to simulate epidemic trajectories over time and weight their likelihood according to both case incidence data and a phylogenetic tree using separate observation models, resulting in the inference of trajectories in agreement with both sets of data. Improvements in our ability to accurately and confidently model pathogen spread help us to respond to infectious disease outbreaks and improve public health.

A Flexible Framework for Local-Level Estimation of the Effective Reproductive Number in Geographic Regions with Sparse Data

A New Framework and Software to Estimate Time-Varying Reproduction Numbers During Epidemics

Generative Bayesian modeling to nowcast the effective reproduction number from line list data with missing symptom onset dates

Early-Phase Local-Area Model for Pandemics Using Limited Data: A SARS-CoV-2 Application

Risk averse reproduction numbers improve resurgence detection

EpiLPS: A fast and flexible Bayesian tool for estimation of the time-varying reproduction number

Refining Estimation of the Instantaneous Reproduction Number During Early Pandemic Stages: Addressing Case-Reporting Variability and Serial Interval Uncertainty

Early real-time estimation of the basic reproduction number of emerging or reemerging infectious diseases in a community with heterogeneous contact pattern: Using data from Hong Kong 2009 H1N1 Pandemic Influenza as an illustrative example

rtestim: Time-varying reproduction number estimation with trend filtering

Likelihood-based Estimation of Spatial Intensity and Variation in Disease Risk from Locations Observed with Error

Adaptive Susceptible-Infectious-Removed Model for Continuous Estimation of the COVID-19 Infection Rate and Reproduction Number in the United States: Modeling Study

The effective reproductive number: Modeling and prediction with application to the multi-wave Covid-19 pandemic

Robust uncertainty quantification in popular estimators of the instantaneous reproduction number

Estimating the effective reproduction number from wastewater (Rt): A methods comparison

A simulation-based approach for estimating the time-dependent reproduction number from temporally aggregated disease incidence time series data

EpiFusion: Joint inference of the effective reproduction number by integrating phylodynamic and epidemiological modelling with particle filtering

Tracking Disease Outbreaks from Sparse Data with Bayesian Inference

Semiparametric inference of effective reproduction number dynamics from wastewater pathogen surveillance data

COVID-19: Nowcasting Reproduction Factors Using Biased Case Testing Data

Risk Estimate under a Nonstationary Autoregressive Model for Data-Driven Reproduction Number Estimation

Temporal evolution of the Covid19 pandemic reproduction number: Estimations from proximal optimization to Monte Carlo sampling