[论文解读] A Hamilton-Jacobi approach to characterize the evolutionary equilibria in heterogeneous environments

Sepideh Mirrahimi|arXiv (Cornell University)|Dec 19, 2016

Evolution and Genetic Dynamics参考文献 32被引用 26

一句话总结

本文提出一种带有校正项的哈密顿-雅可比方法，用于表征在两个异质栖息地中，突变、选择与迁移作用下表型结构化种群的进化均衡。该方法提供了单态或二态的显式条件，并计算了超出高斯假设的高阶近似，为经典数量遗传学模型提供了更严谨的数学基础。

ABSTRACT

In this work, we characterize the solution of a system of elliptic integro-differential equations describing a phenotypically structured population subject to mutation, selection and migration between two habitats. Assuming that the effects of the mutations are small but nonzero, we show that the population's distribution has at most two peaks and we give explicit conditions under which the population will be monomorphic (unimodal distribution) or dimorphic (bimodal distribution). More importantly, we provide a general method to determine the dominant terms of the population's distribution in each case. Our work, which is based on Hamilton-Jacobi equations with constraint, goes further than previous works where such tools were used, for different problems from evolutionary biology, to identify the asymptotic solutions, while the mutations vanish, as a sum of Dirac masses. In order to extend such results to the case with non-vanishing effects of mutations, the main elements are a uniqueness property and the computation of the correctors. This method allows indeed to go further than the Gaussian approximation commonly used by biologists and makes a connection between the theories of adaptive dynamics and quantitative genetics. Our work being motivated by biological questions, the objective of this article is to provide the mathematical details which are necessary for our biological results [16].

研究动机与目标

• 为在突变、选择与迁移作用下，异质环境中表型分布的数学表征提供一个框架。
• 通过证明唯一性并计算下一阶校正项，将哈密顿-雅可比方法扩展至非渐近极限（突变趋于零）的情形。
• 为数量遗传学中常用的高斯近似提供一个严格替代方案，从而实现对进化结果更精确的预测。
• 确定在双栖息地模型中，种群分布呈现单态（单峰）或二态（双峰）的显式条件。
• 通过推导非渐近解，弥合自适应动力学与数量遗传学之间的鸿沟，适用于非零突变效应的情形。

提出的方法

• 建立一个椭圆型积分微分方程组，用于建模两个栖息地之间种群密度、增长率与迁移的动态。
• 应用对数变换（Hopf-Cole变换），将种群密度方程转化为带有约束的哈密顿-雅可比方程。
• 建立黏性解框架，并在非零突变效应下证明哈密顿-雅可比方程黏性解的唯一性。
• 在对数种群密度展开中，计算一阶与二阶校正项（v_i, w_i），以捕捉高阶修正。
• 通过在均衡点附近进行泰勒展开，推导出种群大小修正项（K_i）与校正函数的显式表达式。
• 通过分析极限分布的支集以及哈密顿-雅可比解在临界点的行为，推导出单态与二态的条件。

实验结果

研究问题

• RQ1 在两个异质栖息地中，表型结构化种群在何种条件下会进化为单态或二态分布？
• RQ2 如何将哈密顿-雅可比方法扩展至非零突变极限，以描述非渐近进化均衡？
• RQ3 当突变效应非零时，种群分布中下一阶修正项（校正项）的数学结构是什么？
• RQ4 迁移率与栖息地特异性选择强度如何影响二态性的出现？
• RQ5 能否利用此哈密顿-雅可比框架系统性地校正数量遗传学中的高斯近似？

主要发现

• 种群分布最多具有两个峰，单态或二态的条件由涉及迁移率、选择强度与竞争强度的不等式显式表征。
• 该方法为带有约束的哈密顿-雅可比方程黏性解提供了唯一性结果，从而支持严格的高阶近似。
• 校正项 v_i 与 w_i 被显式计算，使得种群大小修正项 K_i 可精确推导至 O(ε²) 阶。
• 在单态情形下，种群大小修正项 K_2 表示为 K_2 = N_2^{M*} ( (E_2 + 0.5 D_2²)/√g_1 + F_2 )，其中 F_2 通过方程组确定。
• 该方法成功计算了种群分布的主导项，无需假设高斯性，从而修正了进化生物学中先前使用的近似。
• 分析确认，当 ESS 支持两种不同性状时，二态性出现；通过在临界点（如 z = -θ）进行一致性检验，验证了解的结构。