An approximate markov model for the wright–fisher diffusion and its application to time series data

Anna Ferrer-Admetlla, Christoph Leuenberger, Jeffrey D. Jensen, Daniel Wegmann

Research output: Contribution to journalArticlepeer-review

23 Scopus citations

Abstract

The joint and accurate inference of selection and demography from genetic data is considered a particularly challenging question in population genetics, since both process may lead to very similar patterns of genetic diversity. However, additional information for disentangling these effects may be obtained by observing changes in allele frequencies over multiple time points. Such data are common in experimental evolution studies, as well as in the comparison of ancient and contemporary samples. Leveraging this information, however, has been computationally challenging, particularly when considering multilocus data sets. To overcome these issues, we introduce a novel, discrete approximation for diffusion processes, termed mean transition time approximation, which preserves the long-term behavior of the underlying continuous diffusion process. We then derive this approximation for the particular case of inferring selection and demography from time series data under the classic Wright–Fisher model and demonstrate that our approximation is well suited to describe allele trajectories through time, even when only a few states are used. We then develop a Bayesian inference approach to jointly infer the population size and locus-specific selection coefficients with high accuracy and further extend this model to also infer the rates of sequencing errors and mutations. We finally apply our approach to recent experimental data on the evolution of drug resistance in influenza virus, identifying likely targets of selection and finding evidence for much larger viral population sizes than previously reported.

Original languageEnglish (US)
Pages (from-to)831-846
Number of pages16
JournalGenetics
Volume203
Issue number2
DOIs
StatePublished - Jun 2016
Externally publishedYes

Keywords

  • Diffusion approximation
  • Discrete markov model
  • Hidden markov model
  • Time-series data
  • Wright–fisher model

ASJC Scopus subject areas

  • Genetics

Fingerprint Dive into the research topics of 'An approximate markov model for the wright–fisher diffusion and its application to time series data'. Together they form a unique fingerprint.

Cite this