Exact Markov chain and approximate diffusion solution for haploid genetic drift with one-way mutation

Ola Hössjer, Peder A. Tyvand, Touvia Miloh

Research output: Contribution to journalArticlepeer-review

Abstract

The classical Kimura solution of the diffusion equation is investigated for a haploid random mating (Wright-Fisher) model, with one-way mutations and initial-value specified by the founder population. The validity of the transient diffusion solution is checked by exact Markov chain computations, using a Jordan decomposition of the transition matrix. The conclusion is that the one-way diffusion model mostly works well, although the rate of convergence depends on the initial allele frequency and the mutation rate. The diffusion approximation is poor for mutation rates so low that the non-fixation boundary is regular. When this happens we perturb the diffusion solution around the non-fixation boundary and obtain a more accurate approximation that takes quasi-fixation of the mutant allele into account. The main application is to quantify how fast a specific genetic variant of the infinite alleles model is lost. We also discuss extensions of the quasi-fixation approach to other models with small mutation rates.

Original languageEnglish
Pages (from-to)100-112
Number of pages13
JournalMathematical Biosciences
Volume272
DOIs
StatePublished - 2016

Keywords

  • Diffusion
  • Jordan decomposition
  • Markov chain
  • Mutation
  • Quasi-fixation
  • Wright-Fisher model

Fingerprint

Dive into the research topics of 'Exact Markov chain and approximate diffusion solution for haploid genetic drift with one-way mutation'. Together they form a unique fingerprint.

Cite this