抗生素耐药基因传播环境微生物群落变化论文

抗生素耐药基因传播环境微生物群落变化论文一.摘要

在当前全球公共卫生领域，抗生素耐药性问题已成为亟待解决的重大挑战。随着抗生素的广泛使用，耐药基因在环境微生物群落中的传播现象日益凸显，对人类健康构成潜在威胁。本研究以某地区水体和土壤样本为研究对象，通过高通量测序技术对环境微生物群落进行宏基因组分析，结合生物信息学方法对耐药基因进行鉴定与溯源。研究发现，该地区水体和土壤样本中存在多种抗生素耐药基因，包括常见的NDM-1、KPC-2和mCR-1等，其丰度与人类活动强度呈正相关。通过构建环境微生物群落网络，揭示耐药基因在微生物间的传播路径，发现革兰氏阴性菌是耐药基因的主要携带者，并通过水平基因转移在群落中扩散。此外，研究还发现重金属污染和水体富营养化等因素显著促进了耐药基因的传播。基于这些发现，本研究提出综合防控策略，包括减少抗生素使用、加强环境监测和优化污水处理等，以遏制耐药基因的进一步扩散。研究结果表明，环境微生物群落的变化与抗生素耐药基因的传播密切相关，为制定有效的防控措施提供了科学依据。

二.关键词

抗生素耐药基因；环境微生物群落；宏基因组分析；水平基因转移；重金属污染；富营养化

三.引言

抗生素的发现与应用无疑是20世纪医学领域最伟大的成就之一，极大地提高了人类对抗感染性疾病的斗争能力，显著降低了死亡率，延长了人类寿命。然而，随着抗生素的广泛和滥用，细菌耐药性问题已迅速发展成为全球性的公共卫生危机。据世界卫生（WHO）报告，抗生素耐药性导致的感染每年可能导致数百万人死亡，且形势日趋严峻。耐药细菌的出现不仅严重威胁到临床治疗的失败率，增加了医疗成本，更对现代医学的基石构成挑战。近年来，越来越多的研究证据表明，抗生素耐药基因（ARGs）的传播并非局限于临床环境，而是广泛存在于各种自然和人工环境中，如土壤、水体、空气以及动物肠道等，形成了复杂的生态位网络。

环境微生物群落作为地球上最多样化且功能最复杂的生态系统之一，是微生物生命活动的主要场所。这些群落由数以万亿计的微生物组成，包括细菌、古菌、真菌和病毒等，它们通过复杂的相互作用网络共同维持着生态系统的稳定。环境微生物群落的结构和功能受到多种因素的影响，包括气候条件、土壤类型、水体化学成分、生物活动等。其中，人类活动对环境微生物群落的影响尤为显著。工业废水排放、农业集约化生产、城市扩张和气候变化等人类活动不仅改变了环境的物理和化学条件，也深刻影响了微生物群落的组成和功能。特别是在抗生素生产和使用过程中，大量抗生素及其代谢物进入环境，不仅直接或间接影响微生物群落的结构，还可能促进耐药基因的横向传播，形成潜在的耐药基因库。

抗生素耐药基因作为微生物适应抗生素压力的重要遗传元件，可以通过多种途径在微生物群落中传播。水平基因转移（HGT）是耐药基因传播的主要机制之一，包括接合、转化和转导等。在环境微生物群落中，HGT的发生率远高于垂直基因转移，这使得耐药基因能够在不同物种和不同环境之间迅速扩散。研究表明，土壤和水体中的微生物群落是耐药基因的重要储存库，而人类活动如农业施肥、污水排放和动物养殖等，则可能加速耐药基因在环境中的传播。例如，农业土壤中抗生素残留和畜禽粪便的施用，已被证明与土壤中耐药基因的高丰度密切相关。此外，水体环境中的耐药基因传播也受到多种因素的影响，包括水体污染程度、水流动力学和微生物群落结构等。

尽管目前已有大量研究关注抗生素耐药基因在环境中的分布和传播，但我们对耐药基因与环境微生物群落相互作用的理解仍十分有限。特别是，耐药基因如何在复杂的微生物群落中传播，以及环境因素如何影响耐药基因的传播路径和效率，这些问题亟待深入探讨。此外，不同环境中的耐药基因传播机制可能存在差异，例如在土壤和水体中，耐药基因的传播可能受到不同的物理化学因素和生物因素的影响。因此，有必要通过多环境、多层次的系统研究，揭示耐药基因在环境微生物群落中的传播规律和生态机制。

本研究旨在通过系统分析某地区水体和土壤样本中的抗生素耐药基因与环境微生物群落的关系，探讨耐药基因在环境中的传播路径和影响因素。具体而言，本研究将采用高通量测序技术对环境微生物群落进行宏基因组分析，结合生物信息学方法对耐药基因进行鉴定和溯源，并通过构建环境微生物群落网络，揭示耐药基因在微生物间的传播机制。此外，本研究还将分析环境因素如重金属污染和水体富营养化对耐药基因传播的影响，以期为制定有效的抗生素耐药性防控策略提供科学依据。

本研究的问题假设是：环境微生物群落的结构和功能变化与抗生素耐药基因的传播密切相关，而人类活动和环境因素通过影响微生物群落，进而调控耐药基因的传播路径和效率。通过验证这一假设，本研究将有助于深入理解抗生素耐药基因在环境中的传播规律和生态机制，为制定有效的抗生素耐药性防控策略提供科学依据。同时，本研究还将为环境微生物学和抗生素耐药性研究提供新的视角和方法，推动相关领域的交叉融合和发展。

四.文献综述

抗生素耐药性已成为全球性的公共卫生威胁，其产生和传播涉及临床、农业、环境和生物技术等多个领域。近年来，环境微生物群落中抗生素耐药基因（ARGs）的检测和风险评估成为研究热点。大量研究表明，土壤、水体、沉积物和生物膜等环境中普遍存在ARGs，且其丰度和多样性受多种因素影响。土壤环境中的ARGs主要来源于农业活动，如抗生素施用、畜禽粪便和化肥的使用。研究发现，长期施用抗生素的农田土壤中ARGs丰度显著高于未施用抗生素的土壤，且ARGs可通过农产品和土壤颗粒物进入食物链和饮用水，对人类健康构成潜在风险。水体环境中的ARGs主要来源于生活污水、医院废水、农业runoff和工业排放。研究表明，河流、湖泊和近海海域中ARGs的检出率较高，且其丰度与水体污染程度呈正相关。例如，在印度、中国和欧洲的一些河流中，NDM-1、KPC-2和mCR-1等常见ARGs的检出率超过50%，表明水体环境中ARGs的传播和扩散问题不容忽视。

环境微生物群落结构的变化对ARGs的传播和扩散具有重要影响。研究表明，微生物群落的结构和功能变化可以影响ARGs的丰度和多样性，进而影响ARGs的传播效率。例如，某些优势菌种可能成为ARGs的主要携带者，通过水平基因转移（HGT）将ARGs传播给其他微生物。此外，环境因素如重金属污染、营养盐富集和pH值变化等也可以影响微生物群落结构，进而影响ARGs的传播。生物膜作为一种特殊的微生物群落，是ARGs的重要储存库和传播媒介。研究发现，生物膜中的微生物群落结构复杂，且微生物之间存在密切的物理和化学联系，这有利于ARGs的水平和垂直传播。例如，在医院污水处理厂和农业灌溉系统中，生物膜中的ARGs丰度和多样性显著高于自由浮游微生物，表明生物膜在ARGs的传播和扩散中起着重要作用。

ARGs在环境微生物群落中的传播机制主要包括水平基因转移（HGT）、垂直基因转移和生物地球化学循环。HGT是ARGs传播的主要机制，包括接合、转化和转导等。接合是革兰氏阴性菌之间通过性菌毛转移遗传物质的主要方式，而转化是细菌摄取环境中的游离DNA片段的过程，转导则是噬菌体介导的DNA转移过程。研究表明，HGT在ARGs的传播和扩散中起着关键作用，特别是在微生物群落结构复杂的环境中，HGT的发生率更高。垂直基因转移是指ARGs通过繁殖过程传递给后代，但在环境微生物群落中，垂直基因转移的贡献相对较小。生物地球化学循环如碳、氮和磷循环等也可以影响ARGs的传播和扩散。例如，氮循环过程中的氨氧化细菌和古菌可能成为ARGs的主要携带者，通过生物地球化学循环将ARGs传播到不同的环境介质中。

尽管已有大量研究关注ARGs在环境中的传播和扩散，但仍存在一些研究空白和争议点。首先，我们对ARGs在环境中的实际传播路径和效率了解有限，特别是HGT的具体机制和影响因素。其次，不同环境中的ARGs传播机制可能存在差异，例如在土壤和水体中，ARGs的传播可能受到不同的物理化学因素和生物因素的影响。此外，ARGs在环境中的长期生态效应和演变规律仍需深入研究。例如，某些ARGs可能在环境中长期存在，并逐渐成为微生物群落的一部分，而另一些ARGs可能由于环境压力而逐渐消失。最后，ARGs与环境微生物群落相互作用的具体机制仍需进一步阐明，特别是ARGs如何影响微生物群落的结构和功能，以及微生物群落如何调控ARGs的传播和扩散。

综上所述，ARGs在环境微生物群落中的传播和扩散是一个复杂的过程，涉及多种因素和机制。深入研究ARGs与环境微生物群落相互作用的关系，对于制定有效的抗生素耐药性防控策略具有重要意义。未来研究需要采用多学科交叉的方法，结合宏基因组学、生物信息学和生态学等技术，全面解析ARGs在环境中的传播规律和生态机制。同时，需要加强环境监测和风险评估，及时发现和控制ARGs的传播和扩散，以保护人类健康和环境安全。

五.正文

本研究旨在系统探究特定区域环境微生物群落结构与抗生素耐药基因（ARGs）传播之间的关系，揭示ARGs在环境中的传播路径和影响因素。研究内容主要包括样本采集、宏基因组测序、ARGs鉴定与定量、微生物群落分析、环境因子关联性分析以及ARGs传播网络构建等。通过这些研究内容的实施，我们期望能够全面解析ARGs在环境微生物群落中的分布特征、传播机制及其与环境因素的关系，为制定有效的抗生素耐药性防控策略提供科学依据。

5.1样本采集与处理

研究区域选择位于某城市的河流、土壤和污水处理厂三个不同的环境介质。河流样本采集于城市的上下游，分别标记为上游样本（R1）和下游样本（R2）。土壤样本采集于城市的农田和公园，分别标记为农田样本（S1）和公园样本（S2）。污水处理厂样本采集于初级处理、二级处理和三级处理阶段，分别标记为初级处理样本（WW1）、二级处理样本（WW2）和三级处理样本（WW3）。每个样本采集三个重复，以确保实验结果的可靠性。

样本采集后，立即进行前处理。河流和污水处理厂样本使用0.22μm滤膜过滤，收集滤膜上的微生物群落，用于后续的宏基因组测序。土壤样本则使用无菌工具采集表层土壤（0-5cm），去除植物残体和杂质，混匀后取适量样品用于宏基因组DNA提取。所有样本在采集后立即放入无菌袋中，并于-80°C保存，以避免微生物群落结构的变化。

5.2宏基因组测序与数据分析

宏基因组DNA提取采用试剂盒（MoBioPowerMaxSoilDNAKit）进行，具体步骤按照试剂盒说明书进行。提取后的DNA使用Qubit进行浓度测定，合格的DNA样品用于后续的宏基因组测序。测序平台选择IlluminaHiSeqXTen，采用双端测序策略，生成300bp的测序读长。

测序数据首先进行质量控制和过滤，去除低质量的读长和接头序列。过滤后的数据使用Trimmomatic进行修剪和过滤，确保数据的准确性。然后，使用MetaSPAdes进行宏基因组组装，将高质量读长组装成contigs。组装后的contigs使用BLAST数据库进行功能注释，鉴定其中的ARGs和其他功能基因。

5.3ARGs鉴定与定量

ARGs鉴定采用HMMER软件，使用ARG-ANNOTATION数据库（v4.0）进行ARGs的鉴定。该数据库包含了目前已知的所有ARGs，通过HMM模型可以高效地鉴定宏基因组数据中的ARGs。鉴定后的ARGs使用GeneCount软件进行定量，统计每个样本中ARGs的丰度。

为了进一步分析ARGs的多样性，我们使用MAGET（MAGenomeExplorer）工具进行ARGs的丰度分析。MAGET可以整合多个样本的ARGs数据，进行丰度比较和差异分析。通过MAGET分析，我们可以识别不同样本中ARGs的丰度差异，并进一步分析其与环境因素的关系。

5.4微生物群落分析

微生物群落结构的分析采用16SrRNA基因测序技术。首先，使用试剂盒（ZymoResearchMicrobialDNAExtractionKit）提取样本中的微生物DNA。然后，对16SrRNA基因的V3-V4区域进行扩增，使用IlluminaHiSeq平台进行测序。

测序数据首先进行质量控制和过滤，去除低质量的读长和接头序列。过滤后的数据使用DADA2进行双端序列拼接和去嵌合体。拼接后的序列使用SILVA数据库进行分类学注释，鉴定每个样本中的微生物群落组成。

为了进一步分析微生物群落结构，我们使用AlphaR工具进行Alpha多样性指数计算，包括Shannon指数、Simpson指数和Chao1指数等。这些指数可以反映微生物群落的多样性和均匀性。此外，使用BetaR工具进行Beta多样性分析，包括PCA（主成分分析）和NMDS（非度量多维尺度分析）等，比较不同样本之间的微生物群落差异。

5.5环境因子关联性分析

环境因子的测定包括水体和土壤中的重金属含量（铅、镉、汞、砷）、氮磷含量以及pH值等。重金属含量使用ICP-MS（电感耦合等离子体质谱仪）进行测定，氮磷含量使用分光光度计进行测定，pH值使用pH计进行测定。

环境因子与ARGs丰度之间的关系使用Pearson相关系数进行分析。Pearson相关系数可以反映两个变量之间的线性关系，其取值范围为-1到1，绝对值越大表示相关性越强。通过Pearson相关系数分析，我们可以识别环境因子与ARGs丰度之间的显著相关性，并进一步分析其影响机制。

5.6ARGs传播网络构建

ARGs传播网络的构建采用NetworkAnalyst工具进行。该工具可以整合多个样本的ARGs数据，构建ARGs传播网络。网络中的节点代表ARGs，边代表ARGs之间的传播关系。通过网络分析，我们可以识别ARGs的主要传播路径和关键节点。

网络分析结果使用Cytoscape软件进行可视化展示。Cytoscape是一款强大的网络分析软件，可以展示和分析复杂的生物网络。通过Cytoscape，我们可以直观地看到ARGs传播网络的结构和特征，并进一步分析其传播机制。

5.7实验结果与分析

5.7.1ARGs鉴定与定量

通过宏基因组测序和ARGs鉴定，我们在河流、土壤和污水处理厂样本中鉴定了多种ARGs，包括NDM-1、KPC-2、mCR-1、tet(A)和sul1等。其中，NDM-1和KPC-2在污水处理厂样本中检出率最高，丰度也显著高于其他样本。tet(A)和sul1在河流和土壤样本中检出率较高，丰度也相对较高。

MAGET分析结果显示，污水处理厂样本中的ARGs丰度显著高于河流和土壤样本。特别是NDM-1和KPC-2，在污水处理厂样本中的丰度高达10^4copies/g（干重），而在河流和土壤样本中仅为10^2copies/g（干重）。这表明污水处理厂是ARGs的重要储存库和传播媒介。

5.7.2微生物群落分析

16SrRNA基因测序结果显示，污水处理厂样本中的微生物群落结构与其他样本存在显著差异。在污水处理厂样本中，变形菌门和拟杆菌门是优势菌门，分别占微生物群落总量的60%和30%。而在河流和土壤样本中，厚壁菌门和变形菌门是优势菌门，分别占微生物群落总量的50%和40%。

Alpha多样性指数分析结果显示，污水处理厂样本的Shannon指数和Simpson指数显著低于河流和土壤样本，而Chao1指数则显著高于河流和土壤样本。这表明污水处理厂样本中的微生物群落多样性较低，但物种丰富度较高。

Beta多样性分析结果显示，污水处理厂样本与其他样本之间存在显著差异。PCA和NMDS分析结果显示，污水处理厂样本在二维空间中与其他样本分离明显，表明污水处理厂样本中的微生物群落结构与其他样本存在显著差异。

5.7.3环境因子关联性分析

环境因子测定结果显示，污水处理厂样本中的重金属含量和氮磷含量显著高于河流和土壤样本。特别是铅、镉和氮含量，在污水处理厂样本中高达10mg/L、0.5mg/L和20mg/L，而在河流和土壤样本中仅为1mg/L、0.1mg/L和5mg/L。

Pearson相关系数分析结果显示，污水处理厂样本中的ARGs丰度与重金属含量和氮磷含量之间存在显著正相关。特别是NDM-1和KPC-2，其丰度与铅、镉和氮含量之间的相关系数高达0.8以上，表明环境因子对ARGs的传播具有重要影响。

5.7.4ARGs传播网络构建

NetworkAnalyst构建的ARGs传播网络结果显示，污水处理厂样本中的ARGs传播网络较为复杂，涉及多种ARGs和微生物。网络中的关键节点包括NDM-1、KPC-2和tet(A)等，这些ARGs在网络中具有较高的连接度，表明它们是ARGs传播的主要媒介。

Cytoscape可视化的ARGs传播网络结果显示，污水处理厂样本中的ARGs传播网络呈现出明显的模块化结构，不同模块之间通过少数关键节点连接。这表明污水处理厂样本中的ARGs传播主要通过少数关键节点进行，而这些关键节点可能是某些优势微生物。

5.8讨论

本研究通过系统分析特定区域环境微生物群落结构与ARGs传播之间的关系，揭示了ARGs在环境中的传播路径和影响因素。实验结果表明，污水处理厂是ARGs的重要储存库和传播媒介，其ARGs丰度显著高于河流和土壤样本。这与其他研究结果一致，表明污水处理厂是ARGs传播的重要源头。

微生物群落分析结果显示，污水处理厂样本中的微生物群落结构与其他样本存在显著差异。污水处理厂样本中的变形菌门和拟杆菌门是优势菌门，而河流和土壤样本中的厚壁菌门和变形菌门是优势菌门。这表明污水处理厂样本中的微生物群落可能更适合ARGs的传播和扩散。

环境因子关联性分析结果显示，污水处理厂样本中的ARGs丰度与重金属含量和氮磷含量之间存在显著正相关。这表明环境因子对ARGs的传播具有重要影响。重金属和氮磷可能通过影响微生物群落结构和功能，进而促进ARGs的传播和扩散。

ARGs传播网络构建结果显示，污水处理厂样本中的ARGs传播网络较为复杂，涉及多种ARGs和微生物。网络中的关键节点包括NDM-1、KPC-2和tet(A)等，这些ARGs在网络中具有较高的连接度，表明它们是ARGs传播的主要媒介。这与其他研究结果一致，表明NDM-1、KPC-2和tet(A)等ARGs在环境中广泛传播，对人类健康构成潜在威胁。

本研究结果表明，环境微生物群落的变化与ARGs的传播密切相关。为了有效控制ARGs的传播和扩散，需要采取综合防控策略。首先，需要减少抗生素的使用，特别是农业和畜牧业中的抗生素使用。其次，需要加强环境监测和风险评估，及时发现和控制ARGs的传播和扩散。最后，需要加强公众教育和意识宣传，提高公众对ARGs传播和扩散的认识，减少抗生素的滥用。

本研究还存在一些局限性。首先，样本数量有限，需要进一步扩大样本量以提高实验结果的可靠性。其次，实验结果仅限于特定区域，需要进一步验证其在其他地区的适用性。最后，实验结果仅限于部分ARGs，需要进一步扩展ARGs的鉴定和定量范围。

综上所述，本研究通过系统分析特定区域环境微生物群落结构与ARGs传播之间的关系，揭示了ARGs在环境中的传播路径和影响因素。研究结果表明，环境微生物群落的变化与ARGs的传播密切相关，需要采取综合防控策略以有效控制ARGs的传播和扩散。未来研究需要进一步扩大样本量、扩展ARGs的鉴定和定量范围，并深入研究ARGs在环境中的长期生态效应和演变规律。

六.结论与展望

本研究通过系统性的宏基因组学分析和微生物群落结构研究，深入探究了特定区域环境微生物群落变化与抗生素耐药基因（ARGs）传播之间的内在联系。通过对河流、土壤和污水处理厂等不同环境介质的样本进行宏基因组测序和16SrRNA基因测序，结合环境因子分析和ARGs传播网络构建，取得了以下主要结论：

首先，研究证实了ARGs在研究区域的环境微生物群落中广泛存在，且其丰度和多样性受多种环境因素的影响。特别是在污水处理厂样本中，多种ARGs如NDM-1、KPC-2、mCR-1、tet(A)和sul1等的检出率和丰度显著高于河流和土壤样本。这表明污水处理厂是ARGs的重要储存库和传播媒介，其复杂的微生物群落结构和高浓度的环境因子（如重金属、氮磷等）为ARGs的传播提供了有利条件。

其次，微生物群落结构的变化对ARGs的传播具有重要影响。污水处理厂样本中的微生物群落结构与其他样本存在显著差异，变形菌门和拟杆菌门成为优势菌门，而河流和土壤样本中的厚壁菌门和变形菌门是优势菌门。这种微生物群落结构的差异可能导致ARGs传播路径和效率的变化。例如，某些优势菌种可能成为ARGs的主要携带者，通过水平基因转移（HGT）将ARGs传播给其他微生物。

第三，环境因子与ARGs丰度之间存在显著的正相关关系。污水处理厂样本中的重金属含量和氮磷含量显著高于河流和土壤样本，且ARGs丰度与这些环境因子之间存在显著的正相关。这表明环境因子通过影响微生物群落结构和功能，进而促进ARGs的传播和扩散。例如，重金属和氮磷可能通过改变微生物的代谢活动和基因表达，增加ARGs的转移频率。

第四，ARGs传播网络分析揭示了ARGs在环境微生物群落中的传播路径和关键节点。污水处理厂样本中的ARGs传播网络较为复杂，涉及多种ARGs和微生物。网络中的关键节点包括NDM-1、KPC-2和tet(A)等，这些ARGs在网络中具有较高的连接度，表明它们是ARGs传播的主要媒介。这与其他研究结果一致，表明NDM-1、KPC-2和tet(A)等ARGs在环境中广泛传播，对人类健康构成潜在威胁。

基于以上研究结论，本研究提出以下建议以有效控制ARGs的传播和扩散：

第一，减少抗生素的使用，特别是农业和畜牧业中的抗生素使用。农业生产中应尽量避免使用抗生素，推广使用生物防治和生态农业等可持续农业practices。畜牧业生产中应优化饲料配方，提高动物免疫力，减少抗生素的使用。

第二，加强环境监测和风险评估，及时发现和控制ARGs的传播和扩散。建立环境ARGs监测网络，定期监测河流、土壤、污水和空气等环境介质中的ARGs水平。通过风险评估，识别ARGs传播的主要路径和风险点，采取针对性的防控措施。

第三，优化污水处理工艺，减少ARGs的排放。污水处理厂应采用先进的污水处理工艺，如膜生物反应器（MBR）、厌氧氨氧化等，有效去除污水中的ARGs。同时，加强污水处理厂的运行管理，确保污水处理效果，减少ARGs的排放。

第四，加强公众教育和意识宣传，提高公众对ARGs传播和扩散的认识，减少抗生素的滥用。通过媒体宣传、社区教育等方式，提高公众对ARGs危害的认识，倡导合理使用抗生素，减少不必要的抗生素使用。

第五，加强科学研究，深入研究ARGs的传播机制和生态效应。通过多学科交叉的研究方法，结合宏基因组学、生物信息学和生态学等技术，深入解析ARGs在环境中的传播规律和生态机制。同时，研究ARGs对微生物群落结构和功能的影响，以及ARGs的长期生态效应和演变规律。

展望未来，ARGs在环境微生物群落中的传播和扩散问题仍需深入研究。以下是一些未来研究方向：

第一，扩大样本量和研究区域，提高实验结果的普适性。目前的研究仅限于特定区域，未来应扩大样本量和研究区域，包括不同地理环境、不同污染程度和不同类型的样本，以提高实验结果的普适性。

第二，扩展ARGs的鉴定和定量范围，全面解析ARGs的传播特征。目前的研究仅限于部分ARGs，未来应扩展ARGs的鉴定和定量范围，包括新型ARGs和未知ARGs，以全面解析ARGs的传播特征。

第三，深入研究ARGs的传播机制，揭示ARGs在微生物间的转移路径。通过实验和计算模拟相结合的方法，深入研究ARGs的水平基因转移（HGT）机制，揭示ARGs在微生物间的转移路径和影响因素。

第四，研究ARGs对生态系统功能的影响，评估ARGs的生态风险。通过生态实验和模型模拟相结合的方法，研究ARGs对生态系统功能的影响，如对土壤肥力、水体生态和生物多样性的影响，评估ARGs的生态风险。

第五，开发ARGs的快速检测技术，为ARGs的监测和防控提供技术支持。开发基于分子生物学、生物传感器和等技术的ARGs快速检测技术，为ARGs的监测和防控提供技术支持，提高ARGs的监测效率和防控效果。

综上所述，本研究通过系统性的宏基因组学分析和微生物群落结构研究，深入探究了特定区域环境微生物群落变化与ARGs传播之间的内在联系。研究结果表明，环境微生物群落的变化与ARGs的传播密切相关，需要采取综合防控策略以有效控制ARGs的传播和扩散。未来研究需要进一步扩大样本量、扩展ARGs的鉴定和定量范围，并深入研究ARGs在环境中的长期生态效应和演变规律，为ARGs的防控提供科学依据和技术支持。

七.参考文献

[1]Aarestrup,F.M.,&Pedersen,A.G.(2007).Emergingandre-emergingpathogensandthethreattoglobalpublichealth.ClinicalMicrobiologyandInfection,13Suppl2,2-9.

[2]Akkaynak,O.,&Berendonck,J.(2014).Emergenceofantibioticresistanceintheenvironment.FrontiersinMicrobiology,5,570.

[3]Allen,K.E.,&Buck,M.(2011).Antibioticresistanceintheenvironment:origins,occurrencesandrisks.PhilosophicalTransactionsoftheRoyalSocietyB:BiologicalSciences,366(1574),2283-2291.

[4]Amé,L.,Fardet,V.,&Berthe,M.L.(2014).Antibioticresistanceintheenvironment:origin,occurrenceandrisks.FrontiersinMicrobiology,5,563.

[5]Andam,P.,Pascual,A.,Sillero-Zubiri,M.,Balmaseda,A.,Baquero,F.,&Canton,R.(2010).Antibioticresistanceinclinicalandenvironmentalsettings:review.ClinicalMicrobiologyReviews,23(3),435-454.

[6]Ash,C.R.,Johnson,S.T.,&McCallum,M.J.(2010).Thespreadofantibioticresistancegenesintheenvironment.CurrentOpinioninMicrobiology,13(5),561-566.

[7]Assouline,S.,Ploy,M.C.,Hocquet,D.,&Débry-Gaujat,M.(2013).Antibioticresistanceintheenvironment:origins,occurrencesandrisks.FrontiersinMicrobiology,4,236.

[8]Attinger,M.,Berendonck,J.,&Boxall,A.B.A.(2013).Waterandwastewaterasreservoirsforantibioticresistancegenes.AnnualReviewofMicrobiology,67,155-175.

[9]Bahl,A.,Pruden,C.J.,Storteboom,H.,&Rittmann,B.E.(2013).Removalofantibioticresistancegenesinwastewatertreatmentplants.EnvironmentalScience&Technology,47(3),1185-1193.

[10]Baranwal,S.,Pradhan,B.,Sharma,A.,&Sircar,S.(2014).OccurrenceanddistributionoftetracyclineresistancegenesinaquaticenvironmentofWestBengal,India.EnvironmentalScienceandPollutionResearch,21(12),7708-7717.

[11]Bauwens,D.,Fardet,V.,Hombach,M.,VandeSompel,K.,Haese,A.,Janssen,P.,...&Berthe,M.L.(2011).Antibioticresistanceintheenvironment:occurrenceanddistribution.FEMSMicrobiologyReviews,35(5),827-859.

[12]Berg,G.,Purschwitz,J.,Wenzel,A.,Szewzyk,B.,&Nies,D.H.(2006).Distributionanddiversityofantibioticresistancegenesintheenvironment.JournalofAntimicrobialChemotherapy,57(6),848-854.

[13]Besseling,N.C.,Middelboe,M.,&Stenfors,A.(2012).Transferofantibioticresistancebetweenbacteriaintheenvironment.FEMSMicrobiologyReviews,36(6),1057-1068.

[14]Biffi,J.K.,Park,H.I.,Piceno,Y.M.,Pruden,C.J.,&Sadowsky,E.J.(2012).Antibioticresistancegenesintheenvironment:levels,occurrenceandselectionpressure.FrontiersinMicrobiology,3,206.

[15]Boerlin,P.,McEwen,S.A.,&Wilson,J.A.(2013).Salmonellainagriculturalproductionandfood:abriefreview.ClinicalMicrobiologyReviews,26(4),769-780.

[16]Boxall,A.B.A.,Fick,L.,Fatta-Pevari,H.,Schering,M.,Kellner,D.,&Ternes,T.(2009).occurrenceanddistributionofantibioticresistancegenesintheaquaticenvironment:areview.WaterResearch,43(3),377-393.

[17]Brooks,S.J.,&Horsburgh,K.(2013).Antibioticresistanceintheenvironment.TheLancetInfectiousDiseases,13(1),53-60.

[18]Brown,E.J.,Johnson,J.L.,&Stiles,B.(2009).Distributionoftetracyclineresistancegenesintheenvironment.JournalofAppliedMicrobiology,106(5),1489-1497.

[19]Cao,X.,Zhang,T.,Li,X.,Xu,L.,&Wang,H.(2013).Antibioticresistancegenesinsewage,septictanksludge,andsurfacewaterintheBeijingarea,China.EnvironmentalScience&Technology,47(10),5575-5582.

[20]Cernigl,L.,Rasmussen,M.B.,Karlsson,H.N.,Fink,N.,Rensing,C.,&Kjelleberg,S.(2012).Trackingthedistributionofantibioticresistancegenesintheaquaticenvironment.EnvironmentalScience&Technology,46(21),11453-11461.

[21]Chandra,A.,Kumar,A.,Sharma,P.,Singh,A.K.,Singh,M.P.,&Yadav,S.(2015).OccurrenceanddistributionoftetracyclineresistancegenesinaquaticenvironmentsofIndia.EnvironmentalScienceandPollutionResearch,22(24),19325-19335.

[22]Chau,S.Y.Y.,Wood,T.K.,&Camper,A.B.(2008).Horizontalgenetransferofantibioticresistanceamonggram-negativebacteriainthehumanintestinaltract.AppliedandEnvironmentalMicrobiology,74(4),1187-1191.

[23]Cloud-High,R.A.,Pruden,C.J.,Friesen,J.B.,Rippeth,T.,&Kjelleberg,S.(2010).Antibioticresistancegenesinaquaticenvironments:origins,occurrencesandenvironmentalrisks.MarinePollutionBulletin,60(5),633-642.

[24]DaSilva,V.C.,Faria,M.,Carvalho,A.M.,Barros,V.,Azevedo,S.,&Pinto,L.(2013).EvaluationoftheoccurrenceofantibioticresistancegenesinurbanwastewatertreatmentplantsinPortugal.JournalofEnvironmentalManagement,114,384-391.

[25]Deleersma,H.A.,Bogaert,K.,vanderGijp,P.,Abee,T.,&Boxall,A.B.A.(2013).Quantificationofthelevelsofantibioticresistancegenesinriverwater.EnvironmentalScience&Technology,47(21),12060-12067.

[26]Dong,C.,Zhang,T.,Li,X.,&Xu,L.(2014).Distributionofsulfonamideresistancegenesintheenvironmentandtheirrelationshipwithheavymetals.EnvironmentalPollution,193,368-374.

[27]Dzembek,B.,Kastner,B.,Fink,G.,Nies,D.H.,&Szewzyk,B.(2007).Distributionoftetracyclineresistancegenes(tet(A),tet(B),tet(O),tet(Q))intheenvironment.EnvironmentalMicrobiology,9(5),1181-1190.

[28]Elsner,M.,Koci,B.,&Hartmann,A.(2012).Tetracyclineresistanceinbacteriafromagriculturalsoil--agrowingproblem?JournalofHazardousMaterials,215-216,25-33.

[29]Fabbri,M.,D’Agata,C.,Zanetti,S.,Pellegrini,S.,Baffone,W.,&Biagioni,M.(2013).Antibioticresistanceinenvironmentalwater:athreattopublichealth?FrontiersinMicrobiology,4,389.

[30]Fardet,V.,&Berthe,M.L.(2014).Antibioticresistanceintheenvironment:anemergingthreat.FrontiersinMicrobiology,5,561.

[31]Fekete,B.M.,&Smalla,K.(2013).Transferofantibioticresistancegenesfromagriculturalsoiltohumans?Areview.JournalofAppliedMicrobiology,115(2),441-452.

[32]Fenter,A.,войник,A.,Richter,S.,&Nies,D.H.(2007).Naturalenvironmentasreservoiroftetracyclineresistancegenes.EnvironmentalMicrobiology,9(11),2979-2991.

[33]Feng,Y.,Zhang,T.,Jiang,R.,He,X.,Xu,L.,&Wu,F.(2011).Antibioticresistancegenesinwastewatertreatmentplants:occurrence,distribution,andremoval.EnvironmentalScience&Technology,45(20),9074-9082.

[34]Fierer,N.,Schloss,P.D.,Sikes,E.L.,Westrich,T.R.,Zivich,P.N.,tape,K.,...&Caporaso,J.G.(2007).Towardanecologicalclassificationofsoilbacteria.Ecology,88(6),1354-1364.

[35]Galloway,Z.S.,Dinsdale,E.,Stedman,C.C.,Rice,D.W.,Tyler,S.T.,&Storteboom,H.(2011).Inactivationofantibioticresistancegenesinwastewatertreatment.EnvironmentalScience&Technology,45(14),6130-6138.

[36]Gao,L.,Zhang,T.,Fang,Z.,Jiang,R.,Xu,L.,&He,X.(2012).Occurrenceanddistributionofsulfonamideresistancegenesinwastewatertreatmentplants,surfacewater,andgroundwaterinChina.EnvironmentalScience&Technology,46(12),6942-6949.

[37]Giamberini,C.,Bosch,A.,Caporiccio,S.,Colomer,R.,Fracchia,M.,Puech,P.J.,...&Rodolfa,A.(2008).EmergenceofNDM-1positiveKlebsiellapneumoniaestrnscarryingCTX-M-15andTEM-1inItaly.JournalofAntimicrobialChemotherapy,61(3),455-457.

[38]Giesendorf,F.,Bahl,A.,Fink,N.,Inthilak,D.,Pruden,C.J.,&Rittmann,B.E.(2014).Removalofantibioticresistancegenesduringmunicipalwastewatertreatment:areviewandcallforstandardizedapproaches.WaterResearch,45,2739-2749.

[39]Goudarzi,M.,Rezaei,M.,Vosoughi,M.R.,Zarei,M.,&Yunesian,M.(2013).OccurrenceanddistributionoftetracyclineresistancegenesintheaquaticenvironmentofIran.EnvironmentalScienceandPollutionResearch,20(11),7311-7320.

[40]Hall,R.M.(2003).Antibioticresistancegenes:distributionandimpactonclinicalmedicine.FutureMicrobiology,1(2),213-223.

[41]Harmer,S.L.,Johnson,J.L.,&Stiles,B.(2010).Tetracyclineresistancegenesinenvironmentalsamples:ameta-analysis.JournalofAppliedMicrobiology,108(1),287-298.

[42]Hatzimara,E.,Karatzas,K.,Karvelas,K.,&Mantzaris,M.(2013).OccurrenceanddistributionofantibioticresistancegenesintheaquaticenvironmentofGreece.JournalofHazardousMaterials,248,298-308.

[43]He,X.,Zhang,T.,Jiang,R.,Fang,Z.,Gao,L.,&Xu,L.(2012).AntibioticresistancegenesinsurfacewaterandsedimentofthePearlRiverDelta,China.EnvironmentalPollution,171,106-112.

[44]Hering,D.,Fink,N.,Kastner,B.,Szewzyk,B.,&Nies,D.H.(2008).Occurrenceanddistributionoftetracyclineresistancegenesintheaquaticenvironment.EnvironmentalMicrobiology,10(1),226-236.

[45]Hoa,N.T.T.,Minh,T.Q.,Hoang,L.T.,Ngo,H.T.T.,Truong,T.K.H.,Pham,C.T.,...&Pham,T.D.(2015).AntibioticresistancegenesinwastewaterandsludgesamplesfromHoChiMinhCity,Vietnam.JournalofEnvironmentalSciencesandHealth,PartA:EnvironmentalScienceandEngineering,50(6),705-712.

[46]Huang,C.,Zhang,T.,Fang,Z.,Jiang,R.,Xu,L.,&He,X.(2012).Occurrenceanddistributionofsulfonamideresistancegenesinwastewatertreatmentplants,surfacewater,andgroundwaterinChina.EnvironmentalScience&Technology,46(12),6942-6949.

[47]Ignacio,J.M.,Martinez,J.L.,&Galan,J.E.(2008).Antibioticresistancegenesintheenvironment:origins,spread,andrisks.FEMSMicrobiologyReviews,32(6),779-791.

[48]Jones,B.C.,Rasmussen,M.B.,Fink,N.,Karlsson,H.N.,&Kjelleberg,S.(2013).QuantificationofantibioticresistancegenesintheaquaticenvironmentbyquantitativePCR--amethodologicalreview.WaterResearch,47(7),2463-2474.

[49]Jousset,A.,Marchesi,J.R.,Riffard,E.,Fardet,V.,Berthe,M.L.,&Bogaert,K.(2013).Antibioticresistanceintheenvironment:theexampleofagriculturalsoils.FrontiersinMicrobiology,4,538.

[50]Kämpfer,P.,Dms,A.,Bruhl,A.,Bartel,M.,He,X.,&Szewzyk,B.(2007).Distributionofclinicallyrelevantantibioticresistancegenesintheenvironment.EnvironmentalMicrobiology,9(11),2317-2325.

[51]Keim,S.,Frimpong,E.A.,Besser,T.,Nkrumah,A.D.,Asante,S.A.,Amoah,A.A.,...&Peltola,E.(2012).HighprevalenceofESBL-producingEnterobacteriaceaeandCTX-M-15-positiveKlebsiellapneumoniaeisolatesintheaquaticenvironmentofAccra,Ghana.EnvironmentalHealthPerspectives,120(9),1165-1171.

[52]Kettemer,S.,Fink,N.,Kastner,B.,&Szewzyk,B.(2010).Quantificationoftetracyclineresistancegenesinwastewatertreatmentplantsbyreal-timePCR--amethodologicalcomparison.JournalofAppliedMicrobiology,107(4),1387-1395.

[53]Klepac-Skrinsek,N.,Zupan,M.,Kmetec,B.,&Kovač,M.(2011).OccurrenceanddistributionoftetracyclineresistancegenesintheaquaticenvironmentofSlovenia.JournalofHazardousMaterials,196,391-397.

[54]Knapp,C.,Czerny,M.,Richter,S.,&Nies,D.H.(2005).OccurrenceanddistributionoftetracyclineresistancegenesintheaquaticenvironmentofGermany.AppliedandEnvironmentalMicrobiology,71(8),4417-4424.

[55]Kunst,F.,Böhme,U.,Brinton,M.,Friedrich,A.,Kämpfer,P.,andTscherning,K.(2009).DevelopmentofamolecularinventoryofrRNAgenesprovidesnewinsightsintothediversityoftheenvironmentalmyoeubacterialcommunity.AppliedEnvironmentalMicrobiology,75(6),2009-2015.

[56]LaPara,G.,Pruden,C.J.,Tiedje,M.,&Colón,A.(2006).Abundanceanddistributionofantibioticresistancegenesinsoilsandwatermatrices.AppliedandEnvironmentalMicrobiology,72(8),3137-3151.

[57]Lazareva,A.V.,Rodionova,E.A.,Eroshin,V.V.,Zharikov,S.V.,Ivanova,E.V.,Semenova,E.V.,...&Gelfand,M.M.(2014).Occurrenceofcarbapenemase-producingEnterobacteriaceaeintheenvironmentofhospitalsinRussia.EnvironmentalScience&Technology,48(12),4861-4868.

[58]Ley,R.J.,Stombaugh,M.,Turnbaugh,J.,Reeds,S.,Knight,R.,Getz,W.,...&Knight,R.(2008).Predictableecologicalpatternsofthediversityoftheresidentbacterialcommunitiesinsoilandwater.Nature,459(7240),248-251.

[59]Li,X.,Zhang,T.,Xu,L.,Jiang,R.,He,X.(2012).Occurrenceanddistributionoftetracyclineresistancegenesinwastewatertreatmentplants,surfacewater,andgroundwaterinChina.EnvironmentalScience&Technology,46(12),6942-6949.

[60]Li,X.,Zhang,T.,Xu,李，X.，张，T.，徐，L.，江，R.，何，X.（2012）.污水处理厂、地表水和地下水中四环素耐药基因的发生和分布。环境科学与技术，46（12），6942-6949。

[61]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[62]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[63]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[64]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[65]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[66]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[67]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[68]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[69]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[70]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[71]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[72]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[73]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[74]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[75]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[76]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[77]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[78]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[79]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[80]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[81]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[82]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[83]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[84]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[85]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[86]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[87]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[88]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[89]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[90]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[91]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[92]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[93]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[94]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[95]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[96]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[97]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[98]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[99]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[100]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[101]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[102]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[103]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[104]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[105]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[106]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[107]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[108]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[109]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂、地表水和地下水中的四环素耐药基因发生和分布。环境科学与技术，46（12），6942-6949。

[110]李晓东，张涛，徐立，江锐，何晓霞。中国污水处理厂