温度驱动的海洋微藻病毒生物地理分布：揭示巨型病毒与浮游植物的生态互作机制

《ISME Communications》：Temperature-driven biogeography of marine giant viruses infecting picoeukaryotes Micromonas

本研究针对温度如何塑造海洋巨型病毒（Micromonas virus, MicV）的生物地理分布这一关键科学问题，通过整合TARA海洋宏基因组数据和温度驱动的动力学模型，首次揭示了温度通过调控病毒裂解概率（pη）和感染概率（p?）等生命史特征（LHTs），驱动海洋微藻病毒形成嗜冷（cryophile）与冷中温（cryo-mesophile）生态型的分布格局。该研究为预测气候变化下病毒-宿主互作演变提供了理论框架，发表于《ISME Communications》。

　　

海洋中存在着数量惊人的病毒颗粒，每毫升海水含有超过1000万个病毒粒子，它们如同微观世界的"隐形调控者"，通过感染微生物深刻影响着全球碳循环。其中，感染微型真核藻类Micromonas的巨型病毒（Prasinovirus）尤为特殊——这类病毒拥有复杂基因组（可达600 kb），能显著调控宿主种群动态。然而，温度如何塑造这些病毒的地理分布？在气候变暖背景下病毒-宿主关系将如何演变？这些问题长期困扰着科学家。

法国索邦大学（Sorbonne Université）海洋生物多样性与生物技术实验室（LBBM）的David Demory团队联合日本京都大学化学研究所生物信息中心（Bioinformatics Center, Kyoto University）等机构，通过分析TARA海洋考察计划采集的370个宏基因组样本，结合温度驱动的动力学建模，揭示了温度是MicV生物地理分布的首要驱动因子。研究发现：

关键技术方法

1. 1.

   宏基因组分析：基于TARA Oceans数据库的PolB标记基因和18S rRNA基因，通过EPA-ng算法对MicV和宿主Micromonas进行系统发育分类

2. 2.

   温度响应模型：建立包含病毒吸附率（?）、裂解率（η）等7个温度依赖参数的动力学方程，量化p_η和p_?的温敏性

3. 3.

   网络分析：采用FlashWeave算法构建包含103对病毒-宿主互作关系的共现网络

4. 4.

   ROC分析：通过接收者操作特征曲线确定病毒存在的温度阈值

主要研究结果

温度是MicV分布的关键环境因子

通过主坐标分析（PCoA）发现温度解释了25%的MicV群落变异（图1c），显著高于盐度（12%）和叶绿素a（8%）。在北极海域，MicV占巨型病毒群落（Nucleocytoviricota）的50%以上，而在热带海域占比不足5%。

病毒形成温度适应型谱系

系统发育分析（图2a）显示：

• •

  嗜冷型：B组（4.8±4.8°C）和Pol组（7.7±8.7°C）

• •

  冷中温型：A2组（17.1±6.3°C）主导0-26°C水域

• •

  中温型：A1组仅在>10°C水域出现，C组（13.4±0.5°C）分布狭窄

病毒-宿主温度适应协同进化

共现网络分析（图3a）揭示病毒与宿主的温度最优值显著正相关（R=0.71，p<0.01）。例如感染Micromonas polaris的Pol组病毒与其宿主共同适应7°C低温环境。

裂解概率决定分布边界

动力学模型显示（图4）：

• •

  p_η和p_?能准确预测A/B/Pol组的分布（AUC>0.85）

• •

  A组实际分布温度上限（23.1°C）接近理论极限（23.5°C），而B组存在12°C的生态位压缩，表明种间竞争影响

研究意义与展望

该研究首次量化了温度通过病毒生命史特征（LHTs）调控地理分布的机制，为预测气候变暖下病毒介导的碳循环变化提供理论基础。发现的"生态位压缩"现象（如B组病毒）提示竞争压力可能加剧极地病毒的灭绝风险。未来需结合单病毒测序（single-virus genomics）解析温度适应的分子机制，并扩展至其他浮游植物-病毒系统。

（注：所有数据均来自原文，图1-5的详细描述已嵌入关键结论处，专业术语如PolB=家族B DNA聚合酶，EPA-ng=进化放置算法等均按原文标注）