《Environmental Microbiology》：Coexistence of Photosynthetic Marine Microorganisms, Viruses and Grazers: Towards Integration in Ocean Ecosystem Models
编辑推荐：
本研究聚焦于海洋光合微生物、浮游动物捕食者与病毒之间的共存问题，旨在为大型海洋生物地球化学模型整合病毒动态提供理论基础。研究人员开发了一个包含元素配额、适用于海洋生态系统模型的微生物-捕食者-病毒群落模型，探讨了感染类群、对病毒感染易感性的异质性以及捕食者的高阶死亡率对三者共存机制的促进作用。研究发现，病毒潜伏期与毒力之间存在权衡，潜伏期产生的振荡降低了自由病毒的增长速率，促进了共存；易感性异质性通过资源分配支持共存；而高阶死亡率则拓宽了共存范围。模型结果表明，在简化的上层海洋环境中，包含中等到强病毒感染抗性的参数组合能够产生符合生态实际的属性。这项工作强调了在海洋生态系统模型中纳入病毒感染抗性的重要性，为理解和预测病毒在海洋食物网及生物地球化学循环中的作用迈出了关键一步。
在浩瀚的海洋中，生活着一群微小的“太阳能工厂”——光合浮游植物。它们是海洋食物网的基石，通过光合作用产生有机物，支撑着从微小的浮游动物到庞大的鲸鱼等一系列生命，更在全球碳循环和气候调节中扮演着不可或缺的角色。然而，这些“海洋牧草”的繁荣并非没有挑战。它们面临着来自两个方向的致命威胁：一边是浮游动物（如桡足类、鞭毛虫）的捕食（即“ grazing ”），另一边则是数量更为庞大的病毒（病毒体浓度通常超过 107个/毫升）的侵染与裂解。有趣的是，这两种“天敌”虽然都以浮游植物为资源，却在海洋中普遍共存。这引发了一个核心的生态学问题：在竞争同一资源的情况下，病毒和浮游动物是如何实现长期共存的？
理解这种共存机制，不仅关乎基础生态理论，更具有紧迫的实际意义。当前，用于预测全球气候变化和海洋生态变化的大尺度海洋生物地球化学模型日益复杂，包含了多种浮游植物、固氮生物、混合营养生物以及多类捕食者。然而，这些模型普遍缺乏对病毒侵染浮游植物过程的明确表征。这被视为一个重要的未来挑战，因为病毒裂解被认为通过“病毒分流”（viral shunt，将有机物转化为溶解态，促进营养盐再生）和“病毒穿梭”（viral shuttle，促进颗粒有机物聚集与沉降，增强碳输出）等机制，深刻影响着海洋物质流动、营养循环乃至“生物泵”的效率。在将病毒动力学整合进大尺度模型之前，首先需要在更简单的理论框架下，厘清病毒、捕食者与宿主之间的动态关系与共存条件。
为此，由 Jérémy P. P. Laugier 等人领导的研究团队在《Environmental Microbiology》上发表了一项研究。他们从经典的易感-感染-病毒（Susceptible-Infected-Virus, SIV）计数模型出发，构建了一个包含元素（氮）配额、适用于生物地球化学模型的浓度框架。在此基础上，研究人员逐步增加生态复杂性，开发了四个不同复杂度的模型：基础竞争模型（SZV）、包含感染类（I）的模型（SIZV）、包含抗性类型（R）的模型（SZRV）以及同时包含感染类和抗性类型的模型（SIRZV）。他们利用来自 Prochlorococcus（原绿球藻）、Synechococcus（聚球藻）、微小真核生物和小型硅藻及其对应病毒和捕食者的生命史特征数据（如裂解量、潜伏期、吸附率）进行参数化，重点评估了（1）包含明确感染类、（2）宿主对病毒感染易感性的异质性（即抗性类型）、以及（3）捕食者（病毒和浮游动物）的高阶（二次）死亡率，这三种机制在促进病毒与浮游动物共存方面的作用。研究还利用模型代数平衡解，寻找了在简化的寡营养和富营养上层海洋环境中，能产生符合实际观测的生态属性（如浓度、感染百分比）的参数组合。
本研究主要采用了理论建模与数值模拟相结合的方法。首先，基于经典的SIV模型框架，将其从个体计数单位转化为元素浓度单位，构建了包含病毒“分流”和“穿梭”通量的基础模型。其次，通过系统性地增加模型复杂度（SZV, SIZV, SZRV, SIRZV），并利用已有的海洋微生物和病毒生命史特征数据集进行参数化。关键方法包括：1）模型推导与扩展：在元素配额框架下，推导并比较了四个核心模型的微分方程。2）共存理论分析：运用现代共存理论，通过入侵分析和生长率分解，定量评估不同机制对共存的贡献。3）参数优化与验证：针对 SIRZV 模型，在参数空间中进行网格搜索，以最小化模型稳态输出与基于观测数据（如病毒和浮游植物浓度）设定的目标浓度之间的误差，从而找到生态学上合理的参数组合。4）数值模拟：对模型进行长达20年的数值积分，通过傅里叶分析和稳态判别，确定系统最终的动力学状态（共存、排除或振荡）以及各组分浓度、感染百分比等关键生态指标。
结果部分的主要发现如下：
1. 包含感染类或抗性类型有助于拓宽共存区间
研究人员以 Prochlorococcus为例进行了详细分析。在基础的 SZV 模型中，病毒和浮游动物遵循竞争排斥原则，无法稳定共存，总有一方被排除。然而，在包含了明确感染类（I）的 SIZV 模型中，病毒潜伏期的引入带来了动力学振荡。这种内生的资源波动，通过生长率的非线性响应，为两种捕食者的共存创造了条件。同样，在包含抗性类型（R）的 SZRV 模型中，易感性（对病毒吸附和/或细胞内感染的抵抗能力）的异质性将浮游植物种群划分为不同的资源类型（易感S和抗性R），相当于为两种捕食者提供了差异化的资源基础，从而通过资源分配促进了共存。SIRZV 模型兼具了以上两种机制。
2. 病毒潜伏期与毒力存在权衡，并影响共存
在 SIZV 模型中，研究揭示了病毒潜伏期与其“毒力”（这里表现为感染成功率和清除宿主的能力）之间存在着关键的权衡关系。较长的潜伏期会导致感染细胞在裂解释放新病毒子之前，有更多机会被浮游动物捕食（种内捕食）。这虽然削弱了病毒的直接传播效率，但同时也抑制了自由病毒数量的过快增长，从而稳定了系统，促进了病毒与浮游动物的共存。模型分析表明，这种由潜伏期引发的种群振荡，通过降低自由病毒的平均增长率，为浮游动物的入侵和持续存在提供了机会。
3. 高阶（二次）死亡率拓宽共存参数空间
在模型中引入病毒和浮游动物的二次死亡率项（如病毒的非特异性吸附到颗粒物、浮游动物的更高阶捕食），能有效拓宽两者共存的参数空间。这主要是因为二次死亡率对数量占优的捕食者施加了更强的密度制约，防止其过度增长并排除竞争对手，从而起到了稳定群落的作用。
4. 包含病毒抗性的模型能更好地匹配观测数据
通过对 SIRZV 模型进行参数优化，研究人员发现，能够同时匹配寡营养和富营养环境中观测到的浮游植物、病毒浓度，以及合理的感染百分比、病毒致死率比例的参数组合，普遍要求宿主对病毒感染具有中等到强的抗性（包括细胞外抗性，即降低病毒吸附率；和/或细胞内抗性，即降低感染成功概率）。这一发现强烈提示，在未来的海洋生态模型中纳入宿主的病毒抗性机制，对于准确模拟病毒-宿主-捕食者系统的动态及其生物地球化学效应至关重要。
结论与讨论 本研究通过发展一套包含元素配额、可整合进海洋生物地球化学模型的病毒动力学理论框架，系统地探讨了病毒与浮游动物在共享浮游植物资源下的共存机制。研究证实，简单的竞争模型无法实现两者的稳定共存。而引入生物学上合理的复杂性，特别是包含明确的感染类群、宿主对病毒易感性的异质性（抗性）以及捕食者的高阶死亡率，是拓宽和稳定共存区间的关键。其中，病毒的潜伏期通过引发种群振荡，调节了病毒的增长速率；宿主的抗性变异则通过资源分配，为竞争者们提供了差异化的生态位。
这些发现具有重要的理论和应用意义。理论上，它深化了我们对海洋微生物群落中多营养级相互作用，特别是寄生者（病毒）与捕食者（浮游动物）之间复杂竞争关系的理解。应用上，该研究为将病毒过程整合进大尺度海洋生物地球化学和气候模型提供了具体的建模指南和参数化思路。模型优化结果表明，要真实地反映海洋中的病毒丰度、感染水平及其对初级生产力和碳循环的影响，在模型中纳入中等到强水平的宿主病毒抗性机制是必要的。这为未来利用模型量化病毒“分流”和“穿梭”在全球海洋碳循环中的相对重要性，以及预测环境变化下病毒-宿主-捕食者互作网络的响应，奠定了坚实的基础。最终，这项研究朝着构建一个更完整、更能反映海洋生态系统复杂性的数值模型迈出了关键一步，有助于我们更准确地预测海洋在气候变化中的角色和未来。
