[论文解读] Colored Genealogical Trees and Coalescent Theory

Jianjun Paul Tian, Xiao-Song Lin|arXiv (Cornell University)|Oct 24, 2004

Markov Chains and Monte Carlo Methods参考文献 4被引用 1

一句话总结

本文引入了一种着色的共祖过程，该过程生成具有黑色或白色顶点的随机基因谱系树，颜色变化仅在由参数 x 控制的共祖事件期间发生。当 x = 1/2 时，该模型与一个 Wright-Fisher 单倍体种群模型一致，关键结果表明：达到黑色或白色最近共同祖先（MRAC）的期望时间为 3 − 2/n，而达到任意颜色 MRAC（黑色或白色）的期望时间为 2 − 2/n，与标准 Kingman 共祖过程一致。

ABSTRACT

Abstract. We introduce a colored coalescent process which recovers random colored genealogical trees. Here a colored genealogical tree has its vertices colored black or white. Moving backward along the colored genealogical tree, the color of vertices may change only when two vertice coalesce. The rule that governs the change of color involves a parameter x. When x = 1/2, the colored coalescent process can be derived from a variant of the Wright-Fisher model for a haploid population in population genetics. Explicit computations of the expectation and the cumulative distribution function of the coalescent time are carried out. For example, our calculation shows that when x = 1/2, for a sample of n colored individuals, the expected time for the colored coalescent process to reach a black MRAC or a white MRAC, respectively, is 3−2/n. On the other hand, the expected time for the colored coalescent process to reach a MRAC, either black or white, is 2 − 2/n, which is the same as that for the standard Kingman coalescent process. This colored coalescent process with a color mutation process superimposed is also studied in explicit details. 1.

研究动机与目标

• 开发一种随机过程，用于生成具有动态颜色转换的随机着色基因谱系树。
• 对遗传谱系的演化进行建模，其中谱系颜色代表一种可遗传特征，如突变或等位基因。
• 推导在不同颜色规则下共祖时间的显式表达式，特别是针对 x = 1/2 的情况。
• 将着色共祖过程与已知的种群遗传学模型（如 Wright-Fisher 模型）联系起来。
• 分析叠加在着色共祖过程上的颜色突变过程对共祖动力学的影响。

提出的方法

• 着色共祖过程通过仅在共祖事件期间允许颜色变化来定义，颜色转换的概率由参数 x 控制。
• 该过程被构建为在着色基因谱系空间上的连续时间马尔可夫跳跃过程。
• 使用递归和组合技术，推导出共祖时间的期望值和累积分布函数的显式公式。
• 当 x = 1/2 时，该模型与单倍体种群的 Wright-Fisher 模型变体相联系，确保了生物学上的合理性。
• 叠加的颜色突变过程被建模为共祖树上的额外随机转换层。
• 分析包括使用递归矩计算，推导达到特定颜色或任意颜色最近共同祖先（MRAC）的期望时间。

实验结果

研究问题

• RQ1 参数 x 在颜色转换规则中的变化如何影响达到特定颜色（黑色或白色）最近共同祖先（MRAC）的期望时间？
• RQ2 着色共祖过程达到任一颜色 MRAC 的期望时间是多少？与标准 Kingman 共祖过程相比有何差异？
• RQ3 着色共祖过程能否从已知的种群遗传学模型（如 Wright-Fisher 模型）推导得出？在何种条件下成立？
• RQ4 叠加的颜色突变过程如何改变共祖动力学及预期共祖时间？
• RQ5 着色共祖模型中，共祖时间的累积分布函数和期望值的显式表达式是什么？

主要发现

• 当 x = 1/2 时，着色共祖过程达到黑色最近共同祖先（MRAC）的期望时间为 3 − 2/n（样本量为 n 名个体）。
• 在相同参数设置下，达到白色 MRAC 的期望时间同样为 3 − 2/n。
• 达到任意颜色 MRAC（无论颜色）的期望时间为 2 − 2/n，与标准 Kingman 共祖过程的期望共祖时间一致。
• 当 x = 1/2 时，该模型与 Wright-Fisher 单倍体种群模型一致，验证了其生物学相关性。
• 共祖时间的累积分布函数和期望值被显式计算，使得谱系颜色动态的定量分析成为可能。
• 叠加的颜色突变过程得到了详细研究，表明其如何影响共祖过程及共祖时间的分布。