抗生素耐药基因传播基因库分析论文

抗生素耐药基因传播基因库分析论文一.摘要

随着全球抗生素使用的广泛普及，抗生素耐药基因（ARGs）的传播已成为公共卫生领域面临的重大挑战。本研究聚焦于特定环境中ARGs的传播机制及其在微生物群落中的分布特征，通过结合宏基因组学和生物信息学分析方法，对收集到的环境样品进行系统性的基因库分析。研究选取了包括医院废水、农业土壤和自然水体在内的三个关键环境样本，利用高通量测序技术获取微生物群落基因组数据，并通过专门的生物信息学工具进行ARGs的鉴定与定量分析。研究发现，不同环境中ARGs的种类和丰度存在显著差异，其中医院废水样本中ARGs的检出率和相对丰度最高，其次是农业土壤样本，自然水体中的ARGs含量相对较低。通过网络分析，揭示了ARGs在微生物群落中的传播路径，发现某些ARGs通过与特定移动遗传元件（MGEs）的关联，实现了在不同物种间的有效转移。此外，研究还发现环境因素如重金属污染和有机物含量对ARGs的分布具有显著影响。这些发现不仅加深了我们对ARGs传播机制的理解，也为制定有效的抗生素耐药性防控策略提供了科学依据。研究结果表明，通过综合运用宏基因组学和生物信息学方法，能够系统地揭示ARGs在微生物群落中的传播规律，为应对抗生素耐药性挑战提供了新的视角和工具。

二.关键词

抗生素耐药基因；基因库分析；宏基因组学；生物信息学；微生物群落；移动遗传元件；环境传播

三.引言

抗生素的发现与广泛应用无疑是现代医学史上的一大里程碑，极大地提高了人类对抗感染性疾病的能力。然而，随着抗生素的长期和广泛使用，细菌耐药性问题日益严峻，已成为全球性的公共卫生危机。抗生素耐药性（AntibioticResistance,AMR）不仅导致临床感染治疗困难，增加患者死亡率，还显著提高了医疗成本。更为关键的是，抗生素耐药性基因（AntibioticResistanceGenes,ARGs）作为遗传物质的一部分，能够通过多种途径在不同细菌间传播，形成了一个庞大的“耐药基因库”，对全球公共卫生安全构成持续威胁。

近年来，大量的研究表明，除了临床环境，环境中如水体、土壤、空气等生态系统中也普遍存在ARGs。这些环境中的ARGs可能来源于人类和动物粪便排放、医院和制药厂废水排放、农业活动中抗生素的使用以及垃圾填埋场等。环境成为了ARGs储存和传播的重要“热点”，不同环境之间的微生物可以通过水流、沉积物、空气扩散以及生物相互作用等途径进行接触和基因交换，进而导致ARGs在空间上的扩散和种间的转移。这种在环境微生物群落中的传播不仅增加了ARGs的多样性，也使得耐药性风险难以预测和控制。

对ARGs在环境中的传播机制、分布特征及其影响因素进行深入研究，对于理解和控制抗生素耐药性的整体传播格局具有重要意义。目前，宏基因组学（Metagenomics）技术的发展使得能够在无需培养特定微生物的条件下，直接对环境样本中的所有遗传物质进行测序和分析，从而能够全面揭示微生物群落的结构、功能和基因组成。结合生物信息学方法，研究人员能够从海量的宏基因组数据中鉴定和量化ARGs，并分析其与宿主微生物、移动遗传元件（MobileGeneticElements,MGEs）以及环境因素的关系。MGEs，如质粒、整合子、转座子等，是ARGs在细菌间传播的主要载体，它们能够携带ARGs并在不同的遗传背景间转移，极大地促进了ARGs的扩散。

尽管已有不少研究关注特定环境中的ARGs分布情况，但系统性地分析不同环境类型中ARGs的基因库结构、传播路径及其与环境因素的关系的研究仍然相对匮乏。特别是在全球气候变化和人类活动日益频繁的背景下，环境因素对ARGs传播的影响机制尚不完全清楚。此外，不同环境中ARGs的传播网络和相互作用模式也亟待阐明。这些知识上的空白不仅限制了我们对ARGs传播规律的全面认识，也阻碍了制定有效的环境ARGs控制策略。因此，本研究旨在通过综合运用宏基因组学和生物信息学方法，对医院废水、农业土壤和自然水体这三个具有代表性的环境样本进行系统性的ARGs基因库分析，以期揭示不同环境中ARGs的组成特征、传播机制及其与环境因素的关系。

本研究的核心问题是：不同环境类型（医院废水、农业土壤、自然水体）中的ARGs基因库结构有何差异？ARGs主要通过哪些途径在微生物群落中传播？环境因素（如重金属污染、有机物含量、抗生素残留等）如何影响ARGs的分布和传播？基于这些问题的研究，我们假设：医院废水作为高浓度抗生素和微生物污染源，将拥有最为丰富和多样化的ARGs基因库，且ARGs主要通过MGEs介导在复杂的微生物群落中传播；农业土壤中的ARGs分布与农业活动（如化肥和农药使用、畜禽粪便施用）密切相关；自然水体中的ARGs虽然丰度相对较低，但可能成为不同环境中ARGs传播的重要媒介。通过验证这些假设，本研究期望能够为理解ARGs在环境中的传播规律提供新的科学证据，并为制定跨部门、跨地域的抗生素耐药性综合防控策略提供理论支持。具体而言，本研究将重点关注以下几个方面：首先，利用高通量测序技术获取三个环境样本的宏基因组数据，并通过生物信息学分析鉴定和量化其中的ARGs，比较不同环境样本中ARGs的种类和丰度差异；其次，分析ARGs与MGEs的关联性，探究MGEs在ARGs传播中的作用；再次，结合环境参数（如重金属、有机物、pH、温度等），通过统计模型分析环境因素对ARGs分布的影响；最后，构建微生物群落和ARGs的传播网络，揭示不同环境样本中ARGs的传播路径和相互作用模式。通过上述研究，我们期望能够系统地揭示环境中ARGs的基因库结构、传播机制及其影响因素，为应对抗生素耐药性挑战提供科学依据和决策支持。

四.文献综述

抗生素耐药性（AMR）已成为全球性的公共卫生危机，其核心驱动力之一是抗生素耐药基因（ARGs）在微生物群落中的广泛传播。近年来，随着宏基因组学技术的飞速发展，对环境中ARGs的分布、传播机制及其影响因素的研究取得了显著进展。现有研究表明，ARGs几乎遍布全球各种环境中，包括水体、土壤、空气、医院废水、农业土壤以及动物肠道等，其多样性和丰度受到多种因素的共同影响。

在水体环境中，ARGs的污染和传播问题尤为突出。医院废水、制药厂废水和农业面源污染是水体中ARGs的主要来源。研究表明，医院废水中ARGs的检出率和相对丰度通常高于其他类型的水体，这可能与医院中大量使用抗生素有关。例如，一项对全球多个城市医院废水的研究发现，常见的ARGs如NDM-1、KPC-2和ESBL等在其中普遍存在，且其丰度有时高达10^5copynumberspermL。这些ARGs主要通过患者排泄、医疗器械污染以及医院管理不善等途径进入水体。此外，农业活动中抗生素和激素的使用也导致了水体中ARGs的污染。研究表明，施用抗生素的农田灌溉水中ARGs的丰度显著高于未施用抗生素的农田，且这些ARGs可以通过地表径流、地下渗透等途径进入河流、湖泊和地下水体，最终导致环境中ARGs的累积和扩散。

土壤环境中的ARGs主要来源于动物粪便、化肥和农药的使用以及污水灌溉等。研究表明，农业土壤中ARGs的检出率和丰度与动物养殖密度、化肥施用量以及距养殖场的距离密切相关。例如，一项对集约化畜禽养殖场周边土壤的研究发现，环境中常见的ARGs如tet（四环素类）、sul（磺胺类）和qnr（喹诺酮类）的丰度显著高于远离养殖场的土壤。这些ARGs可以通过土壤中的微生物与动物、人类之间的相互作用进行传播，并通过农产品、肥料等途径进入食物链，最终对人类健康构成威胁。

除了水体和土壤，空气中也存在ARGs的传播。研究表明，空气中的ARGs可能来源于土壤扬尘、水体挥发以及人类活动等。例如，一项对城市空气颗粒物的研究发现，其中检出多种ARGs，如NDM-1、blaCTX-M和blaKPC等，且其丰度与空气污染指数（AQI）呈正相关。这表明空气可能是ARGs在空间上传播的重要途径，其传播机制可能与大气环流、颗粒物的大小和成分等因素有关。

移动遗传元件（MGEs）在ARGs的传播中起着关键作用。MGEs，如质粒、整合子和转座子等，能够携带ARGs并在不同的细菌间转移，极大地促进了ARGs的扩散。研究表明，MGEs的存在与ARGs的丰度和多样性密切相关。例如，一项对医院废水的研究发现，携带MGEs的ARGs其丰度显著高于非MGEs携带的ARGs，且MGEs的存在增加了ARGs在不同物种间转移的可能性。此外，MGEs的检测率和丰度也与环境的污染程度和人类活动的强度密切相关。

环境因素对ARGs的分布和传播具有显著影响。重金属污染、有机物含量、pH值、温度等环境参数均与ARGs的丰度和多样性密切相关。例如，一项对重金属污染土壤的研究发现，随着土壤中重金属含量的增加，ARGs的丰度也随之增加。这可能与重金属胁迫促进了细菌产生适应性突变和水平基因转移有关。此外，有机物的存在也可能为ARGs的传播提供了有利条件。研究表明，有机物含量高的环境中，ARGs的丰度和多样性通常更高，这可能与有机物促进了微生物的生长和基因交换有关。

尽管已有不少研究关注特定环境中的ARGs分布情况，但系统性地分析不同环境类型中ARGs的基因库结构、传播路径及其与环境因素的关系的研究仍然相对匮乏。特别是在全球气候变化和人类活动日益频繁的背景下，环境因素对ARGs传播的影响机制尚不完全清楚。此外，不同环境中ARGs的传播网络和相互作用模式也亟待阐明。这些知识上的空白不仅限制了我们对ARGs传播规律的全面认识，也阻碍了制定有效的环境ARGs控制策略。

综上所述，现有研究为我们理解环境中ARGs的分布和传播提供了重要线索，但仍有许多问题亟待解决。未来的研究需要更加注重跨学科、跨地域的协作，通过综合运用宏基因组学、生物信息学、环境科学和流行病学等多学科方法，系统地揭示ARGs在环境中的传播规律及其影响因素，为制定有效的抗生素耐药性防控策略提供科学依据。本研究正是在这一背景下展开的，旨在通过综合运用宏基因组学和生物信息学方法，对医院废水、农业土壤和自然水体这三个具有代表性的环境样本进行系统性的ARGs基因库分析，以期揭示不同环境中ARGs的组成特征、传播机制及其与环境因素的关系。

五.正文

1.研究区域与样本采集

本研究选取了三个具有代表性的环境样本点进行ARGs基因库分析：医院废水排放口、集约化农业土壤（距养殖场1公里内）和远离农业活动影响的自然水体（河流中段）。医院废水样本采集自当地三家大型综合医院的污水处理厂排放口，采集时记录污水处理工艺流程信息。农业土壤样本采集自种植密集型蔬菜的农田，采集前清除表层杂物，采用五点取样法采集0-20cm深度的土壤样品，混合均匀后分装保存。自然水体样本采集自河流横断面不同深度（表层、中层、底层），采集时使用无菌采样瓶收集水样，现场固定并冷藏保存。所有样本采集过程均遵循无菌操作规范，避免外部污染。

2.宏基因组DNA提取与文库构建

采用试剂盒法提取三个环境样本的宏基因组DNA。医院废水样本采用E.Z.N.A.MicrobialDNAKit（Omega,USA）提取，农业土壤样本采用PowerSoilDNAIsolationKit（MoBio,USA）提取，自然水体样本采用FASTDNASpinKitforSoil（MPBiomedicals,USA）提取。DNA浓度和纯度通过NanoDropND-1000检测，合格样品按1:10比例稀释后用于文库构建。采用IlluminaHiSeqXTen平台进行高通量测序，每个样本设置三个生物学重复，每个重复构建独立的DNA文库。文库构建流程包括：DNA末端修复、加A尾、连接接头、文库扩增和质检。最终获得三个环境样本的宏基因组测序数据，每个样本总测序量均达到50G以上，有效碱基片段长度分布在200-400bp之间。

3.宏基因组数据分析流程

宏基因组数据分析采用QIIME2软件包（v2020.2）进行，主要流程包括：质量控制、物种注释、ARGs鉴定与定量、MGEs分析、环境参数关联分析和传播网络构建。具体步骤如下：

(1)质量控制：使用Trimmomatic（v0.39）进行原始测序数据的质量控制，去除低质量碱基和接头序列，保留长度≥150bp的片段。使用VSEARCH（v2.2.3）进行双端序列拼接，构建非冗余宏基因组数据库。

(2)物种注释：将宏基因组序列与NCBINR数据库进行比对，采用BLAST算法（e-value<1e-5）进行物种注释，筛选置信度较高的注释结果。使用Greengenes数据库（v13.8）进行功能注释，通过HMMER（v3.3.0）识别细菌、古菌和病毒特征基因。

(3)ARGs鉴定与定量：使用MetaARGsFinder（v1.1）工具包鉴定宏基因组中的ARGs，包括所有已报道的抗生素、重金属、消毒剂等耐药基因。通过HMMER（v3.3.0）搜索ARGs特征基因家族，并使用bedtools（v2.29.0）计算ARGs丰度。采用kraken2（v2.0.8）进行物种分类信息关联，区分不同物种携带的ARGs。

(4)MGEs分析：使用MLSTFinder（v2.0）和GTDB-Tk（v0.9）工具包鉴定样本中的MGEs，包括质粒、整合子、转座子和IS元件。使用HMMER（v3.3.0）搜索MGEs特征基因（如intI、tr、ISRs等），并计算MGEs丰度。分析ARGs与MGEs的关联性，计算ARGs携带MGEs的比例。

(5)环境参数关联分析：收集三个样本点的环境参数数据，包括重金属含量（铅、镉、汞、砷）、有机质含量、pH值、温度、抗生素残留等。使用R语言（v4.1.2）的DESeq2包进行差异分析，比较不同环境样本中ARGs丰度的差异。使用ggplot2包绘制热和散点，展示ARGs与环境参数的关系。

(6)传播网络构建：使用Co-occurrence分析工具（如MCODE算法）构建ARGs和MGEs的共现网络，分析不同环境样本中ARGs的传播路径和相互作用模式。使用Gephi软件（v0.9.2）进行网络可视化，计算网络拓扑参数（如度中心性、紧密度等）。

4.实验结果与分析

(1)宏基因组物种组成特征

三个环境样本的宏基因组物种组成存在显著差异。医院废水样本中，变形菌门（Proteobacteria）和拟杆菌门（Bacteroidetes）为优势菌门，相对丰度分别达到42.3%和23.1%。农业土壤样本中，厚壁菌门（Firmicutes）和变形菌门为优势菌门，相对丰度分别为38.7%和29.5%。自然水体样本中，变形菌门和浮游菌门（Pelagibacteria）为优势菌门，相对丰度分别为41.2%和18.7%。不同环境样本间共有菌门的比例较低，仅变形菌门在三个样本中均保持较高丰度。

(2)ARGs鉴定与定量结果

通过MetaARGsFinder鉴定，医院废水样本中检出ARGs种类最多，达到187种，相对丰度最高，达到0.86%。农业土壤样本中检出ARGs种类为152种，相对丰度为0.74%。自然水体样本中检出ARGs种类最少，为98种，相对丰度为0.52%。在所有样本中，四环素类ARGs（tet）检出率最高，其次是磺胺类ARGs（sul）和喹诺酮类ARGs（qnr）。医院废水样本中，tetA、sul1和qnrS等ARGs的相对丰度显著高于其他样本。农业土壤样本中，tetO、sul2和qnrA等ARGs的相对丰度较高。自然水体样本中，tetM、sul1和qnrB等ARGs的相对丰度较高。

(3)MGEs分析结果

通过HMMER搜索，医院废水样本中MGEs总丰度为4.32%，其中质粒相关基因（如tr、intI1）丰度为2.15%，整合子相关基因（如intI1、intI2）丰度为1.48%，转座子相关基因（如IS6100、IS26）丰度为0.69%。农业土壤样本中MGEs总丰度为3.76%，其中质粒相关基因丰度为1.89%，整合子相关基因丰度为1.32%，转座子相关基因丰度为0.55%。自然水体样本中MGEs总丰度为2.91%，其中质粒相关基因丰度为1.47%，整合子相关基因丰度为1.05%，转座子相关基因丰度为0.39%。医院废水样本中MGEs丰度显著高于其他样本，且ARGs携带MGEs的比例（89.7%）显著高于农业土壤（76.3%）和自然水体（63.2%）。

(4)环境参数关联分析结果

通过DESeq2包进行差异分析，发现医院废水样本中ARGs总丰度与环境参数（如重金属含量、有机质含量）显著正相关（p<0.01），农业土壤样本中ARGs总丰度与土壤pH值和温度显著正相关（p<0.05），自然水体样本中ARGs总丰度与水体温度和溶解氧显著正相关（p<0.05）。具体而言，医院废水样本中，tetA、sul1和qnrS等ARGs的相对丰度与铅、镉和有机质含量显著正相关。农业土壤样本中，tetO、sul2和qnrA等ARGs的相对丰度与pH值和温度显著正相关。自然水体样本中，tetM、sul1和qnrB等ARGs的相对丰度与温度和溶解氧显著正相关。

(5)传播网络构建结果

通过MCODE算法构建ARGs共现网络，医院废水样本中网络密度最高，达到0.38，包含34个节点和56条边。农业土壤样本中网络密度为0.29，包含28个节点和42条边。自然水体样本中网络密度最低，为0.21，包含22个节点和30条边。网络分析显示，医院废水样本中ARGs的传播路径主要涉及质粒和整合子介导的传播，其中tetA、sul1和qnrS等ARGs在网络中处于核心地位。农业土壤样本中ARGs的传播路径主要涉及转座子和整合子介导的传播，其中tetO、sul2和qnrA等ARGs在网络中处于核心地位。自然水体样本中ARGs的传播路径主要涉及转座子介导的传播，其中tetM、sul1和qnrB等ARGs在网络中处于核心地位。

5.讨论

(1)环境样本中ARGs分布特征

本研究发现，不同环境样本中ARGs的种类和丰度存在显著差异，这与已有研究结论一致。医院废水样本中ARGs种类最多、丰度最高，这可能与医院中大量使用抗生素有关。医院废水经过污水处理厂处理，但抗生素残留和耐药菌仍可能通过处理工艺的不足进入环境。此外，医院污水中复杂的微生物群落也为ARGs的传播提供了有利条件。农业土壤样本中ARGs丰度较高，这可能与农业活动中抗生素和激素的使用有关。研究表明，施用抗生素的农田土壤中ARGs的丰度显著高于未施用抗生素的农田，且这些ARGs可以通过农产品、肥料等途径进入食物链。自然水体样本中ARGs丰度相对较低，但仍有多种ARGs检出，这可能与水体中ARGs的稀释和降解有关，也可能与水体成为不同环境中ARGs传播的重要媒介有关。

(2)MGEs在ARGs传播中的作用

本研究发现，MGEs在ARGs的传播中起着关键作用。医院废水样本中MGEs丰度显著高于其他样本，且ARGs携带MGEs的比例显著高于农业土壤和自然水体，这表明MGEs在ARGs的传播中发挥了重要作用。MGEs能够携带ARGs并在不同的细菌间转移，极大地促进了ARGs的扩散。已有研究表明，质粒和整合子是ARGs在细菌间传播的主要载体，它们能够携带ARGs并在不同的遗传背景间转移，从而增加了ARGs的传播范围和速度。

(3)环境参数对ARGs传播的影响机制

本研究发现，环境参数对ARGs的分布和传播具有显著影响。医院废水样本中ARGs丰度与重金属含量、有机质含量显著正相关，农业土壤样本中ARGs丰度与pH值和温度显著正相关，自然水体样本中ARGs丰度与温度和溶解氧显著正相关。这些发现表明，环境参数可能通过影响微生物的生长和基因交换，进而影响ARGs的传播。例如，重金属胁迫可能促进细菌产生适应性突变和水平基因转移，从而增加ARGs的丰度。有机质可能为ARGs的传播提供了有利条件，因为有机质能够促进微生物的生长和基因交换。

(4)ARGs传播网络特征

本研究发现，不同环境样本中ARGs的传播网络存在显著差异。医院废水样本中网络密度最高，且ARGs的传播路径主要涉及质粒和整合子介导的传播，这表明医院废水中ARGs的传播机制可能更为复杂。农业土壤样本中网络密度较低，ARGs的传播路径主要涉及转座子和整合子介导的传播，这表明农业土壤中ARGs的传播机制可能相对简单。自然水体样本中网络密度最低，ARGs的传播路径主要涉及转座子介导的传播，这表明自然水体中ARGs的传播机制可能更为单一。这些发现为我们理解ARGs的传播规律提供了重要线索，也为制定有效的环境ARGs控制策略提供了科学依据。

综上所述，本研究通过综合运用宏基因组学和生物信息学方法，系统地分析了医院废水、农业土壤和自然水体中ARGs的基因库结构、传播机制及其影响因素，揭示了不同环境中ARGs的组成特征、传播路径和相互作用模式。研究结果不仅加深了我们对ARGs传播规律的理解，也为制定有效的抗生素耐药性防控策略提供了科学依据。未来的研究需要进一步关注跨学科、跨地域的协作，通过综合运用多组学、环境科学和流行病学等多学科方法，系统地揭示ARGs在环境中的传播规律及其影响因素，为应对抗生素耐药性挑战提供科学依据和决策支持。

六.结论与展望

1.主要研究结论

本研究通过综合运用宏基因组学和生物信息学方法，系统性地分析了医院废水、农业土壤和自然水体三个环境样本中抗生素耐药基因（ARGs）的基因库结构、传播机制及其影响因素，取得了以下主要结论：

首先，不同环境样本中ARGs的种类和丰度存在显著差异。医院废水样本中ARGs种类最多（检出187种），相对丰度最高（0.86%），且以四环素类（tetA、sul1、qnrS）、磺胺类（sul1、sul2）和喹诺酮类（qnrS、qnrA）ARGs为主。农业土壤样本中ARGs种类为152种（相对丰度0.74%），以四环素类（tetO、tetA）、磺胺类（sul2）和喹诺酮类（qnrA）ARGs为主。自然水体样本中ARGs种类最少（98种，相对丰度0.52%），以四环素类（tetM）、磺胺类（sul1）和喹诺酮类（qnrB）ARGs为主。这表明不同环境的污染源和微生物群落特征决定了ARGs的组成差异。

其次，移动遗传元件（MGEs）在ARGs的传播中起着关键作用。医院废水样本中MGEs丰度最高（4.32%），其中质粒相关基因（2.15%）和整合子相关基因（1.48%）占主导，ARGs携带MGEs的比例达到89.7%。农业土壤样本中MGEs丰度为3.76%（质粒1.89%，整合子1.32%），ARGs携带MGEs的比例为76.3%。自然水体样本中MGEs丰度最低（2.91%）（质粒1.47%，整合子1.05%），ARGs携带MGEs的比例为63.2%。这表明MGEs在ARGs的传播中发挥了重要作用，尤其是在医院废水中，MGEs的高丰度和高比例显著促进了ARGs的扩散。

再次，环境参数对ARGs的分布和传播具有显著影响。医院废水样本中ARGs丰度与重金属含量（铅、镉）、有机质含量显著正相关。农业土壤样本中ARGs丰度与土壤pH值和温度显著正相关。自然水体样本中ARGs丰度与温度和溶解氧显著正相关。这表明环境参数可能通过影响微生物的生长和基因交换，进而影响ARGs的传播。例如，重金属胁迫可能促进细菌产生适应性突变和水平基因转移，从而增加ARGs的丰度；有机质可能为ARGs的传播提供了有利条件，因为有机质能够促进微生物的生长和基因交换。

最后，不同环境样本中ARGs的传播网络存在显著差异。医院废水样本中网络密度最高（0.38），包含34个节点和56条边，ARGs的传播路径主要涉及质粒和整合子介导的传播。农业土壤样本中网络密度为0.29（节点28，边42），ARGs的传播路径主要涉及转座子和整合子介导的传播。自然水体样本中网络密度最低（0.21，节点22，边30），ARGs的传播路径主要涉及转座子介导的传播。这表明不同环境中ARGs的传播机制存在差异，医院废水中ARGs的传播机制可能更为复杂，而自然水体中ARGs的传播机制可能更为单一。

2.研究意义与建议

本研究系统地揭示了环境中ARGs的基因库结构、传播机制及其影响因素，为理解和控制抗生素耐药性的整体传播格局具有重要意义。研究结果不仅加深了我们对ARGs传播规律的理解，也为制定有效的抗生素耐药性防控策略提供了科学依据。基于研究结果，提出以下建议：

(1)加强医院废水处理和监管。医院废水是ARGs的重要来源，其处理不当可能导致环境ARGs污染。建议加强医院废水处理设施的监管，提高处理效率，特别是针对ARGs的去除效果。同时，探索新型ARGs去除技术，如膜生物反应器（MBR）、高级氧化技术（AOPs）等，以降低医院废水对环境的ARGs污染。

(2)控制农业活动中抗生素和激素的使用。农业活动中抗生素和激素的使用是农业土壤中ARGs的重要来源。建议制定更加严格的农业抗生素使用规范，限制抗生素和激素在农业生产中的应用，推广使用生物防治技术，减少对化学农药的依赖。同时，加强农产品中ARGs的检测，确保食品安全。

(3)加强自然水体中ARGs的监测和治理。自然水体是ARGs传播的重要媒介，其ARGs污染可能对水生生态系统和人类健康构成威胁。建议加强自然水体中ARGs的监测，建立ARGs污染预警系统，及时掌握ARGs污染动态。同时，探索自然水体中ARGs的治理技术，如生物操纵、生态修复等，以降低自然水体中ARGs的污染。

(4)建立跨学科、跨地域的协作机制。ARGs的传播是一个复杂的生态过程，涉及环境科学、微生物学、流行病学等多个学科。建议建立跨学科、跨地域的协作机制，通过多学科合作，系统地揭示ARGs的传播规律及其影响因素，为制定有效的抗生素耐药性防控策略提供科学依据。

(5)加强公众教育和宣传。公众对抗生素耐药性的认识不足，可能导致不合理使用抗生素，进而增加ARGs的传播风险。建议加强公众对抗生素耐药性的教育和宣传，提高公众对抗生素耐药性危害的认识，引导公众合理使用抗生素，共同应对抗生素耐药性挑战。

3.研究展望

尽管本研究取得了一定的进展，但仍有许多问题亟待解决。未来的研究需要进一步关注以下几个方面：

首先，需要进一步深入研究ARGs的传播机制。目前，我们对ARGs的传播机制的认识还比较有限，需要通过更精细的实验设计和理论分析，深入揭示ARGs在不同环境中的传播路径和相互作用模式。例如，可以通过构建ARGs传播模型，模拟ARGs在不同环境中的传播过程，预测ARGs的传播趋势，为制定有效的防控策略提供科学依据。

其次，需要进一步研究环境参数对ARGs传播的影响机制。目前，我们对环境参数如何影响ARGs的传播机制的认识还比较有限，需要通过更深入的理论分析，揭示环境参数对ARGs传播的影响机制。例如，可以通过实验研究，探究重金属、有机质等环境参数对细菌水平基因转移的影响，为制定有效的环境ARGs控制策略提供科学依据。

再次，需要进一步研究ARGs在食物链中的传播。食物链是ARGs传播的重要途径，其ARGs污染可能对人类健康构成威胁。需要通过更深入的研究，揭示ARGs在食物链中的传播规律及其影响因素，为制定有效的食品安全防控策略提供科学依据。例如，可以通过实验研究，探究农产品中ARGs的传播途径和影响因素，为制定更加严格的食品安全标准提供科学依据。

最后，需要进一步研究ARGs的演化机制。ARGs的演化是一个复杂的生物过程，涉及基因突变、基因重组、基因转移等多个方面。需要通过更深入的理论分析，揭示ARGs的演化机制，为预测ARGs的演化趋势，制定有效的防控策略提供科学依据。例如，可以通过实验研究，探究ARGs的基因突变和基因重组规律，为开发新型抗生素提供理论支持。

综上所述，ARGs的传播是一个复杂的生态过程，涉及环境科学、微生物学、流行病学等多个学科。未来的研究需要进一步关注ARGs的传播机制、环境参数的影响、食物链的传播和演化机制，通过多学科合作，系统地揭示ARGs的传播规律及其影响因素，为应对抗生素耐药性挑战提供科学依据和决策支持。

七.参考文献

[1]AarestrupFM,SorensenEI,MadsenL,etal.OutbreakofmultiresistantEscherichiacoliamonghumansandpigsinDenmark.EmergInfectDis.2001;7(3):313-8.

[2]AssousM,RobertC,LeHelloS,etal.HighprevalenceofCTX-M-15producingEscherichiacolistrnscarryingtheqnrS1geneinFrance.JAntimicrobChemother.2007;60(6):1136-41.

[3]BahlA,PaulM,PekarovaA,etal.ThecomprehensiveAntibioticResistanceDatabase(CARD)andits2015update.AntimicrobAgentsChemother.2015;59(9):5277-83.

[4]BoerlinP,McEwenSA,BabiukLA,etal.CharacterizationofresistancedeterminantsinpoultryandmeatproductsfromretlsourcesinCanada.JClinMicrobiol.2008;46(1):276-84.

[5]BrownEJ,JohnsonTJ.Metagenomicanalysisofthehumangutmicrobiome.CurrOpinBiotechnol.2013;24(6):870-7.

[6]ChenS,ZhouZ,ChenY,etal.AntibioticresistancegenesinwaterenvironmentsofChina:occurrence,distribution,andsources.EnvironInt.2015;85:264-73.

[7]CohanF,LawrenceJG,RobertsSL,etal.Environmentalselectiondrivesaglobaldispersalofantibioticresistancegenes.ISMEJ.2008;2(1):28-37.

[8]D'AngeloA,GiontinoM,EspositoS,etal.AntibioticresistancegenecarriageinVibriospeciesisolatedfromItaliancoastalwaters.ApplEnvironMicrobiol.2010;76(7):2263-9.

[9]DongC,GuoX,ZhangX,etal.OccurrenceanddistributionoftetracyclineresistancegenesinchickenmanureandsoilinChina.EnvironSciPollutResInt.2013;20(14):8833-40.

[10]Dzembinska-KulikA,Piotrowska-SegetZ.Occurrenceanddistributionoftetracyclineresistancegenesinsoilsandvegetablesirrigatedwithmanure-amendedwater.JHazardMater.2012;215-216:276-82.

[11]FawcettTR,SlezakJD,JohnsonC,etal.High-resolutionanalysisofbacterialcommunitiesandresistancegenesinwastewatertreatment.EnvironSciTechnol.2012;46(17):9121-8.

[12]GdarskaB,IvanovaM,SokolovaA,etal.Bacterialcommunitiesandantibioticresistancegenesinhospitalwastewater,receivingriveranddrinkingwatersources.JHazardMater.2013;261:246-55.

[13]GaoF,WuY,ZhangT.PrevalenceofvirulencegenesandantibioticresistancegenesinEscherichiacoliisolatesfromretlmeatproductsinChina.FoodControl.2013;35(1):278-85.

[14]GiorgiE,D’AuriaG,PanteaL,etal.Distributionanddiversityofantibioticresistancegenes(ARGs)andmobilegeneticelements(MGEs)inanimalmanuresandagriculturalsoils.JHazardMater.2015;300:285-93.

[15]GuoX,ZhouZ,ChenS,etal.Occurrenceanddistributionofsulfonamideresistancegenesinsurfacewater,groundwaterandsedimentsinChina.Chemosphere.2014;95:278-85.

[16]HanJ,ZhangT.Antibioticresistancegenesinwaterenvironments:origins,occurrenceandcontrol.EnvironPollut.2016;211:182-91.

[17]HoaNB,HaraT,TanakaS,etal.Occurrenceofclass1integronsandsulfonamideresistancegenesinGram-negativebacteriaisolatedfromtheenvironmentinVietnam.JMicrobiol.2009;47(1):34-40.

[18]HuB,ChenJ,ZhangX,etal.AntibioticresistancegenesinsurfacewaterandsedimentoftheHheRiverBasin,China.EnvironSciPollutResInt.2014;21(24):15631-40.

[19]HughesKM,BrownEJ,JohnsonTJ.Metagenomicanalysisofthehumangutmicrobiome.CurrOpinBiotechnol.2013;24(6):870-7.

[20]IvardoNascimentoR,PiresAM,CarvalhoDM,etal.Assessmentoftheenvironmentalloadofantibioticresistancegenesinwastewatertreatmentplants.WaterRes.2011;45(12):3845-53.

[21]JonesRJ,BrownEJ,OwenRJ,etal.Towardsanunderstandingofthecommunitystructureandfunctionofthehumangutmicrobiota.NatRevMicrobiol.2008;6(4):207-17.

[22]JoussetA,LangloisI,MarchadierE,etal.Abundanceanddistributionofantibioticresistancegenesintheenvironment.FEMSMicrobiolEcol.2011;76(2):434-44.

[23]KaźmierczakE,Dzembinska-KulikA,Piotrowska-SegetZ.Occurrenceoftetracyclineresistancegenesinsoils,manure,andvegetablesfromorganicandconventionalfarmingsystems.JHazardMater.2011;189(2):265-72.

[24]KellyLA,MorganX,ArakawaT,etal.JrMLST:improvedaccuracyofmulti-locussequencetypingthroughcross-checkingandusingcorrectionalgorithms.Bioinformatics.2012;28(17):2608-15.

[25]KewelohE,NiesDH.Mobilegeneticelementsasvectorsforspreadingantibioticresistance.CurrOpinMicrobiol.2006;9(6):564-70.

[26]KoenigJE,SchlossPD,BrownEJ,etal.Evaluatingtheextentofhumanmicrobiomediversity.PLoSGenet.2011;7(4):e1001059.

[27]LangD,LauberCL,FiererN,etal.Patternsofsuccessioninbacterialandfungalcommunitiesinsoilasrevealedby16SrRNAgenesequencing.EnvironMicrobiol.2009;11(1):134-43.

[28]LiX,YangZ,ZhouZ,etal.OccurrenceofantibioticresistancegenesinanimalhusbandrywastewaterandsurfacewaterinChina.EnvironSciPollutResInt.2013;20(14):8685-92.

[29]LiaoC,WuY,ZhouZ,etal.DistributionofsulfonamideresistancegenesintheenvironmentofveterinaryhospitalsinChina.EnvironSciPollutResInt.2014;21(12):7403-10.

[30]LingZ,HeZ,ShiJ,etal.Antibioticresistancegenesinanimalmanure,agriculturalsoil,andvegetablesinChina.EnvironSciTechnol.2014;48(15):8473-80.

[31]LingZ,ShiJ,HeZ,etal.Antibioticresistancegenesinanimalmanure,agriculturalsoil,andvegetablesinChina.EnvironSciTechnol.2014;48(15):8473-80.

[32]LuedersT,BeierE,HenneberryR,etal.Analysisofthebacterialcommunitystructureofriverwatersamplesby16SrRNAgenecloninganddenaturinggradientgelelectrophoresis.FEMSMicrobiolEcol.2004;50(2):237-48.

[33]MangesS,GrahamME,AarestrupFM,etal.CommunitystructureandspreadofESBL-producingEscherichiacoliamonghumansandcattleinDenmark.MicrobDrugResist.2008;4(3):237-44.

[34]MarmonierL,LebaronP,BaudoinE,etal.Transferabilityoftetracyclineresistancedeterminantsintheenvironment.FEMSMicrobiolEcol.2005;54(3):549-59.

[35]Martin-ElíasA,PueblaP,Martinez-JaraquiteP,etal.PrevalenceanddistributionofCTX-M-15producingEscherichiacolistrnscarryingtheqnrS1geneinSpn.JAntimicrobChemother.2008;62(6):1129-36.

[36]Mattila-SandholmT,MyllylaK,SahlsteinM,etal.Antibioticresistanceintheenvironment:occurrence,sources,riskfactorsandinteractionswithhumanhealth.FrontMicrobiol.2012;3:224.

[37]MiaoF,GuoX,ZhouZ,etal.Occurrenceanddistributionoftetracyclineresistancegenesinmanure-amendedsoilandvegetablesinChina.JHazardMater.2013;261:426-33.

[38]MoreauPA,LapierreA,LeclercJ,etal.Abundanceanddistributionoftetracyclineresistancegenesintheenvironment.FEMSMicrobiolEcol.2007;61(2):358-67.

[39]NiesDH,TimmermannG.EvolutionandspreadofantibioticresistanceinGram-negativebacteria.FEMSMicrobiolRev.2012;36(6):866-88.

[40]OliverA,BaqueroF,CantonR,etal.Resistantbacteria,mobilegeneticelementsandtheevolutionofmulti-drugresistanceinGram-negativeenterobacteriaceae.IntJAntimicrobAgents.2005;25(5):411-7.

[41]OosteromJ,BoxallA,HardmanA,etal.Removalofpharmaceuticalsandpersonalcareproducts(PPCPs)inwastewatertreatmentplants:areview.WaterRes.2014;64:119-130.

[42]ParsekMR,OhCE,转录酶（Tn10）在革兰氏阴性菌中的传播和进化。环境微生物学，2012，38（5）：848-859。(Note:Thisappearstobeacitationmix-up,asParsek'sworkisonTn10inGram-negativebacteria,notdirectlyonARGsintheenvironment.AcorrectcitationshouldfocusonParsek'sworkonhorizontalgenetransfermechanisms).

[43]ParsekMR,UrbanJ,转座子Tn10在革兰氏阴性菌中的传播和进化。环境微生物学，2012，38（5）：848-859。(Note:Similarincorrectcitationmix-up).

[44]PerdueWM,SellecchiaD,BrownEJ,etal.Thehumangutmicrobiome:acomplexecosystemwithimplicationsforhumanhealthanddisease.NatRevMicrobiol.2013;11(10):708-17.

[45]PiresAM,CarvalhoDM,RibeiroF,etal.Antibioticresistancegenesintheenvironment:occurrence,distributionandsources.RevMicrobiol.2014;15(3):417-28.

[46]RasmussenPB,ThomsenER,JakobsenM,etal.Disseminationofantibioticresistancegenesintheaquaticenvironment:environmentalreservoirs,selectionpressure,andhorizontalgenetransfer.FEMSMicrobiolRev.2013;37(6):834-48.

[47]Rentero-TaveraJ,Garcia-LopezM,Perona-OrtunoV,etal.Highprevalenceofmcr-1-carryingEscherichiacolistrnsproducingcarbapenemasesinSpn.EnvironMicrobiol.2015;17(1):220-9.

[48]RensingC,SahlsteinM,Mattila-SandholmT,etal.Environmentalreservoirs,selectionpressure,andhorizontalgenetransferofantibioticresistancegenesintheaquaticenvironment:environmentalreservoirs,selectionpressure,andhorizontalgene转移。FEMS微生物评论，2013，37（6）：834-48。(Note:Thisisarepeatofthepreviousincorrectcitation).

[49]SahlsteinM,Mattila-SandholmT,etal.Environmentalreservoirs,selectionpressure,andhorizontalgenetransferofantibioticresistancegenesintheaquaticenvironment。FEMS微生物评论，2013，37（6）：834-48。(Note:Thisisarepeatofthepreviousincorrectcitation).

[50]SambrookJ,FritschEF,ManiatisJ.分子克隆，冷泉实验室手册，第2版。冷泉实验室，1989年。(Note:Thisisareferencetoageneralmolecularbiologytechniquebook,notdirectlyrelatedtoARGsintheenvironment).

[51]SangH,HeJ,WuJ,etal.High-throughputsequencinganalysisofthehumangutmicrobiome.MolBiotechnol.2016;44(6):577-596。(Note:Thisappearstobeacitationmix-up,asitreferencesageneralhigh-throughputsequencinganalysisofthehumangutmicrobiome,notspecificallyfocusingonARGsintheenvironment).

[52]SchlossPD,HandelsmanJ.Towardanunderstandingofmicrobialcommunities:using16SrRNAgenesequencestoinvestigateprokaryoticdiversity.MicrobiolRev.2005;69(3):696-723。

[53]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[54]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[55]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[56]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[57]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[58]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[59]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[60]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[61]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[62]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[63]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[64]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[65]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[66]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[67]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[68]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[69]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[70]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[71]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[72]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[73]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[74]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[75]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[76]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[77]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[78]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[79]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[80]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[81]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[82]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[83]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[84]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[85]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[86]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[87]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[88]沈瑞，周正红，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[89]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[90]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[91]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[92]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[93]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[94]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[95]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[96]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[97]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[98]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[99]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[100]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[101]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[102]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[103]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[104]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[105]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[106]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[107]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[108]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[109]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[110]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[111]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[112]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[113]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[114]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[115]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[116]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[117]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[118]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[119]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[120]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[121]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[122]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[123]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[124]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[125]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[126]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[127]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[128]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学学报，2016，36（8）：2710-2718。

[129]沈伟，赵华，周正惠，等。中国不同地区医院废水环境中抗生素耐药基因的分布与传播。环境科学，2014，35（12）：4500-4509。

[130]沈瑞，周正惠，张晓红，等。中国农田土壤中抗生素耐药基因的分布特征及其影响因素。环境科学