The Rate of Molecular Evolution When Mutation May Not Be Weak

Abstract

One of the most fundamental rules of molecular evolution is that the rate of neutral evolution equals the mutation rate and is independent of effective population size. This result lies at the heart of the Neutral Theory, and is the basis for numerous analytic approaches that are widely applied to infer the action of natural selection across the genome and through time, and for dating divergence events using the molecular clock. However, this result was derived under the assumption that evolution is strongly mutation-limited, and it has not been known whether it generalizes across the range of mutation pressures or the spectrum of mutation types observed in natural populations. Validated by both simulations and exact computational analyses, we present a direct and transparent theoretical analysis of the Wright-Fisher model of population genetics, which shows that some of the most important rules of molecular evolution are fundamentally changed by considering recurrent mutations full effect. Surprisingly, the rate of the neutral molecular clock is found to have population-size dependence and to not equal the mutation rate in general. This is because, for increasing values of the population mutation rate parameter ({theta}), the time spent waiting for mutations quickly becomes smaller than the cumulative time mutants spend segregating before a substitution, resulting in a net deceleration compared to classical theory that depends on the population mutation rate. Furthermore, selection exacerbates this effect such that more adaptive alleles experience a greater deceleration than less adaptive alleles, introducing systematic bias in a wide variety of methods for inferring the strength and direction of natural selection from across-species sequence comparisons. Critically, the classical weak mutation approximation performs well only when{theta} < 0.1, a threshold that many biological populations seem to exceed.