Uncovering epidemiological dynamics in heterogeneous host populations using phylogenetic methods

Host population structure has a major influence on epidemiological dynamics. However, in particular for sexually transmitted diseases, quantitative data on population contact structure are hard to obtain. Here, we introduce a new method that quantifies host population structure based on phylogenetic trees, which are obtained from pathogen genetic sequence data. Our method is based on a maximum-likelihood framework and uses a multi-type branching process, under which each host is assigned to a type (subpopulation). In a simulation study, we show that our method produces accurate parameter estimates for phylogenetic trees in which each tip is assigned to a type, as well for phylogenetic trees in which the type of the tip is unknown. We apply the method to a Latvian HIV-1 dataset, quantifying the impact of the intravenous drug user epidemic on the heterosexual epidemic (known tip states), and identifying superspreader dynamics within the men-having-sex-with-men epidemic (unknown tip states).

[1]  Richard G FitzJohn,et al.  Estimating trait-dependent speciation and extinction rates from incompletely resolved phylogenies. , 2009, Systematic biology.

[2]  R. Ribeiro,et al.  Agent-based and phylogenetic analyses reveal how HIV-1 moves between risk groups: injecting drug users sustain the heterosexual epidemic in Latvia. , 2012, Epidemics.

[3]  A. Phillimore,et al.  Diversity-dependence brings molecular phylogenies closer to agreement with the fossil record , 2012, Proceedings of the Royal Society B: Biological Sciences.

[4]  P. J. Hudson,et al.  Large Shifts in Pathogen Virulence Relate to Host Population Structure , 2004, Science.

[5]  O. Pybus,et al.  The Epidemic Behavior of the Hepatitis C Virus , 2001, Science.

[6]  T. Stadler.  Sampling-through-time in birth-death trees. , 2010, Journal of theoretical biology.

[7]  E. Goldberg,et al.  TEMPO AND MODE IN PLANT BREEDING SYSTEM EVOLUTION , 2012, Evolution; international journal of organic evolution.

[8]  M. Keeling,et al.  Modeling Infectious Diseases in Humans and Animals , 2007 .

[9]  Matthias Cavassini,et al.  Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland. , 2010, The Journal of infectious diseases.

[10]  M. Suchard,et al.  Phylogeography takes a relaxed random walk in continuous space and time. , 2010, Molecular biology and evolution.

[11]  David A. Rasmussen,et al.  Inference for Nonlinear Epidemiological Models Using Genealogies and Time Series , 2011, PLoS Comput. Biol..

[12]  Alexei J. Drummond,et al.  Bayesian Phylogeography Finds Its Roots , 2009, PLoS Comput. Biol..

[13]  A. Rambaut,et al.  BEAST: Bayesian evolutionary analysis by sampling trees , 2007, BMC Evolutionary Biology.

[14]  Peter E Midford,et al.  Estimating a binary character's effect on speciation and extinction. , 2007, Systematic biology.

[15]  Beda Joos,et al.  Estimating the basic reproductive number from viral sequence data. , 2012, Molecular biology and evolution.

[16]  Tanja Stadler,et al.  Mammalian phylogeny reveals recent diversification rate shifts , 2011, Proceedings of the National Academy of Sciences.

[17]  M. Pascual,et al.  Global Migration Dynamics Underlie Evolution and Persistence of Human Influenza A (H3N2) , 2010, PLoS pathogens.

[18]  Daniel L Rabosky,et al.  LIKELIHOOD METHODS FOR DETECTING TEMPORAL SHIFTS IN DIVERSIFICATION RATES , 2006, Evolution; international journal of organic evolution.

[19]  T. Leitner,et al.  Phylogenetic analysis of the Latvian HIV-1 epidemic. , 2012, AIDS research and human retroviruses.

[20]  S. Otto,et al.  Linking the Investigations of Character Evolution and Species Diversification , 2012, The American Naturalist.

[21]  Kazutaka Katoh,et al.  Recent developments in the MAFFT multiple sequence alignment program , 2008, Briefings Bioinform..

[22]  R. May,et al.  Infectious Diseases of Humans: Dynamics and Control , 1991, Annals of Internal Medicine.

[23]  Erik M. Volz,et al.  Complex Population Dynamics and the Coalescent Under Neutrality , 2012, Genetics.

[24]  Edward C. Holmes,et al.  Endemic Dengue Associated with the Co-Circulation of Multiple Viral Lineages and Localized Density-Dependent Transmission , 2011, PLoS pathogens.

[25]  S. Bonhoeffer,et al.  Birth–death skyline plot reveals temporal changes of epidemic spread in HIV and hepatitis C virus (HCV) , 2012, Proceedings of the National Academy of Sciences.

[26]  Tanja Stadler,et al.  Simulating trees with a fixed number of extant species. , 2011, Systematic biology.

[27]  T. Leitner,et al.  Rapid epidemic spread of HIV type 1 subtype A1 among intravenous drug users in Latvia and slower spread of subtype B among other risk groups. , 2004, AIDS research and human retroviruses.

[28]  T. Stadler.  How can we improve accuracy of macroevolutionary rate estimates? , 2013, Systematic biology.

[29]  M. Nowak,et al.  HOST POPULATION STRUCTURE AND THE EVOLUTION OF VIRULENCE: A “LAW OF DIMINISHING RETURNS” , 1995, Evolution; international journal of organic evolution.

[30]  Gavin J. D. Smith,et al.  Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic , 2009, Nature.

[31]  Sergei L. Kosakovsky Pond,et al.  Phylodynamic Reconstruction Reveals Norovirus GII.4 Epidemic Expansions and their Molecular Determinants , 2010, PLoS pathogens.

[32]  M. Keeling.  The implications of network structure for epidemic dynamics. , 2005, Theoretical population biology.

[33]  Sergei L. Kosakovsky Pond,et al.  Phylodynamics of Infectious Disease Epidemics , 2009, Genetics.

[34]  Huldrych F. Günthard,et al.  Inferring Epidemic Contact Structure from Phylogenetic Trees , 2012, PLoS Comput. Biol..