A study on a nearly neutral mutation model in finite populations

Abstract

As a nearly neutral mutation model, the house-of-cards model is studied in finite populations using computer simulations. The distribution of the mutant effect is assumed to be normal. The behavior is mainly determined by the product of the population size, N, and the standard deviation, sigma, of the distribution of the mutant effect. If 4N sigma is large compared to one, a few advantageous mutants are quickly fixed in early generations. Then most mutation becomes deleterious and very slow increase of the average selection coefficient follows. It takes very long for the population to reach the equilibrium state. Substitutions of alleles occur very infrequently in the later stage. If 4N sigma is the order of one or less, the behavior is qualitatively similar to that of the strict neutral case. Gradual increase of the average selection coefficient occurs and in generations of several times the inverse of the mutation rate the population almost reaches the equilibrium state. Both advantageous and neutral (including slightly deleterious) mutations are fixed. Except in the early stage, an increase of the standard deviation of the distribution of the mutant effect decreases the average heterozygosity. The substitution rate is reduced as 4N sigma is increased. Three tests of neutrality, one using the relationship between the average and the variance of heterozygosity, another using the relationship between the average heterozygosity and the average number of substitutions and Watterson's homozygosity test are applied to the consequences of the present model. It is found that deviation from the neutral expectation becomes apparent only when 4N sigma is more than two. Also a simple approximation for the model is developed which works well when the mutation rate is very small.