近日，来自浙江大学等单位科研人员，以水稻根际为模型生态系统，围绕三个关键科学问题展开：1）根际环境是否系统性地改变了与砷转化相关病毒的组成与生活史策略。2）这些病毒是否通过AMGs与HGT等机制，重编程宿主代谢网络，从而增强砷氧化能力。3）病毒在群落尺度上对砷氧化的真实功能贡献比例究竟有多大？研究者首次利用添加spike-in内标的宏基因组绝对定量方法，与病毒—宿主关联分析、基因组尺度代谢模型（GEMs）模拟以及体外诱导实验验证，构建了一个从“基因—代谢—群落功能”逐层递进的证据体系，文章发表在了国际知名杂志Nature Communications上。

中文题目：根际触发的病毒溶原性介导了微生物代谢重编程以增强砷氧化

发表时间：2025-4-30

研究者在稻田土壤中构建了砷污染微宇宙系统，针对种植水稻与未种植对照组开展了为期100天的高密度时间序列采样，频次为每48小时一次。该采样方案完整覆盖了水稻从营养生长期到生殖生长期的整个生育周期，系统性地获取了根际土壤与非根际土壤的配对样本，为后续多组学分析提供了高质量的时间序列数据基础。

在此基础上，研究者进行了以下检测分析：砷形态（As(III)/As(V)）的动态监测、宏基因组绝对定量检测以解析微生物与病毒群落结构、病毒生活史（裂解型/溶原型）的判定、病毒辅助代谢基因（AMGs）与水平基因转移（HGT）事件的鉴定、基于宏基因组组装基因组（MAGs）的基因组规模代谢模型（GEMs）定量模拟，以及溶原性病毒诱导的体外砷氧化功能验证实验。这一综合性设计确保了从相关性分析到机制解析再到因果验证，为揭示根际病毒-宿主互作驱动砷转化的生态机制提供了完整的实验结果。

1. 根际显著增强砷氧化，并与“携带病毒的砷氧化微生物”高度相关（图1）

研究首先确认，在整个水稻生长季中，根际土壤中无机砷几乎完全以As(III)与As(V)形式存在，占比超过97%。与非根际土壤相比，根际环境持续降低As(III)的相对比例，在拔节期差异为显著，高可降低约20%。

功能基因层面分析显示，砷氧化相关基因（aioA、aioB、arxA、arxS、arsH）在根际显著富集，而As(V)还原基因arrB在根际反而持续低于非根际。这表明根际并非简单增强“砷耐受”，而是方向性地推动砷氧化通路。进一步关键的是，当研究者将砷氧化微生物划分为“是否可识别到关联病毒”两类后发现：As(III) 降低幅度与“携带病毒的砷氧化微生物”丰度呈显著负相关，而与全部砷氧化微生物丰度并无显著关系。这一结果首次在群落尺度上暗示：病毒可能是驱动根际砷氧化效率的关键调控因子，而非被动旁观者。

图1. 根际与非根际土壤中砷生物转化、砷转化功能基因、砷转化微生物及其相关病毒绝对定量整体概览。（a）实验设计与生物信息学分析流程。（b）根际土壤相对于非根际土壤中 As(III) 比例的相对变化。（c）根际土壤中砷转化相关基因相对于非根际土壤的富集比例。富集比例按照如下公式计算：拷贝数富集比=根际绝对拷贝数/非根际绝对拷贝数。b、c 中虚线表示数值为 0 的参考线。所示曲线为Loess平滑拟合回归结果（n = 50）。阴影区域表示95%置信区间，反映拟合结果随时间变化的波动。（d）根际与非根际土壤中携带已识别病毒的砷氧化微生物（上）和砷还原微生物（下）的绝对定量丰度。圆点表示各时间点对应的微生物绝对定量丰度。虚线回归曲线为 Loess 平滑拟合（n = 50）。阴影区域表示 95% 置信区间，反映拟合结果随时间变化的波动。根际与非根际之间的统计显著性差异分别采用双侧 Student’s t 检验（上）和双侧 Welch’s t 检验（下）进行评估。（e）与砷氧化微生物和砷还原微生物相关病毒的分布（上）与分类（下）。右侧的双向网络图表示砷转化微生物属（前30）与病毒操作分类单元之间的对应关系，其中已注释分类的病毒操作分类单元以红色圆圈标出，并标注其对应的病毒科。（f）总体砷氧化微生物及其携带病毒的砷氧化微生物绝对定量丰度与根际相对于非根际 As(III) 比例变化之间的相关性。黄色与红色曲线表示线性回归模型，阴影带表示95%置信区间。回归系数的统计显著性采用双侧t检验进行评估（n = 50）。

2. 根际选择性促进砷氧化相关病毒的“溶原化”策略（图2）

在明确病毒可能参与砷氧化后，研究者进一步解析病毒生活史。结果显示，根际显著提高了与砷氧化微生物相关的溶原性病毒绝对丰度，同时降低了裂解型病毒比例。对应地，整合酶（integrase）基因在根际显著富集，病毒—宿主比（VHR）明显降低，均指向溶原作用增强。更重要的是，As(III) 降低幅度与“携带溶原性病毒的砷氧化微生物”绝对丰度之间呈现高度显著相关性，而与裂解型病毒无关。在时间尺度上，这种溶原性优势在水稻营养生长期达到峰值，在生殖期逐渐回落至非根际水平，显示出强烈的根系活动依赖性。

图2. 根际与非根际土壤中砷转化相关病毒的生活史特征。分别展示参与砷氧化（a）与砷还原（b）的溶原型与裂解型病毒丰度分布、病毒整合酶基因（c），以及砷氧化与砷还原相关病毒的病毒与宿主丰度比（d）。a–c（右）及 d（左）中，根际与非根际之间的统计显著性采用双侧 Welch’s t 检验评估（n = 50）。c（左）及 d（右）中，根际与非根际之间的统计显著性采用双侧 Student’s t检验评估（n = 50）。

3. 根际溶原性病毒系统性富集砷氧化与磷代谢AMGs绝对定量结果（图3）

在机制层面，研究者重点分析了病毒携带的AMGs的绝对丰度。结果发现，根际显著富集了两大类病毒AMGs：一类直接参与砷氧化，包括 aioB、arsH 和 arxS；另一类参与磷获取、磷饥饿响应及有机磷代谢，包括pstS、phnK/L/M、ugpC、RegX3、phoP等。时间序列显示，这些AMGs在水稻拔节期绝对丰度同步达到峰值，与溶原性病毒富集阶段高度一致。AlphaFold结构预测进一步证实，这些病毒AMGs具备完整、可信的蛋白结构，具备实际催化潜力。在宿主层面，携带上述AMGs的砷氧化微生物在根际显著富集，显示病毒赋予宿主明确的生态适应优势。

图3. 根际与非根际土壤中病毒辅助代谢基因及其对应宿主微生物特征。（a）根际中与砷氧化微生物相关病毒基因的富集分析。统计显著性采用Wald验评估。每个圆点代表一个基因，颜色表示其所属代谢通路。右侧图示显示不同代谢通路中富集基因的数量。（b）根际相对于非根际土壤中辅助代谢基因通路的富集比例。虚线回归曲线为Loess平滑拟合（n = 50）。阴影区域表示95%置信区间，反映拟合结果随时间变化的波动。灰色虚线表示数值为0的参考线。（c）根际中与砷氧化微生物相关的富集辅助代谢基因的蛋白结构。pTM得分用于衡量整体蛋白结构预测的准确性。pLDDT以0–100 的原子级置信度表示，数值越高表示预测置信度越高。（d）根际与非根际土壤中砷氧化辅助代谢基因的绝对丰度。（e）根际与非根际土壤中携带砷氧化辅助代谢基因的宿主微生物在属水平上的绝对丰度与组成。（f）根际与非根际土壤中磷代谢相关辅助代谢基因的绝对丰度。（g）根际与非根际土壤中携带磷代谢辅助代谢基因的宿主微生物在属水平上的绝对丰度与组成。

4. 病毒介导的砷氧化基因HGT在根际显著增强（图4）

除AMGs外，研究还系统鉴定了病毒介导的水平基因转移事件。研究者共识别出40起病毒向宿主转移aioB、arsH、arxS的明确证据，这些事件在系统发育树上表现为典型的HGT拓扑结构。HGT事件在时间上高度集中于水稻拔节期，且转移发生后，受体微生物在根际的绝对丰度显著上升。这一结果表明，根际不仅是病毒溶原化的热点，也是病毒驱动功能进化的“加速器”。

图4. 绝对定量结果解析根际中病毒介导的砷氧化基因水平转移特征。（a）病毒通过水平基因转移将砷氧化基因传递至宿主的系统发育来源。砷氧化基因的表型相关性以颜色标注，统计支持度采用标准自助法评估（n = 1000）。（b）根际中病毒介导水平基因转移的砷氧化基因相对于非根际土壤的富集比例。富集比例按照如下公式计算：读取深度富集比 = 根际读取深度 / 非根际读取深度。实线回归曲线为Loess平滑拟合（n = 50）。阴影区域表示95%置信区间，反映拟合随时间变化的波动。根际与非根际之间的统计显著性采用双侧 Welch’s t 检验评估（n = 44）。（c）根际中水平基因转移受体微生物相对于非根际土壤的拷贝数富集比例。实线回归曲线为Loess平滑拟合（n = 50）。阴影区域表示95%置信区间。根际与非根际之间的统计显著性采用双侧 Welch’s t 检验评估（n = 50）。（d）病毒介导的水平基因转移对受体微生物富集比例的影响。小提琴图展示了水平基因转移发生前后受体微生物的富集比例变化，实心点表示均值，竖线表示标准差。统计显著性采用双侧 Student’s t检验评估。在所示基因组图中，绿色基因表示宿主自身基因，其他颜色标记的基因表示通过病毒介导水平基因转移引入的砷氧化基因位置。

5. 基于GEMs的定量模拟：病毒贡献约21%的砷氧化通量（图5）

为突破“相关性”限制，研究者构建了693个物种水平的砷氧化微生物基因组尺度代谢模型，分别模拟其“有溶原病毒”和“无溶原病毒”两种状态。模拟结果显示：约 16.8% 的砷氧化反应完全依赖病毒基因组才能发生（virus-enabled）；约 4.2% 的砷氧化效率由病毒显著增强（virus-aided）。综合计算，溶原性病毒对根际微生物砷氧化的理论贡献约为 21%。进一步引入根际环境因子后发现，氧气是放大病毒作用的首要因素，其次为 pH、碳源和磷供应，明确揭示了“根系活动—环境梯度—病毒功能放大”的因果链条。

图5. 溶原性病毒对根际微生物砷氧化贡献的计算机模拟验证。（a）基于标志基因构建了693个物种水平砷氧化宏基因组组装基因组的系统发育树。树枝颜色表示门水平分类。内侧两圈颜色分别表示与砷氧化和磷代谢通路相关的辅助代谢基因。外侧下一圈表示宏基因组组装基因组在根际中的平均相对丰度。外圈表示病毒对砷氧化贡献的比例。（b）用于量化病毒对微生物砷氧化贡献的示意图。（c）根际环境变化对溶原状态下病毒促进砷氧化贡献的影响。d 在根际环境因子增强下，溶原性病毒促进砷氧化作用增强的微生物在门水平上的分类组成。

6. 体外诱导实验验证：病毒高贡献 25.63% 的砷氧化（图6）

在终验证环节，研究者通过丝裂霉素 C 诱导溶原性病毒释放，并在体外构建“仅细菌”与“细菌+病毒”体系。结果显示，来源于根际的溶原性病毒可使微生物砷氧化能力提高 25.63%，而非根际病毒则更倾向促进 As(V) 还原。这一实验结果与 GEMs 模拟高度一致，完成了从模型到现实的关键闭环。

图6. 根际溶原性病毒促进微生物砷氧化作用的实验验证。（a）土壤来源细菌中溶原性病毒诱导流程及体外培养实验设计。（b）丝裂霉素C诱导的溶原性病毒的透射电子显微镜图像。每个实验均独立重复三次，各重复之间结果一致。所示为具有代表性的实验图像。图中比例尺为100 nm。(c)根际与非根际溶原性病毒对细菌微生物砷氧化作用的影响。在每个时间点，病毒+细菌处理组与仅细菌处理组均设置三个生物学重复。图中点表示基于平均值计算的病毒+细菌处理中 As(V) 比例的相对变化。实线表示回归曲线（Logistic模型，n = 5）。根际与非根际之间的统计显著性采用双侧 Welch’s t检验评估。

本研究首次利用添加spike-in内标的宏基因组绝对定量检测方法，以及与病毒生态学、系统代谢建模与实验验证的深度整合，在群落尺度上定量证明：根际溶原性病毒可贡献约1/4的微生物砷氧化功能。研究系统揭示了一个多层级机制：根系活动塑造环境梯度 → 选择性促进病毒溶原化 → 病毒通过AMGs与HGT重编程宿主代谢 → 在群落尺度放大砷氧化功能。这一发现不仅颠覆了传统“病毒仅是裂解者”的认知，也将病毒提升为根际砷生物地球化学循环中的关键生态调控因子，为未来通过调控根际—病毒—微生物互作来降低作物砷风险提供了全新的理论基础。