Fixation of a Deleterious Allele under Mutation Pressure and Finite Selection Intensity

TL;DR

This paper develops a Markov chain-based method to accurately calculate the fixation time of deleterious alleles under mutation and finite selection, surpassing traditional diffusion approximations especially for larger selection coefficients.

Contribution

It introduces a direct Markov chain approach for fixation time analysis, improving accuracy over diffusion theory for finite selection intensities.

Findings

The method accurately predicts fixation times across dominance scenarios.

Results agree well with numerical simulations.

Outperforms diffusion approximation when N_e s^2 is not small.

Abstract

The mean fixation time of a deleterious mutant allele is studied beyond the diffusion approximation. As in Kimura's classical work [M. Kimura, Proc. Natl. Acad. Sci. U.S.A. Vol.77, 522 (1980)], that was motivated by the problem of fixation in the presence of amorphic or hypermorphic mutations, we consider a diallelic model at a single locus comprising a wild-type A and a mutant allele A' produced irreversibly from A at small uniform rate v. The relative fitnesses of the mutant homozygotes A'A', mutant heterozygotes A'A and wild-type homozygotes AA are 1-s, 1-h and 1, respectively, where it is assumed that v<< s. Here, we adopt an approach based on the direct treatment of the underlying Markov chain (birth-death process) obeyed by the allele frequency (whose dynamics is prescribed by the Moran model), which allows to accurately account for the effects of large fluctuations. After a general description of the theory, we focus on the case of a deleterious mutant allele (i.e. s>0) and discuss three situations: when the mutant is (i) completely dominant (s=h); (ii) completely recessive (h=0), and (iii) semi-dominant (h=s/2). Our theoretical predictions for the mean fixation time and the quasi-stationary distribution of the mutant population in the coexistence state, are shown to be in excellent agreement with numerical simulations. Furthermore, when s is finite, we demonstrate that our results are superior to those of the diffusion theory that is shown to be an accurate approximation only when N_e s^2 << 1, where N_e is the effective population size.