高分子量溶解性有机质通过增强微生物竞争和病毒碳代谢驱动土壤耐药基因组增殖

主要作者与研究机构

本研究由华东师范大学生态与环境科学学院、浙江舟山岛屿观测与研究站、河口海岸学国家重点实验室的张思宇（Si-Yu Zhang）研究员团队主导完成，第一作者为刘子腾（Zi-teng Liu）。合作者包括赵欣迪（Xin-di Zhao）、李嘉琦（Jia-qi Li）、李舒欣（Shu-xin Li）以及浙江大学环境与资源科学学院的唐先进（Xianjin Tang）。该研究成果发表于 *The ISME Journal*，于2026年8月18日正式刊出，论文编号为wrag212。

研究背景与科学问题

抗生素抗性基因（antibiotic resistance genes，ARGs）的环境扩散已被世界卫生组织列为全球公共卫生的重大威胁。土壤微生物群落是ARGs的重要储存库，而土壤有机碳作为微生物生长和群落结构的主要调控因子，被认为是影响ARGs丰度与多样性的关键因素。然而，以往研究多聚焦于总有机碳（total organic carbon，TOC）的含量效应，忽视了溶解性有机质（dissolved organic matter，DOM）这一最具生物活性的碳组分在分子水平上的异质性。

DOM的生物可利用性主要由其分子量（molecular weight，MW）决定。低分子量（low molecular weight，LMW）组分如糖类和氨基酸可被微生物快速利用，而高分子量（high molecular weight，HMW）组分如木质素和多糖则需要特定酶系和特定类群才能降解。这种碳源利用差异可能对微生物群落施加不同的选择压力，进而影响ARGs的增殖。此外，病毒通过编码辅助代谢基因（auxiliary metabolic genes，AMGs）优化宿主碳代谢的能力，以及病毒裂解和转导作用在ARGs传播中的潜在角色，使得DOM分子特征、微生物竞争、病毒活动与ARGs动态之间的耦合关系成为一个亟待厘清的科学问题。

基于此，本研究提出核心假设：不同分子量的DOM差异性地影响土壤ARGs的增殖，且病毒可能在调节宿主竞争力和碳源适应中发挥关键作用。研究旨在实现三个目标：（i）探究不同DOM分子量组与ARG富集的关联，（ii）阐明DOM分子量影响ARG增殖的机制，（iii）明确病毒在增强宿主竞争力和碳源适应中的作用。

研究设计与方法流程

本研究采用“大陆尺度野外调查+微观控制实验”的整合设计，分为三个主要阶段。

第一阶段：野外采样与分组。 研究团队于2022年8月至2023年7月，在中国18个省份的四大水稻主产区（华北、华东、华中、华南）采集了141份稻田表层土壤（0–20 cm）样本，涵盖47个采样点。所有采样田块仅施用化肥或符合国家NY/T 525-2021标准的有机肥，排除了污泥和粪肥直接施用的干扰。基于傅里叶变换离子回旋共振质谱（FT-ICR MS）分析获得的DOM分子组成，以400 Da为阈值将样本分为LMW组（n=87）和HMW组（n=54）。为验证分组稳健性，还以390、395、405、410和415 Da为替代阈值进行了敏感性分析，结果一致。

第二阶段：多组学数据获取与生物信息学分析。 对141份土壤样本进行了鸟枪法宏基因组测序（Illumina NovaSeq 6000，2×150 bp）。ARGs使用ARGs-OAP v3.2和ARG_ranker v3.0进行注释，碳代谢基因（carbon metabolism genes，CMGs）和碳水化合物活性酶（carbohydrate-active enzymes，CAZymes）基因分别基于KEGG和CAZy数据库通过DIAMOND比对注释。通过共组装（MEGAHIT）和分箱（metaWRAP和BASALT）策略，共回收132个高质量宏基因组组装基因组（metagenome-assembled genomes，MAGs，完整度≥90%且污染度≤5%），并利用AntiSMASH v7.1.0和BiG-SCAPE进行生物合成基因簇（biosynthetic gene clusters，BGCs）预测。病毒序列通过GenoMad、VirSorter2和CheckV三重鉴定流程从大于5000 bp的contigs中识别，经过严格质控后聚类为病毒操作分类单元（viral operational taxonomic units，vOTUs），利用iPHoP进行宿主预测。统计分析采用线性混合效应模型（linear mixed-effects models，LMMs），将采样点作为随机截距以处理样本非独立性，并使用偏最小二乘结构方程模型（PLs-SEM）解析多变量间的直接与间接效应。

第三阶段：微观实验验证。 选取湖南怀化（HN，LMW代表）和福建红塘（FJ，HMW代表）土壤，设置对照、木质素（lignin，MW=505.01，HMW DOM代表）添加和葡萄糖（glucose，MW=180.16，LMW DOM代表）添加三个处理，每个处理三个重复。在添加后第0、1、7和15天采集共44份样本，同时进行宏基因组和宏转录组测序，从DNA和RNA两个层面验证DOM添加对ARGs、CMGs、CAZymes及病毒AMGs的因果效应。

主要研究结果

DOM分子特征与ARG富集格局。 HMW土壤的DOM化学多样性（Chao1指数）、分子稳定性边界（MLB_r）和双键当量（DBE_w）均显著高于LMW土壤（β=0.59–1.59，p<0.05）。HMW土壤中ARGs的相对丰度（0.13 vs. 0.11 copies/cell）和丰富度（135 vs. 116）均显著高于LMW土壤（β=0.48–0.79，p<0.05），且ARG组成在两组间显著分化（PERMANOVA，F=6.9，p=0.001）。在18种核心ARG类型中，8种在HMW土壤中显著富集，包括novobiocin、MLS、tetracycline和β-lactam等。高风险ARGs（Rank I和Rank II）在HMW土壤中的相对丰度也显著更高（β=0.48，p<0.05），其中tetracycline抗性基因占高风险ARGs总量的20.5%。约2.4%的ARG宿主被鉴定为潜在人类病原菌，包括鲍曼不动杆菌（*Acinetobacter baumannii*）、星形诺卡菌（*Nocardia asteroides*）和表皮葡萄球菌（*Staphylococcus epidermidis*），且放线菌门（Actinomycetota）的ARG携带病原菌在HMW土壤中显著富集。

碳代谢与微生物竞争机制。 HMW土壤中CMGs和CAZyme基因的相对丰度显著高于LMW土壤（β=0.56–0.59，p<0.05），特别是降解复杂有机质的CAZyme基因（多糖GH家族、几丁质、木质素和纤维素降解酶）显著富集（β=0.55–0.74，p<0.05）。ARGs与CMGs/CAZyme基因的相对丰度呈显著正相关（r²=0.48–0.63，p<0.001）。在回收的132个MAGs中，68个（51.5%）被鉴定为共宿主（co-hosts），同时携带ARGs、CMGs和CAZyme基因。共宿主的总相对丰度在HMW土壤中显著高于LMW土壤（0.0068 vs. 0.0033 TAD80/GEQ）。共宿主含有的CMGs和CAZyme基因数量显著多于非宿主（147 vs. 115和91 vs. 66，p<0.001），且拥有更多复杂有机质降解酶。BGC分析显示共宿主携带更多核糖体合成和翻译后修饰肽（RiPPs）、非核糖体肽合成酶（NRPS）和其他聚酮合酶（PKS）的基因簇，能产生更多种类的抗菌天然产物。HMW土壤中微生物生态位宽度显著低于LMW土壤（β=−0.56，p=0.041），表明HMW DOM加剧了微生物群落间的资源竞争。

PLs-SEM揭示了驱动机制。 结构方程模型解释了ARG相对丰度60%的方差。HMW DOM相对丰度和CMGs/CAZyme基因对ARG丰度具有直接正效应（路径系数分别为0.10和0.49，p<0.05），而微生物β多样性具有直接负效应（路径系数=−0.27，p<0.001）。DOM化学多样性通过正向影响CMGs/CAZyme基因丰度间接促进ARG富集。HMW DOM丰度的总效应最大（0.30），其次是DOM化学多样性（0.24）。金属共选择的协变量调整模型证实，HMW DOM丰度与ARG丰度的关联在控制砷等金属浓度后仍保持显著（β=0.246，95% CI=0.009–0.483，p<0.05），而砷浓度的条件效应不显著。

微观实验的因果验证。 HMW土壤中添加木质素和葡萄糖均显著增加了ARGs在DNA和RNA水平的丰度（β_g=1.06–1.60，β_l=1.12–1.17，p<0.05），以及ARGs与持家基因的转录比率。CMGs和CAZyme基因的DNA和RNA丰度也同步显著升高。共宿主相对丰度在DOM添加后增加了约2.3倍（β_g=1.35，β_l=1.39，p<0.05），BGCs转录丰度亦显著增加。然而，在LMW土壤中，DOM添加仅在DNA水平增加了ARG丰度，转录水平无显著变化。在ARG-CAZyme共定位contigs（ACCCs）上，多药抗性基因*ranA*与糖基转移酶*gt2*、多粘菌素抗性基因*rosA*与*gt41*的共转录水平在HMW土壤木质素处理中最高。

病毒介导的碳代谢适应。 HMW土壤中病毒总丰度是LMW土壤的1.3倍（β=1.53，p<0.001）。HMW土壤中，DOM添加显著增加了AMGs在DNA和RNA水平的丰度（β=1.39–1.67，p<0.05），尤其是GH家族的复杂碳水化合物降解基因。关键的是，只有木质素添加特异性增加了多糖裂解酶（polysaccharide lyase，PL）家族AMGs的丰度。一个典型实例是vOTU3756（Caudoviricetes），其编码PL7基因，仅在木质素添加后在DNA（β=1.67，p<0.001）和RNA（β=1.27，p<0.05）水平显著富集，其预测宿主为Pseudomonas属的MAG 13，同样仅在木质素处理中显著富集（β=0.95，p<0.05）。这表明病毒AMGs在辅助宿主降解复杂有机质方面发挥了关键的“互助”作用。

研究结论与意义

本研究首次建立了土壤DOM分子量特征与ARG增殖之间的因果关系，证实HMW DOM是驱动土壤耐药基因组发展的关键但此前被忽视的因素。这一发现挑战了传统仅关注总碳含量的研究范式，将研究精度推进到分子水平。机制上，HMW DOM通过三条路径促进ARG富集：（1）低生物可利用性加剧微生物资源竞争，使同时携带ARGs和碳代谢基因的共宿主获得竞争优势；（2）共宿主通过BGCs产抗菌化合物抑制竞争者，进一步强化其生态优势；（3）病毒通过AMGs增强宿主的复杂碳降解能力，形成病毒-宿主协同适应。

该研究的应用价值在于提出了分子级有机质管理的概念。对于农业土壤，通过调控DOM的分子组成（例如减少木质素类难降解组分的积累、促进易利用碳源的周转），可能成为减缓ARGs环境扩散的新策略。此外，高风险ARGs和潜在病原菌在HMW土壤中的选择性富集，为农田生态系统的人体健康风险评估提供了新的预警指标。研究的方法学框架——整合大陆尺度调查、多组学分析和微观因果验证——也为微生物生态学中复杂环境驱动因子的解析提供了可借鉴的范式。