[论文解读] Lotka-Volterra population model of genetic evolution

Mirosław R. Dudek|ArXiv.org|Jan 20, 2007

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

一句话总结

本文提出了一种确定性的年龄结构Lotka-Volterra模型，通过将缺陷基因比例编码为年龄特异性的人口增长参数，复现了Penna模型的遗传演化动态。结果表明，两个模型中均出现相同的年龄特异性人口分布，从而可在无需蒙特卡洛模拟的情况下，对多物种生态系统进行解析研究，关键发现显示，生育年龄和缺陷基因分布决定了平衡状态的稳定性或振荡行为。

ABSTRACT

A deterministic model of an age-structured population with genetics analogous to the discrete time Penna model of genetic evolution is constructed on the basis of the Lotka-Volterra scheme. It is shown that if, as in the Penna model, genetic information is represented by the fraction of defective genes in the population, the population numbers for each specific individual's age are represented by exactly the same functions of age in both models. This gives us a new possibility to consider multi-species evolution without using detailed microscopic Penna model. We discuss a particular case of the predator-prey system representing an ecosystem consisting of a limited amount of energy resources consumed by the age-structured species living in this ecosystem. Then, the increase in number of the individuals in the population under consideration depends on the available energy resources, the shape of the distribution function of defective genes in the population and the fertility age. We show that these parameters determine the trend toward equilibrium of the whole ecosystem.

研究动机与目标

• 开发一种确定性的、年龄结构化的人口模型，以捕捉与随机Penna模型相似的遗传演化动态。
• 用可解析求解的Lotka-Volterra方程替代Penna模型计算成本高昂的蒙特卡洛模拟。
• 利用微分方程研究具有遗传约束的多物种生态系统。
• 探究生育年龄和缺陷基因分布如何影响人口平衡与稳定性。

提出的方法

• 该模型使用基于Lotka-Volterra框架推导出的年龄特异性人口增长方程，其中出生率受各年龄阶段缺陷基因比例的调制。
• 缺陷基因比例εₐ在年龄≤a_F时建模为类似幂函数的形式，在年龄>a_F时建模为类似指数函数的形式，以匹配经验性的Penna模型输出。
• 人口动态由微分方程控制，其增长率取决于前一个年龄组、存活概率以及能量可用性S(t)。
• 通过引入物种特异性参数（a_F, a_D, εₐ）并引入对共享能量资源的竞争，将系统扩展至多物种情形。
• 利用所得微分方程的数值解分析稳定性、平衡状态及振荡行为。
• 模型中引入扰动（如人口灾难）以研究生态系统的恢复力与恢复动力学。

实验结果

研究问题

• RQ1 能否通过将遗传信息编码为年龄特异性缺陷基因比例的确定性Lotka-Volterra框架，精确复现随机Penna模型的年龄特异性人口分布？
• RQ2 生育年龄a_F和缺陷基因分布εₐ的形状变化如何影响人口动态的长期稳定性？
• RQ3 在资源受限的生态系统中，是什么决定了人口收敛至固定点还是表现出振荡行为？
• RQ4 突发的人口扰动（如灾难）如何影响不同物种中人口波动的振幅与稳定性？
• RQ5 竞争物种的年龄结构在决定多物种生态系统中灭绝或共存方面起什么作用？

主要发现

• 当使用相同的缺陷基因比例εₐ时，确定性Lotka-Volterra模型中的年龄特异性人口分布与Penna模型的平衡态完全一致。
• 该系统根据生育年龄a_F和εₐ的形状，表现出收敛至稳定不动点或持续振荡的特性。
• 当生育年龄从a_F = 15提升至a_F = 20时，系统发生定性转变，从稳定平衡转为振荡动力学，可能引发灭绝。
• 在某些情况下，人口灾难可降低振荡振幅，表明生态系统具有自调节行为。
• 最大寿命a_D的微小变化（如从23增至24）可决定两个竞争物种中哪一个能够存活，凸显年龄结构对竞争结果的高度敏感性。
• 该模型表明，生育年龄a_F较低的物种表现出更大的人口波动，使其在扰动下面临更高的灭绝风险。