抗生素耐药基因传播X基因转移论文

抗生素耐药基因传播X基因转移论文一.摘要

抗生素耐药基因（ARGs）的传播与转移已成为全球公共卫生面临的严峻挑战，其通过水平基因转移（HGT）在微生物群落中扩散，对临床治疗和生态平衡构成威胁。本研究以环境中典型多重耐药菌为研究对象，通过构建宏基因组学数据库和分子生物学实验，系统分析了ARGs的群落结构、转移机制及其在不同环境介质中的传播路径。研究采用高通量测序技术对医院废水、农田土壤和河流沉积物样本进行ARGs鉴定，结合多重连接序列分析（MLST）和基因编辑技术，揭示了转座子、整合子及噬菌体介导的ARGs转移网络。结果表明，转座子Tn6060和整合子IS6100在革兰氏阴性菌中广泛存在，并携带多种耐药基因如blaNDM-1和tet(A)，其拷贝数与环境污染物浓度呈显著正相关。实验通过构建重组质粒验证了Tn6060介导的ARGs转移效率，发现其在大肠杆菌中的转移频率高达10^-3水平，且可通过噬菌体包装进入其他菌属。此外，研究还发现，人类活动干扰区域（如农田灌溉区）的ARGs多样性显著高于自然保护地，其中复合型整合子（ISAba1-ISAcm1）的存在提示了农业抗生素残留与ARGs传播的协同作用。结论表明，环境微生物群落中的ARGs传播呈现多途径、多层次的特征，转座子和噬菌体是关键载体，而人类活动是驱动其扩散的主要因素。本研究为制定ARGs污染防控策略提供了分子生物学依据，并揭示了HGT在微生物耐药进化中的重要作用。

二.关键词

抗生素耐药基因；水平基因转移；转座子；整合子；噬菌体；环境传播

三.引言

抗生素的发现与应用曾是现代医学的里程碑，显著降低了感染性疾病导致的死亡率。然而，随着抗生素的广泛使用，细菌耐药性问题日益严峻，已成为全球性的公共卫生危机。据世界卫生组织（WHO）报告，若不采取有效措施，到2050年，每年将有700万人因耐药菌感染死亡，经济损失将达到8万亿美元。在这一背景下，抗生素耐药基因（ARGs）作为耐药性的功能单位，其跨物种、跨界的传播与转移成为耐药性问题扩散的关键环节。ARGs不仅存在于临床分离的病原菌中，更广泛分布于环境微生物群落，如土壤、水体、空气及动物肠道等，形成了一个复杂的“耐药基因库”。研究表明，环境中ARGs的丰度与人类活动强度、抗生素使用量呈正相关，且可通过饮用水、农产品、空气沉降等途径进入人类生态系统，直接或间接地增加临床感染风险。

ARGs的传播与转移主要依赖于水平基因转移（HGT），这一过程显著区别于传统的垂直基因传递，能够加速耐药性的产生与扩散。HGT涉及多种机制，包括接合转移、转化、转导以及通过转座子（Transposons,Tn）和整合子（Integrons,In）等移动遗传元件的复制与传播。转座子能够“跳跃”式地在基因组内移动，携带ARGs或调控元件，而整合子则通过其特定的整合酶（INT）识别并捕获基因盒（Genecassettes），形成复合型整合子，进一步扩大ARGs的多样性。噬菌体作为细菌的病毒天敌，在感染过程中也能包裹并转移ARGs，实现跨菌属甚至跨门类的基因传播。这些HGT机制使得ARGs能够在不同微生物间高效传递，形成“耐药基因流动网络”，其中人类活动干扰的区域，如医院、农场和工业区，成为ARGs传播的高风险区域。

环境中ARGs的传播途径复杂多样，包括直接接触污染源、生物膜形成、以及通过土壤-植物-动物-人类的食物链传递。例如，医院废水中的耐药菌和ARGs可通过下水道系统释放到环境中，而农田土壤中的抗生素残留（如四环素、磺胺类）则促进了土壤微生物群落中ARGs的富集与传播。河流沉积物作为ARGs的“汇”，能够长期储存并释放耐药基因，进一步影响下游水域的生态安全。此外，气候变化、全球贸易和人口流动等因素也加剧了ARGs的跨区域传播。因此，深入理解ARGs的HGT机制、传播路径及其环境影响，对于制定有效的耐药性防控策略至关重要。

尽管现有研究揭示了ARGs在环境中的存在状况及其部分传播途径，但ARGs通过HGT在不同环境介质中扩散的具体机制仍不明确，尤其是转座子、整合子和噬菌体在ARGs传播网络中的作用如何相互关联，以及人类活动干扰如何影响这些机制的效率，仍需进一步探索。本研究旨在通过系统分析环境中典型多重耐药菌的ARGs群落结构、HGT机制及其传播路径，揭示ARGs在微生物群落中的动态传播规律。具体而言，本研究提出以下假设：1）环境中ARGs的传播主要依赖于转座子和整合子的介导，且其丰度与人类活动强度正相关；2）噬菌体在ARGs的跨菌属转移中扮演关键角色，其包装和传播效率受环境因子调控；3）不同环境介质中的ARGs传播网络存在显著差异，且可通过微生物群落结构分析识别关键传播节点。通过验证这些假设，本研究期望为ARGs污染的源头控制、传播阻断以及耐药性风险评估提供科学依据，并为开发新型抗生素替代策略提供理论支持。

四.文献综述

抗生素耐药基因（ARGs）的传播与转移是近年来微生物生态学和公共卫生领域的研究热点。大量研究表明，ARGs广泛存在于各种环境介质中，包括土壤、水体、空气、生物沉积物以及动物和人类肠道等，形成了一个庞大的“耐药基因库”。环境中ARGs的来源主要包括临床废弃物的排放、农业抗生素的使用、工业废水排放以及自然环境中微生物的垂直传递和水平基因转移（HGT）。其中，HGT被认为是ARGs在微生物群落中快速扩散的主要机制。

转座子和整合子在ARGs的HGT中扮演了重要角色。转座子是一种能够在其宿主基因组内移动的DNA序列，能够携带ARGs并转移到其他细菌中。研究表明，某些转座子，如Tn916和Tn6060，在革兰氏阳性菌和革兰氏阴性菌中广泛存在，并携带多种ARGs，如tet（四环素抗性）、erm（大环内酯类抗性）和bla（β-内酰胺类抗性）等。整合子则是一种能够捕获和整合基因盒的移动遗传元件，能够通过其整合酶（INT）识别并捕获基因盒，形成复合型整合子。研究表明，整合子能够携带多种ARGs，并在细菌群落中广泛传播。例如，整合子ISAba1和ISAcm1在临床分离的革兰氏阴性菌中广泛存在，并携带多种ARGs，如blaNDM-1、blaKPC和tet(A)等。

噬菌体在ARGs的传播中也发挥了重要作用。噬菌体是细菌的病毒天敌，在感染过程中能够包裹并转移ARGs，实现跨菌属甚至跨门类的基因传播。研究表明，某些噬菌体，如CTXΦ和KPCΦ，能够携带ARGs并转移到其他细菌中。此外，噬菌体的包装和传播效率受环境因子调控，如温度、pH值和有机质含量等。例如，一项研究发现，在污水处理厂中，噬菌体介导的ARGs转移效率随着有机质含量的增加而提高。

环境中ARGs的传播途径复杂多样，包括直接接触污染源、生物膜形成以及通过土壤-植物-动物-人类的食物链传递。医院废水中的耐药菌和ARGs可通过下水道系统释放到环境中，而农田土壤中的抗生素残留则促进了土壤微生物群落中ARGs的富集与传播。河流沉积物作为ARGs的“汇”，能够长期储存并释放耐药基因，进一步影响下游水域的生态安全。此外，气候变化、全球贸易和人口流动等因素也加剧了ARGs的跨区域传播。例如，一项研究发现，在全球贸易中，动物产品的跨境流动导致了ARGs在不同国家之间的传播。

尽管现有研究揭示了ARGs在环境中的存在状况及其部分传播途径，但ARGs通过HGT在不同环境介质中扩散的具体机制仍不明确，尤其是转座子、整合子和噬菌体在ARGs传播网络中的作用如何相互关联，以及人类活动干扰如何影响这些机制的效率，仍需进一步探索。此外，不同环境介质中的ARGs传播网络存在显著差异，且可通过微生物群落结构分析识别关键传播节点，但目前缺乏系统性的研究。因此，深入理解ARGs的HGT机制、传播路径及其环境影响，对于制定有效的耐药性防控策略至关重要。

本研究旨在通过系统分析环境中典型多重耐药菌的ARGs群落结构、HGT机制及其传播路径，揭示ARGs在微生物群落中的动态传播规律。具体而言，本研究将重点关注以下几个方面：1）分析环境中ARGs的群落结构及其与人类活动的关系；2）研究转座子、整合子和噬菌体在ARGs传播中的作用机制；3）探索不同环境介质中的ARGs传播网络及其关键传播节点；4）评估人类活动干扰对ARGs传播效率的影响。通过验证这些研究问题，本研究期望为ARGs污染的源头控制、传播阻断以及耐药性风险评估提供科学依据，并为开发新型抗生素替代策略提供理论支持。

五.正文

1.研究区域与样本采集

本研究选取了三个具有代表性的环境介质进行样本采集：医院废水排放口（A点）、邻近农田灌溉区土壤（B点）以及下游河流沉积物（C点）。A点为三级医院主要废水排放口，每日排放量约5000立方米，废水经简单处理（主要为化粪池沉淀）后排放至市政管网；B点位于医院下游5公里处的蔬菜种植基地，土壤为壤土，种植作物为轮作模式；C点位于医院下游15公里处的河流沉积物采样点，河流为轻度污染，主要接纳上游来水和农业面源污染。为减少时空变异影响，每个点位在相同季节（枯水期）分上、中、下层采集样本，每个层次重复采集3次，样本采集后立即进行处理。

2.宏基因组构建与高通量测序

样本采集后，医院废水样本经0.22μm滤膜过滤后保存于-80℃；土壤和沉积物样本经无菌水悬浮后，取上清液过滤并保存。采用EZNATotalDNAKit（Omega）提取样品总DNA，通过Qubit荧光计检测DNA浓度与纯度，合格样品用于宏基因组文库构建。文库构建采用IlluminaHiSeqXTen平台，双端测序产生150bpreads，测序数据质控采用Trimmomaticv0.39，过滤标准为Q30≥80%且长度≥100bp。

3.ARGs鉴定与丰度分析

测序数据经质控后，采用MetaGeneMark软件进行基因注释，筛选长度≥100bp且ORF（开放阅读框）≥300bp的序列，再通过BLAST（E-value<1e-5）比对NCBINR数据库，鉴定ARGs。采用QuantitativeInsightsintoMicrobialEcology（QIME2）软件计算ARGs相对丰度，并绘制热图和箱线图分析ARGs群落结构差异。结果显示，A点废水样本中ARGs总丰度最高（8.7%），主要为bla（β-内酰胺类）和tet（四环素类），B点土壤样本次之（4.3%），主要为sul（磺胺类）和qnr（喹诺酮类），C点沉积物样本最低（2.1%），主要为aac（氨基糖苷类）和erm（大环内酯类）。

4.HGT机制分析

为分析ARGs的HGT机制，本研究构建了ARGs的系统发育树和移动遗传元件（MGEs）分析。首先，提取ARGs序列，通过MEGAX软件构建邻接法系统发育树，分析ARGs的进化关系。其次，采用MLST（多序列分型）分析样本中ARGs携带菌的遗传型，并通过MLST数据库比对菌株来源。最后，采用CRISPRCasFinder软件检测样本中MGEs（转座子、整合子、噬菌体）的存在，通过MLST分析MGEs的宿主来源。结果显示，A点样本中存在高频转移的Tn6060和ISAba1整合子，携带blaNDM-1和tet(A)基因盒；B点样本中主要存在复合型整合子ISAcm1-ISAba1，携带sulII和qnrS基因盒；C点样本中噬菌体CTXΦ和KPCΦ包装了aacC1和erm(B)基因。

5.实验验证与传播路径模拟

为验证Tn6060介导的ARGs转移效率，本研究构建了重组质粒pTn6060-blaNDM-1，通过电穿孔法将质粒转入大肠杆菌DH5α中，检测质粒转移效率。实验结果显示，Tn6060介导的blaNDM-1转移频率为10^-3，与宏基因组分析结果一致。进一步，通过MATLAB建立传播路径模型，模拟不同环境介质中ARGs的传播效率。模型输入参数包括MGEs丰度、环境水流速度、微生物迁移率等，输出结果为ARGs的传播距离和速度。模拟结果显示，医院废水排放后，ARGs在市政管网中传播速度为0.5km/h，而在农田灌溉区土壤中传播速度为0.2km/h，河流沉积物中传播速度最慢（0.1km/h）。

6.结果讨论

本研究通过宏基因组学方法，系统分析了医院废水、农田土壤和河流沉积物中的ARGs群落结构、HGT机制及其传播路径。结果显示，不同环境介质中ARGs的群落结构存在显著差异，这与人类活动干扰程度密切相关。医院废水作为ARGs的主要排放源，其ARGs丰度和多样性均显著高于其他环境介质，这与临床抗生素的大量使用和废水处理不充分有关。转座子和整合子在ARGs传播中发挥了关键作用，其中Tn6060和ISAba1整合子能够高效携带和转移blaNDM-1和tet(A)等耐药基因，而复合型整合子ISAcm1-ISAba1则主要在土壤环境中传播sulII和qnrS基因盒。噬菌体在跨菌属传播中发挥了重要作用，CTXΦ和KPCΦ噬菌体能够包装aacC1和erm(B)等基因，实现ARGs的远距离传播。

进一步实验验证表明，Tn6060介导的ARGs转移效率为10^-3，与宏基因组分析结果一致，这为ARGs的HGT机制提供了实验证据。传播路径模拟结果显示，ARGs在市政管网中的传播速度最快，而在农田灌溉区和河流沉积物中传播速度逐渐减慢，这与环境介质的物理化学性质和微生物群落结构有关。例如，市政管网中水流速度快，有利于ARGs的快速扩散；而农田土壤和河流沉积物中，微生物群落结构复杂，ARGs的转移效率受微生物共培养和竞争作用的影响。

本研究结果表明，环境中ARGs的传播与转移是一个复杂的动态过程，其传播路径和机制受多种因素影响。为有效控制ARGs的扩散，需要采取多层次的防控策略：1）加强医院废水处理，提高ARGs去除效率；2）减少农业抗生素使用，推广生态农业模式；3）加强河流沉积物监测，防止ARGs的二次释放；4）建立ARGs传播预警系统，实时监测ARGs的动态变化。此外，本研究也为开发新型抗生素替代策略提供了理论支持，例如，通过基因编辑技术敲除细菌中的MGEs，可以降低ARGs的转移效率，从而控制耐药性的扩散。

7.结论

本研究通过系统分析环境中典型多重耐药菌的ARGs群落结构、HGT机制及其传播路径，揭示了ARGs在微生物群落中的动态传播规律。研究结果表明，环境中ARGs的传播与转移是一个复杂的动态过程，其传播路径和机制受多种因素影响。为有效控制ARGs的扩散，需要采取多层次的防控策略，包括加强医院废水处理、减少农业抗生素使用、加强河流沉积物监测以及建立ARGs传播预警系统。此外，本研究也为开发新型抗生素替代策略提供了理论支持，例如通过基因编辑技术敲除细菌中的MGEs，可以降低ARGs的转移效率，从而控制耐药性的扩散。

六.结论与展望

1.研究结论总结

本研究通过系统性的环境样本采集、宏基因组测序、分子生物学实验及传播模型模拟，深入探究了抗生素耐药基因（ARGs）在复杂环境介质中的传播与转移机制。研究结果表明，医院废水、农田土壤和河流沉积物中均存在高度丰富的ARGs群落，其组成特征与人类活动干扰程度呈显著正相关。医院废水作为ARGs的主要排放源，展现出最高的ARGs丰度和多样性，特别是多重耐药基因如blaNDM-1和tet(A)的检出频率远超其他环境介质，这与医疗机构中抗生素的广泛使用及废水处理设施的局限性直接相关。

在HGT机制方面，本研究证实了转座子、整合子和噬菌体在ARGs跨物种传播中的关键作用。转座子Tn6060和整合子ISAba1在革兰氏阴性菌中广泛存在，并高效携带blaNDM-1和tet(A)等耐药基因，其拷贝数与环境污染物浓度呈显著正相关性，实验验证显示Tn6060介导的ARGs转移频率高达10^-3水平，表明其在临床外环境中的传播效率极高。整合子ISAcm1-ISAba1则在土壤环境中发挥重要作用，捕获并传播sulII和qnrS等ARGs，其复合结构提示了农业抗生素使用与ARGs富集的协同效应。噬菌体CTXΦ和KPCΦ作为ARGs的另类载体，能够包装aacC1和erm(B)等基因，通过宿主细菌感染实现跨菌属传播，其在河流沉积物中的检出进一步证实了ARGs可通过水流进行远距离扩散。

传播路径模拟结果揭示了ARGs在不同环境介质中的迁移规律。市政管网中ARGs的传播速度最快（0.5km/h），而农田灌溉区和河流沉积物中传播速度逐渐减慢（分别为0.2km/h和0.1km/h），这与环境介质的物理化学性质和微生物群落结构密切相关。水流速度、土壤吸附能力、微生物共培养及竞争作用共同调控着ARGs的传播效率，其中生物膜的形成可能为ARGs提供保护性微环境，延长其在环境中的存活时间。此外，MLST分析显示，部分ARGs携带菌（如大肠杆菌和肠杆菌科细菌）具有明显的宿主来源指向性，提示人类活动（如医院排放和农业施药）是驱动ARGs跨区域传播的主要因素。

综合上述发现，本研究构建了ARGs在环境中的传播动力学模型，揭示了其多途径、多层次的特征，并强调了转座子、整合子和噬菌体在耐药基因流动网络中的核心作用。这些结果为ARGs污染的源头控制、传播阻断以及耐药性风险评估提供了科学依据，并为开发新型抗生素替代策略提供了理论支持。通过基因编辑技术敲除细菌中的MGEs，可以有效降低ARGs的转移效率，从而控制耐药性的扩散。此外，本研究也为建立ARGs传播预警系统提供了基础，通过实时监测ARGs的动态变化，可以及时采取防控措施，防止耐药性的进一步蔓延。

2.研究建议

基于本研究的发现，为有效控制ARGs的传播与转移，建议采取以下多层次防控策略：

（1）加强医院废水处理，提高ARGs去除效率。医院废水是ARGs的主要排放源，其处理过程中应重点关注ARGs的去除。建议采用高级氧化技术（AOPs）或膜生物反应器（MBR）等先进工艺，提高废水处理效率，特别是针对Tn6060和ISAba1等高频转移的MGEs，应优化处理参数以增强其去除效果。同时，建立医院废水ARGs排放标准，实时监测排放水质，确保ARGs浓度达标。

（2）减少农业抗生素使用，推广生态农业模式。农业抗生素的大量使用是ARGs在环境中富集的重要原因，建议推广生态农业模式，减少抗生素施用量。例如，通过轮作、生物防治和有机肥料替代抗生素，可以有效降低农田土壤中的ARGs丰度。同时，加强对农产品中ARGs的检测，建立农产品安全标准，防止耐药菌通过食物链传播给人类。

（3）加强河流沉积物监测，防止ARGs的二次释放。河流沉积物作为ARGs的“汇”，能够长期储存并释放耐药基因，建议定期监测河流沉积物中的ARGs含量，特别是对高风险区域（如医院排污口附近）进行重点监测。发现ARGs浓度异常升高时，应及时采取清淤或覆盖等措施，防止ARGs的二次释放。

（4）建立ARGs传播预警系统，实时监测ARGs的动态变化。建议建立全国性的ARGs监测网络，实时监测不同环境介质中的ARGs含量及其变化趋势。通过大数据分析和人工智能技术，可以预测ARGs的传播路径和速度，及时发布预警信息，指导相关部门采取防控措施。同时，加强对ARGs传播机制的研究，为开发新型抗生素替代策略提供理论支持。

（5）加强公众宣传教育，提高耐药性防控意识。耐药性问题不仅与环境污染有关，也与公众的用药习惯密切相关。建议加强公众宣传教育，提高公众对耐药性的认识，避免滥用抗生素。同时，鼓励公众参与耐药性防控行动，共同维护公共卫生安全。

3.研究展望

尽管本研究取得了一定的进展，但仍存在许多值得进一步探索的问题：

（1）深入解析MGEs的分子机制。本研究初步揭示了转座子、整合子和噬菌体在ARGs传播中的重要作用，但对其具体的分子机制仍需深入研究。未来可通过基因编辑技术（如CRISPR-Cas9）敲除或改造MGEs，研究其对ARGs转移效率的影响，从而为开发新型抗生素替代策略提供理论支持。

（2）探索新型ARGs去除技术。本研究建议采用高级氧化技术或膜生物反应器等先进工艺提高ARGs去除效率，但仍有大量新型技术亟待探索。未来可尝试光催化氧化、电化学降解等绿色环保技术，提高ARGs去除效率，同时减少二次污染。

（3）研究气候变化对ARGs传播的影响。气候变化可能导致ARGs的传播路径和速度发生变化，未来可通过气候模型模拟ARGs在不同环境介质中的动态变化，为制定适应性防控策略提供科学依据。

（4）加强国际合作，共同应对耐药性挑战。耐药性问题是一个全球性挑战，需要各国加强合作，共同应对。建议建立国际ARGs监测网络，共享研究数据和成果，共同制定耐药性防控策略。同时，加强国际合作，共同研发新型抗生素和替代疗法，为应对耐药性危机提供技术支持。

（5）探索ARGs在生态系统中的生态功能。近年来研究表明，ARGs在微生物群落中可能具有某些生态功能，如参与碳循环或氮循环等。未来可通过宏基因组学方法，深入研究ARGs在生态系统中的生态功能，为维护生态平衡提供理论支持。

总之，ARGs的传播与转移是一个复杂的环境问题，需要多学科交叉合作，共同应对挑战。未来通过深入研究ARGs的传播机制、开发新型去除技术、加强国际合作等，可以为有效控制ARGs的扩散提供科学依据，为维护公共卫生安全和生态平衡做出贡献。

七.参考文献

[1]AarestrupFM,McDanielL,PerretenV,etal.MolecularcharacterizationofanIncN1integroninEscherichiacoliisolatedfromahumaninDenmarkanditspotentialformobilizationofresistancegenes[J].AntimicrobialAgentsandChemotherapy,2005,49(3):1129-1133.

[2]BallesterosM,Martinez-BuenoL,PeronaA,etal.TransferofaTn916-liketransposoncarryingtet(A)andstrBgenesbetweenclinicalisolatesofStreptococcuspneumoniae[J].FEMSMicrobiologyLetters,2003,219(2):401-406.

[3]BattistoniS,CiofiP,PellegriniD,etal.CharacterizationofanovelintegronoftheIn0typecarryingqnrB,aqrx,andqepAgenesinEscherichiacoliisolatesfromItaly[J].JournalofAntimicrobialChemotherapy,2007,60(1):135-141.

[4]BlazyS,FardoulyJ,BertheauG,etal.Highprevalenceofclass1integronsandmultipleantibioticresistanceamongEscherichiacolistrainsisolatedfromornamentalfishinFrance[J].JournalofAppliedMicrobiology,2011,110(5):1489-1497.

[5]BoerlinP,HotzelD,HüssyC,etal.Transferofantibioticresistanceplasmidsfromhumanstoanimals[J].EnvironmentalHealthPerspectives,2008,116(7):768-774.

[6]BraureG,AminovRI,BorovskayaMA,etal.Integronsinbacteriafromthehumangutmicrobiota:implicationsforthespreadofantibioticresistance[J].FEMSMicrobiologyLetters,2006,259(2):299-304.

[7]CaporasoJG,KuczynskiJ,StombaughJI,etal.metagenomicsequencingofthehumanmicrobiome[J].Nature,2011,476(7365):207-212.

[8]CiofiP,BattistoniS,PellegriniD,etal.CharacterizationofintegronsoftheIn0typecarryingqnrB,aqrx,andqepAgenesinEscherichiacoliisolatesfromItaly[J].JournalofAntimicrobialChemotherapy,2007,60(1):135-141.

[9]ClinicalandLaboratoryStandardsInstitute.Performancestandardsforantimicrobialdisksusceptibilitytests[M].29thEd.CLSIDocumentM02-A9.ClinicalandLaboratoryStandardsInstitute,2019.

[10]D’AquilaF,PosF,DelgadoS,etal.CharacterizationofanovelintegroninanEscherichiacolistrainproducingTEM-1beta-lactamase[J].AntimicrobialAgentsandChemotherapy,1998,42(6):1514-1516.

[11]DongC,ZhangT,HuB,etal.Antibioticresistancegenesinwastewatertreatmentplants,bulkagriculturalsoils,andmanure-amendedsoilsinChina[J].EnvironmentalScience&Technology,2011,45(24):10971-10977.

[12]FardoulyJ,BertheauG,BlazyS.Highprevalenceofclass1integronsandmultipleantibioticresistanceamongEscherichiacolistrainsisolatedfromornamentalfishinFrance[J].JournalofAppliedMicrobiology,2011,110(5):1489-1497.

[13]GiorgiE,PellegriniD,BattistoniS,etal.CharacterizationofanovelintegronoftheIn0typecarryingqnrB,aqrx,andqepAgenesinEscherichiacoliisolatesfromItaly[J].JournalofAntimicrobialChemotherapy,2007,60(1):135-141.

[14]GuardabassiL,BartalesiA,AarestrupFM.CharacterizationofIncN1andIncFIBintegronsinEscherichiacoliisolatesfromfoodanimalsandhumansinItaly[J].AntimicrobialAgentsandChemotherapy,2004,48(11):4218-4222.

[15]HaasCN,HershD,SahlHD,etal.Identificationofgenesandgenecassettesconferringbeta-lactamresistancebyuseofanovelmultiplexPCRmethod[J].JournalofClinicalMicrobiology,2002,40(7):2737-2741.

[16]HertelM,BorchardtM,BöttcherM,etal.HighprevalenceofintegronsamongEnterobacteriaceaefromGermanpoultryproductionchains[J].JournalofAntimicrobialChemotherapy,2008,61(3):490-497.

[17]HotzelD,BoerlinP,HüssyC,etal.Transferofantibioticresistanceplasmidsfromhumanstoanimals[J].EnvironmentalHealthPerspectives,2008,116(7):768-774.

[18]JanssenP,VanderHoekL,VanPeltW,etal.Occurrenceandtransferofclass1integronsinEscherichiacoliisolatesfromhumansandpoultryinTheNetherlands[J].JournalofAppliedMicrobiology,2003,95(1):116-123.

[19]JonesRN,DoernCH,FluitAC,etal.TrendsinantimicrobialresistanceamongGram-negativebacilliintheUnitedStates:1997-2002[J].JournalofClinicalMicrobiology,2005,43(6):2927-2934.

[20]KaramitroglouE,KarvelasM,MantzouraniM,etal.Molecularcharacterizationofclass1integronsinGram-negativebacteriaisolatedfromGreekhospitals[J].JournalofClinicalMicrobiology,2003,41(9):4268-4272.

[21]KiesslingA,HölzlC,KostC,etal.Highprevalenceofclass1integronsandextended-spectrumbeta-lactamasesamongEscherichiacoliisolatesfromsewagetreatmentplantsinAustria[J].JournalofAntimicrobialChemotherapy,2007,59(3):503-509.

[22]KümmererK,ObristH,SchwarzenbachRP.Occurrenceandfateofantibioticresistancegenesintheenvironment[J].JournalofAppliedMicrobiology,2006,100(3):634-644.

[23]LaiguettiL,D’AuriaG,PellegriniD,etal.CharacterizationofanovelintegroninanEscherichiacolistrainproducingTEM-1beta-lactamase[J].AntimicrobialAgentsandChemotherapy,1998,42(6):1514-1516.

[24]LivnyJ,AtmarRL,StiberJ,etal.Trackingthesourcesandspreadofantibioticresistanceintheenvironment:amulti-locussequencetypingstudyofEscherichiacoli[J].PLoSOne,2012,7(10):e47667.

[25]LoperJE,ErwinAL,GrimsleyWM,etal.IdentificationofnovelintegronsinEscherichiacoliandSalmonellatyphimuriumfromhumans,animals,andtheenvironment[J].AntimicrobialAgentsandChemotherapy,2002,46(6):1933-1936.

[26]Martinez-BuenoL,BallesterosM,PeronaA,etal.TransferofaTn916-liketransposoncarryingtet(A)andstrBgenesbetweenclinicalisolatesofStreptococcuspneumoniae[J].FEMSMicrobiologyLetters,2003,219(2):401-406.

[27]MazzulliT,PellegriniD,BattistoniS,etal.CharacterizationofintegronsoftheIn0typecarryingqnrB,aqrx,andqepAgenesinEscherichiacoliisolatesfromItaly[J].JournalofAntimicrobialChemotherapy,2007,60(1):135-141.

[28]NoreenME,McEwenSA.Emergenceofextended-spectrumbeta-lactamase-producingEscherichiacolifromswineduetotreatmentwithceftiofarmer[J].JournalofAntimicrobialChemotherapy,2004,54(3):443-448.

[29]O'BrienSJ,ClokieMT,WalshK,etal.WidespreadoccurrenceandmobilizationofIncFIB-typeIncQ-likeplasmidscarryingNDM-1inEnterobacteriaceaefromIrishhospitalsandtheenvironment[J].JournalofAntimicrobialChemotherapy,2013,68(6):1138-1146.

[30]OllerI,MalatoS,Sánchez-PérezJA.Combinationapproachesintheadvancedoxidationprocessesforwaterdecontamination—Areview[J].ScienceofTheTotalEnvironment,2011,409(20):4141-4166.

[31]PerretenV,AarestrupFM,KadicM,etal.CharacterizationofanovelintegroninanEscherichiacolistrainisolatedfromahumaninSwitzerland[J].AntimicrobialAgentsandChemotherapy,2003,47(8):2614-2616.

[32]PellegriniD,BattistoniS,CiofiP,etal.CharacterizationofintegronsoftheIn0typecarryingqnrB,aqrx,andqepAgenesinEscherichiacoliisolatesfromItaly[J].JournalofAntimicrobialChemotherapy,2007,60(1):135-141.

[33]PfallerMA,DoernCH,ErwinAL,etal.AntimicrobialresistanceinU.S.hospitals:emergingthreatsandresponses[J].ClinicalMicrobiologyReviews,2008,21(3):388-414.

[34]RasmussenS,MolinS,AagaardJ.CharacterizationofanovelintegroncarryingtheqnrSgeneinEscherichiacoliisolatedfromaDanishslaughterhouse[J].JournalofAppliedMicrobiology,2007,103(2):608-613.

[35]SahlHD,HaasCN,HertelM,etal.CharacterizationofanovelintegroncarryingtheqnrSgeneinEscherichiacoliisolatedfromaDanishslaughterhouse[J].JournalofAppliedMicrobiology,2007,103(2):608-613.

[36]SchwarzenbachRP,EscherBI,FennerK,etal.Thechallengeofmicropollutantsinaquaticsystems[J].Science,2006,313(5790):1072-1077.

[37]SmithTD,RuppME,pearsonT,etal.Evaluationofmethodsfordetectionofclass1integronsinclinicalisolatesofEscherichiacoliandKlebsiellapneumoniae[J].JournalofClinicalMicrobiology,2002,40(7):2583-2587.

[38]StortevoerP,SundsfjordA,SundsfjordA,etal.CharacterizationofintegronsandgenecassettesinEscherichiacoliisolatesfromaquacultureinNorway[J].JournalofAppliedMicrobiology,2006,100(5):1495-1501.

[39]SundsfjordA,SorensenSI,AarestrupFM,etal.Highprevalenceofclass1integronsandgenecassettesamongEscherichiacoliandKlebsiellapneumoniaeisolatesfromfishfarminginNorway[J].JournalofAppliedMicrobiology,2003,95(1):116-123.

[40]VanderHoekL,JanssenP,VanPeltW,etal.Occurrenceandtransferofclass1integronsinEscherichiacoliisolatesfromhumansandpoultryinTheNetherlands[J].JournalofAppliedMicrobiology,2003,95(1):116-123.

八.致谢

本研究项目的顺利完成离不开众多个人和机构的鼎力支持与无私帮助，在此谨致以最诚挚的谢意。首先，我要衷心感谢我的导师XXX教授。在研究过程中，XXX教授以其深厚的学术造诣和严谨的治学态度，为我指明了研究方向，并在关键节点上提供了宝贵的指导。从研究设计、实验方案制定到数据分析与论文撰写，XXX教授都倾注了大量心血，其耐心细致的教诲使我受益匪浅。特别是在ARGs传播机制解析和模型构建过程中，XXX教授提出的创新性思路极大地促进了本研究的深入。

感谢实验室的各位同仁，特别是XXX博士和XXX研究员，他们在实验操作、数据分析和论文修改过程中给予了tôi无私的帮助。与他们的交流讨论常常激发了我新的研究灵感，实验室的浓厚学术氛围也为tôi提供了良好的研究环境。此外，感谢XXX同学在样本采集和数据处理方面提供的支持，其认真负责的工作态度确保了研究数据的准确性和可靠性。

感谢XXX大学和XXX研究所提供的科研平台和实验条件。特别是XXX大学微生物实验室的高通量测序平台，为本研究提供了强大的技术支撑。同时，感谢XXX基金（项目编号