[论文解读] Mathematical Modeling on Obligate Mutualism: Interactions between leaf-cutter ants and their fungus garden

Yun Kang, Michael Makiyama+2 more|arXiv (Cornell University)|Feb 11, 2011

Insect and Arachnid Ecology and Behavior被引用 4

一句话总结

本研究基于米氏- Menten动力学构建了一种新颖的数学模型，用于分析早期蚁巢发育阶段切叶蚁与其菌圃之间专性共生关系的机制。研究发现，分工行为与初始种群规模对共存至关重要，敏感性分析表明，真菌生长速率和蚂蚁种群参数的影响大于转化率或死亡率。

ABSTRACT

We propose a simple mathematical model by applying Michaelis-Menton equations of enzyme kinetics to study the mutualistic interaction between the leaf cutter ant and its fungus garden at the early stage of colony expansion. We derive the sufficient conditions on the extinction and coexistence of these two species. In addition, we give a region of initial condition that leads to the extinction of two species when the model has an interior attractor. Our global analysis indicates that the division of labor by workers ants and initial conditions are two important {factors} that determine whether leaf cutter ants colonies and their fungus garden survive and grow can exist or not. We validate the model by doing the comparing between model simulations and data on fungal and ant colony growth rates under laboratory conditions. We perform sensitive analysis and parameter estimation of the model based on the experimental data to gain more biological insights on the ecological interactions between leaf cutter ants and their fungus garden. Finally, we give conclusions and {discuss} potential future {work}.

研究动机与目标

• 理解切叶蚁与其菌圃在早期适应性生长阶段的种群动态。
• 解决专性共生关系缺乏理论模型的问题，特别是涉及复杂行为特征（如分工）的共生关系。
• 利用实验室蚁巢的实证数据验证模型，并估算难以直接测量的参数。
• 通过全局动力系统分析，探究互惠系统的稳定性和共存条件。
• 探讨生命周期特异性行为与初始条件对蚁巢生存与生长的影响。

提出的方法

• 将米氏- Menten动力学适配用于建模蚂蚁-真菌相互作用中的功能响应与数值响应。
• 基于Holland与DeAngelis的消费者-资源框架，构建一个包含密度依赖性生长与捕食项的双物种系统（蚂蚁与真菌）。
• 通过将工蚁在真菌栽培与叶片采集中的角色建模为不同的功能响应，引入分工机制。
• 进行全局稳定性分析，推导出两种物种灭绝与共存的充分条件。
• 利用参数估计与敏感性分析，结合实验数据对模型进行校准，以确定种群增长率。
• 将模型预测与实验室数据进行对比，特别关注第6至第9周的偏差，以揭示模型的局限性。

实验结果

研究问题

• RQ1 分工行为与初始种群规模如何影响切叶蚁与其菌圃的共存或灭绝？
• RQ2 该专性共生关系中，稳定共存与相互灭绝的充分条件是什么？
• RQ3 真菌生长速率、蚂蚁种群增长速率与死亡率等关键参数如何影响蚁巢发育结果？
• RQ4 为何模型在第6至第9周与数据拟合较差？这反映了该阶段的哪些生物过程？
• RQ5 哪些参数对模型输出影响最大？这些参数如何从实证数据中估计？

主要发现

• 模型识别出蚂蚁与真菌共存所需的初始种群规模阈值，若未达到该阈值则导致相互灭绝。
• 敏感性分析显示，半饱和常数 $ b $ 对模型输出的影响最大，其次为真菌生长速率 $ r_f $ 与蚂蚁生长速率 $ r_a $。
• 参数估计结果表明，转化率 $ r_c $ 与蚂蚁死亡率 $ d_a $ 极小且对初始值高度敏感，暗示其对早期动态影响微弱。
• 模型在早期蚁巢生长阶段与实验数据拟合良好，但在第6至第9周出现不一致，表明模型可能难以准确捕捉后期发育阶段的动态。
• 模型中存在多个吸引子，表明早期蚁巢发育阶段的种群动态可能对初始条件极为敏感，表现出高度不稳定性。
• 本研究结论认为，需构建更详细的模型，例如引入发育阶段（卵、幼虫、蛹）或随机性，以更好地捕捉早期蚁巢动态的复杂性。