Greenbaum Gili
Department of Solar Energy and Environmental Physics and Mitrani Department of Desert Ecology, The Jacob Blaustein Institutes for Desert Research, Ben-Gurion University of the Negev, Midreshet Ben-Gurion, 84990, Israel.
J Theor Biol. 2015 Sep 7;380:98-102. doi: 10.1016/j.jtbi.2015.05.019. Epub 2015 May 21.
Evaluation of the time scale of the fixation of neutral mutations is crucial to the theoretical understanding of the role of neutral mutations in evolution. Diffusion approximations of the Wright-Fisher model are most often used to derive analytic formulations of genetic drift, as well as for the time scales of the fixation of neutral mutations. These approximations require a set of assumptions, most notably that genetic drift is a stochastic process in a continuous allele-frequency space, an assumption appropriate for large populations. Here equivalent approximations are derived using a coalescent theory approach which relies on a different set of assumptions than the diffusion approach, and adopts a discrete allele-frequency space. Solutions for the mean and variance of the time to fixation of a neutral mutation derived from the two approaches converge for large populations but slightly differ for small populations. A Markov chain analysis of the Wright-Fisher model for small populations is used to evaluate the solutions obtained, showing that both the mean and the variance are better approximated by the coalescent approach. The coalescence approximation represents a tighter upper-bound for the mean time to fixation than the diffusion approximation, while the diffusion approximation and coalescence approximation form an upper and lower bound, respectively, for the variance. The converging solutions and the small deviations of the two approaches strongly validate the use of diffusion approximations, but suggest that coalescent theory can provide more accurate approximations for small populations.
评估中性突变固定的时间尺度对于从理论上理解中性突变在进化中的作用至关重要。赖特 - 费希尔模型的扩散近似最常用于推导遗传漂变的解析公式,以及中性突变固定的时间尺度。这些近似需要一组假设,最显著的是遗传漂变是连续等位基因频率空间中的一个随机过程,这一假设适用于大种群。在此,使用一种溯祖理论方法推导出等效的近似,该方法依赖于与扩散方法不同的一组假设,并采用离散的等位基因频率空间。从这两种方法得出的中性突变固定时间的均值和方差的解对于大种群是收敛的,但对于小种群略有不同。对小种群的赖特 - 费希尔模型进行马尔可夫链分析,以评估所获得的解,结果表明,溯祖方法对均值和方差的近似都更好。与扩散近似相比,溯祖近似对于固定的平均时间代表了一个更紧的上界,而扩散近似和溯祖近似分别对方差形成了一个上界和一个下界。两种方法的收敛解和小偏差有力地验证了扩散近似的使用,但表明溯祖理论可以为小种群提供更准确的近似。