抗生素耐药基因传播宏基因组学研究论文

抗生素耐药基因传播宏基因组学研究论文一.摘要

抗生素耐药性已成为全球公共卫生面临的严峻挑战，其核心驱动因素之一是耐药基因的广泛传播。本研究聚焦于某区域环境样本的宏基因组学分析，旨在揭示环境中抗生素耐药基因（ARGs）的群落结构、传播途径及其潜在风险。研究选取了包括土壤、水体和沉积物在内的多个环境样本，采用高通量测序技术对样本中的ARGs进行深度测序和生物信息学分析。通过构建ARGs的丰度图谱，结合环境参数（如pH值、有机质含量和重金属浓度）进行相关性分析，我们发现特定ARGs（如NDM-1、mcr-1和blaCTX-M）在重金属污染区域呈现显著富集现象。此外，通过网络分析揭示了ARGs之间及与移动遗传元件（MGEs）的相互作用，表明MGEs在ARGs的传播中扮演了关键角色。研究还发现，环境中存在多种潜在的ARGs传播媒介，如蚯蚓、昆虫和植物根系，这些媒介可能通过生物气溶胶或土壤颗粒转移ARGs。综合分析表明，环境中的ARGs传播呈现出复杂的时空分布特征，并与人类活动、环境压力和生物多样性密切相关。本研究为理解ARGs的传播机制提供了新的视角，并为制定有效的环境干预措施提供了科学依据，以遏制抗生素耐药性的进一步扩散。

二.关键词

抗生素耐药基因；宏基因组学；环境样本；移动遗传元件；传播途径；公共卫生

三.引言

抗生素的发现与应用曾是现代医学的里程碑，极大地降低了感染性疾病导致的死亡率，显著提升了人类平均寿命。然而，随着抗生素的广泛和滥用，抗生素耐药性（AntibioticResistance,AMR）问题日益严峻，已成为全球性的公共卫生危机。据世界卫生组织（WHO）报告，每年约有70万人死于耐药菌感染，如果不采取有效措施，到2050年，这一数字可能增至1000万。抗生素耐药性的产生涉及多种因素，包括细菌自发突变、水平基因转移（HorizontalGeneTransfer,HGT）以及抗生素的不合理使用。其中，抗生素耐药基因（AntibioticResistanceGenes,ARGs）作为耐药性的功能单元，通过HGT在细菌间快速传播，是导致耐药性蔓延的关键因素。

ARGs的传播途径复杂多样，主要包括直接接触传播、医疗器械污染、人类排泄物排放以及环境介导的传播。近年来，越来越多的研究表明，环境（如土壤、水体、沉积物和空气）已成为ARGs的重要储存库和传播媒介。环境样本中的ARGs可能来源于农业活动（如动物粪便和抗生素残留）、城市污水排放、工业废水泄漏以及医疗废弃物处理不当等。这些ARGs可通过多种途径进入人类和动物体内，例如饮用水污染、食物链富集以及直接接触受污染环境。此外，环境中存在的移动遗传元件（MobileGeneticElements,MGEs），如质粒、转座子和整合子，能够携带ARGs在不同物种间转移，进一步加速了ARGs的扩散。

尽管已有大量研究关注临床样本中的ARGs，但环境中ARGs的群落结构、传播机制及其对公共卫生的影响仍需深入探究。特别是，不同环境样本中ARGs的丰度、多样性及其与环境参数（如重金属含量、pH值和有机质含量）的关系，以及MGEs在ARGs传播中的作用，尚缺乏系统性的研究。此外，环境中ARGs的传播媒介（如昆虫、植物根系和土壤生物）及其介导的传播途径也亟待阐明。这些问题的解决不仅有助于揭示ARGs的生态学特征，还为制定有效的环境干预措施提供了科学依据。

基于上述背景，本研究旨在通过宏基因组学技术，系统分析某区域环境样本中的ARGs群落结构、传播途径及其与环境参数的关系，并探究MGEs在ARGs传播中的作用。具体而言，本研究提出以下假设：1）环境中存在多种ARGs，其丰度和多样性与环境参数密切相关；2）MGEs是ARGs传播的关键媒介；3）环境中存在多种潜在的ARGs传播媒介，如昆虫、植物根系和土壤生物。通过验证这些假设，本研究将有助于揭示ARGs的传播机制，并为制定有效的环境管理策略提供科学支持。

本研究选取了包括土壤、水体和沉积物在内的多个环境样本，采用高通量测序技术对样本中的ARGs进行深度测序和生物信息学分析。通过构建ARGs的丰度图谱，结合环境参数进行相关性分析，我们发现特定ARGs在重金属污染区域呈现显著富集现象。此外，通过网络分析揭示了ARGs之间及与MGEs的相互作用，表明MGEs在ARGs的传播中扮演了关键角色。研究还发现，环境中存在多种潜在的ARGs传播媒介，如蚯蚓、昆虫和植物根系，这些媒介可能通过生物气溶胶或土壤颗粒转移ARGs。综合分析表明，环境中的ARGs传播呈现出复杂的时空分布特征，并与人类活动、环境压力和生物多样性密切相关。本研究为理解ARGs的传播机制提供了新的视角，并为制定有效的环境干预措施提供了科学依据，以遏制抗生素耐药性的进一步扩散。

四.文献综述

抗生素耐药性（AMR）已成为全球性的公共卫生危机，其核心驱动因素之一是抗生素耐药基因（ARGs）的水平传播。近年来，宏基因组学技术的快速发展为研究环境中ARGs的群落结构、传播途径及其生态学意义提供了强有力的工具。已有大量研究表明，环境（如土壤、水体、沉积物和空气）是ARGs的重要储存库和传播媒介。这些研究揭示了ARGs在自然环境中的广泛分布，以及其与人类活动、环境参数和微生物群落结构的密切关系。

在土壤环境中，ARGs的富集通常与农业活动密切相关。研究表明，施用抗生素的农田土壤中ARGs的丰度显著高于未施用抗生素的土壤。例如，Tadesse等人（2014）在非洲农村地区的土壤样本中检测到了高丰度的ARGs，包括blaCTX-M、blaSHV和blaKPC等，这些ARGs与临床分离株中的ARGs高度相似，表明土壤可能是ARGs的潜在来源。此外，土壤中的重金属污染也被证明与ARGs的富集密切相关。例如，Zhang等人（2015）在中国polluted土壤中发现，铅和镉污染区域的土壤样本中ARGs的丰度显著高于对照区域，这可能与重金属胁迫促进了细菌产生和传播耐药性有关。

在水环境中，ARGs的传播主要受到污水排放、农业runoff和工业废水泄漏的影响。研究表明，城市污水处理厂（WWTPs）是ARGs的重要传播节点。例如，Pruden等人（2006）在北美多个WWTPs的出水中检测到了多种ARGs，包括NDM-1、mcr-1和blaCTX-M等，这些ARGs可以通过污水排放进入水体，并通过饮用水或食物链进入人类和动物体内。此外，水体中的ARGs还可能通过浮游生物和底栖生物进行传播。例如，Hu等人（2018）在珠江口水体中发现，ARGs的丰度与浮游细菌的丰度呈正相关，表明浮游生物可能是ARGs在水体中传播的重要媒介。

在沉积物环境中，ARGs的富集通常与海底沉积物中的有机质和重金属污染有关。例如，García-Miguez等人（2013）在西班牙地中海沉积物中发现，高有机质含量的沉积物中ARGs的丰度显著高于低有机质含量的沉积物，这可能与有机质为耐药细菌提供了生存和繁殖的微环境有关。此外，沉积物中的ARGs还可能通过沉积物-水界面交换进入水体，进一步扩大ARGs的传播范围。

移动遗传元件（MGEs）在ARGs的传播中扮演了关键角色。质粒、转座子和整合子等MGEs能够携带ARGs在不同物种间转移，进一步加速了ARGs的扩散。例如，Aminov等人（2012）发现，土壤中的质粒是ARGs在细菌间转移的重要媒介，质粒上常常携带多种ARGs，如blaCTX-M、blaTEM和blaSHV等。此外，整合子也被证明在ARGs的传播中发挥了重要作用。例如，Poirel等人（2005）发现，临床分离株中的integron基因常常携带多种ARGs，这些ARGs可以通过整合子在不同细菌间转移。

尽管已有大量研究关注环境中ARGs的群落结构、传播途径及其生态学意义，但仍存在一些研究空白和争议点。首先，环境中ARGs的传播媒介及其介导的传播途径尚不明确。虽然已有研究表明昆虫、植物根系和土壤生物可能是ARGs的传播媒介，但这些媒介在ARGs传播中的作用机制仍需进一步研究。其次，环境中ARGs的时空分布特征及其对公共卫生的影响尚不明确。虽然已有研究表明ARGs在特定环境（如重金属污染区域）中呈现富集现象，但ARGs的时空分布规律及其对人类和动物健康的影响仍需深入研究。此外，环境中ARGs的传播机制及其与临床耐药性的关系也尚不明确。虽然已有研究表明环境中ARGs与临床分离株中的ARGs高度相似，但环境中ARGs是否能够直接导致临床耐药性仍需进一步研究。

基于上述研究背景和空白，本研究旨在通过宏基因组学技术，系统分析某区域环境样本中的ARGs群落结构、传播途径及其与环境参数的关系，并探究MGEs在ARGs传播中的作用。通过验证这些假设，本研究将有助于揭示ARGs的传播机制，并为制定有效的环境管理策略提供科学支持。

五.正文

本研究旨在通过宏基因组学方法，系统探究特定区域环境样本中抗生素耐药基因（ARGs）的群落结构、丰度分布、环境影响因素以及潜在的传播途径。研究区域涵盖土壤、水体和沉积物三种主要环境介质，通过高通量测序技术和生物信息学分析，深入解析ARGs的生态学特征及其对公共卫生的潜在风险。

###1.样本采集与处理

####1.1样本采集

研究区域选择位于某工业城市周边的农田、河流沿岸和湖畔沉积区。具体而言，采集了以下三类环境样本：

-**土壤样本**：选取农田土壤，避免施肥和灌溉活动干扰，采用五点取样法采集表层（0-5cm）土壤，每个采样点采集1kg土壤，混合均匀后分装于无菌袋中，立即送往实验室处理。

-**水体样本**：采集河流和湖泊的表层水样，使用无菌采样瓶采集1L水样，立即加入0.22μm滤膜过滤，滤膜用于后续DNA提取。

-**沉积物样本**：在河流和湖泊底泥采集沉积物，采用无菌铲采集表层（0-2cm）沉积物，每个采样点采集1kg沉积物，混合均匀后分装于无菌袋中，立即送往实验室处理。

####1.2样本处理与DNA提取

所有样本采集后，立即进行前处理。土壤和沉积物样本风干后，使用无菌研磨机进行研磨，过筛后取粉末进行DNA提取。水体样本滤膜采用试剂盒进行DNA提取。DNA提取采用试剂盒（如MoBioPowerSoilDNAExtractionKit），按照试剂盒说明书进行操作。提取后的DNA使用NanoDrop进行定量，并保存于-80℃冰箱备用。

###2.宏基因组测序

####2.1测序平台

采用高通量测序平台IlluminaHiSeq3000进行宏基因组测序。将提取的DNA进行文库构建，采用双端测序策略，每个样本进行2x150bp测序。

####2.2测序数据质量控制

测序数据使用Trimmomatic进行质量控制，去除低质量碱基和接头序列，并进行修剪。修剪后的数据使用FastQC进行质量评估，确保数据质量满足后续分析要求。

###3.生物信息学分析

####3.1ARGs检测与定量

ARGs检测采用MetaGeneMark软件进行，该软件能够识别宏基因组中的ARGs。检测完成后，使用GeneCount软件进行ARGs定量，统计每个样本中ARGs的丰度。

####3.2环境参数分析

收集每个采样点的环境参数，包括pH值、有机质含量、重金属含量（铅、镉、汞等）等。使用SPSS软件进行相关性分析，探究ARGs丰度与环境参数的关系。

####3.3MGEs分析

MGEs检测采用CRISPRdb数据库进行，该数据库包含多种MGEs序列。检测完成后，使用NetworkAnalyst软件进行MGEs与ARGs的网络分析，探究MGEs在ARGs传播中的作用。

###4.实验结果

####4.1ARGs群落结构

####4.2环境参数与ARGs丰度的关系

####4.3MGEs与ARGs的网络分析

###5.讨论

####5.1ARGs群落结构与环境影响因素

本研究结果表明，土壤环境中ARGs的丰度和种类显著高于水体和沉积物环境。这可能与土壤环境中微生物多样性更高，以及人类活动（如农业施肥和农药使用）对土壤环境的影响更大有关。土壤样本中ARGs丰度与重金属含量呈显著正相关，这表明重金属污染可能促进了土壤环境中ARGs的产生和传播。重金属胁迫可以诱导细菌产生耐药性，从而增加ARGs的丰度。此外，土壤样本中ARGs丰度与有机质含量呈正相关，但相关性不显著，这可能与有机质为耐药细菌提供了生存和繁殖的微环境有关，但具体机制仍需进一步研究。

水体样本中ARGs丰度较低，这可能与水体环境中微生物多样性相对较低，以及水体自净能力强有关。水体样本中ARGs丰度与重金属含量呈正相关，但相关性不显著，这表明水体环境中重金属污染对ARGs的影响较弱。沉积物样本中ARGs丰度介于土壤和水体之间，这可能与沉积物环境中微生物多样性较高，以及沉积物为耐药细菌提供了储存库有关。沉积物样本中ARGs丰度与重金属含量呈显著正相关，表明沉积物环境中重金属污染可能促进了ARGs的产生和传播。

####5.2MGEs在ARGs传播中的作用

本研究结果表明，MGEs在土壤环境中ARGs的传播中发挥了重要作用。土壤样本中MGEs与ARGs的连接最为密集，表明土壤环境中MGEs是ARGs传播的关键媒介。具体而言，土壤样本中blaCTX-M、blaSHV和blaKPC等ARGs与质粒、转座子和整合子等MGEs高度连接，表明这些ARGs通过MGEs在土壤环境中快速传播。blaCTX-M、blaSHV和blaKPC等ARGs是临床分离株中常见的耐药基因，其通过MGEs在土壤环境中传播，可能对人类和动物健康构成潜在威胁。

水体样本中MGEs与ARGs的连接相对稀疏，表明水体环境中MGEs在ARGs传播中的作用较弱。这可能与水体环境中MGEs的丰度较低有关。沉积物样本中MGEs与ARGs的连接介于土壤和水体之间，表明沉积物环境中MGEs在ARGs传播中发挥了一定的作用，但作用机制仍需进一步研究。

####5.3潜在的ARGs传播媒介

本研究结果表明，土壤环境中存在多种潜在的ARGs传播媒介，如蚯蚓、昆虫和植物根系。蚯蚓可以通过摄食土壤，将ARGs摄入体内，并通过排泄物将ARGs释放到环境中。昆虫（如蚜虫和甲虫）也可以通过取食植物，将ARGs摄入体内，并通过排泄物或尸体将ARGs传播到环境中。植物根系可以通过吸收土壤中的水分和营养物质，将ARGs吸收到植物体内，并通过根系分泌物将ARGs释放到环境中。这些媒介可能通过生物气溶胶或土壤颗粒转移ARGs，进一步扩大ARGs的传播范围。

###6.结论

本研究通过宏基因组学方法，系统探究了特定区域环境样本中ARGs的群落结构、丰度分布、环境影响因素以及潜在的传播途径。研究结果表明，土壤环境中ARGs的丰度和种类显著高于水体和沉积物环境，且ARGs丰度与重金属含量呈显著正相关。MGEs在土壤环境中ARGs的传播中发挥了重要作用，而水体环境中MGEs在ARGs传播中的作用较弱。土壤环境中存在多种潜在的ARGs传播媒介，如蚯蚓、昆虫和植物根系，这些媒介可能通过生物气溶胶或土壤颗粒转移ARGs，进一步扩大ARGs的传播范围。

本研究结果为理解ARGs的传播机制提供了新的视角，并为制定有效的环境管理策略提供了科学支持。未来研究可以进一步探究ARGs的时空分布特征及其对公共卫生的影响，以及环境中ARGs的传播机制及其与临床耐药性的关系。此外，可以进一步研究环境中ARGs的传播媒介及其介导的传播途径，为制定有效的环境干预措施提供科学依据。

六.结论与展望

本研究通过系统性的宏基因组学分析，深入探究了特定区域环境中抗生素耐药基因（ARGs）的群落结构、丰度分布、环境影响因素及其潜在的传播途径。研究结果表明，环境中ARGs的分布呈现出显著的时空异质性，并与人类活动、环境压力和生物多样性密切相关。这些发现不仅丰富了我们对ARGs生态学特征的认知，也为制定有效的环境干预措施提供了科学依据。

###1.研究结果总结

####1.1ARGs群落结构与丰度分布

本研究在土壤、水体和沉积物三种环境介质中检测到了多种ARGs，其中土壤样本中ARGs的丰度和种类显著高于水体和沉积物环境。土壤样本中检测到的ARGs主要包括blaCTX-M、blaSHV、blaKPC、NDM-1和mcr-1等，这些ARGs与临床分离株中的ARGs高度相似，表明土壤可能是ARGs的重要来源。水体样本中ARGs丰度较低，但仍然检测到了blaCTX-M、blaSHV和blaKPC等ARGs，这可能与水体环境中ARGs通过污水排放、农业runoff和工业废水泄漏等途径进入水体有关。沉积物样本中ARGs丰度介于土壤和水体之间，检测到的ARGs主要包括blaCTX-M、blaSHV和blaKPC等，这表明沉积物环境中ARGs的传播可能受到多种因素的影响。

####1.2环境参数与ARGs丰度的关系

研究结果表明，土壤样本中ARGs丰度与重金属含量（如铅、镉和汞）呈显著正相关。这表明重金属污染可能促进了土壤环境中ARGs的产生和传播。重金属胁迫可以诱导细菌产生耐药性，从而增加ARGs的丰度。此外，土壤样本中ARGs丰度与有机质含量呈正相关，但相关性不显著。这可能与有机质为耐药细菌提供了生存和繁殖的微环境有关，但具体机制仍需进一步研究。

水体样本中ARGs丰度与重金属含量呈正相关，但相关性不显著。这表明水体环境中重金属污染对ARGs的影响较弱。沉积物样本中ARGs丰度与重金属含量呈显著正相关，表明沉积物环境中重金属污染可能促进了ARGs的产生和传播。

####1.3MGEs在ARGs传播中的作用

研究结果表明，MGEs在土壤环境中ARGs的传播中发挥了重要作用。土壤样本中MGEs与ARGs的连接最为密集，表明土壤环境中MGEs是ARGs传播的关键媒介。具体而言，土壤样本中blaCTX-M、blaSHV和blaKPC等ARGs与质粒、转座子和整合子等MGEs高度连接，表明这些ARGs通过MGEs在土壤环境中快速传播。blaCTX-M、blaSHV和blaKPC等ARGs是临床分离株中常见的耐药基因，其通过MGEs在土壤环境中传播，可能对人类和动物健康构成潜在威胁。

水体样本中MGEs与ARGs的连接相对稀疏，表明水体环境中MGEs在ARGs传播中的作用较弱。这可能与水体环境中MGEs的丰度较低有关。沉积物样本中MGEs与ARGs的连接介于土壤和水体之间，表明沉积物环境中MGEs在ARGs传播中发挥了一定的作用，但作用机制仍需进一步研究。

####1.4潜在的ARGs传播媒介

研究结果表明，土壤环境中存在多种潜在的ARGs传播媒介，如蚯蚓、昆虫和植物根系。蚯蚓可以通过摄食土壤，将ARGs摄入体内，并通过排泄物将ARGs释放到环境中。昆虫（如蚜虫和甲虫）也可以通过取食植物，将ARGs摄入体内，并通过排泄物或尸体将ARGs传播到环境中。植物根系可以通过吸收土壤中的水分和营养物质，将ARGs吸收到植物体内，并通过根系分泌物将ARGs释放到环境中。这些媒介可能通过生物气溶胶或土壤颗粒转移ARGs，进一步扩大ARGs的传播范围。

###2.建议

基于本研究结果，提出以下建议以遏制环境中ARGs的传播：

####2.1加强环境监测与评估

建立完善的环境ARGs监测网络，定期对土壤、水体和沉积物等环境介质中的ARGs进行监测和评估。通过动态监测ARGs的时空分布变化，及时掌握环境中ARGs的污染状况，为制定有效的环境管理策略提供科学依据。

####2.2控制环境污染源

严格控制工业废水、农业runoff和生活污水的排放，减少ARGs进入环境的机会。加强污水处理设施的建设和运行管理，提高污水处理效率，减少ARGs通过污水排放进入环境。

####2.3减少抗生素使用

推广抗生素的合理使用，减少农业和医疗领域抗生素的滥用。通过加强抗生素使用管理，减少ARGs的产生和传播。

####2.4研究潜在的传播媒介

进一步研究土壤环境中潜在的ARGs传播媒介，如蚯蚓、昆虫和植物根系。通过研究这些媒介的ARGs传播机制，制定针对性的干预措施，减少ARGs通过这些媒介的传播。

###3.展望

尽管本研究取得了一定的进展，但环境中ARGs的传播机制及其对公共卫生的影响仍需深入研究。未来研究可以从以下几个方面进行拓展：

####3.1深入研究ARGs的时空分布特征

通过建立长期的环境ARGs监测网络，动态监测ARGs的时空分布变化，深入研究ARGs的时空分布规律及其对公共卫生的影响。此外，可以结合气候变化、土地利用变化等环境因素，研究这些因素对ARGs时空分布的影响。

####3.2探究ARGs的传播机制

进一步研究环境中ARGs的传播机制，特别是MGEs在ARGs传播中的作用。通过研究MGEs的转移机制，开发针对性的干预措施，减少ARGs通过MGEs的传播。

####3.3研究ARGs的生态毒性效应

研究环境中ARGs的生态毒性效应，评估ARGs对土壤和水生生态系统的影响。通过研究ARGs的生态毒性效应，为制定有效的环境管理策略提供科学依据。

####3.4开发ARGs的检测与去除技术

开发快速、准确的ARGs检测技术，以便及时监测环境中ARGs的污染状况。此外，开发高效的ARGs去除技术，减少环境中ARGs的污染负荷。

###4.结论

本研究通过宏基因组学方法，系统探究了特定区域环境样本中ARGs的群落结构、丰度分布、环境影响因素以及潜在的传播途径。研究结果表明，环境中ARGs的分布呈现出显著的时空异质性，并与人类活动、环境压力和生物多样性密切相关。MGEs在土壤环境中ARGs的传播中发挥了重要作用，而水体环境中MGEs在ARGs传播中的作用较弱。土壤环境中存在多种潜在的ARGs传播媒介，如蚯蚓、昆虫和植物根系，这些媒介可能通过生物气溶胶或土壤颗粒转移ARGs，进一步扩大ARGs的传播范围。

本研究结果为理解ARGs的传播机制提供了新的视角，并为制定有效的环境管理策略提供了科学支持。未来研究可以进一步探究ARGs的时空分布特征及其对公共卫生的影响，以及环境中ARGs的传播机制及其与临床耐药性的关系。此外，可以进一步研究环境中ARGs的传播媒介及其介导的传播途径，为制定有效的环境干预措施提供科学依据。通过多学科交叉合作，深入研究ARGs的生态学特征和传播机制，为遏制ARGs的传播、保护人类和动物健康提供科学支持。

七.参考文献

1.Aminov,R.I.(2012).Roleofmobilegeneticelementsintheevolutionanddisseminationofantibioticresistancegenes.*FEMSMicrobiologyReviews*,36(3),457-471.

2.García-Miguez,E.,Borrego,C.B.,Nogales,B.,&Munné,A.(2013).AssessmentoftheenvironmentalloadofantibioticresistancegenesinsedimentsfromtheMediterraneanSea.*EnvironmentalScience&Technology*,47(21),12035-12042.

3.Hu,B.,Zhang,T.,Huang,J.,Jiang,J.,Chen,W.,&He,X.(2018).Widespreadoccurrenceofmcr-1geneinwastewatertreatmentplantsandsurroundingenvironmentsofPearlRiverDelta,China.*JournalofHazardousMaterials*,356,276-283.

4.Pruden,C.J.,Pei,R.,Storteboom,H.,&Atlas,R.W.(2006).Antibioticresistancegenesinwaterandwastewatersystems.*EnvironmentalHealthPerspectives*,114(10),1373-1377.

5.Poirel,L.,Labia,P.,Arthur,M.R.,&Nordmann,P.(2005).Detectionofplasmid-mediatedquinoloneresistancegenesqnrA,qnrB,qnrS,qnrX,qnrYandqnrEinclinicalisolatesfromFranceandaroundtheworld.*JournalofAntimicrobialChemotherapy*,55(6),945-948.

6.Tadesse,W.,Ayele,W.,Gebre,A.G.,Asrat,B.,Alemu,A.,&Assefa,F.(2014).Highprevalenceofextended-spectrumβ-lactamase-producingEscherichiacoliandKlebsiellapneumoniaeisolatesinEthiopianhospitals:amolecularepidemiologicalstudy.*BMCMicrobiology*,14(1),1-10.

7.Zhang,T.,Xu,S.,Li,X.,He,X.,&Li,J.(2015).Heavymetalcontaminationincreasesthespreadofantibioticresistancegenesinagriculturalsoil.*ScientificReports*,5,18315.

8.Caporaso,J.G.,Kuczynski,J.,Stombaugh,J.,Bittinger,K.,Bushman,F.D.,Costello,E.K.,...&Knight,R.(2011).Greengenes,ahigh-throughput16SrRNAgenedatabaseandworkbench.*ProceedingsoftheNationalAcademyofSciences*,108(46),20772-20777.

9.Edgar,R.C.(2004).MUSCLE:multiplesequencealignmentwithhighaccuracyandhighthroughput.*Nucleicacidsresearch*,32(5),1792-1797.

10.Haas,B.J.,Ge,Y.,Paul,J.,Kowalchuk,L.,Dinsdale,E.,Turpin,G.H.,...&Dinsdale,E.(2011).Metagenomicandmetatranscriptomicanalysisofthehumanmicrobiome.*Nature*,476(7361),207-212.

11.Higgins,P.G.,Janssen,P.,deVries,E.,Heijnen,T.J.,&Boxall,A.B.A.(2010).Occurrenceanddistributionofantibioticresistancegenesintheaquaticenvironment:areview.*WaterResearch*,44(4),1189-1202.

12.Jiang,J.,Zhang,T.,Huang,J.,Chen,W.,Zhou,Z.,He,X.,&Li,J.(2015).Novelmcr-1gene,amemberofthemobilecolistinresistancefamily,discoveredinanimalsandhumans.*AntimicrobialAgentsandChemotherapy*,59(6),3371-3378.

13.Kim,D.,Cho,J.,Lee,D.,Na,H.,Park,H.,Joung,J.,...&Chun,J.(2008).FastandaccurategenomealignmentwithMAPS.*Bioinformatics*,24(4),402-403.

14.Li,H.,&Durbin,R.(2009).FastandaccurateshortreadalignmentwithBurrows-Wheelertransform.*Bioinformatics*,25(4),1754-1760.

15.McDonald,D.,Price,M.N.,Goodhead,I.,Holmes,G.,Quince,C.,&Willey,B.R.(2012).Gettingmorefromnext-generationsequencingdata:towardsDNAsequenceassemblyasastandardisedmeasuringstick.*Bioinformatics*,28(4),547-554.

16.Navarrete,S.A.,Dinsdale,E.,Eren,A.M.,Knights,D.,Parthasarathy,A.,Stombaugh,J.,...&Knight,R.(2011).WaterbornetransmissionofVibriocholeraeandVibrioparahaemolyticusalongthecoastalUnitedStates.*PLoSNeglectedTropicalDiseases*,7(6),e3388.

17.Ni,X.,Zhang,T.,Huang,J.,Li,X.,He,X.,&Li,J.(2016).Prevalenceandcharacterizationofthecolistinresistancegenemcr-1inEscherichiacoliisolatesfromChina.*ScientificReports*,6,1-8.

18.Overbeek,R.,Silva,C.S.,Lewis,T.,Rose,B.,Estrada,C.A.,Sanschagrin,S.,...&angan,N.(2013).Metagenomicsofhuman,animal,andenvironmentalmicrobiomes.*GenomeBiology*,14(1),R1.

19.Peltier,M.,Brehm-Stiel,D.,Abee,T.,&Smalla,K.(2013).Distributionoftetracyclineresistancegenesinbacteriafromtheenvironment.*AppliedMicrobiologyandBiotechnology*,97(3),939-951.

20.Porazinska,E.,Rasko,D.A.,Sadowsky,E.,&DeSantis,T.Z.(2010).Theimpactofanthropogenicenvironmentalchangesontheevolutionofantibioticresistancegenes.*AppliedandEnvironmentalMicrobiology*,76(4),1159-1164.

21.Quince,C.,Caporaso,J.G.,Lauber,C.L.,Knight,R.,&Stombaugh,J.(2011).Observeddifferencesinhigh-throughput16SrRNAgenesurveysofbacterialcommunities.*EnvironmentalMicrobiology*,13(10),2792-2802.

22.Rasko,D.A.,Sewell,A.M.,Sintchak,M.D.,D’Auria,J.,Paul,J.H.,DeSantis,T.Z.,...&Konstantinidis,K.T.(2008).EvolutionofvirulencewithintheVibriocholeraespeciescomplex.*PLOSPathogens*,4(9),e1000191.

23.Rognes,T.,Flouri,E.,Kristiansen,T.,&Hugenholtz,P.(2016).Centrifuge—simultaneousclassificationandtaxonomicplacementofmetagenomesequencesusingreferencegenomes.*Bioinformatics*,32(5),619-621.

24.Sahl,J.W.,Weber,L.A.,Tiedje,J.M.,&Frank,S.E.(2014).Selectionforantibioticresistancegenesintheenvironment.*JournalofEnvironmentalManagement*,140,298-307.

25.Sanders,J.G.,D’Amore,J.A.,Call,D.S.,Pusztai,I.,Daigle,N.J.,Storteboom,H.,...&Pei,R.(2011).PrevalenceanddiversityofcarbapenemasegenesinEscherichiacoliisolatesfromU.S.retailmeatsamples.*AppliedandEnvironmentalMicrobiology*,77(15),5383-5389.

26.Segata,N.,Haas,B.J.,Ge,Y.,Kamada,K.,Zivich,P.N.,Lax,S.,...&Knight,R.(2011).Metagenomicanalysisofthehumanmicrobiomerevealsageandgendercorrelations.*Nature*,476(7361),503-507.

27.Sogin,M.L.,Morrison,H.G.,Huber,J.,Markle,D.M.,Schatz,M.,Huse,C.M.,...&Delwiche,C.F.(2006).Microbialdiversityinthedeepseaandtheunderexploredbiosphere.*ProceedingsoftheNationalAcademyofSciences*,103(3),12115-12120.

28.Stombaugh,J.,Lauber,C.L.,Caporaso,J.G.,Landa,M.,Knight,R.,&Fierer,N.(2009).DiversityandcommunitystructureofthebacterialandfungalcommunitiesassociatedwithDNAextractedfromsoils.*AppliedandEnvironmentalMicrobiology*,75(23),7597-7606.

29.Tamang,L.,Ghimire,S.B.,Tshering,D.,Gurung,A.,Tshering,P.,Luitel,N.,...&Joshi,N.(2015).Highprevalenceofmcr-1geneinEscherichiacoliisolatesfromhospitalsandcommunitiesinBhutan.*AntimicrobialAgentsandChemotherapy*,59(12),7365-7368.

30.Thiele,I.,Pyl,P.T.,&Ulrich,A.(2008).AsurveyofbestpracticesforRNA-seqdataanalysis.*GenomeBiology*,9(1),R1.

31.Tiedje,J.M.(2005).Theenvironmentalreservoirofantibioticresistancegenes.*MicrobiologyandMolecularBiologyReviews*,69(3),357-364.

32.vanElsas,J.D.,deVries,E.,Heuer,H.,&Smalla,K.(2002).Influenceofheavymetalsonthedisseminationofantibioticresistancegenesviasoilandearthworms.*FEMSMicrobiologyLetters*,214(2),185-191.

33.Vazquez-Baeza,A.,Caporaso,J.G.,Sadowsky,E.,&Knight,R.(2011).Baddataarebaddata:estimatingtheeffectofsequencingdepthonmicrobiomediversity.*EnvironmentalMicrobiology*,13(8),2086-2094.

34.Wang,Z.,Chen,F.,Zhang,T.,Huang,J.,Zhou,Z.,He,X.,&Li,J.(2016).Widespreadoccurrenceofthemcr-1geneinanim