中国近海海底泥里的病毒到底有多丰富，它们可能在海洋物质循环里做什么

中国近海海底泥里藏着一个非常庞大、绝大多数还没被认识过的病毒世界。这些病毒不只是杀细菌，它们还可能帮宿主改写碳、氮、硫、磷的代谢方式，甚至把人类活动带来的抗药基因也卷进了海洋沉积物生态系统。

大陆架虽然只占海洋表面积的一小部分，却承担了非常高比例的海洋初级生产、有机质沉降和生物地球化学循环。海底沉积物中的微生物本来就对碳、氮、硫、磷循环非常关键，而病毒作为海洋中数量最多的生命形式之一，不只是会杀死宿主、调节宿主数量，还会通过辅助代谢基因改变宿主代谢，甚至通过水平基因转移影响宿主进化。问题在于，过去关于海洋沉积物病毒的研究更多集中在深海热液、冷泉、海沟、海山这些特殊环境，对大陆架表层沉积物病毒了解得非常少。作者这篇文章想回答的，就是中国东部大陆架海域的沉积物中到底有多少病毒、这些病毒新不新、它们在空间上怎样分布、可能感染哪些宿主、是否携带与碳氮硫磷循环有关的基因，以及人类活动会不会把抗生素抗性基因也带进这些病毒基因组里。因为中国东部大陆架海域本身就是一个很典型、也很复杂的研究区域。这里包括渤海、黄海和东海，既有黑潮、台湾暖流、沿岸流这些复杂海流控制的温盐梯度，也有黄河、长江等河流不断带来的营养盐、有机质和陆源污染输入。换句话说，这个区域既能代表大陆架作为海洋物质循环热点的一面，也能代表近岸人类活动影响很强的一面。如果这里的沉积物病毒真的很丰富、很独特，还携带很多调控宿主代谢甚至抗性传播的基因，那它们就不能再被看成“附属角色”，而应该被看成大陆架生态系统中的重要参与者。

在方法上，这篇文章做的是一次大范围的沉积物宏基因组调查。作者一共采集了 48 个表层沉积物样品，覆盖渤海、黄海和东海大陆架，采样水深范围大约 12 到 130 米，只保留最上面的 0 到 1 厘米 沉积物做分析。所有样品做了宏基因组测序，总共得到大约 1.6 TB 的高质量 Illumina 数据。之后作者用多个病毒识别工具联合筛选病毒序列，包括 VirSorter2、DeepVirFinder、VIBRANT，再用 Cobra 做延伸，最后把序列按 95% 平均核苷酸一致性聚类成 vOTUs，也就是近似于“病毒种”层面的单元。与此同时，他们还从同一批数据里拼装出了 80 个 MAGs，也就是宿主微生物基因组草图，用来分析宿主组成和进行病毒宿主预测。宿主预测本身又结合了多种方法，包括序列同源性、tRNA 匹配、寡核苷酸频率和 CRISPR spacer 匹配。之后作者再用 DRAM-v 和 VIBRANT 去找病毒辅助代谢基因，用 CARD-RGI 去找抗生素抗性基因，并结合环境因子、海流分区、RDA、VPA、Mantel 检验等方法分析病毒群落的空间格局和驱动因素。

文章最先给出的一个核心结果是，这些大陆架沉积物病毒非常丰富，而且绝大多数都很新。作者从所有样品中一共识别出 12,540 个 vOTUs。其中只有 81 个完整基因组、94 个高质量基因组、240 个中等质量基因组，绝大多数仍然是低质量或不完整片段，这本身就说明沉积物病毒组装难度很高。但更重要的是，它们的“新颖性”非常惊人。按家族水平做分类时，93.6% 的 vOTUs 无法被分到已知病毒家族。作者又对长度超过 10 kb 的 2,908 条 vOTUs 做了基因共享网络分析，结果发现 60.4% 都是单独的 singleton，也就是和参考数据库里任何已知病毒都聚不到一起。就连能形成病毒簇的那部分里，也有不少是中国大陆架独有的独立簇。简单说，这片海底泥里藏着大量此前从未被认识过的病毒“暗物质”。作者还特别指出，真正能和 NCBI RefSeq 中已知病毒直接聚到一类的只有 14 个 vOTUs，这进一步说明目前公共数据库对大陆架沉积物病毒的覆盖仍然非常有限。

尽管大多数病毒都很新，作者还是尽量做了分类学归属。结果显示，在能被注释的部分中，最多的是 Uroviricota，其次是 Nucleocytoviricota 和 Artverviricota。其中 Uroviricota 里的大部分都属于尾噬菌体类群。值得注意的是，文章还发现大陆架沉积物里 NCLDV 相关病毒比例比很多其他沉积环境更高，这说明这里不只是细菌病毒多，感染真核生物的大型 DNA 病毒也不容忽视。能分到的家族中，较常见的包括 Kyanoviridae、Autographiviridae、Herelleviridae、Zobellviridae，而在 Nucleocytoviricota 里，较多的是 Phycodnaviridae、Mesomimiviridae、Mimiviridae。不过作者始终强调，这些已知家族只占很小一部分，真正的主体仍然是未知病毒库。

在空间分布上，作者发现，季节差异对病毒群落的影响没有空间差异大，也就是说，同一海区不同季节之间的差异，不如不同海区之间大。为了更好理解这种空间格局，作者不是简单按行政海区分组，而是按海流和近岸、外海特征把样品重新分成几个区域，包括 BSNYS、SYS、ECSA、ECSB 和 ECSC。这样分组之后，病毒群落结构差异就更清楚了。总体趋势是，近岸组的病毒多样性更高，离岸越远、多样性越低。例如 ECSA、ECSB、ECSC 之间就表现出随着离岸距离增大和暖流影响增强，Shannon 多样性、Simpson 多样性和 Chao1 丰富度逐渐下降的趋势。作者认为，这说明大陆架沉积物病毒并不是随机分布的，而是受海流、近岸输入和环境梯度共同组织起来的。进一步分析驱动因素后，作者发现，温度和水深 是最关键的环境变量，同时宿主分布也非常重要。变异分解分析显示，单独看宿主因素可以解释 29.2% 的病毒群落变异，单独看环境因素能解释 12.9%，两者还有 19.2% 的共享解释量。这说明环境并不是直接单独决定病毒，而是往往先改变宿主分布，再通过宿主去影响病毒。RDA 也显示环境和宿主一起能解释约 46.7% 的病毒群落变异。Mantel 分析进一步表明，不同区域受控因素还不一样。比如渤海和北黄海更受温度、氮浓度和 pH 影响，南黄海更受纬度和浊度影响，东海近岸更受盐度和叶绿素 a 影响，而东海外海最远的那组则更受宿主 MAG 多样性影响。作者因此得出一个非常重要的结论：大陆架沉积物病毒的空间格局，本质上是病毒、宿主和环境三者耦合的结果，而海流通过塑造温盐梯度和宿主更替，间接塑造了病毒群落。

这篇文章并没有停在病毒有多少、分布在哪，而是进一步去看这些病毒可能感染谁。作者一共为 568 个 vOTUs 预测到了潜在宿主，大约占全部 vOTUs 的 4.5%。这个比例不算高，但在环境病毒组里已经是常见水平。被预测的宿主覆盖了 16 个细菌门和 2 个古菌门，其中最主要的是 Proteobacteria，占所有宿主预测配对的 65.2%。后面比较多的是 Myxococcota 和 Desulfobacterota。作者特别强调，Proteobacteria 在这些沉积物里本来就很重要，因为它们代谢类型非常多样，和碳、氮、硫循环密切相关。Desulfobacterota 则是典型的硫酸盐还原相关类群，在缺氧沉积物中非常关键。这意味着这些病毒并不是随机感染一些无关紧要的宿主，而很可能是瞄准了沉积物里那些最有代谢影响力的微生物。更进一步，作者分析了这些宿主 MAGs 的代谢潜力，发现它们常常带有复杂多糖降解、硝酸盐还原、硫代谢等多种代谢途径。比如一些预测宿主带有纤维素酶、几丁质酶、pullulanase 等复杂碳源利用酶，也带有 napAB、nrfADH、nirBD、nirKS 这类硝酸盐和亚硝酸盐转化基因，以及 dsrABD、asrABC、sdo、sor 等硫代谢相关基因。作者据此认为，这些病毒如果感染这些宿主，就有可能通过改变宿主代谢间接影响沉积物中的碳、氮和硫循环。换句话说，这里的病毒生态意义不只是“裂解宿主”，而是可能通过宿主去调节整个沉积物代谢网络。

文章还分析了病毒生活方式。结果显示，这些大陆架沉积物中的病毒以裂解型为主，在各区域都占超过 75%。温和型病毒比例从近岸向外海逐渐降低，而“未定型”生活方式的病毒比例则在外海增加。作者认为，这说明沿着大陆架从近岸到外海，病毒感染策略也在变化，可能与环境梯度和宿主群落组成共同有关。不过他们也很谨慎地指出，这种空间变化背后的具体机制，目前还不能完全解释清楚。在生态功能上，是病毒携带的 AMGs。作者一共识别出 3,321 个候选 AMGs，其中有一部分来自前噬菌体。这个数量非常大，说明大陆架沉积物病毒在代谢干预方面的潜力非常强。KEGG 注释显示，这些 AMGs 最多和 碳代谢、辅因子和维生素代谢、氨基酸代谢 有关。作者特别关注了与 碳、氮、磷、硫循环 直接相关的基因，因为这些最能反映大陆架沉积物生态过程。

在碳代谢方面，作者发现病毒 AMGs 不只是包含常见的糖酵解、糖异生和 TCA 循环相关基因，还带有一些能帮助宿主利用“非常规碳源”的基因，比如 gatZ/kbaZ 这类与 D-塔格糖利用相关的基因，以及 prpB 这种参与丙酸代谢的基因。更关键的是，他们还发现了 16 类与淀粉和蔗糖代谢相关的 AMGs，包括编码 alpha-trehalase、cellobiosidase、alpha-amylase 等复杂多糖水解酶的基因。作者的解释很直接，这意味着病毒可以通过把这些基因带给宿主，扩展宿主可利用的碳源范围，从而改变沉积物中的碳降解路径。在硫循环方面，结果尤其突出。作者发现了一套非常完整、很像“从硫酸盐一路走到硫化物”的病毒编码酶，包括 sat、cysC、cysH、sir 等。这组基因基本覆盖了硫酸盐还原通路中的关键步骤。作者认为，这在缺氧沉积物里特别重要，因为硫循环本来就是这类环境中的核心过程之一。更有意思的是，其中的 cysC 甚至来自前噬菌体，说明不只是游离裂解型病毒，整合到宿主里的病毒也可能参与调控关键硫代谢步骤。在氮循环方面，作者找到了 12 类、29 个 与氮代谢相关的 AMGs，包括 nrfH、nirA、nirK、norB 等。这些基因涉及硝酸盐还原、亚硝酸盐还原和反硝化的不同步骤。作者特别指出，硫还原和氨氧化在缺氧沉积物中往往是耦联发生的，因此病毒携带这些氮、硫代谢基因，很可能会共同影响沉积物中的关键氮硫耦合过程。在磷循环方面，作者发现了 phoD 和 phoN 这类酸性或碱性磷酸酶相关基因，能够帮助宿主把有机磷分解释放为无机磷。值得注意的是，phoN 出现在更受河流输送影响的区域，而 phoD 主要出现在渤海和黄海偏碱性的环境中。作者认为，这说明病毒 AMGs 不只是“有”，而且在不同环境条件下表现出很明确的地理分布特征，反映了它们可能与环境磷获取策略相匹配。

除了 AMGs，这篇文章另一个很受关注的结果是 ARGs。作者在这些大陆架沉积物病毒中一共检测到了 33 个抗生素抗性基因，可归为 8 大类抗生素 和 28 个亚型。最丰富的几类分别是 β-内酰胺类、糖肽类、大环内酯类、喹诺酮类 和 四环素类。更特别的是，作者还发现了一个与季铵盐类消毒剂苯扎氯铵耐受相关的 qacG，并指出这是首次在病毒基因组中报道到这个基因。空间格局上，这些 ARGs 呈现出很清楚的近岸高、外海低的趋势，尤其是 β-内酰胺抗性基因 主要集中在更靠近陆源输入的 BSNYS、SYS 和 ECSA，而离岸更远的 ECSB 和 ECSC 基本没有这类基因。作者把这一现象和中国近年医疗体系中抗生素使用结构做了比较，发现沉积物病毒中的 ARG 类型与人类使用抗生素的主要类别高度一致，因此认为这些 ARGs 很可能和陆源污染、尤其是污水和河流输入密切相关。也就是说，病毒在这里不仅是自然生态参与者，也可能是人类活动影响海洋沉积物的遗传记录者和潜在传播者。

这篇文章的主要结果可以归纳成这样一条非常清楚的逻辑链：中国东部大陆架表层沉积物中存在一个极其丰富、极其新颖的病毒库。它们的空间分布不是随机的，而是受 温度、水深、海流以及宿主分布 共同控制。它们主要感染那些与碳氮硫磷循环关系密切的宿主，且大多以裂解型生活方式存在。与此同时，它们还携带大量与碳、氮、磷、硫代谢有关的辅助代谢基因，说明它们有潜力通过改写宿主代谢来参与沉积物中的关键生物地球化学过程。最后，这些病毒还携带多种抗生素抗性基因，而且在近岸更强，提示人类活动已经明显影响到了大陆架沉积物病毒库的遗传组成。

结论

大陆架沉积物病毒是大陆架生态系统中不可忽视的组成部分，它们的作用远远不只是感染和裂解宿主。它们还可能通过 AMGs 调节宿主代谢，通过 ARGs 参与抗性基因传播，并通过与环境和宿主的耦合关系塑造整个大陆架沉积物生态过程。作者同时也很诚实地指出，目前仍有几个明显局限。比如大部分病毒无法分类，宿主预测仍然主要是计算推断，恢复到的 MAG 数量相对整个微生物多样性还偏少，而且关于 ARG 来源和海流输送的推断仍然偏间接。因此，这篇文章更像是把大陆架沉积物病毒这个巨大黑箱打开了一个很重要的入口，而不是把所有问题都已经解决完了。今后研究大陆架和近海沉积物时，不能只关注细菌和古菌，因为病毒不仅数量巨大，而且在代谢调控和遗传传播上可能承担了被长期低估的角色。未来如果想真正理解大陆架在全球碳中和、养分循环和近岸污染扩散中的作用，就需要把病毒也纳入核心框架。更进一步说，这篇文章也提示，今后的重点不仅是发现更多新病毒，还要结合长读长测序、转录组、实验宿主验证和海流粒子输运模型，去真正弄清这些病毒什么时候活跃、感染谁、如何调控宿主，以及它们在多大程度上参与了大陆架生态系统对人类活动的响应。