Aaron A. King, Ph.D.

Nelson G. Hairston Collegiate Professor of Ecology, Evolutionary Biology,
Complex Systems, and Mathematics, University of Michigan
External Professor, Santa Fe Institute
Fellow of the American Association for the Advancement of Science

Infectious disease dynamics inferred from genetic data via sequential Monte Carlo

R. A. Smith, E. L. Ionides, and A. A. King
Molecular Biology and Evolution 34(8): 2065–2084, 2017.
A general, statistically efficient, plug-and-play method to jointly estimate both disease transmission and phylogeny using genetic data and, if desired, other epidemiological observations.

Genetic sequences from pathogens can provide information about infectious disease dynamics that may supplement or replace information from other epidemiological observations. Most currently available methods first estimate phylogenetic trees from sequence data, then estimate a transmission model conditional on these phylogenies. Outside limited classes of models, existing methods are unable to enforce logical consistency between the model of transmission and that underlying the phylogenetic reconstruction. Such conflicts in assumptions can lead to bias in the resulting inferences. Here, we develop a general, statistically efficient, plug-and-play method to jointly estimate both disease transmission and phylogeny using genetic data and, if desired, other epidemiological observations. This method explicitly connects the model of transmission and the model of phylogeny so as to avoid the aforementioned inconsistency. We demonstrate the feasibility of our approach through simulation and apply it to estimate stage-specific infectiousness in a subepidemic of human immunodeficiency virus in Detroit, Michigan. In a supplement, we prove that our approach is a valid sequential Monte Carlo algorithm. While we focus on how these methods may be applied to population-level models of infectious disease, their scope is more general. These methods may be applied in other biological systems where one seeks to infer population dynamics from genetic sequences, and they may also find application for evolutionary models with phenotypic rather than genotypic data.


The official version of the paper is here.   Please contact Prof. King if you'd like a reprint.

© 2024 Aaron A. King
3038 Biological Sciences Building
1105 North University Avenue
Ann Arbor MI 48109-1085 USA