[论文解读] Fixation probability and fixation time in structured populations

Josef Tkadlec, Andreas Pavlogiannis+2 more|arXiv (Cornell University)|Sep 27, 2018

Evolution and Genetic Dynamics参考文献 44被引用 4

一句话总结

本文提出了一类结构化种群模型，可在保持高固定概率的同时缩短固定时间，证明了某些二分图可在固定时间仅适度增加的情况下放大选择效应。其主要贡献在于提出了一类种群结构家族，当突变率较低时，通过优化固定概率与固定时间之间的权衡，实现了比充分混合种群更高的有效进化速率。

ABSTRACT

The rate of biological evolution depends on the fixation probability and on the fixation time of new mutants. Intensive research has focused on identifying population structures that augment the fixation probability of advantageous mutants. But these `amplifiers of natural selection' typically increase fixation time. Here we study population structures that achieve a trade-off between high fixation probability and short fixation time. First, we show that no amplifiers can have asymptotically lower absorption time than the well-mixed population. Then we design population structures that substantially augment the fixation probability with just a minor increase in fixation time. Finally, we show that those structures enable higher effective rate of evolution than the well-mixed population provided that the rate of generating advantageous mutants is relatively low. Our work sheds light on how population structure affects the rate of evolution. Moreover, our structures could be useful for lab-based, medical or industrial applications of evolutionary optimization.

研究动机与目标

• 识别在不显著增加固定时间的前提下提升固定概率的种群结构。
• 解决在进化动力学中放大选择与减缓固定之间的权衡问题。
• 设计通过平衡固定概率与固定时间以最大化有效进化速率的种群结构。
• 利用图论模型，为结构化种群中的吸收时间与固定时间提供解析边界。
• 证明在低突变率下，某些结构化种群可在进化优化中优于充分混合种群。

提出的方法

• 使用进化图论将种群结构建模为加权图，其中顶点代表个体，边代表繁殖连通性。
• 引入α-平衡与α-加权二分图，系统性地调节固定概率与固定时间。
• 在图上应用出生-死亡过程，建模具有依赖适应度的繁殖与死亡率的莫兰过程。
• 利用递推方程与调和数近似，推导固定概率、吸收时间与条件固定时间的解析表达式。
• 采用渐近分析与极限计算（如调和级数、对数项），推导大N下固定动力学的行为。
• 通过结合固定概率与固定时间来比较有效进化速率，使用固定概率与固定时间倒数的乘积作为度量指标。

实验结果

研究问题

• RQ1 是否存在既能提高固定概率又能最小化固定时间的种群结构？
• RQ2 在结构化种群中，固定概率与固定时间之间存在何种根本性权衡？
• RQ3 是否存在有效进化速率高于充分混合种群的种群结构？
• RQ4 能否设计出仅使固定时间略有增加的选择放大器？
• RQ5 不同的初始化方案（均匀分布与温度分布）如何影响结构化种群中的固定动力学？

主要发现

• 没有任何选择放大器的吸收时间能渐近低于充分混合种群，确立了固定时间的根本下界。
• α-加权二分图可实现固定概率趋近于 1 - 1/r²，且与充分混合情况相比，固定时间仅呈对数增长。
• 在低突变率下，由于固定概率与时间之间有利的权衡，所提结构的有效进化速率超过充分混合种群。
• 所提结构中的固定时间呈 O(N log N) 标度，与充分混合种群一致，而固定概率则显著增强。
• 所提模型中的灭绝时间与 N 呈线性关系，与中性及有利突变体的理论预期一致。
• 解析推导证实，α-加权二分图的固定时间约为 (r+1)/(r-1) · (N-1)H_{N-1} + o(N)，显示出接近最优的时间标度。