抗生素耐药基因传播X移动性遗传元件论文

抗生素耐药基因传播X移动性遗传元件论文一.摘要

抗生素耐药性问题已成为全球公共卫生领域面临的主要挑战之一，其中耐药基因的传播机制对于理解和控制其扩散至关重要。本研究聚焦于抗生素耐药基因（ARGs）通过移动性遗传元件（MGEs）的传播现象，通过整合宏基因组学和生物信息学分析方法，系统探究了ARGs与MGEs在环境样品中的共现模式及其传播途径。研究选取了典型的高污染农业土壤和医院废水作为样本来源，利用高通量测序技术获取环境微生物组的基因组数据，并通过特定的生物信息学工具对ARGs和MGEs进行鉴定与分类。研究发现，多种常见的ARGs，如ntermicin-resistantgenes和vancomycin-resistantgenes，与转座子、整合子等MGEs存在显著的正相关关系，表明MGEs在ARGs的传播过程中起到了关键的媒介作用。进一步的功能预测分析揭示了这些MGEs携带的ARGs可能赋予微生物多重耐药性，从而增强了其在复杂环境中的生存竞争力。此外，研究还发现了ARGs与MGEs的共现热点区域，这些区域通常与人类活动干扰频繁的环境（如农田灌溉区、医院排污口）高度重合。这些发现不仅揭示了ARGs通过MGEs传播的具体机制，也为制定有效的干预策略提供了科学依据，例如通过控制MGEs的传播来遏制ARGs的扩散，从而缓解抗生素耐药性对公共卫生构成的威胁。本研究强调了MGEs在ARGs传播中的核心作用，为抗生素耐药性治理提供了新的视角和思路。

二.关键词

抗生素耐药基因；移动性遗传元件；宏基因组学；转座子；整合子；环境传播

三.引言

抗生素的发现与应用无疑是现代医学史上最重大的突破之一，极大地提高了人类对抗感染性疾病的斗争能力。然而，随着抗生素的广泛和反复使用，细菌耐药性问题日益严峻，已成为全球性的公共卫生危机。据世界卫生（WHO）报告，抗生素耐药性导致的感染每年可能造成数百万人死亡，且形势正在不断恶化。在这一背景下，理解抗生素耐药基因（ARGs）的传播机制，特别是其如何在不同微生物种群乃至不同环境中扩散，对于制定有效的控制策略至关重要。

近年来，大量的研究证据表明，ARGs的传播并不仅仅是通过传统的水平基因转移（HGT）途径，移动性遗传元件（MGEs）在其中扮演了关键角色。MGEs是一类能够在不同基因组之间移动的DNA片段，包括转座子（Transposons,Tn）、整合子（Integrons,In）、质粒（Plasmids）和噬菌体（Phages）等。它们不仅能够携带基因本身，还能介导基因在不同宿主间的转移，从而极大地促进了遗传物质的交换和演化。特别是在微生物群落中，MGEs如同“基因剪刀”和“基因粘贴”工具，能够捕获、移动和重组ARGs，使其能够在不同的细菌物种间传播，形成复杂的耐药性基因库。

环境水体和土壤作为微生物基因库的重要储存库，是ARGs和MGEs传播的关键媒介。农业活动、医院废水排放、动物粪便等人类活动产生的废弃物中，往往含有高浓度的抗生素和微生物，这些环境介质为ARGs和MGEs的富集与传播提供了有利条件。例如，在集约化农业系统中，大量使用抗生素不仅直接诱导细菌产生耐药性，还可能通过选择性压力促进携带ARGs的MGEs在微生物群落中的传播。医院废水处理厂（WWTPs）则被认为是ARGs和MGEs传播的重要节点，由于污水中含有大量的抗生素、抗菌药物以及多样化的微生物群落，这些地方可能成为耐药基因的“熔炉”，通过初级处理、二级处理和污泥处置等环节，将耐药基因释放到环境中，进而扩散到更广泛的生态系统。

尽管已有研究揭示了某些特定ARGs或MGEs的存在及其潜在风险，但系统性地研究ARGs与MGEs在环境样品中的共现模式、传播特征及其环境影响因素仍面临诸多挑战。首先，环境中微生物组的复杂性使得精确鉴定和定量ARGs和MGEs成为一项难题。其次，MGEs与ARGs的关联性研究通常需要整合多组学数据（如基因组学、转录组学和蛋白质组学），并结合环境参数进行综合分析，这对研究技术和数据处理能力提出了较高要求。此外，目前对于ARGs通过MGEs传播的具体路径和动力学过程尚缺乏深入理解，特别是在不同环境介质中的传播机制可能存在显著差异。

鉴于上述背景，本研究旨在系统探究抗生素耐药基因（ARGs）通过移动性遗传元件（MGEs）在环境样品中的传播现象。具体而言，本研究将重点关注以下几个方面：（1）利用高通量宏基因组测序技术，在农业土壤和医院废水等典型环境样品中鉴定和定量ARGs与MGEs；（2）分析ARGs与MGEs的共现关系，揭示MGEs在ARGs传播中的介导作用；（3）结合环境参数（如抗生素残留浓度、微生物群落结构等），探讨影响ARGs与MGEs共现模式的关键因素；（4）评估ARGs通过MGEs传播的环境风险，为制定有效的抗生素耐药性控制策略提供科学依据。本研究的意义在于，通过揭示ARGs与MGEs在环境中的传播机制，不仅能够深化对细菌耐药性演化规律的理解，还能为开发新的干预措施提供理论支持，例如通过靶向抑制MGEs的活性来阻断ARGs的传播，从而缓解抗生素耐药性对公共卫生构成的威胁。同时，本研究结果也将为环境微生物基因组的深入研究提供新的思路和方法。

四.文献综述

抗生素耐药性（AntibioticResistance,AR）已成为全球性的公共卫生危机，威胁着现代医学的基石。在众多导致AR的因素中，抗生素耐药基因（AntibioticResistanceGenes,ARGs）的传播机制是理解其扩散规律和控制策略制定的关键。近年来，移动性遗传元件（MobileGeneticElements,MGEs）在ARGs传播中的作用日益受到关注。MGEs是一类能够自我复制并在不同基因组间转移的DNA片段，包括转座子（Transposons,Tn）、整合子（Integrons,In）、质粒（Plasmids）和噬菌体（Phages）等。它们不仅能够携带基因本身，还能介导基因在不同宿主间的转移，从而极大地促进了遗传物质的交换和演化。

早期关于MGEs的研究主要集中在病原菌的基因组中，发现MGEs是ARGs传播的重要载体。例如，转座子Tn21被证明能够携带多个ARGs，并在不同的革兰氏阴性菌中转移，导致多重耐药性的出现。整合子被认为是ARGs捕获和传播的关键工具，它们能够通过重组作用将ARGs整合到自己的结构中，并通过转座酶的作用将ARGs转移到其他基因组中。质粒则是ARGs传播的另一重要媒介，许多临床分离的多重耐药菌株携带的质粒上就整合了多种ARGs，使其能够在不同的临床和环境样品中传播。

在环境微生物组中，MGEs和ARGs的共现现象也得到了广泛报道。研究表明，在农业土壤、医院废水、污水处理厂（WWTPs）和海洋等环境中，MGEs和ARGs的存在率往往高于在纯净环境中。例如，一项在农业土壤中的研究发现，携带ARGs的MGEs（主要是整合子和转座子）的存在率高达80%，且与环境中抗生素的使用密切相关。在WWTPs中，由于高浓度的抗生素和多样化的微生物群落，MGEs和ARGs的富集和传播更为显著。研究表明，WWTPs不仅是ARGs的“汇”，也是其向周围环境释放的“源”。通过初级处理、二级处理和污泥处置等环节，WWTPs可能将耐药基因释放到环境中，进而扩散到更广泛的生态系统。

尽管已有研究揭示了MGEs在ARGs传播中的重要作用，但关于ARGs与MGEs在环境样品中共现模式及其传播机制的研究仍存在诸多空白和争议。首先，关于ARGs与MGEs的共现关系，目前的研究结果并不完全一致。一些研究发现ARGs与MGEs存在显著的正相关关系，表明MGEs在ARGs的传播中起到了关键的媒介作用；而另一些研究则发现ARGs与MGEs的共现关系并不显著，甚至存在负相关关系。这种不一致性可能是由于研究区域、环境介质和微生物群落结构的差异导致的。

其次，关于MGEs在ARGs传播中的具体作用机制，目前的研究仍存在一些争议。例如，转座子和整合子在ARGs传播中的作用机制尚不完全清楚。转座子主要通过“复制-粘贴”机制将ARGs转移到其他基因组中，而整合子则通过“捕获-转移”机制将ARGs整合到自己的结构中，并通过转座酶的作用将ARGs转移到其他基因组中。然而，这些机制在环境样品中的具体表现和效率仍需要进一步研究。

此外，关于ARGs通过MGEs传播的环境风险，目前的研究结果也存在一些争议。一些研究表明，MGEs和ARGs的共现热点区域（如农田灌溉区、医院排污口）可能成为耐药基因的“熔炉”，通过人类活动的影响，将耐药基因释放到环境中，进而扩散到更广泛的生态系统。然而，另一些研究则认为，MGEs和ARGs的传播可能受到环境因素的影响，如抗生素残留浓度、微生物群落结构等，这些因素可能抑制或促进耐药基因的传播。

综上所述，尽管已有研究揭示了MGEs在ARGs传播中的重要作用，但关于ARGs与MGEs在环境样品中共现模式及其传播机制的研究仍存在诸多空白和争议。未来的研究需要进一步整合多组学数据（如基因组学、转录组学和蛋白质组学），并结合环境参数进行综合分析，以揭示ARGs通过MGEs传播的具体路径和动力学过程。同时，还需要加强对MGEs在ARGs传播中的具体作用机制和环境风险的研究，为制定有效的抗生素耐药性控制策略提供科学依据。

五.正文

1.研究设计与环境样品采集

本研究旨在系统探究抗生素耐药基因（ARGs）通过移动性遗传元件（MGEs）在特定环境中的传播现象。研究区域选取了两个具有代表性且存在显著人类活动干扰的环境：一个集约化农业示范区（包括蔬菜种植区和使用相同抗生素类别的养殖场）及其周边农田灌溉渠，另一个综合性医院及其配套的污水处理厂（WWTP）。选择这两个区域是因为它们分别代表了农业环境和医疗环境，这两个环境都是ARGs和MGEs潜在的高丰度区域。

样品采集遵循严格的无污染操作规程。在农业示范区，采集了表层（0-10cm）土壤样品，同时采集了灌溉渠的水样和沉积物样品。每个采样点设置3个重复，样品采集于不同季节（春季、夏季、秋季）以捕捉环境变化的动态。在医院环境，采集了医院污水的进水口和出水口样品，以及医院废水处理厂不同处理单元（初沉池、生化池、二沉池、消毒池）的样品。同时，还采集了处理厂产生的污泥样品。所有样品采集后迅速冷藏（4°C），并在4小时内运至实验室进行前处理。所有样品的采集和保存过程均遵循标准操作程序，以避免ARGs和MGEs的二次污染或降解。

2.宏基因组DNA提取与质量控制

根据样品类型（土壤、水、沉积物、污泥），采用不同的宏基因组DNA提取试剂盒进行DNA提取。对于土壤和沉积物，采用基于试剂盒的方法，通常包括细胞裂解、DNA纯化和过滤等步骤。对于水样，采用过滤浓缩法富集微生物，然后对滤膜进行DNA提取。污泥样品由于有机质含量高，采用更高效的试剂盒进行DNA提取。提取后的DNA样品通过0.8%琼脂糖凝胶电泳检测DNA条带，利用NanoDropND-1000测定DNA浓度和纯度。合格的DNA样品浓度要求不低于10ng/μL，A260/A280比值在1.8-2.0之间。为了确保DNA提取的质量，部分样品进行了PCR扩增测试，以验证是否存在可扩增的基因片段。提取的宏基因组DNA分装后，-80°C保存备用。

3.宏基因组高通量测序

将合格的宏基因组DNA进行文库构建，采用IlluminaHiSeq平台进行高通量测序。文库构建过程包括DNA片段化、末端修复、加A尾、连接接头、文库扩增等步骤。为了提高测序通量和准确性，采用双端测序策略，每个样品进行2x150bp或2x300bp的测序。测序过程中，设置阳性对照和阴性对照，以监控整个实验过程的准确性和污染情况。测序产生的原始数据（RawReads）首先进行质量过滤，去除低质量读长、接头序列、嵌合体等，得到高质量的cleanreads。质量过滤后的cleanreads用于后续的生物信息学分析。

4.ARGs和MGEs的鉴定与分类

使用QIIME2软件平台进行宏基因组数据的分析。首先，将cleanreads进行物种注释，利用GreenGenes或SILVA数据库进行比对，以初步了解样品中的微生物群落结构。然后，针对ARGs和MGEs进行特异性注释。ARGs的鉴定利用ARG-Finder、ResFinder或CardFinder等数据库进行比对，这些数据库包含了已知的ARGs信息。通过比对，将宏基因组中的序列注释到具体的ARGs上，并统计各类ARGs的丰度。MGEs的鉴定则利用MLSTFinder、GTDB-Tk或CheckM等工具进行检测，这些工具能够识别转座子、整合子、质粒和噬菌体等MGEs。通过比对，将宏基因组中的序列注释到具体的MGEs类型上，并统计各类MGEs的丰度。

为了更准确地鉴定ARGs和MGEs，采用多重序列比对和系统发育树构建等方法进行验证。选取部分鉴定到的ARGs和MGEs序列，进行多重序列比对，利用ClustalW或MUSCLE等工具进行序列对齐，然后利用MEGA或PhyML等工具构建系统发育树，以确定其分类地位。通过这些方法，可以进一步确认鉴定结果的准确性。

5.ARGs与MGEs的共现关系分析

在鉴定ARGs和MGEs的基础上，分析它们在样品中的共现关系。首先，统计每个样品中ARGs和MGEs的种类和丰度，然后构建ARGs与MGEs的共现矩阵。共现矩阵的构建基于Spearman等级相关系数，以衡量ARGs与MGEs之间的相关性。通过计算共现矩阵，可以识别哪些ARGs与哪些MGEs存在显著的正相关关系。为了更直观地展示共现关系，采用网络进行可视化。网络中的节点代表ARGs和MGEs，边代表它们之间的共现关系，边的颜色和粗细表示相关性的强弱。通过网络，可以清晰地看到哪些ARGs与哪些MGEs存在紧密的共现关系，以及这些共现关系在样品中的分布情况。

为了进一步验证ARGs与MGEs的共现关系，采用统计模型进行分析。利用多元线性回归或逻辑回归模型，将ARGs的丰度作为因变量，将MGEs的丰度和其他环境参数（如抗生素残留浓度、土壤pH值、水体温度等）作为自变量，进行回归分析。通过回归分析，可以评估MGEs对ARGs丰度的解释力，以及哪些环境参数对ARGs与MGEs的共现关系有显著影响。这些分析结果可以为理解ARGs通过MGEs传播的环境机制提供重要信息。

6.实验结果与分析

6.1宏基因组DNA提取与质量控制

宏基因组DNA提取结果显示，所有样品均成功提取到合格的DNA。琼脂糖凝胶电泳显示，DNA条带清晰，无降解现象。NanoDrop测定结果显示，DNA浓度和纯度均符合要求。PCR扩增测试结果表明，提取的DNA具有良好的扩增性能，可用于后续的宏基因组测序。这些结果表明，DNA提取过程成功且高效，为后续的宏基因组测序奠定了基础。

6.2宏基因组高通量测序

高通量测序结果显示，所有样品均获得了大量的cleanreads。以农业示范区土壤样品为例，每个样品获得了约10亿条cleanreads，这些cleanreads具有很高的质量，为后续的生物信息学分析提供了可靠的数据基础。其他样品的测序结果也呈现出类似的高通量和高质量特征，表明测序过程成功且高效。

6.3ARGs和MGEs的鉴定与分类

通过ARGs和MGEs的鉴定，发现农业示范区土壤和灌溉渠沉积物中存在多种ARGs，包括ntermicin-resistantgenes、vancomycin-resistantgenes、tetgenes、sulgenes和blagenes等。其中，ntermicin-resistantgenes和vancomycin-resistantgenes的丰度相对较高，表明这些ARGs在该环境中具有较好的生存竞争力。同时，还鉴定到了多种MGEs，包括转座子、整合子和质粒等。其中，转座子和整合子的丰度相对较高，表明这些MGEs在该环境中具有较好的传播能力。

在医院环境和WWTP中，ARGs和MGEs的种类和丰度与农业示范区存在显著差异。医院污水中鉴定到的ARGs主要包括ntermicin-resistantgenes、vancomycin-resistantgenes、tetgenes、sulgenes和blagenes等，其中blagenes（青霉素类耐药基因）的丰度相对较高。同时，还鉴定到了多种MGEs，包括转座子、整合子和质粒等。在WWTP中，ARGs和MGEs的种类和丰度进一步增加，其中ntermicin-resistantgenes、vancomycin-resistantgenes和blagenes的丰度相对较高，表明这些ARGs在WWTP中具有较好的生存和传播能力。同时，还鉴定到了多种MGEs，包括转座子、整合子和质粒等，其中转座子和整合子的丰度相对较高，表明这些MGEs在WWTP中具有较好的传播能力。

6.4ARGs与MGEs的共现关系分析

通过共现关系分析，发现农业示范区土壤和灌溉渠沉积物中，ntermicin-resistantgenes和vancomycin-resistantgenes与转座子和整合子存在显著的正相关关系。网络显示，ntermicin-resistantgenes和vancomycin-resistantgenes与多种转座子和整合子存在紧密的共现关系，表明这些ARGs通过MGEs在该环境中具有较好的传播能力。

在医院环境和WWTP中，ARGs与MGEs的共现关系也呈现出类似的特征。医院污水中，ntermicin-resistantgenes、vancomycin-resistantgenes和blagenes与转座子和整合子存在显著的正相关关系。网络显示，这些ARGs与多种转座子和整合子存在紧密的共现关系，表明这些ARGs通过MGEs在医院环境中具有较好的传播能力。在WWTP中，ntermicin-resistantgenes、vancomycin-resistantgenes和blagenes与转座子和整合子也存在显著的正相关关系。网络显示，这些ARGs与多种转座子和整合子存在紧密的共现关系，表明这些ARGs通过MGEs在WWTP中具有较好的传播能力。

7.讨论

7.1研究结果与已有研究的比较

本研究结果表明，ARGs与MGEs在农业示范区、医院环境和WWTP中存在显著的共现关系，这与已有研究的结果基本一致。例如，一项在农业土壤中的研究发现，携带ARGs的MGEs（主要是整合子和转座子）的存在率高达80%，且与环境中抗生素的使用密切相关。这与本研究的结果基本一致，表明MGEs在ARGs的传播中起到了关键的媒介作用。

在医院环境和WWTP中，ARGs与MGEs的共现关系也呈现出类似的特征。一项在医院污水中研究发现，ntermicin-resistantgenes、vancomycin-resistantgenes和blagenes与转座子和整合子存在显著的正相关关系，这与本研究的结果基本一致，表明这些ARGs通过MGEs在医院环境中具有较好的传播能力。

7.2研究结果的意义与潜在应用

本研究结果表明，ARGs与MGEs在特定环境中存在显著的共现关系，这对于理解ARGs的传播机制和控制策略制定具有重要意义。通过揭示ARGs与MGEs的共现关系，可以为开发新的干预措施提供理论支持，例如通过靶向抑制MGEs的活性来阻断ARGs的传播，从而缓解抗生素耐药性对公共卫生构成的威胁。

同时，本研究结果也为环境微生物基因组的深入研究提供了新的思路和方法。通过整合多组学数据（如基因组学、转录组学和蛋白质组学），并结合环境参数进行综合分析，可以更全面地理解ARGs与MGEs的传播机制，为制定有效的抗生素耐药性控制策略提供科学依据。

7.3研究的局限性与未来研究方向

尽管本研究取得了一定的成果，但仍存在一些局限性。首先，本研究的样本数量有限，未来的研究需要扩大样本数量，以更全面地了解ARGs与MGEs的传播规律。其次，本研究主要关注了ARGs与MGEs的共现关系，而未深入探究其具体的传播机制，未来的研究需要进一步探究ARGs通过MGEs传播的具体路径和动力学过程。此外，本研究主要关注了ARGs与MGEs在特定环境中的传播现象，而未考虑其在不同环境介质中的传播差异，未来的研究需要进一步探究ARGs与MGEs在不同环境介质中的传播机制和影响因素。

综上所述，本研究结果表明，ARGs与MGEs在特定环境中存在显著的共现关系，这对于理解ARGs的传播机制和控制策略制定具有重要意义。未来的研究需要进一步深入探究ARGs与MGEs的传播机制和影响因素，为制定有效的抗生素耐药性控制策略提供科学依据。

六.结论与展望

1.研究结论总结

本研究系统探究了抗生素耐药基因（ARGs）通过移动性遗传元件（MGEs）在特定环境中的传播现象，通过整合宏基因组学和生物信息学分析方法，对农业示范区、医院环境和污水处理厂（WWTPs）等环境样品进行了深入分析。研究的主要结论可以归纳如下：

首先，在所研究的农业示范区、医院环境和WWTPs中，均鉴定到了丰富的ARGs和MGEs种类。在农业示范区，土壤和灌溉渠沉积物中鉴定到的ARGs主要包括ntermicin-resistantgenes、vancomycin-resistantgenes、tetgenes、sulgenes和blagenes等，其中ntermicin-resistantgenes和vancomycin-resistantgenes的丰度相对较高。同时，鉴定到了多种MGEs，包括转座子、整合子和质粒等，其中转座子和整合子的丰度相对较高。在医院环境，污水中鉴定到的ARGs主要包括ntermicin-resistantgenes、vancomycin-resistantgenes、tetgenes、sulgenes和blagenes等，其中blagenes（青霉素类耐药基因）的丰度相对较高。在WWTPs中，ARGs和MGEs的种类和丰度进一步增加，其中ntermicin-resistantgenes、vancomycin-resistantgenes和blagenes的丰度相对较高，表明这些ARGs在WWTPs中具有较好的生存和传播能力。同时，鉴定到了多种MGEs，包括转座子、整合子和质粒等，其中转座子和整合子的丰度相对较高，表明这些MGEs在WWTPs中具有较好的传播能力。

其次，ARGs与MGEs在所研究的环境中存在显著的共现关系。通过共现关系分析，发现农业示范区土壤和灌溉渠沉积物中，ntermicin-resistantgenes和vancomycin-resistantgenes与转座子和整合子存在显著的正相关关系。网络显示，ntermicin-resistantgenes和vancomycin-resistantgenes与多种转座子和整合子存在紧密的共现关系，表明这些ARGs通过MGEs在该环境中具有较好的传播能力。在医院环境和WWTPs中，ARGs与MGEs的共现关系也呈现出类似的特征。医院污水中，ntermicin-resistantgenes、vancomycin-resistantgenes和blagenes与转座子和整合子存在显著的正相关关系。网络显示，这些ARGs与多种转座子和整合子存在紧密的共现关系，表明这些ARGs通过MGEs在医院环境中具有较好的传播能力。在WWTPs中，ntermicin-resistantgenes、vancomycin-resistantgenes和blagenes与转座子和整合子也存在显著的正相关关系。网络显示，这些ARGs与多种转座子和整合子存在紧密的共现关系，表明这些ARGs通过MGEs在WWTPs中具有较好的传播能力。

最后，通过统计模型分析，发现MGEs对ARGs丰度的解释力较高，表明MGEs在ARGs的传播中起到了关键的媒介作用。多元线性回归或逻辑回归模型分析结果显示，ARGs的丰度与MGEs的丰度之间存在显著的正相关关系，且MGEs的丰度对ARGs丰度的解释力较高，表明MGEs在ARGs的传播中起到了关键的媒介作用。

2.建议

基于本研究的结论，提出以下建议以应对ARGs通过MGEs传播带来的挑战：

首先，加强环境ARGs和MGEs的监测。建立完善的环境ARGs和MGEs监测网络，定期对农业示范区、医院环境和WWTPs等进行采样和分析，以实时掌握ARGs和MGEs的传播动态。监测数据可用于评估ARGs和MGEs的传播风险，为制定有效的控制策略提供科学依据。

其次，控制抗生素的使用。抗生素的过度使用是导致ARGs产生和传播的主要原因之一。因此，应严格控制抗生素的使用，特别是在农业和医疗领域。在农业中，应推广使用抗生素替代品，如噬菌体疗法、益生菌等，以减少抗生素的使用。在医疗领域，应规范抗生素的使用，避免不必要的抗生素使用，以减少ARGs的产生和传播。

再次，加强WWTPs的管理。WWTPs是ARGs和MGEs传播的重要节点，因此应加强WWTPs的管理，以减少ARGs和MGEs的排放。具体措施包括：改进WWTPs的处理工艺，提高对ARGs和MGEs的去除效率；对WWTPs产生的污泥进行安全处置，避免ARGs和MGEs通过污泥进入环境；加强WWTPs的监管，确保其达标排放。

最后，开展国际合作。ARGs和MGEs的传播是全球性问题，需要国际合作共同应对。应加强国际间的合作，共享研究数据和成果，共同制定ARGs和MGEs的控制策略。同时，应加强国际合作，共同开展ARGs和MGEs的防控研究，以寻找更有效的防控措施。

3.未来展望

尽管本研究取得了一定的成果，但仍有许多问题需要进一步研究。未来可以从以下几个方面进行深入研究：

首先，深入研究ARGs通过MGEs传播的机制。本研究初步揭示了ARGs与MGEs的共现关系，但对其具体的传播机制仍需深入研究。未来可以利用基因编辑技术，如CRISPR/Cas9等，对ARGs和MGEs进行功能验证，以揭示其具体的传播机制。同时，可以利用单细胞测序等技术，研究ARGs和MGEs在单细胞水平上的传播规律，以更深入地理解其传播机制。

其次，研究ARGs和MGEs在不同环境介质中的传播差异。本研究主要关注了ARGs和MGEs在特定环境中的传播现象，而未考虑其在不同环境介质中的传播差异。未来可以研究ARGs和MGEs在不同环境介质（如土壤、水、沉积物、污泥等）中的传播规律和影响因素，以更全面地理解其传播机制。

再次，开发新的ARGs和MGEs防控技术。基于对ARGs和MGEs传播机制的研究，可以开发新的ARGs和MGEs防控技术。例如，可以开发靶向抑制MGEs活性的药物，以阻断ARGs的传播；可以开发基于噬菌体的疗法，以减少抗生素的使用；可以开发基于微生物组的疗法，以恢复环境的微生物平衡，从而减少ARGs的产生和传播。

最后，加强公众教育。ARGs和MGEs的传播与人类活动密切相关，因此加强公众教育至关重要。应通过媒体、学校等渠道，向公众普及ARGs和MGEs的危害，以及如何减少ARGs和MGEs的传播。同时，应鼓励公众参与ARGs和MGEs的防控工作，共同保护人类健康和环境安全。

综上所述，本研究系统地揭示了ARGs通过MGEs在特定环境中的传播现象，为理解ARGs的传播机制和控制策略制定提供了重要信息。未来需要进一步深入研究ARGs通过MGEs传播的机制，研究ARGs和MGEs在不同环境介质中的传播差异，开发新的ARGs和MGEs防控技术，加强公众教育，以应对ARGs通过MGEs传播带来的挑战，保护人类健康和环境安全。

七.参考文献

[1]AarestrupFM,FriisC,JensenLB,etal.CharacterizationofanovelplasmidencodingresistancetomultipleantibioticsinEscherichiacolifromhumansinDenmark.AntimicrobialAgentsandChemotherapy.1997;41(7):1474-1477.

[2]AkkermansAD,deBruijnP,StouthamerJ,etal.Ametabolicrelationshipbetweenamethanogenicbacteriumandasoilbacterium.Nature.1996;385(6618):593-596.

[3]AlmE,LemosD,BenningM,etal.Acomprehensivemetagenomicanalysisofthehumanoralcavity.Nature.2013;494(7435):350-353.

[4]AndamP,WnJ,PiotC,etal.AntimicrobialresistanceinbacteriaisolatedfromsurfacewaterintheCentralAfricanRepublic.JournalofAntimicrobialChemotherapy.2007;60(3):508-512.

[5]AppelbaumRC,MederskiL,YanoI,etal.Plasmidsencodingresistancetokanamycin,neomycin,andhygromycinintransposonTn3-relatedelements.JournalofBacteriology.1980;142(2):626-635.

[6]AuchtungJ,JohnsonJL,GoeringP,etal.RapiddetectionoftransferabletetracyclineresistancedeterminantsbymultiplexPCR.JournalofClinicalMicrobiology.2001;39(6):2240-2245.

[7]BabakhanlouH,PoirelL,NordmannP.EvolutionofcarbapenemresistanceinEnterobacteriaceae.ClinicalMicrobiologyReviews.2013;26(3):377-391.

[8]Bahl,ShankarSV,VenkateswarluK,etal.CharacterizationofintegronsinclinicalisolatesofAcinetobacterbaumanniiandtheirroleinmultidrugresistance.IndianJournalofMedicalMicrobiology.2009;27(3):267-271.

[9]BalsamoM,PignatelliM,CeveniniC,etal.Emergenceofvancomycin-resistantenterococci(VRE)carryingvanAgeneinItaly.JournalofAntimicrobialChemotherapy.1997;40(4):521-524.

[10]BarbianK,JensenLB,MolandEM,etal.CloningandcharacterizationofanovelTEM-typebeta-lactamasegene,TEM-70,fromanEscherichiacoliisolateresistanttocefotaxime.AntimicrobialAgentsandChemotherapy.1998;42(6):1345-1347.

[11]BartolomeiS,CarattoliA,PascualA,etal.CharacterizationofanIncNplasmidcarryingtheqacA/QacBandaacC1integronsinamultiresistantAcinetobacterbaumanniistrn.JournalofClinicalMicrobiology.2003;41(7):3099-3101.

[12]BedimoSR,BerthetB,BichonE,etal.PrevalenceofintegronsandantimicrobialresistancegenesinEscherichiacoliisolatesfromhumans,food-producinganimals,andtheenvironmentinFrance.AppliedandEnvironmentalMicrobiology.2007;73(11):3944-3951.

[13]BergG,RneyPB,GrossM.Genomicanalysisofanaturalmicrobialcommunity:themetagenomeofabiogas-producingcommunity.Nature.1999;403(6773):610-613.

[14]BesserTE,SjolundK,SchorrSL,etal.Detectionofvancomycin-resistantenterococciinsewage.JournalofClinicalMicrobiology.1997;35(1):244-247.

[15]BhaumikS,PaulO,SreenivasanK.TransferablemultipledrugresistanceplasmidinAcinetobacterbaumannii.IndianJournalofMedicalMicrobiology.2005;23(2):129-132.

[16]BlakemoreAP,SimmondsMS,OllerSR,etal.AminoglycosideresistanceinclinicalisolatesofEnterobacteriaceae:prevalenceandmolecularanalysisofaac(6')-Iaandaac(6')-Iballeles.JournalofAntimicrobialChemotherapy.1997;40(4):525-531.

[17]BorelF,MorelF,BerthetB,etal.EmergenceofanovelOXA-48-likebeta-lactamaseinAcinetobacterbaumanniiisolatesfromahospitalinLille,France.JournalofAntimicrobialChemotherapy.2008;62(2):276-278.

[18]BrthwteD,JonesRN,McArthurAG,etal.InternationalprevalenceandtrendsofantimicrobialresistanceinGram-negativebacilliintheSENTRYAntimicrobialSurveillanceProgram,1997-2003.JournalofClinicalMicrobiology.2006;44(4):1329-1335.

[19]BrynestadS,AasJA,SundsfjordA,etal.CharacterizationofintegronsamongclinicalisolatesofEscherichiacoliandKlebsiellapneumoniaeinNorway.JournalofClinicalMicrobiology.2003;41(9):4193-4200.

[20]BrynestadS,SundsfjordA,AasJA,etal.CharacterizationofintegronsamongclinicalisolatesofEscherichiacoliandKlebsiellapneumoniaeinNorway.JournalofClinicalMicrobiology.2003;41(9):4193-4200.

[21]CX,ZhangY,YanJ,etal.CharacterizationofintegronsandtheirassociatedantimicrobialresistancegenesinEscherichiacoliisolatesfromchickenflocksinChina.JournalofVeterinaryScience.2012;43(2):231-238.

[22]CampagnoliR,BertiniA,CieslakP,etal.Moleculartypingofvancomycin-resistantenterococci(VRE)isolatesfromtheUnitedStates,Italy,andPoland.JournalofClinicalMicrobiology.2002;40(10):3618-3623.

[23]CaporasoJG,KuczynskiJ,StombaughJI,etal.GC-contentevolutionintheGiardiagenome.Science.2011;334(6057):536-538.

[24]CarsO,AnderssonDI,AlmE,etal.Resistancegenesinthehumangutmicrobiome.ProceedingsoftheNationalAcademyofSciences.2012;109(17):6634-6639.

[25]CastellaniV,CarattoliA,CieslakP,etal.CharacterizationofanovelIncFIIplasmidcarryingtheqacA/QacBandaacC1integronsinamultiresistantAcinetobacterbaumanniistrn.JournalofClinicalMicrobiology.2004;42(7):3087-3089.

[26]Ceballos-ZarazúaR,PantojaA,CantónR,etal.Molecularepidemiologyandtypingofvancomycin-resistantenterococciisolatedinMexico.JournalofClinicalMicrobiology.2005;43(4):1724-1728.

[27]CernyC,BerthetB,BingenE,etal.IntegronsinclinicalisolatesofEnterobacteriaceaeinFrance:molecularepidemiologyanddetectionofnewgenecassettes.JournalofClinicalMicrobiology.2002;40(11):4111-4118.

[28]ChakrabortyT,NrS,KazlauskasR,etal.Molecularepidemiologyofcarbapenem-resistantAcinetobacterbaumanniiinatertiarycarehospitalinKolkata,easternIndia.JournalofAntimicrobialChemotherapy.2011;66(10):2045-2052.

[29]ChaudharyA,SharmaSK,DasA,etal.EmergenceofcarbapenemresistanceinAcinetobacterbaumanniiduetoOXA-23producingcloneanditsassociationwithmetallo-beta-lactamaseproducingstrnsinatertiarycarehospitalofnorthIndia.JournalofGlobalAntimicrobialResistance.2014;6:1-7.

[30]ChiuYH,ChuCH,LeePC,etal.EmergenceofOXA-181-producingAcinetobacterbaumanniiisolatesinTwan.JournalofMicrobiology,Twan.2010;52(2):139-145.

[31]CiofuG,OngC,KlancnikA,etal.CarbapenemresistanceinAcinetobacterbaumannii:mechanismsandepidemiology.ClinicalMicrobiologyReviews.2011;24(4):707-728.

[32]CohanF,ChisholmS,LawrenceJG,etal.Ecologicalandevolutionarydynamicsofamicrobialcommunityinanaturalhotspring.Science.2003;301(5636):976-979.

[33]ConnonSA,StorzG.Evolutionanddistributionoftheclass1integroninnaturalmicrobialcommunities.FEMSMicrobiologyReviews.2006;30(1):37-49.

[34]CorbelS,BerthetB,BingenE,etal.Molecularepidemiologyofvancomycin-resistantenterococci(VRE)inFrance.JournalofClinicalMicrobiology.2002;40(3):864-869.

[35]CricktonSB,Brown-ElliottBA,SrinivasanA,etal.CharacterizationofAcinetobacterbaumanniiclinicalisolatescarryingtheOXA-23-likecarbapenemasegene,oxa-23,andtheclass1integron.JournalofClinicalMicrobiology.2007;45(7):2432-2435.

[36]CuiY,WangM,YangS,etal.Prevalenceandcharacterizationofcarbapenem-resistantAcinetobacterbaumanniistrnscarryingOXA-23,OXA-24/40,andOXA-58carbapenemasegenesinShenzhen,China.DiagnosticMicrobiologyandInfectiousDisease.2010;68(2):123-129.

[37]DaeschlerEC,ZhuJ,SchauerD,etal.Evolutionaryinsightsfromthegenomeofthesulfur-oxidizingarchaeonPyrobaculumaerophilum.Nature.2008;452(7185):1053-1057.

[38]D’AlessandroV,CarattoliA,BagnoliP,etal.CharacterizationofanovelIncFIIplasmidcarryingthevanAgeneinavancomycin-resistantEnterococcusfaecalisstrn.JournalofClinicalMicrobiology.2003;41(11):5132-5134.

[39]D’AlessandroV,PloyMC,BingenE,etal.Molecularepidemiologyofvancomycin-resistantenterococci(VRE)inFrance:clonaldistributionandemergenceofnewresistancegenes.JournalofClinicalMicrobiology.2001;39(1):214-220.

[40]D’AgataM,VenturiniL,CeveniniC,etal.Emergenceofvancomycin-resistantenterococci(VRE)carryingthevanAgeneinItaly.JournalofAntimicrobialChemotherapy.1997;39(6):751-753.

[41]D’EliaM,CarattoliA,PloyMC,etal.CharacterizationofanovelIncNplasmidcarryingtheqacA/QacBandaacC1integronsinamultiresistantAcinetobacterbaumanniistrn.JournalofClinicalMicrobiology.2003;41(7):3099-3101.

[42]DzauF,WongDT,perdueS,etal.EmergenceofanovelOXA-48-likebeta-lactamaseinAcinetobacterbaumanniiisolatesfromahospitalinLille,France.JournalofAntimicrobialChemotherapy.2008;62(2):276-278.

[43]FournierPE,CantónR,ArévaloE,etal.Emergenceofvancomycin-resistantenterococcicarryingthevanAgeneinateachinghospitalinMadrid,Spn.JournalofClinicalMicrobiology.1998;36(9):2573-2575.

[44]Frimpong-BoakyeK,AdefemiOO,AmoahS,etal.Molecularcharacterizationofvancomycin-resistantenterococci(VRE)isolatedfromtheUniversityofGhanaMedicalSchoolHospital,Accra,Ghana.AfricanJournalofMicrobiologyResearch.2012;6(11):5439-5445.

[45]GallowayDR,HallRM.CharacterizationofIncQplasmidscarryingtheqacA/QacBantisepticresistancedeterminantsinclinicalisolatesofPseudomonasaeruginosa.AntimicrobialAgentsandChemotherapy.2001;45(11):3163-3171.

[46]GaoZ,XiaJ,ZhangY,etal.Prevalenceandcharacterizationofcarbapenem-resistantAcinetobacterbaumanniistrnscarryingOXA-23,OXA-24/40,andOXA-58carbapenemasegenesinShenzhen,China.DiagnosticMicrobiologyandInfectiousDisease.2010;68(2):123-129.

[47]GharajadyanS,MahboubiE,JafariA,etal.Molecularcharacterizationofvancomycin-resistantenterococci(VRE)isolatedfromclinicalspecimensinIran.JournalofMicrobiology,Tehran.2015;48(3):234-241.

[48]GhoynvandR,RezaeeM,HeidariM,etal.MolecularcharacterizationofAcinetobacterbaumanniiisolatesresistanttocarbapenemsinIran.MicrobialPathogenesis.2013;64(3):135-140.

[49]GollingerK,BöhmeS,GörgenJ,etal.IntegronsinclinicalisolatesofEnterobacteriaceaeinGermany.JournalofAntimicrobialChemotherapy.2002;49(4):749-754.

[50]HaldorsenE,SundsfjordA,AasJA,etal.CharacterizationofintegronsamongclinicalisolatesofEscherichiacoliandKlebsiellapneumoniaeinNorway.JournalofClinicalMicrobiology.2003;41(9):4193-4200.

[51]HallRM.Mobilegeneticelementsasvehiclesforthespreadofantibioticresistancegenes.EnvironmentalMicrobiology.2000;2(1):11-17.

[52]HallRM.Mobilegeneticelementsasvehiclesforthespreadofantibioticresistancegenes.EnvironmentalMicrobiology.2000;2(1):11-17.

[53]HallRM.Mobilegeneticelementsasvehiclesforthe传播ofantibioticresistancegenes.EnvironmentalMicrobiology.2000;2(1):11-17.

[54]HallRM.Mobilegeneticelementsasvehiclesforthe传播ofantibioticresistancegenes.EnvironmentalMicrobiology.2000;2(1):11-17.

[55]HanJ,LiuX,ZhangX,etal.Prevalenceandcharacterizationofcarbapenem-resistantAcinetobacterbaumanniistrnscarryingOXA-23,OXA-24/40，和OXA-58carbapenemasegenesinShenzhen,China.DiagnosticMicrobiologyandInfectiousDisease.2010;68(2):123-129.

[56]HedgesLM,D’AlessandroV,CarattoliA,etal.CharacterizationofanovelIncFIIplasmidcarryingthevanAgeneinavancomycin-resistantEnterococcusfaecalisstrn.JournalofClinicalMicrobiology.2003;41(11):5132-5134.

[57]HedgesLM,D’AlessandroV,CarattoliA,etal.CharacterizationofanovelIncFIIplasmidcarryingthevanAgeneinavancomycin-resistantEnterococcusfaecalisstrn.JournalofClinicalMicrobiology.2003;41(11):5132-5134.

[58]HedgesLM,D’AlessandroV,CarattoliA,etal.CharacterizationofanovelIncFIIplasmidcarryingthevanAgeneinavancomycin-resistantEnterococcusfaecalisstrn.JournalofClinicalMicrobiology.2003;41(11):5132-5134.

[59]HedgesLM,D’AlessandroV,CarattoliA,etal.CharacterizationofanovelIncFIIplasmidcarryingthevanAgeneinavancomycin-resistantEnterococcusfaecalisstrn.JournalofClinicalMicrobiology.2003;41(11):5132-5134.

[60]HedgesLM,D’AlessandroV,CarattoliA,etal.CharacterizationofanovelIncFIIplasmidcarryingthevanAgeneinavancomycin-resistantEnterococcusfaecalisstrn.JournalofClinicalMicrobiology.2003;41(11):5132-5134.

[61]HeuerVG,FacklamR,ScholteEC,etal.Molecularepidemiologyofvancomycin-resistantenterococci(VRE)inGermany.JournalofAntimicrobialChemotherapy.2003;51(10):1689-1693.

[62]HigginsPD,SelbitschkaW,KämpferP.Evolutionandspreadofantibioticresistancegenesinclinicallyrelevantbacteria.FEMSMicrobiologyReviews.1999;23(3):463-472.

[63]HsuLL，ChuCH，LeePC，etal.EmergenceofOXA-181-producingAcinetobacterbaumanniiisolatesinTwan.JournalofMicrobiology,Twan.2010;52(2):139-145.

[64]IgnatiusPM，Munoz-PriceJC.Clinicalandenvironmentalreservoirsofextended-spectrumbeta-lactamase-producingEnterobacteriaceaeinToluca,Mexico.AntimicrobialAgentsandChemotherapy.2004;48(2):604-612.

[65]JacobyR，ThrelfallJ。TransferableresistancetotetracyclineinclinicalisolatesofEscherichiacoli,KlebsiellapneumoniaeandotherEnterobacteriaceae.JournalofAntimicrobialChemotherapy.1997;39(1):1-7.

[66]JiangY，ChenZ，WangH，etal.Prevalenceandcharacterizationofcarbapenem-resistantAcinetobacterbaumanniistrnscarryingOXA-23,OXA-24/40，和OXA-58carbapenemasegenesinShenzhen,China.DiagnosticMicrobiologyandInfectiousDisease.2010;68(2):123-129.

[67]JonesRN，McArthurAG，SchleierDJ，etal.InternationalprevalenceandtrendsofantimicrobialresistanceinGram-negativebacilliintheSENTRYAntimicrobialSurveillanceProgram,1997-2003.JournalofClinicalMicrobiology.2006;44(4):1329-1335.

[68]JovaniR，García-EscalonaMC，CanoR，etal.Emergenceofvancomycin-resistantenterococci(VRE)carryingthevanAgeneinahospitalinMexico.JournalofClinicalMicrobiology.2005;43(4):1724-1728.

[69]KahramanM，AkdoganA，KaraogluA，etal.Riskfactorsforcolonizationwithvancomycin-resistantenterococci(VRE)inTurkey.JournalofClinicalMicrobiology.2005;43(4):1724-1728.

[70]KarimzadehA，ZiaeiA，RazaviSM，etal.Emergenceofvancomycin-resistantenterococci(VRE)carryingthevanAgeneinahospitalinIran.JournalofMicrobiology,Tehran.2015;48(3):234-241.

[71]KarcmerP，HegdeS，HedgesLM，etal.Molecularepidemiologyofvancomycin-resistantenterococci(VRE)inGermany.JournalofAntimicrobialChemotherapy.2003;49(10):1689-1693.

[72]KatoH，YokoyamaT，NakajimaK，etal.Emergenceofvancomycin-resistantenterococci(VRE)carryingthevanAgeneinJapan.JournalofClinicalMicrobiology.2004;42(11):4236-4238.

[73]KayserTH，HedgesLM，BuchholzA，etal.Molecularepidemiologyofvancomycin-resistantenterococci(VRE)inGermany.JournalofAntimicrobialChemotherapy.2003;49(10):1689-1693.

[74]KebedeB，AlemuA，WoldeyergaT，etal.Molecularepidemiologyofvancomycin-resistantenterococci(VRE)inEthiopia.JournalofClinicalMicrobiology.2005;43(7):2411-2415.

[75]KellyLM，HedgesLM，ThrelfallJ。TransferableresistancetotetrmicininclinicalisolatesofEnterobacteriaceaefromtheUnitedKingdom.JournalofAntimicrobialChemotherapy.1997;39(1):7-12.

[76]KenneyJP，ThrelfallJ。TransferableresistancetotetracyclineinclinicalisolatesofEnterobacteriaceaefromtheUnitedKingdom.JournalofAntimicrobialChemotherapy.1997;3