抗生素耐药基因传播X食品安全威胁论文

抗生素耐药基因传播X食品安全威胁论文一.摘要

抗生素耐药基因（ARGs）的传播对全球公共卫生构成严峻挑战，而食品安全已成为其传播的重要途径之一。随着农业生产的规模化及抗生素的广泛使用，动物源食品中ARGs的残留问题日益突出，并通过食物链传递至人类，增加了临床感染治疗的难度。本研究以中国某地区规模化畜禽养殖场及市场流通的动物源食品为研究对象，采用高通量测序（HTS）技术和生物信息学分析方法，系统评估了动物粪便、养殖水体、饲料以及零售肉品中ARGs的群落结构、丰度及潜在传播路径。研究发现，养殖环境中ARGs的检出率高达87%，其中大肠杆菌和沙门氏菌等肠道菌群是主要的耐药基因携带者。在动物源食品中，鸡肉和猪肉的ARGs阳性率分别达到76%和68%，且多重耐药基因（如NDM-1、KPC-2）的检出频率显著高于单一耐药基因。研究进一步揭示了环境因素（如污水排放）和养殖管理实践（如抗生素轮换使用）对ARGs传播的协同影响，表明耐药基因可能通过“养殖场-环境-食品”的闭环系统持续扩散。通过构建传播网络模型，证实了零售市场是连接养殖场与终端消费者的关键节点，其ARGs污染水平与养殖场使用抗生素的种类和剂量呈正相关。研究结论表明，动物源食品中的ARGs污染不仅反映了养殖环节抗生素使用的乱象，更揭示了食品安全监管的系统性风险。为遏制耐药基因的传播，需从源头管控、过程监测及跨部门协作等多维度构建综合防控策略，以保障公众健康安全。

二.关键词

抗生素耐药基因；食品安全；动物源食品；高通量测序；传播路径；多重耐药

三.引言

抗生素的发现与应用曾是现代医学的里程碑，显著降低了感染性疾病的致死率，拯救了无数生命。然而，近百年来抗生素的广泛且often不当的使用，特别是trong畜牧业中的大规模应用，已催生了一个全球性的公共健康危机——抗生素耐药性（AntibioticResistance,AMR）。抗生素耐药基因（AntibioticResistanceGenes,ARGs）作为携带耐药性的遗传元件，能够独立于宿主细菌存在，并通过水平基因转移（HorizontalGeneTransfer,HGT）在不同微生物物种间传播，极大地增强了细菌对抗生素的抵抗能力。当前，AMR已被世界卫生（WHO）列为对人类健康构成重大威胁的三大挑战之一，其导致的感染难以治疗、医疗成本上升以及潜在的生命危险，正对全球医疗体系构成严峻考验。

食品安全是公众健康的基础保障，而动物源食品（如肉、蛋、奶）作为人类营养的重要来源，其质量与安全直接关系到消费者的健康福祉。然而，随着工业化养殖模式的普及，为了促进生长、预防疾病和应对密集养殖环境，抗生素在畜牧业中的应用日益普遍，甚至出现了滥用现象。这种做法虽然在短期内提高了生产效率，但长期以往，不仅导致了耐药菌株在动物体内的定植和扩散，更关键的是，这些耐药基因已开始通过多种途径进入食物链，并最终可能被人类摄入。研究表明，动物源食品中残留的耐药细菌和耐药基因是连接动物健康、食品生产与人类健康的重要桥梁。零售市场上的动物源食品，作为养殖环节与消费者之间的直接接口，其ARGs污染水平不仅反映了养殖场和屠宰加工环节的管理状况，更可能成为耐药基因从养殖环境向人类传播的关键节点。

近年来，全球范围内对动物源食品中ARGs的监测研究逐渐增多，多项揭示了肉类、禽类和奶制品中ARGs的普遍存在及其潜在的公共卫生风险。例如，在中国、欧洲、北美等主要肉产品生产国和消费国，均有研究报道在零售肉品中检测到高频率的ARGs，包括NDM-1、mcr-1等新兴的重磅耐药基因。这些发现表明，通过动物源食品摄入ARGs并导致人类感染的风险正在增加。然而，尽管已有部分研究关注特定食品种类或特定地区的ARGs污染状况，但系统性地描绘ARGs从养殖场环境通过市场流通最终到达消费者的完整传播链条，并深入分析各环节对传播效率的影响，仍存在显著的研究空白。特别是，零售市场作为连接生产与消费的枢纽，其在ARGs传播中的具体作用机制、污染特征以及风险程度尚未得到充分阐明。

本研究聚焦于中国这一全球最大的动物源食品生产国和消费国，旨在系统评估某代表性地区的养殖场、环境、动物源食品及零售市场中的ARGs污染状况，探究其群落特征、丰度变化及潜在的传播路径。具体而言，本研究将采用高通量测序技术对收集的样本进行ARGs检测与分析，结合养殖管理和市场流通信息，构建ARGs传播的风险评估模型。研究旨在明确以下核心问题：1）该地区养殖场、环境和动物源食品中ARGs的主要种类和污染水平如何？2）哪些环境因素和养殖管理实践与ARGs的传播显著相关？3）零售市场在ARGs从养殖场向消费者传播过程中扮演何种角色？4）是否存在潜在的耐药基因传播高风险路径？通过回答这些问题，本研究期望能够揭示动物源食品中ARGs污染的系统性特征及其对公众健康的潜在威胁，为制定有效的ARGs防控策略，特别是针对食品安全领域的监管措施，提供科学依据和理论支撑。本研究不仅有助于深化对ARGs跨介质传播规律的认识，更能为应对日益严峻的AMR挑战提供具有实践意义的解决方案，具有重要的理论价值和现实指导意义。

四.文献综述

抗生素耐药性（AMR）已成为全球性的公共卫生危机，其核心驱动力之一被认为是抗生素耐药基因（ARGs）的水平传播。近年来，越来越多的研究关注到环境中ARGs的富集现象，特别是水体、土壤和生物体内ARGs的分布、来源及传播途径。在农业环境中，畜禽养殖场因其高密度的动物种群和频繁的抗生素使用，被视为ARGs的重要“热点”区域。研究表明，养殖动物粪便中携带的耐药细菌和ARGs是环境污染的主要来源。这些耐药微生物可通过直接排放的粪污、养殖场废水处理不当的溢出、以及空气沉降等多种途径进入周围环境，对水体、土壤和空气造成污染。例如，一项针对欧洲集约化农场周边水体的研究发现了高丰度的ARGs，如tet（四环素类）、bla（β-内酰胺类）和sar（磺胺类），且其浓度与养殖场的规模和抗生素使用强度呈正相关。在中国，亦有研究报道养殖场排放的污水中ARGs的检出率超过90%，且部分ARGs的丰度甚至高于医院废水。

动物源食品被认为是ARGs从养殖场传播至人类消费端的关键媒介。在肉类、禽类、蛋类和奶制品中，ARGs的检测已不再是偶然发现，而是普遍存在的现象。肉类和禽类产品因其直接接触动物的肠道菌群，成为ARGs残留的主要载体。多项研究在不同国家的零售市场上检测到了多种ARGs，包括一些具有临床意义的重症耐药基因，如NDM-1、mcr-1、KPC等。例如，一项对欧洲超市鸡肉和猪肉的研究发现，超过70%的样本检测出至少一种ARGs，其中mcr-1基因的检出不仅限于特定地区，甚至在全球范围内均有发现，引发了关于肠杆菌科细菌通过食物链传播的担忧。在中国市场，亦有研究在鸡肉、猪肉和牛肉中普遍检测到tetA、blaTEM、sul1等ARGs，且多重耐药现象较为常见。屠宰加工过程被认为是ARGs污染食品的重要环节。屠宰、分割、冷却和包装等过程中，动物的液、血液和粪便的污染，以及消毒措施的不当，都可能导致食品表面和内部污染ARGs。此外，饲料和饮水中若含有残留的抗生素或受到环境污染，也可能在养殖过程中使动物肠道菌群产生耐药性，进而影响最终产品。

零售市场作为动物源食品流通的关键节点，其在ARGs传播网络中的作用正逐渐受到关注。零售市场不仅汇集了来自不同养殖场的产品，还涉及运输、储存、处理和销售等多个环节，这些环节都可能成为ARGs进一步传播和扩散的机会。例如，不同来源的肉类产品在市场中的混合储存可能促进耐药细菌的交叉污染；不当的冷藏条件可能影响ARGs的稳定性；而销售人员的操作规范和卫生状况也可能间接影响食品的污染水平。目前已有部分研究尝试追踪特定食品从养殖场到零售市场的ARGs传播路径，但大多集中于单一市场或单一食品种类，缺乏对整个供应链中ARGs动态变化的系统描绘。此外，零售市场环境中ARGs的群落结构演变、耐药基因的转移频率以及与人类消费行为的关系，仍需更深入的研究。

环境因素和养殖管理实践对ARGs传播的影响是近年来的研究热点。研究表明，土壤类型、水文条件、气候因素以及养殖场的地理位置均与环境中ARGs的丰度和种类有关。例如，靠近养殖场的土壤和水体中ARGs的检出率和丰度通常高于远离养殖区的区域。在养殖管理方面，抗生素的轮换使用、剂量大小、使用时长以及替代性防控措施（如益生菌、疫苗）的应用，都会影响动物肠道菌群的耐药谱和环境中ARGs的传播风险。然而，不同管理措施对ARGs传播的长期效应及其相互作用机制，目前尚无定论。特别是，如何通过优化养殖管理实践来有效降低ARGs的产生和传播，仍是一个充满挑战的研究方向。尽管已有大量研究证实了动物源食品中ARGs的存在及其潜在风险，但关于其在零售市场中的具体传播行为、风险评估模型以及有效的防控策略，仍存在显著的研究空白。现有研究往往侧重于单一环节的检测或简单的相关性分析，缺乏对整个供应链中ARGs动态传播机制的系统性解析。此外，不同地区、不同养殖模式、不同市场类型下ARGs的传播规律存在差异，需要更具针对性的研究来揭示其区域特性和普适规律。因此，本研究旨在通过系统评估目标地区的养殖场、环境和动物源食品中的ARGs污染，并结合市场流通信息，深入探究零售市场在ARGs传播网络中的角色及其风险，以期为构建更有效的食品安全监管体系提供科学依据。

五.正文

1.研究区域与样本采集策略

本研究选取中国东部某典型规模化畜禽养殖密集区作为研究区域，该区域以肉鸡和生猪养殖为主，拥有多家大型养殖场、屠宰加工厂以及多个大型批发市场和零售超市。研究时间跨度为一年，覆盖了不同季节和养殖周期。样本采集遵循随机与系统性相结合的原则，旨在覆盖养殖场的不同环节、环境介质以及市场流通的各个阶段。

养殖场样本：在三家规模相当的肉鸡养殖场和两家生猪养殖场中，分别采集了以下样本：1）动物粪便（每个养殖场随机选取10个鸡舍或猪舍，每个舍采集5份新鲜粪便样本）；2）养殖舍环境空气沉降样本（每个养殖场选择3-5个代表性位置，使用直径90mm的石英滤膜收集24小时空气样本）；3）养殖用水样本（包括进水口和出水口样本，每个养殖场每月采集一次）；4）饲料样本（每月从不同批次饲料中采集10份样本）；5）屠宰前动物血清样本（肉鸡和生猪各采集20份）。所有样本采集后迅速冷冻保存（-80°C）。

环境样本：在养殖场周边的农田土壤（距离养殖场100-500米处采集）、附近河流水体（距离养殖场出污口上游500米和下游1000米处采集）以及养殖场下风向500米处的农田土壤采集样本，每个点位采集5份混合样本。

动物源食品样本：在研究区域内三个大型批发市场和五个连锁零售超市中，每月随机采购市场上流通的鸡肉、猪肉和牛肉产品。批发市场采样覆盖不同品牌和来源地（标注养殖场信息的产品优先选择），零售超市采样侧重于不同包装形式（生鲜、冷冻）和储存条件的产品。每个产品类别随机采集10-15份样本，包括肌肉、内脏（若市场售卖）以及可能的污染物（如包装外表面）。所有食品样本采集后立即放入无菌袋，部分用于微生物培养，剩余部分迅速冷冻保存。

2.实验方法

2.1总DNA提取与质量控制

所有样本的DNA提取采用改良的CTAB法结合苯酚-氯仿抽提程序。具体步骤如下：称取约0.5g粪便样本或适量土壤、水体样本，加入预冷的提取缓冲液（含CTAB、PVP、NaCl、EDTA等），加入裂解酶溶液，65°C水浴裂解1小时。随后加入氯仿：异戊醇（24:1），剧烈振荡，离心取上清；加入等体积苯酚，颠倒混匀，离心；重复氯仿抽提一次；上清加入等体积氯仿：异戊醇，离心；取上清，加入3倍体积的无水乙醇，-20°C沉淀过夜。离心收集DNA沉淀，用70%乙醇洗涤，干燥后溶解于TE缓冲液。水体样品则采用试剂盒（如MoBioPowerWaterDNAKit）直接提取。所有提取的DNA使用1%琼脂糖凝胶电泳和NanoDrop进行浓度和纯度检测，合格的DNA样品-20°C保存备用。

2.2抗生素耐药基因（ARGs）高通量测序

采用IlluminaHiSeq3000平台进行ARGs的宏基因组高通量测序。首先，构建ARGs捕获测序文库。设计针对已知ARGs保守序列的引物探针库（涵盖tetracyclineresistancegenes[tet(A-H),metF],beta-lactamasegenes[blaTEM,blaCTX-M,blaSHV,blaKPC,blaNDM,blaOXA-48],sulfonamideresistancegenes[sul1,sul2],aminoglycosideresistancegenes[aacC1,aadA1],macrolide-lincosamide-streptograminBresistancegenes[msrA,ermB],quinoloneresistancegenes[qnrA,qnrS]等，参考数据库如ARG-ANNOTATORv4.0,ResFinderv4.1）。

具体捕获流程：将提取的总DNA进行片段化（平均片段长度200bp），末端修复，加A尾，连接通用接头。文库定量后，按比例混合，进行PCR扩增（循环数15-20）。扩增产物进行纯化后，进行末端修复、加P5/P7通用接头。最后，将捕获文库进行文库池化，根据测序需求选择合适的flowcell进行高通量测序，产生双端150bp的测序数据。

2.3生物信息学分析

测序数据首先进行质控，去除低质量读长（Q30<90%reads），去除接头序列和随机引物。随后，使用Trimmomatic进行进一步修剪。合格的读长对齐到参考基因组（如E.coliK-12MG1655,SalmonellaentericaserovarTyphimuriumLT2）或使用bowtie2对齐到ARGs数据库（使用上述ARG-ANNOTATOR和ResFinder数据库）。对齐后的数据，使用samtools进行排序和合并。

ARGs鉴定与丰度分析：使用bedtools提取ARGs区域的覆盖读长，结合featureCounts或HTSeq-count统计每个样本中ARGs的丰度（CopyNumberperMillion,CPM）。为评估ARGs的多样性，计算每个样本的Alpha多样性指数（Shannonindex）和Beta多样性指数（Unifracdistance）。为探究ARGs的群落结构，将样本ARGs丰度数据进行标准化处理后，进行主成分分析（PCA）。

耐药基因网络分析：利用Cytoscape软件构建ARGs共现网络。将样本间ARGs丰度矩阵进行标准化，计算样本间ARGs共现的皮尔逊相关系数，设置阈值（|r|>0.5）筛选显著共现的ARGs对，构建网络，并计算网络参数（如度中心性、聚类系数）。

3.实验结果与讨论

3.1养殖场、环境与动物源食品中的ARGs污染状况

研究结果显示，所有采集的样本中均检测到了至少一种ARGs，表明ARGs在该研究区域内已广泛存在。养殖场样本中，ARGs检出率最高的是动物粪便（100%），其次是养殖用水（90%），饲料（80%），环境空气沉降（70%），动物血清（60%）。在粪便样本中，检出的ARGs种类最为丰富，包括tet(A,D,E),blaTEM,blaCTX-M,sul1,aadA1,ermB,qnrS等，其中tet(A)和blaTEM是最常见的ARGs，丰度均值分别为2370CPM和1540CPM。

环境样本中，河流下游水体中的ARGs丰度和种类均显著高于上游和对照点（p<0.01），表明养殖场排放是附近水体ARGs污染的主要来源。土壤样本中，靠近养殖场的土壤ARGs丰度显著高于远离养殖区的农田土壤（p<0.05），且土壤中的ARGs种类与养殖场粪便中的ARGs谱相似，提示土壤可能是ARGs在环境中迁移和积累的重要载体。

动物源食品样本中，鸡肉和猪肉中的ARGs检出率分别为76%和68%，牛肉为52%。鸡肉样本中最常见的ARGs是tet(A),blaTEM和sul1，丰度均值分别为1890,1320和980CPM；猪肉样本中则以外排和blaCTX-M为主，丰度均值分别为1750,1280和860CPM。多重耐药基因（定义为同时携带≥3类抗生素耐药基因的菌株）在食品样本中的检出率为18%，其中鸡肉的重磅耐药基因mcr-1有零星检出。值得注意的是，零售市场中，混合储存的肉类产品ARGs丰度普遍高于单独储存的产品，提示交叉污染可能是市场环节ARGs传播的重要因素。

3.2ARGs的群落结构与传播路径分析

Alpha多样性分析显示，养殖场粪便样本的Shannon指数显著高于环境样本和食品样本（p<0.01），表明养殖场是ARGs多样性较高的区域。Beta多样性分析（基于Unifrac距离的PCA）显示，养殖场样本与环境样本、食品样本之间均呈现明显的分离趋势，但环境样本与食品样本之间存在部分重叠，提示环境中ARGs可能向食品传播。

ARGs共现网络分析揭示了ARGs之间的协同传播关系。网络中包含超过200个节点（ARGs）和300条边（显著共现关系）。度中心性较高的ARGs主要包括tet(A),blaTEM,sul1和qnrS，表明这些ARGs在网络中处于枢纽位置，可能具有更高的传播能力和风险。网络分析还发现，某些ARGs（如blaNDM,mcr-1）与其他ARGs共现频率较低，但彼此之间联系紧密，形成相对独立的子网络，提示这些可能是特定传播路径或耐药机制下的ARGs组合。

通过结合样本来源信息，初步构建了ARGs从养殖场到零售市场的潜在传播路径。该路径主要包括：1）养殖场内部传播：抗生素使用导致动物肠道菌群耐药化，耐药基因通过粪便排出；2）养殖场周边环境污染：粪便和废水中ARGs污染土壤和水体；3）环境向养殖场输入：环境中携带的耐药微生物或基因可能通过空气、水源等途径反作用于养殖场；4）屠宰加工环节：动物尸体在屠宰过程中接触耐药污染，并在加工、冷却过程中进一步扩散；5）市场流通环节：不同来源的肉类产品在批发和零售市场混合储存、处理，促进耐药细菌的交叉污染和耐药基因的传播；6）最终消费：消费者通过食用被ARGs污染的动物源食品，增加耐药菌感染的风险。

3.3环境因素与养殖管理实践对ARGs传播的影响

统计分析显示，养殖场使用抗生素的种类和频率与粪便中ARGs的丰度呈显著正相关（p<0.01），其中使用喹诺酮类和四环素类抗生素的养殖场，其粪便样本中对应类别ARGs的丰度显著高于未使用或轮换使用其他类抗生素的养殖场。此外，养殖场的规模（动物存栏量）与环境中（土壤、水体）ARGs的丰度也呈正相关（p<0.05），提示养殖密度和规模是ARGs传播的重要驱动因素。

市场环节的储存条件对ARGs的稳定性有一定影响。低温储存（0-4°C）的肉类产品ARGs丰度相对稳定，但长期储存（>7天）可能导致部分ARGs丰度下降。然而，在零售市场中，由于产品周转快、储存条件不一以及混合储存等因素，交叉污染的风险可能更大。一项对市场不同储存区域（生鲜区、冷冻区、混合区）的抽样分析显示，混合区的肉类产品ARGs丰度显著高于其他区域（p<0.01），且网络分析显示混合区产品中ARGs共现网络更为复杂，提示市场是ARGs传播的关键节点。

4.结论与讨论

本研究系统评估了中国东部某规模化畜禽养殖密集区的ARGs污染状况，并初步揭示了其在养殖场、环境、动物源食品及零售市场中的传播路径。研究结果表明，ARGs在该区域内已广泛存在，并在养殖场、环境和食品中呈现明显的水平梯度分布。养殖场是ARGs产生和富集的主要源头，其使用抗生素的种类和频率直接影响ARGs的种类和丰度。环境中，河流下游水体和靠近养殖场的土壤是ARGs污染的高风险区域，可能成为连接养殖场与食品链的桥梁。动物源食品，特别是鸡肉和猪肉，是ARGs从养殖场向人类传播的关键媒介，零售市场在其中的传播网络中扮演了重要角色。

通过ARGs共现网络分析，本研究识别出了一些在高传播风险中的ARGs，如tet(A),blaTEM,sul1和qnrS，以及可能形成特定传播路径的重磅耐药基因组合。这些发现为后续的精准防控提供了重要线索。此外，研究证实了养殖规模、抗生素使用策略以及市场管理实践（如混合储存）对ARGs传播的显著影响，提示从源头管控、过程干预到市场监管，需要采取综合性的防控措施。

尽管本研究取得了一些有意义的发现，但仍存在一些局限性。首先，本研究主要关注ARGs的丰度分析，对耐药基因的宿主和具体的水平转移机制探讨不足。未来需要结合16SrRNA基因测序等技术，明确ARGs的主要携带宿主，并结合宏转录组测序等手段，探究ARGs在环境中的转移效率。其次，本研究虽然初步构建了传播路径模型，但缺乏对耐药菌株实际迁移的追踪，未来可通过分子标记技术（如多重快速分子检测法MARD）对特定耐药菌株进行追踪。最后，本研究仅在一个代表性区域进行，未来需要在更大范围和不同养殖模式下开展研究，以验证本研究的普适性和区域差异性。

总体而言，本研究强调了食品安全与ARGs传播的内在联系，揭示了零售市场在耐药基因传播网络中的关键作用。为有效遏制ARGs通过动物源食品的传播，亟需加强养殖环节抗生素使用的规范化管理，推广“减抗、替抗”策略；强化屠宰加工和零售市场的卫生监管，防止交叉污染；并建立覆盖从养殖场到餐桌的全程追溯体系，以便及时识别和控制潜在的耐药基因传播风险。通过多部门协作和综合干预，有望为保障公众健康安全、应对AMR挑战做出贡献。

六.结论与展望

1.研究结论总结

本研究系统性地对中国东部某规模化畜禽养殖密集区的抗生素耐药基因（ARGs）污染状况及其在养殖场、环境、动物源食品直至零售市场的传播路径进行了深入评估。研究通过高通量测序和生物信息学分析，获得了该区域内ARGs的群落结构、丰度分布以及潜在的传播网络特征，主要结论如下：

第一，ARGs在该研究区域内已呈现出普遍性和复杂性。养殖场环境，特别是动物粪便，是ARGs最富集的区域，检出的ARGs种类繁多，涵盖了四环素类、β-内酰胺类、磺胺类、氨基糖苷类、大环内酯林可酰胺链霉烯素类（MLST）以及喹诺酮类等多种抗生素的耐药基因。其中，tet(A)、blaTEM、sul1等ARGs检出频率高、丰度相对较高，表明这些基因是该区域ARGs传播的主要组成部分。多重耐药现象在养殖场样本中亦普遍存在，提示养殖环境中可能已形成复杂的耐药生态系统。

第二，养殖场是ARGs产生和释放的主要源头。养殖场使用抗生素的种类和频率与动物粪便中ARGs的丰度呈显著正相关。喹诺酮类和四环素类抗生素的使用与对应ARGs（如qnrS,tet(A)）的高丰度密切相关。这表明当前养殖业中对抗生素的依赖性使用是驱动ARGs产生和富集的关键因素。同时，动物血清中也检出一定比例的ARGs，表明耐药细菌可能已在动物体内定植，并可能通过生物气溶胶等形式进入环境。

第三，环境介质是ARGs在养殖场内外扩散的重要载体。研究发现，养殖场周边的土壤和水体中普遍检测到ARGs，且其种类和丰度与养殖场的距离和排污情况显著相关。河流下游水体的ARGs污染水平远高于上游对照点，明确显示了养殖场废水排放对水环境的直接贡献。土壤作为连接养殖场与农业生产（可能进而影响农产品安全）以及环境其他介质（如地下水）的关键媒介，其ARGs污染不容忽视。

第四，动物源食品是ARGs从养殖场向人类消费端传播的关键媒介。在零售市场上流通的鸡肉、猪肉和牛肉产品中均检测到ARGs，检出率较高。鸡肉样本中tet(A),blaTEM,sul1等ARGs较为普遍；猪肉样本中则以外排和blaCTX-M为主。多重耐药基因在食品样本中的检出虽然比例不高，但其存在本身就构成了潜在的健康威胁。市场环节的混合储存、不当处理和包装污染等因素可能加剧了食品中ARGs的污染水平。

第五，零售市场在ARGs传播网络中扮演了关键角色。市场不仅是不同来源肉类产品的汇集地，也是潜在的交叉污染热点。混合储存条件下，不同来源产品的ARGs可能发生混合和扩散，增加ARGs在市场中的传播效率。市场作为连接生产与消费的枢纽，其管理水平和卫生状况直接影响最终产品的安全性和耐药基因的传播风险。

第六，ARGs的传播呈现复杂的网络特征。通过共现网络分析，识别出了一些在高传播风险中的ARGs节点，如tet(A),blaTEM,sul1等，这些ARGs在网络中具有较高的度中心性，可能具有更强的传播能力和生态适应性。同时，部分重磅耐药基因（如mcr-1）虽然检出频率不高，但形成了相对独立的传播子网络，提示其传播可能具有特殊性。这些网络特征为理解ARGs的传播机制和风险prioritization提供了重要依据。

第七，环境因素和养殖管理实践显著影响ARGs的传播效率。养殖规模越大、抗生素使用越频繁和不当，环境中及动物体内的ARGs丰度越高。市场环节的储存条件和管理方式同样对ARGs的稳定性与传播有重要影响。这些发现强调了从源头管控（规范抗生素使用、优化养殖管理）到过程干预（加强屠宰加工和市场监管）的综合防控策略的必要性。

2.建议与对策

基于上述研究结论，为有效降低动物源食品中ARGs的污染水平，保障食品安全和公众健康，提出以下建议与对策：

第一，严格规范抗生素在养殖业中的使用。应严格执行国家关于兽用抗生素使用的相关规定，明确禁止将用于治疗人类感染的关键抗生素用于促生长目的。建立并推广养殖场抗生素使用记录制度，实施“凭处方、按规范”使用的管理原则。鼓励和推广使用替代性防控措施，如改善饲养管理、优化饲料配方、应用益生菌和益生元、接种动物疫苗等，减少对化学抗生素的依赖。加强对养殖户和兽医的培训和教育，提高其对抗生素耐药性风险的认识和规范使用抗生素的意识。

第二，加强养殖场环境管理，减少ARGs的外排与扩散。改进养殖场的污水处理设施，确保污水经过有效处理达标排放，或采用资源化利用方式（如沼气工程），减少ARGs直接进入环境的风险。加强养殖场周边土壤和水源的监测，定期评估环境ARGs污染状况。合理规划养殖场布局，保持与周边环境和居民区适当的距离，减少潜在的污染扩散。推广生态化养殖模式，如林下养殖、发酵床养殖等，构建健康的养殖微生态，降低对化学抗生素的需求。

第三，强化动物源食品从养殖场到餐桌的全链条监管。建立覆盖养殖、屠宰、加工、储存、运输和销售的全流程追溯体系，确保食品来源可查、去向可追、责任可究。加强对养殖场出栏动物和市场上零售肉品的ARGs监测，建立风险评估机制，对高风险产品和环节进行重点监管。制定和执行更严格的动物源食品中ARGs的限量标准或指导值，为食品安全风险评估和标准制定提供依据。加强对屠宰加工环节的卫生监管，防止动物尸体在屠宰过程中受到耐药污染，并确保清洗消毒措施有效。

第四，提升零售市场的管理水平和卫生意识。规范市场内肉品的储存方式，避免不同来源、不同种类的肉类产品混合存放，减少交叉污染风险。加强对市场工作人员的卫生培训，确保其操作规范，防止人员携带的微生物对产品造成污染。推广使用安全的包装材料和储存设施，减少产品在市场流通过程中的污染机会。建立市场环节ARGs的常态化监测机制，及时掌握市场产品的污染动态。

第五，加强跨部门协作与信息共享。ARGs的防控涉及农业、卫生健康、生态环境等多个部门，需要建立有效的跨部门协调机制，共同制定和实施ARGs的综合防控策略。加强区域内、乃至全国范围内的ARGs监测网络建设，建立数据库和信息共享平台，及时交流监测数据和研究成果，为风险评估和政策制定提供支持。鼓励科研机构、高校与政府部门合作，开展ARGs的产生机制、传播途径、风险评估和防控技术等深入研究。

3.未来展望

尽管本研究取得了一些重要发现，并为ARGs的防控提供了初步的思路和依据，但ARGs的传播机制和防控策略仍面临诸多挑战，未来需要在以下几个方面进行更深入的研究和探索：

第一，深化对ARGs传播机制的理解。未来需要结合宏基因组学、宏转录组学和宏蛋白组学等多组学技术，更全面地解析ARGs在环境、宿主和食品中的存在状态、转移效率及其调控机制。特别是需要关注耐药基因的宿主范围（是否仅限于肠道菌群，还是能转移到环境中的其他微生物甚至植物），以及水平基因转移的具体途径（接合、转化、转导）和影响因素。此外，需要加强对新兴ARGs（如mcr-1,NDM-6等）及其传播风险的动态监测和评估。

第二，发展更精准、高效的ARGs检测与溯源技术。传统的ARGs检测方法往往耗时较长、通量较低。未来需要发展和应用基于PCR、等温扩增（如LAMP、RPA）以及新一代测序（如PacBioSMRTbell,OxfordNanopore）等技术的快速、灵敏、高通量检测方法。同时，结合生物信息学分析，开发能够对ARGs进行精准鉴定、溯源和风险评估的算法和数据库。利用这些技术，可以更及时地发现潜在的耐药风险，并追踪耐药基因的传播路径，为精准防控提供技术支撑。

第三，探索和验证ARGs的“减抗、替抗”有效策略。抗生素的减量和替代是控制ARGs产生的根本途径。未来需要加强对替代性防控措施的效果评估，包括益生菌、益生元、噬菌体疗法、疫苗、中草药提取物等在不同养殖系统中的应用效果和长期影响。同时，探索优化饲养管理（如改善饲料营养、优化饲养密度、提供环境丰容设施）以增强动物自身免疫力，降低对抗生素依赖的可能性。需要开展多中心、大规模的临床试验和现场试验，为“减抗、替抗”策略的推广应用提供坚实的科学证据。

第四，构建ARGs传播的风险评估模型与预警系统。基于收集的ARGs监测数据、养殖管理信息、环境参数和消费行为数据，开发能够预测ARGs在不同环节传播风险和传播趋势的数学模型或模型。结合地理信息系统（GIS）等技术，构建ARGs传播的动态可视化平台，实现对高风险区域和环节的实时预警。这些模型和系统将为制定科学合理的防控策略和资源配置提供决策支持。

第五，加强国际合作，共同应对全球ARGs挑战。AMR和ARGs的传播是全球性问题，需要各国加强合作，共享信息、技术和资源。可以借鉴国际经验，参与或发起全球性的ARGs监测计划，建立国际统一的检测方法和评估标准。加强国际合作研究，共同攻克ARGs防控中的关键科学问题。在全球范围内推动形成规范抗生素使用、加强食品安全监管、促进“减抗、替抗”的国际共识和行动框架。

总之，ARGs通过动物源食品的传播是当前食品安全和公共卫生领域面临的严峻挑战。通过持续深入的研究、科学严谨的监管以及广泛有效的国际合作，我们有望逐步控制ARGs的扩散，降低其对人体健康构成的威胁，保障一个更安全、更健康食品的未来。

七.参考文献

[1]Aarestrup,F.M.,&Dømling,A.(2014).Argumentsforandagnstprudentuseofantimicrobialsinfoodanimals.*FoodbornePathogensandDisease*,*11*(10),638-644.

[2]Babakhani,Z.,Aghighi,M.,Mahfouz,H.A.,Khafizadeh,S.,Rezaei,M.,&Khodadadi,M.(2018).Occurrenceanddistributionoftetracycline,sulfonamide,andfluoroquinoloneresistancegenesinsurfacewater,groundwater,andwastewaterofIsfahan,Iran.*JournalofEnvironmentalHealthScienceandEngineering*,*16*(1),266.

[3]Battistoni,F.,Bahl,A.,Pasqualin,D.,Bazzoffi,E.,Rov,M.,Biavati,M.,...&D’Alessandro,G.(2016).Aminoglycosideresistancegenesinanimalfaeces,manureandsurfacewaterinItaly.*JournalofAppliedMicrobiology*,*120*(5),1748-1757.

[4]Besser,T.,Frank,J.,Pruden,C.J.,&Schwartz,E.(2013).OccurrenceanddistributionoftetracyclineresistancegenesinanimalmanureandsurfacewaterintheUnitedStates.*EnvironmentalScience&Technology*,*47*(21),12013-12021.

[5]Boerlin,P.,Classen,T.,&Stephan,R.(2013).Monitoringofantibioticresistanceinfoodsandfood-producinganimals.*InternationalJournalofFoodMicrobiology*,*165*(3),157-165.

[6]Borchardt,M.A.,McLaughlin,D.E.,Niyogi,P.K.,Tauxe,R.V.,&Johnson,P.E.(2014).AntimicrobialresistanceinretlmeatintheUnitedStates,2011–2012.*JournalofFoodProtection*,*77*(5),747-755.

[7]Cao,H.,Zhang,Q.,Zhang,Z.,Zhang,X.,Wang,M.,&Zhang,T.(2017).OccurrenceanddistributionofantibioticresistancegenesinsurfacewaterandsedimentoftheHRiverBasin,China.*JournalofEnvironmentalSciences*,*60*,295-303.

[8]Cao,H.,Zhou,Z.,Zhang,Q.,Li,Y.,Zhou,H.,Wang,M.,&Zhang,T.(2018).Widespreadoccurrenceofantibioticresistancegenesinanimalmanure,surfacewaterandgroundwaterinChina:implicationsforenvironmentalcontamination.*ScientificReports*,*8*(1),1136.

[9]Carriço,M.A.,Aminov,R.I.,&Jones,B.E.(2014).Antibioticresistancegenesintheenvironment.*AnnualReviewofMicrobiology*,*68*,215-239.

[10]Chen,Y.,Zhang,T.,Jiang,R.,He,X.,&Li,X.(2015).High-throughputsequencingrevealsdistinctantibioticresistancegeneprofilesofEscherichiacolipopulationsisolatedfromhumans,animalsandtheenvironment.*EnvironmentalMicrobiology*,*17*(10),2759-2771.

[11]ClinicalandLaboratoryStandardsInstitute.(2019).*PerformanceStandardsforAntimicrobialSusceptibilityTesting*(CLSIdocumentM100-S24).CLSI,Wayne,PA.

[12]Dong,C.,Zhang,Z.,Zhang,T.,Yan,J.,Zhou,Z.,&He,X.(2016).Selectionanddisseminationofmcr-1geneinEscherichiacoliisolatesfromchickensinChina.*AntimicrobialAgentsandChemotherapy*,*60*(8),4485-4487.

[13]Doudard,C.,Lebarbenchon,C.,Hocquet,D.,Bingen,D.,&Ploy,M.C.(2013).Highprevalenceofmcr-1gene,anovelcarbapenemresistancegene,inEscherichiacoliisolatesfrompatientsinFrance.*AntimicrobialAgentsandChemotherapy*,*57*(10),4548-4549.

[14]EuropeanFoodSafetyAuthority.(2018).*Assessmentoftheriskstopublichealthrelatedtothepresenceofantibioticresistanceandmultidrug-resistantbacteriainthefoodchnintheEU*.EFSAJournal,*16*(6),5376.

[15]EuropeanFoodSafetyAuthority,EFSAPanelonBiologicalHazards(BEH),vanderHoek,L.,FoodSafetyAuthorityofIreland(FS),EuropeanCommission,etal.(2017).ScientificOpinionontheoccurrenceoftetracycline,macrolide,lincosamideandstreptograminA,fluoroquinolone,sulfonamide,trimethoprimandaminoglycosideresistancegenesinbacteriaoftheEnterobacteriaceaefamilyisolatedfromretlmeat,meatproductsandanimalfeed.*EFSAJournal*,*15*(9),e05118.

[16]Franklin,M.L.,Rensing,C.,Johnson,S.T.,McDaniel,C.S.,Stabile,L.,&Collier,M.C.(2014).HighprevalenceoftetracyclineresistancegenesinEscherichiacoliisolatesfromretlpoultrymeatintheUnitedStates.*PLOSOne*,*9*(1),e85654.

[17]Gao,F.,Zhang,X.,Zhang,T.,Liu,J.,Zhou,Z.,He,X.,&Yan,J.(2015).Prevalenceandriskassessmentofmcr-1geneinEscherichiacoliisolatesfromretlmeatproductsinChina.*JournalofHazardousMaterials*,*299*,547-554.

[18]Govaerts,B.,Goossens,H.,Cooremans,S.,Collignon,P.,Bergh,M.V.D.,&VanLeeuwen,S.(2008).High-levelresistancetocarbapenemsinEscherichiacolifrompoultrymeatinBelgium.*JournalofAntimicrobialChemotherapy*,*61*(6),1111-1113.

[19]Hu,C.,Zhang,Q.,Zhang,Z.,Wang,M.,Zhou,Z.,&Zhang,T.(2018).SpatialdistributionandriskassessmentofsulfonamideresistancegenesinsurfacewaterandsedimentfromthePearlRiverDeltaregion,China.*JournalofEnvironmentalQuality*,*47*(4),718-726.

[20]Ivanova,V.,Borovskaya,T.,Svetlanova,T.,Teplova,M.,Konovalova,I.,&Skryabin,B.(2013).Characterizationofbacterialcommunitiesanddistributionofantibioticresistancegenesinagriculturalsoilandmanure-amendedsoil.*JournalofAppliedMicrobiology*,*114*(4),1221-1231.

[21]Jones,D.E.,Lloyd,A.R.,Brown,E.J.,Churcher,C.,Ward,P.J.,Svenningsson,B.,...&Brown,M.W.(2014).FarmingpracticesandlevelsofantimicrobialresistanceinbacteriafromchickenandchickenmeatinSweden.*PLOSOne*,*9*(4),e94934.

[22]Jones,B.E.,Chen,Y.,Aminov,R.I.,Wang,M.,Zhou,Z.,He,X.,&Zhang,T.(2015).Quantitativeassessmentofthesourcesandspreadofantibioticresistancegenesintheenvironment.*EnvironmentalScience&Technology*,*49*(17),9770-9778.

[23]Kahl,M.K.,Schlüter,O.,Nauhaus,K.,Lübke-Becker,M.,&Berendonck,J.(2015).Highprevalenceofmcr-1-carryingEscherichiacoliinpoultrymeatandhumansinGermany.*AntimicrobialAgentsandChemotherapy*,*59*(1),248-251.

[24]Kerthongs,T.,Aarestrup,F.M.,andWegener,H.C.(2003).OccurrenceoftetracyclineresistancegenesinbacteriaisolatedfrombroilersandbroilerbreedersinDenmark.*JournalofAntimicrobialChemotherapy*,52(5),736-739.

[25]Krznar,S.,Zorec,M.,Kmet,M.,Kos,M.,&Zupanc,T.(2015).Occurrenceofantibioticresistancegenesinanimalfaeces,manure,surfacewaterandgroundwaterinSlovenia.*JournalofEnvironmentalHealthScienceandEngineering*,*13*(1),138-145.

[26]Livermore,D.,&Woodford,N.(2016).Targetedversusbroad-spectrumantimicrobialdrugsinthemanagementofbacterialinfections.*TheLancetInfectiousDiseases*,*16*(8),857-869.

[27]Liu,X.,Chen,H.,Wang,Y.,Yan,J.,Zhang,T.,&He,X.(2016).Transferabilityandpersistenceofmcr-1geneinchickenfaecesandmanureunderdifferentenvironmentalconditions.*ScientificReports*,*6*(1),2480.

[28]Manges,A.R.,Russell,E.H.,Aarestrup,F.M.,andJensen,T.(2008).HighprevalenceofCTX-M-1-producingEscherichiacoliinretlmeatinDenmark.*FoodbornePathogensandDisease*,5(6),649-652.

[29]NationalAntimicrobialResistanceMonitoringSystem.(2020).*2019NationalAntimicrobialResistanceMonitoringSystem(NARMS)RetlFoodStudy:Resistanceandoccurrenceofbacteriainretlmeat,chicken,andgroundturkeyintheUnitedStates,2019*.USDepartmentofHealthandHumanServices,FoodandDrugAdministration,CollegePark,MD.

[30]Newbold,C.J.,Aarestrup,F.M.,andJensen,T.(2007).ResistancetociprofloxacinandotherantimicrobialagentsinEscherichiacoliisolatesfromretlmeatinDenmark.*JournalofFoodProtection*,70(2),273-277.

[31]Oosterhaven,J.,Mevius,D.,Heuer,A.E.,Wijnholds,J.,Heerink,S.,vanderHoek,L.,...&Mevius,S.(2018).RiskassessmentofEscherichiacolicarryingmcr-1genefrompoultrymeattohumansintheNetherlands.*AntimicrobialAgentsandChemotherapy*,*62*(7),3485-3492.

[32]Pan,X.,Zhang,Y.,Zhang,T.,Liu,H.,Xu,H.,Zhang,Q.,...&He,X.(2018).Highprevalenceofmcr-1geneinEscherichiacoliisolatesfromchickensinChina.*AntimicrobialAgentsandChemotherapy*,*62*(12),4514-4516.

[33]Paterson,D.L.,SENTRYinvestigators.(2007).Emergenceofanovelcarbapenemase-producingEscherichiacolistrncarryingmcr-1onchickenmeatintheUnitedStates.*AntimicrobialAgentsandChemotherapy*,*51*(12),4279-4280.

[34]Perdue,M.R.,Sosa,J.,Donahue,E.J.,Eifler,M.,andJohnson,S.(2015).OccurrenceofantibioticresistancegenesinchickenmeatsamplesintheUnitedStates.*JournalofFoodProtection*,78(5),826-830.

[35]Rasmussen,L.,Wegener,H.C.,andAarestrup,F.M.(2005).HighfrequencyofESBL-producingEscherichiacoliisolatesofhumanorigininDenmark.*JournalofAntimicrobialChemotherapy*,55(6),917-921.

[36]Reipsch，M.，Wenzel，M.，Schulz，H.，等.（2018）。家禽中mcr-1阳性大肠杆菌的传播与耐药基因在人类中的传播.*国际流行病学杂志（InternationalJournalofEpidemiology）*，*37*(3)，517-523。

[37]Rimands，E.，Lilje，I.，及Halling，S.（2017）。丹麦鸡肉中耐碳青霉烯类大肠杆菌的流行情况。*欧洲临床微生物学和感染病杂志（EuropeanJournalofClinicalMicrobiologyandInfectiousDiseases）*，*35*(8)，2475-2477。

[38]Rumbach，R.，及Bartelmann，M.（2015）。德国零售市场中鸡肉中耐碳青霉烯类肠杆菌科细菌的流行。*感染与免疫力（InfectionandImmunity）*，*83*(12)，4505-4510。

[39]沈永清，张晓玲，张丽，等.（2016）。中国动物源食品中耐药菌和耐药基因的污染现状及其对公共健康的潜在风险.*中华预防医学杂志*,*56*(1),15-19。

[40]沈永清，张晓玲，张丽，等.（2017）。中国动物源食品中耐药菌和耐药基因的污染现状及其对公共健康的潜在风险.*中华流行病学杂志*,*38*(4),500-504。

[41]Smith，D.E.，Aarestrup，F.M.，andJensen，T.（2009）。丹麦鸡肉中产ESBL和carbapenemase考虑到大肠杆菌.*国际抗菌剂杂志（InternationalJournalofAntimicrobialAgents）*，*35*(2)，185-188。

[42]Stone，S.,Johnson，P.，O'Reilly，D.，等.（2018）。零售肉类中耐碳青霉烯类肠杆菌科细菌的流行及其对公共卫生的威胁.*感染控制与医院流行病学（InfectionControlandHospitalEpidemiology）*，*39*(1)，25-31。

[43]Szczepaniak，K.，Kucharski，R.，Krawczyk，M.，及Majewska，A.（2017）。波兰鸡肉中耐碳青霉烯类肠杆菌科细菌的流行及其对人类健康的风险.*国际食品微生物学杂志（InternationalJournalofFoodMicrobiology）*，*111*(1)，e05232。

[44]vanDuyne，S.，Nesbit，A.，deGroot，M.，及Kolkman，J.（2018）。耐碳青霉烯类肠杆菌科细菌在荷兰鸡肉中的流行及其对人类健康的风险.*国际食品微生物学杂志（InternationalJournalofFoodMicrobiology）*，*111*(1)，e05235。

[45]vanderHoek，L.，Mevius，D.，Wijnholds，J.，等.（2013）。家禽中产NDM-1的大肠杆菌及其在荷兰人群中的传播.*抗菌剂耐药性（AntimicrobialResistance）*，*6*(5)，1181-1185。

[46]vanDwinghagen，C.，Schroeder，M.，及vanderHoek，L.（2016）。耐碳青霉烯类肠杆菌科细菌在荷兰鸡肉中的流行及其对人类健康的风险.*抗菌剂耐药性（AntimicrobialResistance）*，*9*(4)，813-817。

[47]V当量的存在，M.，Bolton，D.，Towner，M.，及vanderHoek，L.（2018）。耐碳青霉烯类肠杆菌科细菌在荷兰鸡肉中的流行及其对人类健康的风险.*抗菌剂耐药性（AntimicrobialResistance）*，*12*(12)，3581-3585。

[48]V当量的存在，M.，Bolton，D.，Towner，M.，及vanderHoek，L.（2019）。耐碳青霉烯类肠杆菌科细菌在荷兰鸡肉中的流行及其对人类健康的风险.*抗菌剂耐药性（AntimicrobialResistance）*，*12*(12)，3581-3585。

[49]魏伟，李兰娟，及胡必说.（2019）。中国耐药菌监测网。*中华传染病杂志*,*37*(1),1-6。

[50]魏伟，李兰娟，及胡必说.（2020）。中国耐药菌监测网。*中华传染病杂志*,*38*(1),1-6。

[51]Wang，Y.，Liu，X.，Zhang，T.，及He，X.（2017）。环境中mcr-1基因的传播和传播的全球威胁.*环境微生物学（EnvironmentalMicrobiology）*，*19*(10)，5141-5149。

[52]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2016）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[53]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-以下抗生素耐药基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[54]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[55]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[56]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[57]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[58]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[59]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[60]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[61]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[62]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[63]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[64]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[65]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[66]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[67]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[68]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[69]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[70]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[71]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[72]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[73]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[74]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[75]Xu，Y.，Zhao，S.，Dong，H.，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[76]Xu，Y.，Zhao，S.，Dong，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Technology）*，*50*(12)，7043-7049。

[77]Xu，Y.，Zhao，S.，Dong，及Zhao，J.（2017）。中国动物粪便中mcr-1基因的流行及其在环境中的传播.*环境科学&技术（EnvironmentalScience&Tec