A computationally tractable birth-death model that combines phylogenetic and epidemiological data

Abstract

AbstractInferring the dynamics of pathogen transmission during an outbreak is an important problem in both infectious disease epidemiology. In mathematical epidemiology, estimates are often informed by time series of confirmed cases, while in phylodynamics genetic sequences of the pathogen, sampled through time, are the primary data source. Each data type provides different, and potentially complementary, insight; recent studies have recognised that combining data sources can improve estimates of the transmission rate and number of infected individuals. However, inference methods are typically highly specialised and field-specific and are either computationally prohibitive or require intensive simulation, limiting their real-time utility.We present a novel birth-death phylogenetic model and derive a tractable analytic approximation of its likelihood, the computational complexity of which is linear in the size of the dataset. This approach combines epidemiological and phylodynamic data to produce estimates of key parameters of transmission dynamics and the number of unreported infections. Using simulated data we show (a) that the approximation agrees well with existing methods, (b) validate the claim of linear complexity and (c) explore robustness to model misspecification. This approximation facilitates inference on large datasets, which is increasingly important as large genomic sequence datasets become commonplace.Author summaryMathematical epidemiologists typically studies time series of cases, ie the epidemic curve, to understand the spread of pathogens. Genetic epidemiologists study similar problems but do so using measurements of the genetic sequence of the pathogen which also contain information about the transmission process. There have been many attempts to unite these approaches so that both data sources can be utilised. However, striking a suitable balance between model flexibility and fidelity, in a way that is computationally tractable, has proven challenging; there are several competing methods but for large datasets they are intractable. As sequencing of pathogen genomes becomes more common, and an increasing amount of epidemiological data is collected, this situation will only be exacerbated. To bridge the gap between the time series and genomic methods we developed an approximation scheme, called TimTam, which can accurately and efficiently estimate key features of an epidemic such as the prevalence of the infection and the effective reproduction number, ie how many people are currently infected and the degree to which the infection is spreading.