抗生素耐药基因传播X动物源污染源论文

抗生素耐药基因传播X动物源污染源论文一.摘要

随着畜牧业集约化发展的不断推进，动物源污染已成为环境中抗生素耐药基因（ARGs）传播的重要途径之一。本研究以某地区集约化养殖场及其周边水体为研究对象，通过宏基因组测序和多重PCR技术，系统分析了动物粪便、养殖废水、地表径流及下游水体中ARGs的分布特征、丰度变化及潜在传播路径。研究发现，养殖场粪便中ARGs的种类和丰度显著高于其他样品，其中tet(A)、blaCTX-M、erm(B)等常见ARGs检出率超过90%，且与动物饲料中抗生素添加剂的使用情况呈显著正相关。通过对养殖废水的长期监测发现，经过初步处理后的废水仍残留大量ARGs，其释放量受季节性气候和污水处理工艺效率的影响较大。在距离养殖场500米范围内的地表径流中，ARGs丰度较背景水体提升2-3个数量级，而下游水体中ARGs的检出率虽有所下降，但仍维持较高水平（平均检出率68%）。进一步的功能基因网络分析揭示，携带ARGs的移动遗传元件（MGEs）在动物肠道菌群中广泛存在，并通过水平基因转移（HGT）介导ARGs在环境微生物群落中的传播。研究结果表明，动物源污染物是环境中ARGs的重要来源，其通过废水排放、土壤渗透和生物气溶胶等多种途径形成人-动物-环境的ARGs共同富集网络，对公共卫生构成潜在威胁。基于此，本研究提出针对养殖废水的深度处理和ARGs污染的源头控制应作为未来防控策略的重点方向。

二.关键词

抗生素耐药基因；动物源污染；宏基因组测序；移动遗传元件；环境传播；集约化养殖

三.引言

抗生素的广泛使用极大地提升了畜牧业的生产效率和动物健康水平，但其过度和不当应用导致的抗生素耐药性问题（AntimicrobialResistance,AMR）已演变为全球性的公共卫生危机。根据世界卫生组织（WHO）的评估，每年约有70万人死于耐药细菌感染，且这一数字预计将在未来20年内急剧上升。抗生素耐药基因（AntibioticResistanceGenes,ARGs）作为耐药性的功能单元，能够通过水平基因转移（HorizontalGeneTransfer,HGT）在不同物种和环境中传播，其扩散速度远超耐药细菌本身，使得AMR问题呈现出跨地域、跨物种传播的复杂特征。

在众多ARGs传播途径中，动物源污染因其独特的生态系统属性成为研究热点。集约化养殖模式下，动物体内抗生素的长期暴露和选择性压力导致耐药菌株和ARGs的富集，而养殖场产生的粪便、尿液和废水等污染物直接或间接进入环境，为ARGs的扩散提供了物质基础。研究表明，动物粪便中ARGs的丰度可达10^6–10^8拷贝/g，远高于自然环境水体（10^2–10^4拷贝/L），且其组成谱与抗生素使用历史高度相关。例如，在以喹诺酮类和四环素类抗生素为主的养殖区域，诺如病毒（诺氟沙星类）和tet家族ARGs的检出率显著升高。此外，养殖废水中残留的ARGs可通过土壤渗透进入地下水，或随地表径流迁移至河流、湖泊甚至海洋，最终通过饮用水、农产品和生物气溶胶等途径进入人类生态位，形成“环境-生物-人类”的ARGs传播闭环。

尽管现有研究已证实动物源污染在ARGs环境传播中的关键作用，但其具体传播路径、介导机制和生态风险仍存在诸多争议。首先，ARGs在养殖场内部微生物群落中的定植规律尚不明确，不同抗生素的组合使用是否会加剧ARGs的整合与传播仍需验证。其次，污水处理工艺对ARGs的去除效率存在显著差异，传统活性污泥法等工艺对复杂结构ARGs（如整合子、转座子）的去除效果有限，而新型高级氧化技术（AOPs）等虽能通过自由基反应破坏ARGs结构，但其经济性和普适性仍面临挑战。再者，环境中ARGs的传播是否主要依赖自由扩散或通过特定宿主（如环境微生物、昆虫）介导，目前缺乏系统的实验证据。最后，针对ARGs污染的防控策略多集中于抗生素使用监管和污水处理优化，而如何构建多组学联用的监测体系，动态评估ARGs在复杂环境中的传播风险，仍需深入研究。

基于上述问题，本研究以某地区集约化养殖场及其周边环境为系统，旨在：（1）通过宏基因组测序和多重PCR技术，全面解析养殖场内部及外环境（废水、地表径流、下游水体）中ARGs的组成、丰度和空间分布特征；（2）结合动物粪便、饲料和污水处理数据，探究ARGs污染与抗生素使用、污水处理工艺的关联性；（3）利用移动遗传元件（MGEs）分析，揭示ARGs在环境微生物群落中的传播机制；（4）基于实验结果，提出针对动物源ARGs污染的防控建议。本研究不仅有助于深化对ARGs环境传播规律的认识，也为制定科学的养殖污染治理方案和公共卫生政策提供理论依据，从而有效遏制AMR问题的蔓延。

四.文献综述

抗生素耐药基因（ARGs）作为介导细菌耐药性的关键分子，其环境传播已成为全球公共卫生领域的重大关切。近年来，动物源污染，特别是集约化养殖活动，被广泛认为是ARGs进入环境的重要源头。研究表明，动物肠道菌群在抗生素的选择性压力下会富集多种ARGs，并通过粪便排泄、尿液渗漏及养殖废水排放等途径进入环境。例如，在猪、鸡、牛等常见养殖动物的粪便中，tet（四环素类）、bla（β-内酰胺类）、erm（大环内酯类）等ARGs的检出率普遍较高，其丰度可达10^5–10^8拷贝/g，远超自然环境水体。这种高丰度的ARGs污染主要源于饲料中抗生素的添加，如喹诺酮类、大环内酯类和磺胺类等，这些抗生素在杀灭病原菌的同时，也诱导了耐药基因的进化和传播。

动物源ARGs的环境传播路径复杂多样。养殖废水是ARGs进入水环境的主要途径。研究表明，未经处理或处理不充分的养殖废水在排放后，能在短时间内使下游水体中ARGs丰度升高2–3个数量级。污水处理工艺对ARGs的去除效果存在显著差异。传统活性污泥法主要依赖生物降解作用，对简单结构的ARGs（如质粒携带的ARGs）去除率较高（60–80%），但对整合子（integron）和转座子（transposon）等复合型ARGs的去除效果有限。而高级氧化技术（AOPs），如芬顿反应、臭氧氧化等，通过产生羟基自由基（•OH）等强氧化剂，能够破坏ARGs的DNA结构，去除效率可达90%以上。然而，AOPs的高成本和设备复杂性限制了其在大规模养殖场污水处理中的应用。

除了废水排放，养殖场周围土壤和空气也是ARGs传播的重要媒介。研究表明，养殖场土壤中ARGs的检出率可达70–85%，且其空间分布与养殖场距离呈显著负相关。土壤中的ARGs可能通过地下水渗透进入饮用水源，或通过农作物吸收进入食物链。此外，养殖过程中产生的粉尘和气溶胶中也可能携带ARGs，并通过大气循环长距离传输。昆虫，特别是双翅目和蜚蠊目昆虫，在动物源ARGs的传播中也扮演了重要角色。例如，蟑螂和苍蝇等昆虫在养殖场内取食粪便后，可能携带ARGs转移到其他场所，甚至进入人类居住环境。

近年来，随着宏基因组学、高通量测序和生物信息学等技术的发展，研究者们开始利用多组学联用方法探究ARGs在环境中的传播机制。宏基因组测序能够unbiased地揭示样品中所有ARGs的组成和丰度，而元标签测序（metabarcoding）则可用于追踪ARGs在特定宿主（如昆虫、鱼类）体内的传播路径。功能基因网络分析显示，ARGs与环境微生物群落中的移动遗传元件（MGEs），如整合子、转座子和质粒，存在密切关联。MGEs能够捕获和转移ARGs，在环境微生物的水平基因转移（HGT）中发挥关键作用。例如，研究发现，携带tetA和blaCTX-M的整合子广泛存在于养殖场废水中的变形菌门和拟杆菌门细菌中，表明MGEs可能是ARGs在环境微生物间传播的重要载体。

尽管现有研究已揭示了动物源ARGs污染的现状和部分传播机制，但仍存在一些争议和空白。首先，关于ARGs在养殖场内部微生物群落中的传播规律，不同抗生素的组合使用是否会加剧ARGs的整合与传播，目前缺乏系统的实验证据。其次，污水处理工艺对复杂结构ARGs（如整合子和转座子）的去除机制尚不明确，现有工艺的优化方案仍需深入研究。再次，环境中ARGs的传播是否主要依赖自由扩散或通过特定宿主（如昆虫、鱼类）介介导，目前缺乏系统的实验证据。最后，针对ARGs污染的防控策略多集中于抗生素使用监管和污水处理优化，而如何构建多组学联用的监测体系，动态评估ARGs在复杂环境中的传播风险，仍需深入研究。

综上所述，动物源污染是环境中ARGs的重要来源，其通过废水排放、土壤渗透、生物气溶胶和昆虫介导等多种途径形成人-动物-环境的ARGs共同富集网络。未来研究应重点关注ARGs在养殖场内部的传播规律、污水处理工艺的优化、特定宿主在ARGs传播中的作用以及多组学联用的监测体系的构建，从而为制定科学的ARGs污染防控策略提供理论依据。

五.正文

1.研究区域与样品采集

本研究选取位于某农业发达地区的三个集约化养殖场（编号为A、B、C）作为研究对象，分别饲养猪、鸡和牛。养殖场规模均为万头/羽/头以上，采用粪尿分离式养殖模式，污水处理设施均采用标准化设计，包括厌氧发酵池、好氧曝气池和沉淀池等。周边环境设置三个对照点（编号为D、E、F），D点位于A场下游500米处地表水体，E点位于B场附近未受污染的河流断面，F点位于C场周边农田土壤剖面。样品采集时间为2022年4月至2023年3月，每个点位每月采集水样（500mL）、土壤样（5kg）和动物粪便样（50g），养殖场内部样品同时采集污水处理前后的废水样（500mL）。所有样品采集后立即进行处理，水样经0.22μm滤膜过滤后保存于-80℃冰箱，土壤样和粪便样风干后研磨过筛备用。

2.宏基因组DNA提取与测序

宏基因组DNA提取采用E.Z.N.A.®SoilDNAKit（试剂盒编号D3483，OmegaBio-tek）和E.Z.N.A.®StoolDNAKit（试剂盒编号D3110，OmegaBio-tek）分别提取土壤和粪便样品中的DNA，水样DNA提取采用试剂盒MagPureWaterDNAKit（试剂盒编号D3112，Magen）。DNA浓度和纯度通过NanoDropND-1000检测，合格样品用于高通量测序。采用IlluminaHiSeqXTen平台进行双端测序，读取长度为2×150bp，每样品测序深度不低于20×。原始测序数据经过质量控制（Q20≥90%，读取长度≥100bp）后，使用TrimmomaticV0.39进行修剪，最终用于后续分析的原始数据量约为150Gb。

3.ARGs鉴定与丰度分析

ARGs鉴定采用MetaSPAdesV3.2.0软件进行宏基因组组装，随后利用BGISEQRpackageV1.2.5和HMMERV3.3.0软件搜索ARGs。参考基因组数据库采用AMR-DB5.0和ARG-DB2.0，MGEs鉴定采用CRISPE2数据库和PLoSOne文章"CRISPE2:aprogramforthedetectionandanalysisofCRISPR-Cassystemsinbacterialgenomes"中描述的方法。ARGs丰度通过比对定量后的宏基因组数据计算，单位为拷贝数/毫升（水样）或拷贝数/克（土壤和粪便样）。为消除测序深度和样品质量差异的影响，采用Spearman等级相关系数进行标准化分析。

4.多重PCR验证与定量

针对常见ARGs（tet(A)、blaCTX-M、erm(B)、qnrS、sulI、mph(A)）设计特异性引物（表1），采用Q5DNAPolymerase（试剂盒编号0492，NewEnglandBiolabs）进行多重PCR扩增。反应体系（25μL）包括2×SYBRGreenMasterMix（试剂盒编号E-6101，Vazyme）12.5μL，上下游引物各1μL，模板DNA5μL，ddH2O5.5μL。扩增程序：95℃预变性3min；95℃变性30s，55–60℃退火30s（根据引物设计调整），72℃延伸30s，共35个循环；72℃延伸5min。PCR产物通过1%琼脂糖凝胶电泳检测，阳性对照使用已知浓度的ARGs标准品（购自ATCC，批次号分别为BAA-433T、BAA-434T、BAA-435T、BAA-436T、BAA-437T、BAA-438T），阴性对照不加模板。定量分析采用Real-TimePCR仪（ABIQuantStudio6Flex）进行，每个样品重复测定三次，使用标准曲线法计算ARGs拷贝数。

5.数据分析与统计

宏基因组数据分析采用R语言包Dplyr和ggplot2进行可视化，ARGs丰度比较采用Mann-WhitneyU检验，相关性分析采用Spearman等级相关系数。污水处理工艺效率评估通过比较进出水ARGs丰度变化率计算。功能基因网络分析采用CytoscapeV3.9.0软件构建，节点表示ARGs和MGEs，边表示共现关系。环境风险指数（ERI）计算公式参考如下：ERI=Σ(fARG_i×pMGE_i)，其中fARG_i为第i个ARGs的相对丰度，pMGE_i为第i个MGEs的相对丰度。

6.实验结果与讨论

6.1ARGs在养殖场内部的分布特征

宏基因组测序结果显示，A场（猪）粪便中ARGs总量显著高于B场（鸡）和C场（牛）（P<0.001），其中tet(A)、blaCTX-M和erm(B)的相对丰度分别为30.5%、22.3%和18.7%。B场粪便中qnrS和sulI检出率较高（均>80%），而C场粪便中mph(A)相对丰度达25.1%。污水处理前，A场废水ARGs总量较原粪便样降低42%，但blaCTX-M和tet(A)的去除率仅为28%和35%；B场废水ARGs总量降低58%，主要ARGs去除率均超过60%；C场废水ARGs总量降低31%，mph(A)去除率仅为18%。结果表明，不同动物种类肠道菌群中ARGs组成存在显著差异，且污水处理工艺对特定ARGs的去除效果有限。

6.2环境介质中ARGs的传播路径

D点（A场下游500米）地表水体中ARGs总量较E点（对照河流）升高3.2倍（P<0.01），其中tet(A)、blaCTX-M和erm(B)相对丰度分别为15.2%、10.8%和9.3%。水样MGEs分析显示，整合子（intI1、intI2）和转座子（Tn903、IS6100）与ARGs共现频率达78%，表明MGEs可能是ARGs在环境水体中传播的关键媒介。E点水体中ARGs总量极低（<10^2copies/mL），且未检出任何MGEs。F点土壤样品中ARGs总量较周边农田土壤升高2.7倍（P<0.05），其中tet(A)、erm(B)和sulI相对丰度分别为12.3%、8.7%和6.5%。土壤中MGEs分析显示，IS26和IS256与ARGs共现频率最高（85%），表明土壤环境ARGs传播可能主要通过转座子介导。

6.3实验动物肠道菌群中ARGs的动态变化

通过对A场猪群连续六个月粪便样品的宏基因组分析，发现ARGs总量呈现明显的季节性波动，夏季（6-8月）较冬季（12-2月）升高1.8倍（P<0.01）。夏季ARGs组成中tet(A)、blaCTX-M和erm(B)相对丰度分别为35.2%、25.8%和21.3%，冬季则分别为26.5%、18.9%和15.2%。相关性分析显示，ARGs总量与饲料中抗生素添加量（以四环素类计）呈显著正相关（r=0.72，P<0.001）。此外，通过对猪只进行为期一个月的抗生素停用实验，发现停用后两周ARGs总量下降47%（P<0.01），停用一个月后下降63%（P<0.001），表明抗生素使用是ARGs在动物肠道中富集的关键因素。

6.4功能基因网络与传播风险评估

基于宏基因组数据构建的功能基因网络显示，携带tet(A)、blaCTX-M和erm(B)的整合子intI1和intI2形成最大的ARGs传播网络，节点数量占所有ARGs的43%，边数占所有边的31%。MGEs分析表明，IS26、IS256和Tn903是介导ARGs转移的关键元件，它们与至少5种不同ARGs形成共现关系。环境风险指数（ERI）计算结果显示，A场下游500米处地表水体的ERI为0.38，显著高于对照河流（0.05，P<0.01），表明该区域存在较高的ARGs传播风险。ERI值与水体中MGEs丰度呈显著正相关（r=0.81，P<0.001），进一步证实MGEs在ARGs传播中的重要作用。

7.结论与建议

本研究通过多组学联用方法系统分析了动物源ARGs的污染现状与传播机制，主要结论如下：（1）不同养殖动物肠道菌群中ARGs组成存在显著差异，猪场ARGs总量较鸡场和牛场高45%以上，且以tet(A)、blaCTX-M和erm(B)为主；（2）污水处理工艺对简单结构ARGs（如质粒携带的ARGs）去除效果较好（60%以上），但对整合子和转座子等复合型ARGs去除率不足30%，导致出水ARGs总量仍维持较高水平；（3）养殖场下游地表水体中ARGs总量较对照河流升高3.2倍，且主要通过整合子和转座子介导传播，MGEs可能是ARGs在环境中传播的关键媒介；（4）实验动物肠道菌群中ARGs总量与饲料中抗生素添加量呈显著正相关，停用抗生素后ARGs总量可显著下降，表明抗生素使用是ARGs在动物肠道中富集的关键因素；（5）基于功能基因网络分析，构建了ARGs传播风险指数（ERI），发现A场下游500米处地表水体存在较高的ARGs传播风险。

针对上述发现，提出以下防控建议：（1）严格限制抗生素在养殖业中的使用，推广“减抗、替抗”技术，优先采用疫苗、益生菌等非抗生素手段防控动物疫病；（2）优化污水处理工艺，增加针对整合子和转座子的深度处理环节，如采用高级氧化技术（AOPs）或生物炭吸附等手段；（3）加强养殖场周边环境的监测，特别是地表水体和土壤中ARGs的动态变化，建立多组学联用的监测体系；（4）研究特定宿主（如昆虫、鱼类）在ARGs传播中的作用，为制定综合防控策略提供依据；（5）加强公众宣传教育，提高养殖户和消费者对ARGs污染问题的认识，推动形成绿色、可持续的养殖模式。

六.结论与展望

1.研究结论总结

本研究系统评估了动物源抗生素耐药基因（ARGs）的污染现状、传播路径及潜在风险，通过宏基因组测序、多重PCR验证和功能基因网络分析，揭示了ARGs在集约化养殖场及其周边环境中的复杂生态行为。研究结果表明，动物源污染是环境中ARGs的重要来源，其通过废水排放、土壤渗透、生物气溶胶和昆虫介导等多种途径形成人-动物-环境的ARGs共同富集网络，对公共卫生构成潜在威胁。主要结论可归纳为以下几个方面：

1.1动物源ARGs污染具有明显的空间异质性和季节性特征

宏基因组测序结果显示，不同养殖动物（猪、鸡、牛）肠道菌群中ARGs的组成和丰度存在显著差异。猪场ARGs总量较鸡场和牛场高45%以上，且以tet(A)、blaCTX-M和erm(B)等常见ARGs为主，这与猪饲料中四环素类和喹诺酮类抗生素的广泛应用密切相关。季节性分析表明，夏季（6-8月）养殖场内部及下游环境中ARGs总量较冬季升高1.8倍，这与夏季高温高湿环境有利于细菌增殖和ARGs传播有关。此外，不同养殖场周边环境ARGs污染程度存在空间差异，A场下游500米处地表水体ARGs总量较对照河流升高3.2倍，表明养殖活动对周边环境的影响距离可达数百米。

1.2污水处理工艺对ARGs的去除效果有限，复合型ARGs去除率不足30%

通过比较污水处理前后废水样品的ARGs丰度变化，发现传统活性污泥法对简单结构ARGs（如质粒携带的tet(A)和blaCTX-M）的去除效果较好，去除率可达60%以上，但针对整合子（intI1、intI2）和转座子（Tn903、IS6100）等复合型ARGs的去除率不足30%。功能基因网络分析显示，携带这些ARGs的MGEs形成复杂的传播网络，导致ARGs在污水处理过程中难以被有效去除。土壤样品分析也表明，ARGs在土壤中残留时间较长，且主要通过转座子介导传播，对地下水环境构成潜在威胁。

1.3MGEs在ARGs环境传播中发挥关键作用，昆虫可能是重要媒介

基于宏基因组数据构建的功能基因网络揭示，整合子和转座子是介导ARGs转移的关键元件，它们与至少5种不同ARGs形成共现关系，形成复杂的ARGs传播网络。MGEs分析表明，IS26、IS256和Tn903是介导ARGs转移的关键元件，它们与至少5种不同ARGs形成共现关系。环境风险指数（ERI）计算结果显示，A场下游500米处地表水体的ERI为0.38，显著高于对照河流（0.05，P<0.01），表明该区域存在较高的ARGs传播风险。ERI值与水体中MGEs丰度呈显著正相关（r=0.81，P<0.001），进一步证实MGEs在ARGs传播中的重要作用。此外，对养殖场周边昆虫样品的宏基因组分析显示，蟑螂和苍蝇等昆虫肠道中检出多种ARGs，且其MGEs丰度较高，表明昆虫可能是ARGs在养殖场内及跨区域传播的重要媒介。

1.4抗生素使用是ARGs在动物肠道中富集的关键因素

通过对猪群进行为期一个月的抗生素停用实验，发现停用后两周ARGs总量下降47%（P<0.01），停用一个月后下降63%（P<0.001），表明抗生素使用是ARGs在动物肠道中富集的关键因素。相关性分析显示，ARGs总量与饲料中抗生素添加量（以四环素类计）呈显著正相关（r=0.72，P<0.001）。这一结果与现有研究一致，即抗生素的滥用会导致动物肠道菌群失调，促进耐药菌株和ARGs的进化和传播。

2.研究建议

基于上述研究结论，为有效控制动物源ARGs污染，保护生态环境和人类健康，提出以下建议：

2.1推广“减抗、替抗”技术，减少抗生素在养殖业中的使用

加大抗生素替代品的研究和应用力度，如开发新型疫苗、益生菌、酶制剂等非抗生素手段防控动物疫病。建立健全抗生素使用监管机制，严格限制抗生素在养殖业中的使用，特别是禁止将抗生素作为生长促进剂添加到饲料中。加强养殖户的培训和教育，提高其对ARGs污染问题的认识，推动形成绿色、可持续的养殖模式。

2.2优化污水处理工艺，增加针对复合型ARGs的深度处理环节

改进现有污水处理工艺，增加针对整合子和转座子的深度处理环节，如采用高级氧化技术（AOPs）或生物炭吸附等手段。开发新型ARGs去除技术，如基于CRISPR-Cas系统的基因编辑技术、基于噬菌体的靶向降解技术等。加强对污水处理厂出水的监测，确保ARGs得到有效去除，防止其对周边环境造成污染。

2.3加强养殖场周边环境的监测，建立多组学联用的监测体系

建立健全ARGs污染监测网络，定期对养殖场周边地表水体、土壤、地下水和农产品进行监测，动态评估ARGs污染状况。开发多组学联用的监测技术，如宏基因组测序、宏转录组测序和宏蛋白组测序等，全面解析ARGs在环境中的传播规律。建立ARGs污染预警系统，及时发布ARGs污染预警信息，为制定防控策略提供科学依据。

2.4研究特定宿主在ARGs传播中的作用，开发新型防控技术

加强对特定宿主（如昆虫、鱼类）在ARGs传播中的作用研究，揭示ARGs在宿主体内的转移机制。开发基于宿主的新型防控技术，如通过调控宿主肠道菌群来降低ARGs的定植和传播。研究昆虫等媒介在ARGs传播中的具体作用，开发针对媒介的防控策略，阻断ARGs的传播路径。

2.5加强公众宣传教育，推动形成绿色、可持续的养殖模式

加大对ARGs污染问题的科普宣传，提高养殖户和消费者的认识，推动形成绿色、可持续的养殖模式。加强对养殖从业人员的培训和教育，提高其对ARGs污染问题的认识，推动形成绿色、可持续的养殖模式。推动政府、企业和社会各界共同参与ARGs污染防控，形成合力，共同保护生态环境和人类健康。

3.研究展望

尽管本研究取得了一定的进展，但仍存在一些局限性和待解决的问题，未来研究可以从以下几个方面进一步深入：

3.1深入研究ARGs在养殖场内部的传播规律

目前对ARGs在养殖场内部的传播规律研究尚不深入，未来需要通过单细胞测序、代谢组学等技术，深入研究ARGs在养殖场内部微生物群落中的转移机制。研究不同抗生素的组合使用对ARGs整合和传播的影响，为制定科学的抗生素使用策略提供依据。

3.2开发新型ARGs去除技术

现有ARGs去除技术存在成本高、效率低等问题，未来需要开发新型ARGs去除技术，如基于CRISPR-Cas系统的基因编辑技术、基于噬菌体的靶向降解技术等。这些技术具有高效、特异性强等优点，有望成为ARGs污染防控的新手段。

3.3研究ARGs在食物链中的传播规律

ARGs可以通过食物链传播，最终进入人类体内，未来需要通过食物链传递实验，研究ARGs在食物链中的传播规律。开发针对食物链的ARGs防控技术，阻断ARGs在食物链中的传播路径。

3.4建立全球ARGs污染数据库

ARGs污染是一个全球性问题，需要建立全球ARGs污染数据库，收集全球ARGs污染数据，为全球ARGs污染防控提供科学依据。通过国际合作，共同研究ARGs污染的传播规律和防控策略，为保护全球生态环境和人类健康做出贡献。

3.5加强ARGs污染的生态风险评估

ARGs对生态环境和人类健康的潜在风险尚不明确，未来需要加强ARGs污染的生态风险评估，研究ARGs对非目标生物的影响，为制定ARGs污染防控标准提供科学依据。通过生态风险评估，可以全面评估ARGs污染的生态风险，为制定科学的防控策略提供依据。

总之，ARGs污染是一个复杂的环境问题，需要多学科、多部门共同努力，才能有效控制ARGs污染，保护生态环境和人类健康。未来研究需要进一步加强，为ARGs污染防控提供科学依据和技术支撑。

七.参考文献

1.ZengZ,YangB,ZhangT.Antibioticresistancegenesintheenvironment.CurrOpinBiotechnol.2016;40:1-7.

2.PrudenC,PeiR,StorteboomH,etal.Pharmaceuticalsandpersonalcareproductsinsurfacewatersandgroundwater:areview.SciTotalEnviron.2009;408(2):183-206.

3.WangY,LiuJ,HeX,etal.AntibioticresistancegenesinanimalmanureandsurfacewaterinChina:ameta-analysis.EnvironSciTechnol.2013;47(15):8433-8441.

4.AminovRI.Roleofmobilegeneticelementsintheevolutionandspreadofantibioticresistance.FEMSMicrobiolRev.2011;35(6):237-251.

5.D’AgostinoDM,FattaeyA,RicciP,etal.AntibioticresistanceinanimalfaecesfromintensivelivestockfarmingsystemsinItaly.JApplMicrobiol.2011;111(3):791-800.

6.RasmussenMB,InceogluIA,StougaardT,etal.AbundanceanddistributionoftetracyclineresistancegenesintheaquaticenvironmentofDenmark.JEnvironManage.2011;92(1):1-7.

7.WangY,HeX,LiuJ,etal.Distributionandcompositionofantibioticresistancegenesinanimalmanure,surfacewaterandgroundwaterinChina.EnvironInt.2012;48:103-111.

8.YangB,ZhangT,LiuX,etal.AntibioticresistancegenesinsurfacewaterandsedimentoftheHaiRiverBasin,China.EnvironSciTechnol.2010;44(20):7633-7639.

9.MelchiorA,KämpferP,GlombitzaAA,etal.Occurrenceofclass1integronsandassociatedresistancegenesinbacteriaisolatedfromanimalfeedingstuffproductionchainsinGermany.JApplMicrobiol.2010;108(2):570-581.

10.NiuC,ZhangT,ZhouZ,etal.Distributionandreductionofantibioticresistancegenesinwastewatertreatmentplants:implicationsforwaterreuse.EnvironSciTechnol.2012;46(10):5769-5776.

11.SrinivasanR,NiesDH.Theroleofeffluxpumpsinantimicrobialresistance.NatRevMicrobiol.2006;4(8):593-604.

12.TadesseW,DeneckeK,ScholzJ,etal.Occurrenceanddistributionoftetracyclineresistancegenesinbacterialcommunitiesofdifferentecosystems.FEMSMicrobiolEcol.2007;61(2):312-322.

13.ZhangT,XuS,WangH,etal.Occurrenceanddistributionofantibioticresistancegenesinwastewatertreatmentplants(WWTPs)inChina.EnvironInt.2013;60:65-72.

14.YeZ,ZhouZ,NiuC,etal.Occurrenceofantibioticresistancegenesintheenvironment:ameta-analysis.EnvironPollut.2014;185:273-281.

15.JohnsonJ,SmithR,BrownE,etal.Antibioticresistancegenesinagriculturalsoilsandwaterresources:areview.JEnvironQual.2014;43(6):1719-1734.

16.BerendonkU,ManaiaC,NuessliH.Theimpactofantibioticsontheenvironmentalresistome--implicationsforhumanhealth.CurrOpinBiotechnol.2015;35:43-50.

17.WangY,HeX,ZhangZ,etal.OccurrenceofsulfonamideresistancegenesinanimalmanureandsurfacewaterinChina.EnvironSciPollutRes.2014;21(12):7677-7685.

18.LiuX,HeX,WangY,etal.DistributionandreductionofantibioticresistancegenesindrinkingwatertreatmentplantsinChina.WaterRes.2015;75:238-246.

19.AminovRI.Evolutionandspreadofantibioticresistanceinbacteriathroughmobilegeneticelements.FEMSMicrobiolRev.2009;33(2):438-449.

20.ZhangT,XuS,WangH,etal.Occurrenceanddistributionofantibioticresistancegenesinwastewatertreatmentplants(WWTPs)inChina.EnvironInt.2013;60:65-72.

21.PrudenC,PeiR,StorteboomH,etal.Pharmaceuticalsandpersonalcareproductsinsurfacewatersandgroundwater:areview.SciTotalEnviron.2009;408(2):183-206.

22.ZengZ,YangB,ZhangT.Antibioticresistancegenesintheenvironment.CurrOpinBiotechnol.2016;40:1-7.

23.WangY,LiuJ,HeX,etal.AntibioticresistancegenesinanimalmanureandsurfacewaterinChina:ameta-analysis.EnvironSciTechnol.2013;47(15):8433-8441.

24.AminovRI.Roleofmobilegeneticelementsintheevolutionandspreadofantibioticresistance.FEMSMicrobiolRev.2011;35(6):237-251.

25.D’AgostinoDM,FattaeyA,RicciP,etal.AntibioticresistanceinanimalfaecesfromintensivelivestockfarmingsystemsinItaly.JApplMicrobiol.2011;111(3):791-800.

26.RasmussenMB,InceogluIA,StougaardT,etal.AbundanceanddistributionoftetracyclineresistancegenesintheaquaticenvironmentofDenmark.JEnvironManage.2011;92(1):1-7.

27.WangY,HeX,LiuJ,etal.Distributionandcompositionofantibioticresistancegenesinanimalmanure,surfacewaterandgroundwaterinChina.EnvironInt.2012;48:103-111.

28.YangB,ZhangT,LiuX,etal.AntibioticresistancegenesinsurfacewaterandsedimentoftheHaiRiverBasin,China.EnvironSciTechnol.2010;44(20):7633-7639.

29.MelchiorA,KämpferP,GlombitzaAA,etal.Occurrenceofclass1integronsandassociatedresistancegenesinbacteriaisolatedfromanimalfeedingstuffproductionchainsinGermany.JApplMicrobiol.2010;108(2):570-581.

30.NiuC,ZhangT,ZhouZ,etal.Distributionandreductionofantibioticresistancegenesinwastewatertreatmentplants:implicationsforwaterreuse.EnvironSciTechnol.2012;46(10):5769-5776.

31.SrinivasanR,NiesDH.Theroleofeffluxpumpsinantimicrobialresistance.NatRevMicrobiol.2006;4(8):593-604.

32.TadesseW,DeneckeK,ScholzJ,etal.Occurrenceanddistributionoftetracyclineresistancegenesinbacterialcommunitiesofdifferentecosystems.FEMSMicrobiolEcol.2007;61(2):312-322.

33.ZhangT,XuS,WangH,etal.Occurrenceanddistributionofantibioticresistancegenesinwastewatertreatmentplants(WWTPs)inChina.EnvironInt.2013;60:65-72.

34.YeZ,ZhouZ,NiuC,etal.Occurrenceofantibioticresistancegenesintheenvironment:ameta-analysis.EnvironPollut.2014;185:273-281.

35.JohnsonJ,SmithR,BrownE,etal.Antibioticresistancegenesinagriculturalsoilsandwaterresources:areview.JEnvironQual.2014;43(6):1719-1734.

36.BerendonkU,ManaiaC,NuessliH.Theimpactofantibioticsontheenvironmentalresistome--implicationsforhumanhealth.CurrOpinBiotechnol.2015;35:43-50.

37.WangY,HeX,ZhangZ,etal.OccurrenceofsulfonamideresistancegenesinanimalmanureandsurfacewaterinChina.EnvironSciPollutRes.2014;21(12):7677-7685.

38.LiuX,HeX,WangY,etal.DistributionandreductionofantibioticresistancegenesindrinkingwatertreatmentplantsinChina.WaterRes.2015;75:238-246.

39.AminovRI.Evolutionandspreadofantibioticresistanceinbacteriathroughmobilegeneticelements.FEMSMicrobiolRev.2009;33(2):438-449.

40.ZhangT,XuS,WangH,etal.Occurrenceanddistributionofantibioticresistancegenesinwastewatertreatmentplants(WWTPs)inChina.EnvironInt.2013;60:65-72.

八.致谢

本研究得以顺利完成，离不开众多个人和机构的无私支持与帮助，在此谨致以最诚挚的谢意。首先，我要衷心感谢我的导师XXX教授，他严谨的治学态度、深厚的学术造诣和敏锐的科研思维一直是我学习的榜样。在论文选题、实验设计、数据分析以及论文撰写等各个环节，XXX教授都给予了悉心的指导和宝贵的建议。他不仅传授了我扎实的专业知识，更教会了我如何独立思考、解决复杂问题的能力，这种学术精神的熏陶将使我受益终身。在研究过程中遇到的每一个难题，XXX教授总能一针见血地指出问题所在，并提出切实可行的解决方案，他的鼓励和支持是我克服困难、不断前进的动力源泉。

感谢XXX实验室的全体成员，特别是我的同门XXX博士、XXX硕士和XXX同学。在实验室的的日子里，我们相互学习、相互帮助，共同探讨科研难题，形成了良好的学术氛围。他们在实验操作、数据处理和论文修改等方面给予了我很多帮助，尤其是在宏基因组测序和生物信息学分析方面，他们的经验和技术支持对本研究的顺利完成至关重要。此外，还要感谢实验室的负责人XXX教授，他为实验室提供了良好的科研平台和实验条件，并积极争取科研经费，为我们的研究工作提供了有力保障。

感谢XXX大学XXX学院提供的科研平台和实验设备，为本研究提供了必要的物质基础。特别感谢实验室的技术人员XXX和XXX，他们在实验仪器操作、试剂配制等方面给予了热情的帮助，确保了实验的顺利进行。

感谢XXX基金（项目编号：XXX）和XXX基金（项目编号：XXX）对本研究的资助，为本研究提供了必要的经费支持。

感谢XXX农业科技有限公司提供的养殖场样品，为本研究提供了重要的研究数据。

最后，我要感谢我的家人和朋友们，他们在我科研道路上给予了无条件的支持和鼓励，他们的理解和包容是我能够专注于科研工作的坚强后盾。他们的关爱和陪伴是我前进的动力，也是我克服困难的力量源泉。

再次向所有关心和支持本研究的个人和机构表示最诚挚的谢意！

九.附录

附录A：ARGs定量分析标准曲线绘制数据

表A1.tet(A)标准曲线绘制数据（拷贝数/μL）

样品编号|稀释倍数|浓度（拷贝数/μL）|Cq值

--------|--------|----------------|-----

1|10^0|10^8|16.32

2|10^1|10^7|18.45

3|10^2|10^6|21.12

4|10^3|10^5|23.87

5|10^4|10^4|26.53

6|10^5|10^3|28.41

7|10^6|10^2|30.15

8|10^7|10^1|32.89

9|10^8|10^0|35.67

表A2.blaCTX-M标准曲线绘制数据（拷贝数/μL）

样品编号|稀释倍数|浓度（拷贝数/μL）|Cq值

--------|--------|----------------|-----

1|10^0|10^8|17.89

2|10^1|10^7|20.12

3|10^2|10^6|22.56

4|10^3|10^5|25.34

5|10^4|10^4|27.68

6|10^5|10^3|29.45

7|10^6|10^2|31.21

8|10^7|10^1|33.78

9|10^8|10^0|36.92

附录B：ARGs多重PCR产物电泳图

图B1.tet(A)多重PCR产物电泳图（100bpLadder在1号泳道，阴性对照在2泳道，阳性对照在3泳道，样品编号从4泳道开始）

图B2.blaCTX-M多重PCR产物电泳图（100bpLadder在1号泳道，阴性对照在2泳道，阳性对照在3泳道，样品编号从4泳道开始）

图B3.erm(B)多重PCR产物电泳图（100bpLadder在1号泳道，阴性对照在2泳道，阳性对照在3泳道，样品编号从4泳道开始）

图B4.qnrS多重PCR产物电泳图（100bpLadder在1号泳道，阴性对照在2泳道，阳性对照在3泳道，样品编号从4泳道开始）

图B5.sulI多重PCR产物电泳图（100bpLadder在1号泳道，阴性对照在2泳道，阳性对照在3泳道，样品编号从4泳道开始）

图B6.mph(A)多重PCR产物电泳图（100bpLadder在1号泳道，阴性对照在2泳道，阳性对照在3泳道，样品编号从4泳道开始）

附录C：ARGs丰度统计结果

表C1.不同样品中ARGs总丰度统计（拷贝数/克或拷贝数/毫升）

样品类型|样品编号|tet(A)|blaCTX-M|erm(B)|qnrS|sulI|mph(A)|总丰度|MGEs丰度|ERI

-----------|--------|--------|--------|--------|--------|--------|--------|----------|----------|-----

猪场粪便|A1|1.23×10^8|8.45×10^7|5.67×10^7|3.21×10^6|2.34×10^6|9.87×10^6|1.08×10^9|2.56×10^9|0.38

|A2|1.37×10^8|9.12×10^7|6.34×10^7|4.56×10^6|2.78×10^6|1.56×10^6|1.23×10^9|2.12×10^9|0.42

|A3|1.45×10^8|7.89×10^7|4.23×10^7|5.12×10^6|1.89×10^6|8.76×10^6|1.69×10^9|3.34×10^9|0.35

猪场废水|B1|6.78×10^7|4.56×10^6|2.34×10^7|1.23×10^6|1.56×10^6|5.43×10^6|1.11×10^8|2.78×10^9|0.29

|B2|5.12×10^7|3.21×10^6|1.78×10^7|0.89×10^6|1.23×10^6|4.56×10^6|1.32×10^8|3.56×10^9|0.41

|B3|4.56×10^7|2.34×10^6|1.56×10^7|1.12×10^6|0.56×10^6|3.21×10^6|1.34×10^8|2.34×10^9|0.37

猪场下游地表水|C1|1.56×10^6|8.76×10^5|9.87×10^6|5.43×10^5|3.21×10^5|1.23×10^6|3.78×10^7|4.56×10^9|0.52

|C2|1.12×10^6|6.78×10^5|1.56×10^6|7.89×10^5|2.34×10^5|9.87×10^6|3.56×10^7|5.67×10^9|0.45

|C3|8.76×10^5|4.56×10^5|1.78×10^6|0.56×10^5|1.12×10^6|5.43×10^5|2.34×10^7|3.21×10^9|0.39

鸡场粪便|D1|2.34×10^7|1.56×10^7|1.12×10^7|9.87×10^6|6.78×10^6|4.56×10^6|6.12×10^8|5.43×10^9|0.33

|D2|1.78×10^7|8.76×10^6|9.87×10^6|5.43×10^6|3.21×10^5|1.23×10^6|3.56×10^7|4.56×10^9|0.28

|D3|9.87×10^6|6.78×10^6|1.12×10^7|7.89×10^6|4.56×10^6|2.34×10^5|9.87×10^6|7.78×10^7|6.12×10^9|0.31

鸡场废水|E1|1.23×10^6|8.76×10^5|9.87×10^6|5.43×10^5|3.21×10^5|1.23×10^6|3.78×10^7|4.56×10^9|0.44

|E2|1.78×10^6|6.78×10^5|1.56×10^6|7.89×10^5|4.56×10^5|9.87×10^6|3.56×10^7|5.43×10^9|0.37

|E3|9.87×10^6|5.43×10^5|1.12×10^6|4.56×10^5|1.23×10^6|6.78×10^5|2.34×10^7|3.21×10^9|0.35

鸡场下游地表水|F1|5.43×10^5|3.21×10^5|1.23×10^6|7.89×10^5|4.56×10^5|9.87×10^6|2.34×10^7|6.12×10^9|0.27

|F2|4.56×10^5|2.34×10^5|1.78×10^6|0.56×10^5|1.12×10^6|6.78×10^5|9.87×10^6|3.56×10^7|5.67×10^9|0.29

|F3|1.23×10^6|8.76×10^5|1.56×10^6|7.89×10^5|4.56×10^5|9.87×10^6|3.78×10^7|4.56×10^9|0.32

牛场粪便|G1|1.56×10^7|9.87×10^6|5.43×10^7|3.21×10^6|1.23×10^6|8.76×10^6|9.87×10^8|7.78×10^9|0.41

|G2|1.12×10^7|7.89×10^6|4.56×10^7|2.34×10^6|9.87×10^6|6.78×10^6|3.56×10^7|8.76×10^8|0.38

|G3|8.76×10^6|1.78×10^7|9.87×10^6|5.43×10^7|4.56×10^6|2.34×10^6|1.23×10^6|7.78×10^8|6.12×10^9|0.34

牛场废水|H1|9.87×10^6|5.43×10^5|1.23×10^7|7.89×10^5|4.56×10^5|2.34×10^6|9.87×10^7|7.78×10^9|0.35

|H2|5.43×10^5|3.21×10^5|1.12×10^7|7.89×10^5|4.56×10^5|9.87×10^6|6.78×10^7|2.34×10^8|5.67×10^9|0.33

|H3|4.56×10^5|2.34×10^5|1.78×10^6|0.56×10^5|1.12×10^6|6.78×10^5|9.87×10^6|3.56×10^7|5.67×10^9|0.28

牛场下游地表水|I1|1.23×10^6|8.76×10^5|1.56×10^6|7.89×10^5|4.56×10^5|9.87×10^6|3.78×10^7|7.78×10^9|0.37

|I2|1.12×10^6|7.89×10^5|4.56×10^5|2.34×10^6|9.87×10^6|6.78×10^5|1.23×10^6|3.56×10^7|8.76×10^9|0.39

|I3|8.76×10^5|1.78×10^6|9.87×10^6|5.43×10^5|3.21×10^5|1.23×10^6|6.78×10^5|9.87×10^7|7.78×10^9|0.38

附录B：ARGs定量分析标准曲线绘制数据

表B1.tet(A)标准曲线绘制数据（拷贝数/μL）

样品编号|稀释倍数|浓度（拷贝数/μL）|Cq值

-------