基于宏基因组学的限氧条件下PET微塑料诱导宿主关联抗生素抗性基因富集机制研究

《Water》：Metagenomic Insights into Host-Associated Enrichment of Antibiotic Resistance Genes Under Oxygen-Limited Conditions Induced by PET Microplastics Yinhe Zhao, Jun Li, Kunpeng Jiang, Zhaoming Zheng and Zehao Zhang

　　

抗生素抗性基因（antibiotic resistance genes，ARGs）作为废水处理系统中日益受关注的新兴污染物，其在曝气不足导致的溶解氧（dissolved oxygen，DO）限制条件下的响应机制，尤其是微塑料共存场景下的行为特征尚不明晰。研究人员构建了3组序批式反应器（sequencing batch reactors，SBRs），在31 d运行周期内设置渐进式限氧条件，并添加不同浓度聚对苯二甲酸乙二醇酯（polyethylene terephthalate，PET）微塑料，通过宏基因组学技术系统解析了二者联合作用对处理效能、微生物群落、ARGs、移动遗传元件（mobile genetic elements，MGEs）及PET降解相关基因的影响。结果表明，长期限氧条件下有机物去除效能维持相对稳定，但铵态氮去除能力与污泥沉降性能呈渐进式恶化，PET的投加显著加剧了上述负面效应。PET暴露大幅提升了ARGs的绝对丰度，但未显著改变抗性组组成与主导抗性机制，表明其作用是放大而非重构抗性组。相关性分析显示，ARGs富集以宿主关联为主，由有限类群微生物的增殖驱动；多个潜在ARG宿主同时与PET降解相关基因存在强关联，提示PET功能与抗生素抗性由共享微生物类群承载。此外，ARGs与MGEs的强正相关表明基因流动性在限氧条件下抗性组动态中发挥关键作用。综上，限氧与PET微塑料的联合作用促进了废水系统中宿主关联的ARG富集，凸显了其环境与公共健康风险，强调维持运行稳定性对阻断抗性基因扩散的重要性。

本研究针对废水处理厂曝气不足引发的溶解氧（DO）限制与微塑料复合污染场景，系统揭示了聚对苯二甲酸乙二醇酯（PET）微塑料对抗生素抗性基因（ARGs）传播的调控机制，全文发表于《Water》。当前研究存在三方面空白：一是实际工程中普遍存在的非优化型DO限制对ARGs的驱动效应缺乏系统性解析；二是多数微塑料与抗性基因的交互研究集中于好氧短期实验，忽略了限氧这一常见工况的调控作用；三是PET降解功能与抗性组的宿主耦合关系尚未明确。研究人员通过序批式反应器（SBR）模拟实际工程限氧梯度，结合宏基因组学与网络分析，证实PET微塑料通过促进丝状菌宿主增殖，在限氧条件下显著放大ARGs丰度，且抗性组与PET降解功能存在宿主层面的耦合，为废水处理系统的抗性风险防控提供了新的理论依据。

关键技术方法方面，研究人员采集北京高碑店污水处理厂二沉池回流污泥作为接种物，构建3组有效容积1.0 L的SBR，设置DO渐进降低的3个阶段（0.8–1.0 mg/L、0.5–0.6 mg/L、0.2–0.4 mg/L），分别对应节能曝气下限、硝化菌半饱和临界值与严重限氧工况，其中R0为无PET对照组，R1、R2分别投加1 mg/L与5 mg/L的10 μm PET微塑料。运行周期包含进水、厌氧、曝气、沉淀、排水与闲置阶段，水力停留时间（HRT）24 h，固体停留时间（SRT）约20 d。运行结束后采集污泥样品，采用宏基因组鸟枪法测序，通过Prodigal预测开放阅读框（ORFs），注释ARGs、移动遗传元件（MGEs）与PET降解相关基因，并基于Pearson相关系数（r≥0.9，p≤0.01）构建微生物-基因关联网络，识别潜在宿主类群。

研究结果部分，首先在3.1节“限氧条件下的处理性能恶化”中，通过连续31 d的水质监测发现，所有反应器的化学需氧量（COD）去除率维持在80%左右，仅后期略有下降，表明异养菌代谢对限氧的耐受性较强；而铵态氮（NH₄⁺–N）去除效率呈渐进式下降，R0末期降至约50%，R1、R2分别降至30%–40%与低于40%，证实PET以浓度依赖方式抑制硝化过程；污泥体积指数（SVI）从初始约100 mL/g持续上升，R1超过200 mL/g，R2达300 mL/g以上，表明限氧与PET共同诱发丝状菌膨胀，导致污泥沉降性能恶化。

3.2节“限氧与PET暴露下的ARG特征”显示，限氧本身使ARGs总丰度较初始水平提升18%，PET投加进一步将增幅提升至34%（R1）与38%（R2）。抗性组中四环素类与糖肽类抗性基因占比最高，合计超过40%，大环内酯类与氨基糖苷类次之；抗性机制以外排泵（占43.9%–46.7%）与靶点修饰（占36.3%–38.3%）为主导，二者合计占比超80%，PET暴露仅轻微提升外排机制的占比，未重构抗性组结构，证实PET主要发挥丰度放大器作用。

3.3节“微生物群落与MGE响应”表明，限氧驱动了显著的群落重构：初始优势菌属Candidatus_Competibacter相对丰度从15.65%降至R2的1.79%，而Hyphomicrobium从3.13%增至R1的26.22%与R2的25.80%。丝状菌类群呈浓度依赖性富集，其中Thiothrix从0.02%飙升至R2的29.96%，所有丝状菌总丰度从初始7.22%升至R2的37.94%，与SVI上升趋势完全吻合。MGEs以转座酶（tnpA，占50%–65%）与插入序列（IS）为核心，整合子相对丰度在PET暴露组提升至12.6%左右，表明基因水平转移潜力增强。

3.4节“ARGs的潜在微生物宿主”通过网络分析识别出103组高置信度宿主-ARG关联，富集于少数类群：Hyphomicrobium与Caldilinea各与11种ARGs强相关，涵盖mlaF、tetA(58)、TolC等外排相关基因；Filomicrobium、Ignavibacterium与Candidatus_Denitrolinea分别关联10、9、9种ARGs。值得注意的是，Caldilinea与Candidatus_Denitrolinea均为典型丝状菌，其富集与SVI升高直接耦合，证实ARG积累与丝状菌增殖紧密关联。MGE宿主同样集中于Anaerolinea与Caldilinea等丝状菌，且与tnpA等转座元件强相关，提示PET通过促进丝状菌增殖同时提升ARG丰度与转移潜力。

3.5节“PET降解基因的潜在宿主”检测到PET降解相关基因绝对丰度随PET浓度显著上升，R2较R0提升近80%。相关性分析识别出268组微生物-PET基因关联，Ignavibacterium、Candidatus_Denitrolinea、Candidatus_Gottesmanbacteria与Anaerolinea各与12种PET降解基因强相关，其中Anaerolinea与Candidatus_Denitrolinea同时为ARGs与MGEs的核心宿主，证实抗性组、PET降解功能与基因流动性由同一批丝状菌类群承载，形成功能耦合。

3.6节“环境因子与ARG富集的关系”显示，ARGs与出水COD、NH₄⁺–N浓度呈强正相关，与对应去除率呈强负相关，表明处理效能恶化与抗性组扩增紧密耦合；SVI是所有环境因子中与ARGs相关性最强的指标，证实污泥膨胀是ARG积累的关键驱动因素；ARGs与胞外聚合物（EPS）中多糖组分呈正相关，与蛋白组分呈负相关，提示松散的絮体结构为ARG宿主提供了保护性微生境。

讨论与结论部分指出，曝气不足导致的限氧条件本身会诱发污泥膨胀与硝化抑制，PET微塑料以浓度依赖方式加剧上述效应。在功能基因层面，PET主要通过提升ARGs绝对丰度发挥作用，而非重构抗性组结构，该效应由宿主选择驱动——限氧与PET共同促进携带ARGs、MGEs与PET降解基因的丝状菌增殖，实现抗性传播与塑料降解的功能耦合。研究强调，维持废水处理系统的运行稳定性、控制微塑料输入是阻断ARG扩散的关键举措。未来需进一步优化工艺参数、开发微塑料高效去除技术，并深入解析微塑料暴露下ARG水平转移的分子机制，为降低废水处理设施的抗性风险提供支撑。