水文情势主导粪源抗生素抗性基因在稻–蚕豆轮作中的时空归趋及碳–微生物–MGEs级联调控

《Antibiotics》：Hydrological Regime Governs the Fate of Manure-Derived Antibiotic Resistance Genes in Rice–Faba Bean Rotation

　　

背景/目标：农业土壤中抗生素抗性基因（ARGs）的增殖对“同一健康”（One Health）构成日益严重的威胁。然而，稻–旱轮作中的水文波动如何塑造土壤抗性组（resistome）仍不清楚。本研究考察了施用粪肥的水稻–蚕豆轮作系统中土壤ARGs的时空动态及其驱动因素。方法：开展为期一年的田间试验，采用定量PCR（qPCR）和16S rRNA基因测序表征ARGs、可移动遗传元件（MGEs）和细菌群落。结果：水文情势与粪肥施用交互驱动抗性组演替。粪肥引入多种ARGs和MGEs，而氧化还原波动重构了细菌宿主群落及ARG–MGE网络。与湿相相比，干相使intI1富集约两个数量级。具体而言，log₁₀转换的intI1/16S rRNA基因拷贝数比值从?7.11升至?5.14，表明水平基因转移（HGT）潜力更大。休耕过渡期是抗性组重构的关键阶段，具有持久遗留效应。结构方程模型（SEM）识别出以碳–微生物–MGE为主的级联。活性有机碳组分，包括溶解性有机碳（DOC）和易氧化碳（ROC），塑造微生物群落结构，进而调控MGE丰度并最终影响ARG动态。MGEs是ARGs最强的直接驱动因子（路径系数=0.575，p < 0.001）。结论：干湿交替通过碳–微生物–MGE级联调控土壤抗性风险。因此，休耕过渡期为通过管理水文和碳输入以减缓ARG传播、实现可持续农业生态系统提供了关键靶点。

研究背景与问题：抗生素抗性基因（ARGs，赋予微生物抗抗生素能力的基因）已成为持久性环境污染物，在“同一健康”（One Health）框架下对全球公共卫生构成重大威胁。畜禽养殖和作物生产中抗生素的集约使用使农业土壤成为ARGs的重要储库；ARGs不仅可经垂直遗传在微生物群落中持续存在，更可借助可移动遗传元件（MGEs，介导基因水平转移的遗传因子）介导的水平基因转移（HGT，遗传物质跨物种屏障转移）跨物种传播。粪肥施用和污水灌溉是外源ARGs进入表土的主要途径，且农业土壤是ARGs在环境、共生及潜在致病细菌间循环的重要界面。现有土壤ARG污染缓解策略主要依赖物理化学修复和农艺生态管理，但大规模物理修复常经济不可行，通过优化种植模式调控土壤微环境更具可持续性。轮作作为经典保护性耕作措施可提升土壤健康，但调控土壤抗生素抗性组（resistome，土壤中全部抗性基因集合）的具体机制仍不清楚。既有研究多集中于连作障碍或有机肥外源输入，且多来自旱地轮作系统；稻–旱轮作具有独特的氧化还原电位波动、淹水–非淹水循环及好氧–厌氧微生物群落演替等生态位特征，可能通过调控MGE活性和HGT潜力根本性重组ARG动态，但相关机制尚未阐明。水稻–蚕豆轮作是中国典型稻–旱种植系统，整合淹水稻作与旱地豆科栽培，形成独特的水–旱和碳–氮耦合环境，然而不同生育季土壤ARGs的时空动态及其环境因子—抗性组组装的机制通路仍缺乏研究。为此，研究人员开展田间小区试验，以表征干湿水文循环下土壤ARGs的时空变异、识别关键环境驱动因子（尤其活性有机碳组分与季节性水文情势的交互效应），并阐明介导粪肥土壤中ARG持久性和水平传播的微生物共现网络与宿主演替机制。研究结论表明，水文情势主导粪肥施用稻–蚕豆轮作系统中土壤抗性组的时间演替；粪肥引入多样ARGs和MGEs，循环淹水–非淹水转变重塑ARG丰度和组成，休耕过渡期是抗性组重构关键窗口；活性有机碳组分（尤其DOC和ROC）是关键生态驱动因子，通过微生物群落结构和MGE丰度形成碳–微生物–MGE调控级联。其重要意义在于突出水文管理和关键过渡期是减轻稻–旱农业生态系统中粪肥土壤ARG风险的重要靶点。

主要关键技术方法：在云南大理洱海流域云南大理农业生态系统国家野外科学观测研究站开展一年田间试验；设置无粪肥无种植（CU）、无粪肥轮作（CP）、粪肥无种植（MU）、粪肥+轮作（MP）4个处理，3次重复，共12个小区，并于P1–P9共9个时间点重复采样。采用qPCR绝对定量22类ARGs及intI1、intI2、traA三个MGE标记，采用16S rRNA基因V3–V4区扩增子测序分析细菌群落；统计分析采用重复测量线性混合效应模型、主坐标分析（PCoA）、置换多元方差分析（PERMANOVA）、冗余分析（RDA）、Mantel检验、共现网络分析及偏最小二乘结构方程模型（PLS-SEM）。

研究结果：
3.1 水文情势与粪肥施用对土壤抗性组结构的交互影响：通过qPCR、PCoA及重复测量线性混合效应模型发现，22类ARGs和3个MGE相关目标均在稻–蚕豆轮作全周期检出，sul-ARGs最丰富；粪肥在P1使MU和MP总ARG丰度分别较CU高4.1倍和7.7倍，并影响13个ARGs和intI1；时间×粪肥交互影响25个目标中的22个，时间×粪肥×种植三因素交互对20个目标仍显著。总ARG丰度在P1–P2间急剧下降，水稻生长季相对稳定，P5休耕过渡期反弹，蚕豆季再次下降；intI1与总ARG在粪肥后和水稻生长季相关，休耕和蚕豆季解耦。结论为粪肥脉冲并不均匀持续，种植和水文阶段对粪肥富集的抗性组形成过滤，抗性组由粪肥关联的ARG–intI1组合转向受作物阶段和水文条件过滤的后期抗性组。

3.2 水稻–蚕豆轮作驱动典型ARGs的时间演替：通过变异系数分析和核心/可变抗性组框架发现，qnrB和tetL属核心抗性组成员，strA、strB、tetX、tetQ和tetO等为可变抗性组成员；P5休耕过渡期是抗性组重塑关键窗口，排水触发总ARG池同步下降和全面重构，随后好氧蚕豆栽培期ARG丰度逐步达到新平衡；由淹水厌氧向好氧旱地条件转变后，多数ARGs相对丰度明显增加，表明好氧土壤微生境和豆科根际效应促进多种携带ARG的微生物类群增殖，蚕豆旱作期是抗性组回升而非衰减的关键窗口。结论为土壤ARG时间演替并非线性，而是与水分情势转变耦合的阶段性重构，P5是急剧转折点并产生遗留效应。

3.3 循环淹水–非淹水情势下微生物群落组装与演替：通过16S rRNA测序、α多样性、β多样性及PERMANOVA发现，MP处理周年轮作中微生物周转明显；Actinomycetota从P1的42.77%降至P9的31.55%，Chloroflexota、Acidobacteriota和Pseudomonadota相对丰度持续增加；Bacillota在水稻淹水季增加超过221%并于P5急剧下降，Pseudomonadota由P1的16.76%升至P9的19.61%；119个细菌属构成核心微生物组， accessory microbiome呈高时间可塑性；Chao1、Shannon和observed features在P4和P6出现拐点；PERMANOVA显示period、treatment及交互均显著，period的R²=0.767。结论为水文情势而非粪肥或种植是塑造整体土壤微生物组装的主导驱动因子，淹水–非淹水交替作为强生态过滤器顺序选择不同微生物宿主并引导抗性组演替。

3.4 不同轮作阶段和粪肥处理下ARGs、MGEs与细菌群落的共现模式：通过共现网络分析发现，水稻淹水季MP网络较CP扩大，蚕豆季CP网络缩小而MP仍较大；水稻网络模块度较高（0.743–0.814，14–17个模块），蚕豆网络模块度较低（0.571–0.675，7–8个模块），表明厌氧淹水促进区室化ARG–宿主关联和生态位特化，好氧旱地增强跨类群连接并可能促进HGT；中心ARG随季节和粪肥输入变化，水稻季CP中strB和qnrB为枢纽，MP中blaTEM-1、blaOXA-1和tetM突出，blaOXA-1仅在淹水加粪肥条件下与Bacillota和Actinomycetota属形成广泛正相关；蚕豆季CP中tetL、tetO和fexa居中心，MP中blaampC和mcr-1成为枢纽；宿主关联按水文阶段划分，blaNDM-1宿主可塑性高，intI1宿主由淹水期厌氧类群转向旱地好氧Actinomycetota和Pseudomonadota，Candidatus Nitrotoga与intI1跨季稳定相关。结论为水文情势是ARG网络重组的主要驱动因子，并改变潜在ARG宿主组合和HGT潜力。

3.5 从碳到抗性：轮作农业生态系统中驱动ARG动态的生态级联：通过RDA、Mantel检验和PLS-SEM发现，MGEs对ARG丰度具有最强直接路径（路径系数=0.575，p < 0.001）；微生物群落结构与MGE丰度负相关（路径系数=?0.458）；RDA中DOC和ROC沿第一轴载荷分别为0.90和0.79，pH方向相反；Mantel检验确认多种ARG亚型与DOC和ROC显著正相关；水文情势对ARGs总体负效应最强，既直接抑制ARGs，又通过强烈促进微生物群落经微生物–MGE途径间接调控ARGs；β-葡萄糖苷酶和蔗糖酶作为土壤氧化还原状态敏感生物指示物；粪肥施用总体正效应但作用路径多样，可直接抑制某些土壤性质，同时直接促进酶活性和ARGs。结论为水文转变、活性有机碳、微生物群落、MGEs和ARGs构成层级通路，碳–微生物–MGE通路为抗性组演替提供机制解释。

讨论总结：讨论强调MGEs是ARGs最强直接驱动因子，微生物群落周转改变MGE库规模或组成，DOC和ROC等活性有机碳组分是调控微生物群落结构并进而影响抗性组的核心环境过滤器；水文情势具有最强总体负效应，粪肥总体呈正效应但路径更复杂，酶活性变化反映淹水–非淹水转变下的氧化还原状态。讨论还指出，不同作物可能通过根系生长、养分吸收、残体输入、根际招募和根系代谢物释放影响ARG动态，未来需在控制水分和施肥条件下比较多种作物，并整合植物与根际代谢组学、宏基因组学和宏转录组学以阐明作物特异性代谢过程如何调控微生物群落、MGEs和ARG归趋。

研究结论（翻译）：本研究证明，水文情势主导粪肥施用稻–蚕豆轮作系统中土壤抗性组的时间演替。在整个周年轮作周期中，粪肥引入多样ARGs和MGEs，而循环淹水–非淹水转变重塑ARG丰度和组成，休耕过渡期作为抗性组重构的关键窗口。多变量环境分析进一步揭示，活性有机碳组分，尤其DOC和ROC，是连接水文变化与ARG动态的关键生态驱动因子。这些碳介导效应通过微生物群落结构和MGE丰度变化发挥作用，形成碳–微生物–MGE调控级联。网络和宿主关联分析表明，水文转变重组ARG–MGE–细菌关联，改变跨轮作阶段的潜在ARG宿主组合和HGT潜力。总体而言，研究结果突出水文管理和关键过渡期是减轻稻–旱农业生态系统中粪肥土壤ARG风险的重要靶点。