生态进化建模、适合度与Lotka质量猜想：互惠共生进化最大化指标的检验

《Journal of Theoretical Biology》：Eco-evolutionary modelling, fitness and Lotka’s mass conjectures

　　

研究人员考虑一个互惠共生（mutualism）的生态进化模型（eco-evolutionary model）的结果，以检验三个进化指标的作用：总适合度（total fitness）以及Lotka猜想（Lotka's conjectures），即生态系统可能进化以最大化其总生物量（total biomass）和通过系统的物质总通量（total flux of mass）。本文提供了支持适合度和生物量指标的分析论证。结果表明，总适合度和Lotka总质量猜想是互惠共生模型进化结果的有用指标。这些生态系统层面属性的最大化与每个种群在另一种群进化所施加的约束内最大化该指标是一致的。研究人员注意到，这些指标不一定是全局属性，因为独特的生态结果可以通过多个进化结果实现，每个进化结果在其自身域内最大化生物量、循环和适合度。本文提供了一个例子，其中模型的解空间被划分为不同区域，每个区域在其最大值处具有一个稳定的生态进化平衡。每个区域代表兼性（facultative）和/或专性（obligate）互惠共生的不同组合，但都与相同的生态结果相关。在这个简单模型中，实现哪一个解取决于种群的初始进化状态。

论文解读：生态进化建模中的适合度与Lotka质量猜想

一、研究背景与目的

最大化原则在生态学中具有长期研究历史。尽管适合度景观（fitness landscape）是理解进化过程的核心概念，但许多已提出的进化指标并未能真正揭示自然选择机制的运作方式。Lotka（1922）提出早期猜想：生态系统将通过自然选择进化以最大化其总生物量、物质循环速率和能量通量。然而，这些猜想在复杂生态进化系统中的适用性，尤其在互惠共生（mutualism）这一最具挑战性的相互作用形式中，尚未得到充分检验。互惠共生建模长期受困于早期模型产生的“互惠盛宴”（orgy of mutual benefaction）等不现实行为，而采用显式资源核算的方法能够生成符合现实的动态。因此，研究人员利用生态进化建模（eco-evolutionary modelling）框架，构建了一个两个自养种群互惠共生的模型，以检验总生物量（total biomass，BIO）、总物质循环速率（total mass recycling rate，CYC）和总适合度（total fitness，FIT）三个指标对进化结果的预测能力。

二、研究内容与主要结论

研究人员基于Cropp和Norbury建立的生态进化建模框架，构建了一个以有限循环限制性营养盐总量N_T为货币的两种群模型。每个种群具有β分布的性状，其性状均值γ_i和性状分化Q_i随进化动态调整，从而改变互惠共生的义务（obligation）、收益（benefit）和成本（cost）。研究通过解析论证支持总生物量和总适合度指标，并数值求解模型，评估三个指标在稳定生态进化平衡点处的表现。结果表明，总适合度、总生物量和总循环速率的最大化均与稳定平衡点一致，支持Lotka质量猜想和适合度作为进化指南。然而，模型解空间可能被分隔线（separatrices）分割，导致指标可能仅为局部属性；同时，多个不同的进化结果可以对应同一个生态结果。

三、主要技术方法

研究使用了Cropp和Norbury（2022, 2024）的生态进化建模框架，将生态与进化动态耦合为常微分方程组（ODE）。模型以系统中循环的限制性营养盐总量N_T为恒定环境背景，种群生物量表示为该营养盐的质量分数以满足资源守恒。数值求解采用Matlab中的四阶-五阶Runge-Kutta方法（ode45），绝对和相对误差设为机器精度。通过解空间截面、零等斜线、平衡点分析、相关分析以及1,600个不同初始条件的积分试验，研究系统动力学与分隔线作用。参数选取两组互惠成本系数χ，使其分别位于分岔两侧。无样本队列来源。

四、研究结果

5.1 简单生态进化结果

当成本系数取值较高（χ₁=18=χ₂）时，系统仅存在一个稳定内部生态进化平衡点，两个种群均进化为专性互惠（obligate mutualism）。通过向量场和平衡面分析发现，总生物量、总循环速率和总适合度三个指标的最大值均位于同一稳定平衡点{γ₁=0.28, γ₂=0.43}，且三个平衡面高度相关（生物量与循环的相关系数r>0.96，适合度与二者的相关系数r≈0.90）。该结果从数值上支持了Lotka关于总生物量和物质循环的猜想以及适合度最大化的观点。

5.2 多重进化结果

将成本系数降低至χ₁=8、χ₂=13，系统跨越分岔，出现四个稳定内部生态进化平衡点，分别代表专性-专性、兼性-兼性（facultative-facultative）、专性-兼性和兼性-专性互惠组合。1600次积分显示，四个进化结果对应几乎相同的生态种群平衡状态，表明单一生态结果可由多个进化策略实现。分隔线将进化解空间划分为四个独立子空间，各子空间内的初始条件只能到达对应的局部最大点，无法跨越分隔线。因此，进化结果取决于初始平均性状值，指标最大化仅为局部属性。

5.3 分岔细节

通过改变成本系数并观察零等斜线的相切关系，确定分岔临界值约为χ₁≈8.25、χ₂≈14.25，并识别出鞍结分岔（saddle-node bifurcation）。在接近分岔的参数下，适合度景观已开始弯曲为四个区域，但分隔线尚未完全形成。短时间积分（如1,000代）可能产生假的多稳态迹象，而长时间积分（超过9,000代）显示系统仍只存在一个稳定平衡点。该结果提示，在分岔附近进行进化模拟时需足够长的积分时间，以避免错误结论。

五、讨论与结论

讨论指出，本模型所最大化的三个指标均为生态系统层面属性，单个种群不一定最大化其自身生物量、循环或适合度，而是在其他种群进化约束下达到折中。Lotka的总生物量最大化在质量封闭系统中等价于最小化未利用资源，称为N^*理论，与Tilman（1982）的R^*竞争理论一致，架起了生态与进化属性的桥梁。分隔线的存在则警示，生物量、循环或适合度的峰值可能是局部而非全局属性，初始状态决定系统最终到达哪个吸引域。此外，多个进化性状组合可产生同一生态结果，说明仅凭最大化指标难以预测具体进化路径，需结合零等斜线、平衡点和适合度景观进行综合判断。

研究结论可翻译如下：研究人员分析了一个具有解析可处理性的简单生态进化系统，该系统涉及互惠共生这一“经典”生态相互作用，并以限制性营养盐作为恒定的非生物环境。结果表明，一些关键的生态系统层面指标可能被进化最大化：总适合度、总生物量和总循环速率。Lotka猜想中这些质量方面的最大化，或可为更常用的能量方面提供启示。具体地，模拟结果支持以下三个猜想：生态系统进化最大化系统总生物量；生态系统进化最大化物质总循环速率；生态系统进化最大化种群总适合度。这些结果基于现代进化观点，即每个种群在其相互作用种群进化所施加的约束内进化。每个被最大化的指标均为生态系统层面，故没有种群必然最大化其个体生物量、循环或适合度，而是可能进化至较低值，以使整个生态系统在整体上最大化该指标。Lotka的总生物量最大化建议，在质量守恒系统中等价于系统进化最小化未利用资源，这与R^*理论一致。分隔线的存在表明，单一生态结果可与多个进化结果关联，所有互惠组合（专性-专性、兼性-兼性、专性-兼性、兼性-专性）均可达到同一生态结果，因此仅靠指标可能无法可靠预测进化结果。动力系统分析工具与适合度景观结合可有效预测低维生态进化系统，而分岔行为则提示接近分岔的模拟需谨慎解释。此外，关于互惠共生进化，当提供互惠收益的成本较高时，成为专性互惠者可能是有利的；当成本降低时，更多互惠策略变得可行。