[论文解读] Moment closure in a Moran model with recombination

Ellen Baake, Thiemo Hustedt|arXiv (Cornell University)|May 4, 2011

Advanced Thermodynamics and Statistical Mechanics参考文献 14被引用 3

一句话总结

该论文将具有重组的Moran模型扩展至一般重组和突变，证明在无遗传漂变时，矩层次的动力学可在有限步内精确闭合。令人惊讶的是，当引入重采样（遗传漂变）后，这种精确矩闭合被破坏，尽管系统结构仍保持线性。

ABSTRACT

We extend the Moran model with single-crossover recombination to include general recombination and mutation. We show that, in the case without resampling, the expectations of products of marginal processes defined via partitions of sites form a closed hierarchy, which is exhaustively described by a finite system of differential equations. One thus has the exceptional situation of moment closure in a nonlinear system. Surprisingly, this property is lost when resampling (i.e., genetic drift) is included.

研究动机与目标

• 研究在具有广义重组和突变的随机Moran模型中，矩闭合是否会发生。
• 确定矩动力学层次是否在有限步后闭合，从而避免典型非线性系统中出现的无限层次。
• 考察遗传漂变（重采样）对矩闭合的影响，特别是其对连锁不平衡和相关性衰减的影响。
• 探索在随机种群遗传模型中，精确矩闭合得以保持或丧失的结构条件。

提出的方法

• 作者使用连续时间马尔可夫过程对种群进行建模，其状态空间由位点上的等位基因和重组事件构成。
• 通过位点的划分定义子过程，并追踪这些过程乘积的期望。
• 利用马尔可夫过程的生成元，推导出矩动力学的耦合常微分方程组（ODEs）。
• 分析表明，在无重采样时，由于过程的划分结构，矩层次在有限步后闭合。
• 当引入重采样时，推导显示高阶矩持续出现，破坏了闭合性。
• 该方法明确比较了在重组和重采样下相关函数（如连锁不平衡）的动力学，表明其呈指数衰减。

实验结果

研究问题

• RQ1 在Moran模型中引入广义重组是否会导致矩动力学的有限、闭合ODE系统？
• RQ2 即使重组引起非线性相互作用，能否在无遗传漂变时实现精确矩闭合？
• RQ3 引入重采样（遗传漂变）如何影响矩层次的闭合性？
• RQ4 基于划分的子过程在实现精确矩闭合中起什么作用？
• RQ5 为何在加入重采样后矩闭合性丧失，尽管重组和重采样在结构上均为线性？

主要发现

• 在无重采样时，子过程乘积的矩层次在有限步后闭合，从而得到有限ODE系统。
• 当引入突变时，这种精确矩闭合仍被保持，但过程不再能用位点划分描述。
• 引入重采样会破坏矩闭合，因为求导后持续产生高阶矩。
• 连锁不平衡项 $\mathbb{E}[1,2]_tN - \mathbb{E}[1]_t\mathbb{E}[2]_t$ 在重组和重采样下均呈指数衰减。
• 对于三基因位点系统，$\mathbb{E}[1]_t[2]_t[3]_t$ 的导数包含如 $\mathbb{E}[1]^2_t\mathbb{E}[1,2,3]_t$ 的项，这些项生成更高阶矩，证明层次无法闭合。
• 该系统表现出Haldane线性化的随机类比，但仅在无遗传漂变时成立。