[论文解读] The continuous limit of the Moran process and the diffusion of mutant genes in infinite populations

Fabio A. C. C. Chalub, Max O. Souza|ArXiv.org|Feb 23, 2006

Stochastic processes and statistical mechanics参考文献 37被引用 7

一句话总结

本文推导了无限种群中频率依赖性 Moran 过程的连续极限，表明离散随机动力学在漂移-扩散标度下收敛于一个退化的抛物型偏微分方程（PDE）。关键结果是该 PDE 推广了 Kimura 的固定方程，并在频率无关适应度情况下形式上恢复该方程，且已证明在大种群下解的存在性与收敛性。

ABSTRACT

We consider the so called Moran process with frequency dependent fitness given by a certain pay-off matrix. For finite populations, we show that the final state must be homogeneous, and show how to compute the fixation probabilities. Next, we consider the infinite population limit, and discuss the appropriate scalings for the drift-diffusion limit. In this case, a degenerated parabolic PDE is formally obtained that, in the special case of frequency independent fitness, recovers the celebrated Kimura equation in population genetics. We then show that the corresponding initial value problem is well posed and that the discrete model converges to the PDE model as the population size goes to infinity. We also study some game-theoretic aspects of the dynamics and characterize the best strategies, in an appropriate sense.

研究动机与目标

• 建立具有频率依赖适应度的有限种群 Moran 过程的热力学极限。
• 识别适当的标度（漂移与漂移-扩散），以在大种群下获得连续 PDE 近似。
• 证明所得退化抛物型 PDE 的适定性，并证明离散模型在无限种群极限下收敛于 PDE。
• 阐明演化稳定策略（ESS）与大有限种群中固定动力学之间的关系。
• 在频率无关适应度条件下，将推导出的 PDE 作为特例恢复 Kimura 的经典方程。

提出的方法

• 使用带有收益矩阵的 Moran 过程来建模有限种群中大小为 $N$ 的频率依赖性选择。
• 应用两种时间标度方案：$\Delta t = 1/N$（漂移极限）和 $\Delta t = 1/N^2$（漂移-扩散极限），后者采用弱选择。
• 推导出无限种群极限下的退化抛物型 PDE，其中扩散系数在边界 $x=0$ 和 $x=1$ 处消失。
• 采用形式渐近展开与弱形式化方法推导 PDE 并验证边界条件。
• 通过分部积分与测试函数分析，确定边界项的结构并确认 PDE 的一致性。
• 证明初值问题的适定性，并证明当 $N \to \infty$ 时，离散过程收敛于 PDE。

实验结果

研究问题

• RQ1 具有频率依赖适应度的有限种群 Moran 过程在无限种群极限下如何收敛于连续 PDE？
• RQ2 需要何种时间标度才能获得非平凡的扩散极限，其对所得 PDE 有何影响？
• RQ3 在何种条件下，所推导的 PDE 会退化为 Kimura 的经典方程（用于突变基因固定的方程）？
• RQ4 连续复制者动力学中的演化稳定策略（ESS）与大有限种群中的固定概率之间有何关系？
• RQ5 边界扩散在确保固定或丢失中的作用是什么，其在退化 PDE 中如何被捕捉？

主要发现

• 有限种群 Moran 过程在无限种群极限下产生一个退化抛物型 PDE，用于建模突变基因频率随时间的概率分布。
• 在漂移-扩散标度（$\Delta t = 1/N^2$）和弱选择下，PDE 同时捕捉了漂移与扩散效应，且扩散系数在 $x=0$ 和 $x=1$ 处消失。
• 在频率无关适应度情况下，PDE 形式上退化为 Kimura 的方程，恢复了著名的固定概率前向-后向方程对。
• 所推导 PDE 的初值问题适定，确保在适当条件下解的存在性与唯一性。
• 当种群大小 $N \to \infty$ 时，离散 Moran 过程收敛于 PDE 模型，确立了连续近似的有效性。
• 在大有限种群中表现最佳的策略对应于演化稳定策略（ESS），其与复制者方程的动力学行为一致。