抗生素耐药基因传播X环境样品论文

抗生素耐药基因传播X环境样品论文一.摘要

抗生素耐药基因（ARGs）的传播已成为全球公共卫生领域面临的重要挑战，其广泛存在于各类环境样品中，对人类健康和生态系统稳定构成潜在威胁。本研究聚焦于典型环境样品中ARGs的传播特征及其影响因素，选取了工业废水、农业土壤、医院污水及城市污水处理厂出水和河流沉积物等代表性样品进行系统分析。研究采用高通量测序技术结合生物信息学分析，对样品中的ARGs种类、丰度及遗传背景进行深度解析。结果表明，工业废水和医院污水中ARGs的种类和丰度显著高于其他样品，其中多重耐药基因如NDM-1、KPC-2等检出率较高，提示人类活动产生的污染物是ARGs的重要来源。农业土壤样品中，与农业抗生素使用相关的ARGs（如tetA、ermB）检出率显著增加，与化肥和农药的长期施用密切相关。城市污水处理厂出水及河流沉积物中ARGs的多样性呈现梯度变化，靠近排污口区域ARGs丰度明显升高，表明污水处理设施和自然水体的交互作用是ARGs传播的关键环节。遗传背景分析显示，ARGs常与移动遗传元件（MGEs）如整合子、转座子等共现，加速了ARGs在不同环境介质间的转移。研究还揭示了环境因素（如pH值、有机质含量、重金属浓度）对ARGs稳定性和传播效率的调控作用。结论表明，人类活动产生的污染源是ARGs环境传播的主导因素，污水处理系统和自然水体的交互区域是ARGs扩散的关键节点。通过构建ARGs传播风险评估模型，为制定有效的环境保护和感染控制策略提供了科学依据。本研究不仅揭示了ARGs在复杂环境介质中的传播规律，也为未来开发基于环境监测的耐药性预警系统奠定了理论基础。

二.关键词

抗生素耐药基因；环境样品；高通量测序；移动遗传元件；污水处理厂；河流沉积物；多重耐药性；环境传播；风险评估模型

三.引言

抗生素的发现与应用无疑是20世纪医学领域最重大的成就之一，极大地提高了人类对抗感染性疾病的抵抗能力，显著降低了相关疾病的死亡率。然而，随着抗生素的广泛和滥用，抗生素耐药性问题（AntibioticResistance,AMR）已迅速发展成为全球性的公共卫生危机。据世界卫生（WHO）报告，AMR每年导致全球数十万人死亡，且形势日益严峻，部分细菌已对现有几乎所有抗生素产生耐药性，使得感染治疗变得异常困难。抗生素耐药基因（AntibioticResistanceGenes,ARGs）作为耐药性的功能载体，是驱动AMR传播的关键分子。与传统观点认为耐药菌是AMR传播的主要源头，越来越多的证据表明，即使在不具备活性抗生素的环境样品中，如土壤、水体、沉积物等，ARGs也以游离或整合于宿主基因组的形式广泛存在，并通过水平基因转移（HorizontalGeneTransfer,HGT）等机制在不同微生物群落间传播，构成了“环境库”（EnvironmentalReservoir），对人类健康构成潜在威胁。

近年来，环境样品中ARGs的污染状况及其生态风险受到广泛关注。研究表明，医院污水、生活污水、农业runoff（农田排水）、工业废水、养殖场废水等是人类活动产生的ARGs进入环境的主要途径。医院污水中由于含有大量耐药菌和抗生素，其ARGs浓度通常远高于其他类型污水。城市污水处理厂（WastewaterTreatmentPlants,WWTPs）作为处理大量人类废水的场所，不仅是ARGs的巨大汇（sink），同时也是潜在的源（source），处理过程中可能发生微生物群落结构改变、抗生素浓度波动以及生物化学过程（如厌氧消化），这些都可能促进ARGs的释放和水平转移。农业环境中，为了提高作物产量和防治病虫害，抗生素（如四环素类、大环内酯类）被广泛用作兽药和农药，导致土壤和动物粪便中ARGs检出率显著升高，并通过农产品、水体和空气等途径扩散。工业废水中的重金属和有机污染物可能通过影响微生物群落结构和功能，间接调控ARGs的传播。自然水体和沉积物中的ARGs来源复杂，可能包括上游输入、沉积物释放以及水生微生物的自身合成。值得注意的是，河流沉积物不仅吸附和储存来自陆地的污染物和ARGs，还可能成为厌氧条件下微生物水平基因转移的高效场所，因此对沉积物中ARGs的研究对于理解水生生态系统中AMR的长期风险具有重要意义。

尽管对环境样品中ARGs的污染状况已有诸多研究，但现有研究多集中于特定类型环境样品或特定ARGs的检测，对于不同环境样品中ARGs的组成特征、丰度水平、空间分布格局及其与人类活动、环境因子和微生物群落结构的关联性尚缺乏系统性的比较和深入探讨。特别是，如何准确评估ARGs在不同环境介质间的传播路径和风险，以及如何揭示环境因素对ARGs稳定性和转移效率的具体作用机制，仍然是当前研究面临的重要挑战。此外，环境ARGs的“环境库”作用机制，即ARGs如何在环境中维持存在、扩散和在特定条件下重新进入病原菌的生态动力学过程，仍需进一步阐明。这些知识的缺乏限制了我们对AMR传播规律的科学认知，也不利于制定有效的环境保护和感染控制策略。因此，本研究旨在通过对典型环境样品（包括工业废水、农业土壤、医院污水、城市污水处理厂出水和河流沉积物）进行系统性的ARGs分析，明确不同样品中ARGs的群落特征和丰度差异，探究ARGs与环境介质理化性质、微生物群落结构以及人类活动强度之间的关联性，并评估环境样品作为ARGs传播媒介的潜在风险。本研究的意义在于：首先，通过多类型环境样品的比较研究，可以更全面地揭示ARGs在自然环境中的分布格局和传播特征，为识别AMR传播的关键路径和风险源提供科学依据；其次，通过分析ARGs的遗传背景，特别是与移动遗传元件（MGEs）的共现关系，可以深入理解ARGs在环境中的传播机制；再次，研究环境因子对ARGs的影响，有助于揭示ARGs在环境中的生态动力学过程，为预测和控制ARGs的扩散提供理论支持；最后，本研究结果将为制定基于环境监测的AMR预警系统、优化污水处理工艺、加强农业抗生素管理以及保护生态环境健康提供重要的参考数据和科学建议。基于上述背景，本研究提出以下核心研究问题：不同类型环境样品中ARGs的群落组成和丰度是否存在显著差异？这些差异与环境介质特征、微生物群落结构以及人类活动因素之间存在怎样的关联？环境因素如何影响ARGs的稳定性和传播效率？环境样品（特别是污水处理厂和河流沉积物）在ARGs的传播网络中扮演着怎样的角色？通过回答这些问题，本研究期望能够为深入理解环境ARGs的传播规律和风险机制提供新的见解，并为应对全球AMR挑战提供科学支撑。

四.文献综述

抗生素耐药基因（ARGs）作为细菌耐药性的遗传基础，其环境分布和传播问题已成为全球公共卫生和生态安全领域的研究热点。现有研究表明，ARGs广泛存在于各种自然和人工环境中，包括土壤、水体、沉积物、空气以及废弃物等，形成了一个庞大的“环境库”，对人类健康和生态系统稳定构成潜在威胁。对环境ARGs的研究主要集中在其来源、分布、传播途径以及生态风险等方面。

在ARGs的来源方面，人类活动是ARGs进入环境的主要驱动因素。医院污水和污水处理厂（WWTPs）是ARGs的重要输入源。研究发现，未经处理或处理不充分的医院污水中含有高浓度的ARGs和耐药菌，其中检出的ARGs种类多，丰度高，特别是与临床常用抗生素（如β-内酰胺类、大环内酯类）相关的基因。WWTPs作为大量人类废水的处理场所，不仅是ARGs的巨大汇，也是潜在的源。处理过程中，微生物群落结构发生剧烈变化，抗生素浓度波动，以及厌氧消化等生物化学过程，都可能促进ARGs的释放和水平基因转移（HGT）。研究表明，WWTPs的出水、污泥以及生物膜中均检测到丰富的ARGs，且部分WWTPs出水的ARGs浓度甚至高于原污水，表明处理过程可能并未有效去除ARGs，反而可能通过污泥扩散或生物膜结构促进了ARGs的释放和传播。农业活动也是ARGs环境传播的重要途径。为了提高农作物产量和防治动植物疾病，抗生素（如四环素、大环内酯类、喹诺酮类）被广泛用作兽药和农药。研究发现，施用抗生素的农田土壤、动物粪便、农作物以及灌溉水中均检测到相应的ARGs，表明抗生素可通过农业runoff、土壤吸附、农产品以及动物排泄等途径进入环境。此外，工业废水、养殖场废水以及垃圾渗滤液等也是ARGs的重要来源，其中可能含有特定行业使用的抗生素或与抗生素结构类似的化合物，导致ARGs的多样性增加。

在ARGs的分布方面，环境样品中ARGs的丰度和种类存在显著的空间异质性。一般来说，城市化和人类活动干扰强烈的区域，如工业区、农业区以及WWTPs周边，ARGs的丰度通常较高。研究发现，WWTPs出水和沉积物中的ARGs丰度显著高于远离污染源的区域，且在排污口附近形成高浓度梯度。河流沉积物作为ARGs的储存库和潜在的转移场所，其ARGs含量与上游污染输入、沉积物释放以及水生生物活动等因素密切相关。一些研究表明，沉积物中的ARGs可以通过再悬浮过程释放到水体中，或在水生生物体内转移。土壤中的ARGs分布受土壤类型、有机质含量、pH值以及微生物群落结构等多种因素影响。例如，施用抗生素和化肥的农田土壤中，特定ARGs（如tet类、erm类）的丰度显著升高。值得注意的是，即使在看似“清洁”的自然环境中，如偏远地区的土壤和水体中，也检测到了ARGs，表明ARGs具有全球性的分布特征，并通过大气传输、洋流等长距离途径进行扩散。

在ARGs的传播途径方面，HGT被认为是ARGs在微生物群落间传播的主要机制。与传统的垂直基因转移相比，HGT使得ARGs能够在不同物种、不同环境之间快速传播，极大地加速了耐药性的扩散。移动遗传元件（MobileGeneticElements,MGEs），如质粒、整合子、转座子等，是介导ARGsHGT的关键分子工具。研究发现，ARGs常常与MGEs共现于各种环境样品中，且MGEs的存在显著增加了ARGs的拷贝数和传播能力。例如，在WWTPs中，携带ARGs的整合子和大肠杆菌-肠杆菌整合子（ICEs）非常普遍，它们能够捕获和转移多种ARGs，促进了耐药基因的聚合和传播。此外，噬菌体介导的基因转移也是ARGs传播的重要途径之一。噬菌体作为细菌的天敌，在捕食过程中可以携带宿主细菌的ARGs，并将其转移到其他细菌体内。环境因素，如重金属污染、有机污染物以及抗生素压力等，可以影响微生物群落结构和功能，进而影响MGEs的活性，从而调控ARGs的传播效率。例如，研究表明，重金属胁迫可以诱导某些细菌产生更多的MGEs，增加ARGs的传播风险。

在ARGs的生态风险方面，环境ARGs已被证实可以通过多种途径重新进入人类和动物体内，构成潜在的感染风险。饮用水和农产品是ARGs进入人体的主要途径。研究发现，部分地区饮用水源和灌溉水中检测到ARGs，表明通过饮用水摄入是ARGs暴露的重要途径。农产品中的ARGs可能来源于土壤污染、灌溉水污染以及动物粪便污染等，通过食物链传递进入人体。此外，手接触受污染的环境介质（如土壤、水体）后经口摄入，以及与受污染动物接触等，也是ARGs暴露的潜在途径。研究表明，暴露于环境ARGs的人群，其肠道菌群中的耐药性可能增加，甚至可能发生耐药菌感染。虽然目前直接证据表明环境ARGs导致人类感染的案例还比较有限，但其潜在的生态风险不容忽视。此外，ARGs的存在不仅威胁人类健康，也可能对生态环境造成负面影响。例如，ARGs可能通过影响微生物群落结构和功能，改变生态系统的物质循环和能量流动。某些ARGs可能赋予微生物新的生存优势，导致某些耐药菌在环境中过度增殖，从而抑制其他有益微生物的生长，破坏生态平衡。

尽管对环境ARGs的研究取得了显著进展，但仍存在一些研究空白和争议点。首先，关于环境ARGs的“环境库”作用的机制尚不明确。虽然ARGs在环境中广泛存在，但其能否在环境中长期维持、如何维持以及何时重新进入病原菌，这些生态动力学过程仍需深入研究。其次，不同环境介质中ARGs的传播路径和风险评估模型尚不完善。目前的研究多集中于特定类型环境样品或特定ARGs的检测，对于不同环境介质间ARGs的传播网络和风险传递机制缺乏系统性的研究。如何建立基于多介质监测的ARGs传播风险评估模型，仍然是当前研究面临的挑战。再次，关于环境因素对ARGs稳定性和转移效率的具体作用机制，特别是非抗生素压力环境下的ARGs传播机制，需要进一步阐明。此外，针对环境ARGs的防控策略研究相对滞后。虽然已有研究表明某些措施（如优化污水处理工艺、减少抗生素使用、加强环境监测等）可以降低环境ARGs的污染水平，但如何制定基于成本效益和环境影响的综合防控策略，仍需更多实证研究的支持。最后，关于ARGs在极端环境（如高温、高盐、高辐射等）下的生存和传播规律，以及气候变化对ARGs环境分布和传播的影响，也是需要关注的研究方向。总之，深入理解环境ARGs的传播规律和风险机制，对于制定有效的防控策略、保障人类健康和生态安全具有重要意义。未来需要加强多学科交叉研究，整合宏基因组学、代谢组学、环境化学、生态学和公共卫生等领域的知识和方法，以期更全面地揭示环境ARGs的生态行为和风险机制。

五.正文

为系统探究典型环境样品中抗生素耐药基因（ARGs）的传播特征及其影响因素，本研究选取了工业废水、农业土壤、医院污水、城市污水处理厂（WWTP）出水和河流沉积物五种代表性样品类型，采用高通量测序技术结合生物信息学分析，对样品中ARGs的种类、丰度、遗传背景及其与环境因子、微生物群落结构的关联性进行了深入研究。研究旨在揭示不同环境介质中ARGs的组成差异、传播规律及其潜在风险，为理解和控制ARGs的环境传播提供科学依据。

1.研究区域与样品采集

研究区域涵盖工业发达地区、农业密集区以及城市流域，以反映不同人类活动强度下的环境ARGs污染状况。具体样品采集点信息如下：

-工业废水：某化工厂排放口下游500米处，采集未经处理的生产废水。

-农业土壤：某集约化养殖场周边农田，距离养殖场100米，采集表层（0-20cm）土壤样品。

-医院污水：某三甲医院污水处理站进水口，采集未经处理的原污水。

-城市污水处理厂出水：某大型城市WWTP的最终出水口，采集处理后的排放水。

-河流沉积物：某城市主要河流下游沉积物，采集位于排污口下游1公里和5公里处的沉积物样品。

样品采集时间为2023年夏季，每个样品采集至少三个重复，混合后分装于无菌袋中，立即冷藏保存，并尽快进行实验室分析。

2.样品预处理与ARGs提取

所有样品均经过预处理以去除杂质并提取总DNA。具体步骤如下：

-工业废水：通过0.22μm滤膜过滤去除悬浮物，使用试剂盒（如MoBioPowerWaterDNAKit）提取滤液中的总DNA。

-农业土壤：风干后过筛（孔径0.25mm），取适量样品加入裂解缓冲液，使用试剂盒提取土壤总DNA。

-医院污水：同样通过0.22μm滤膜过滤，使用试剂盒提取滤液中的总DNA。

-城市污水处理厂出水：通过0.45μm滤膜过滤，使用试剂盒提取滤液中的总DNA。

-河流沉积物：风干后过筛（孔径0.25mm），取适量样品加入裂解缓冲液，使用试剂盒提取沉积物总DNA。

DNA提取后，使用纳米分光光度计（如NanoDrop）检测DNA浓度和纯度，合格的DNA样品用于后续高通量测序。

3.高通量测序与生物信息学分析

本研究采用IlluminaHiSeq平台进行高通量测序，具体流程如下：

-文库构建：将提取的总DNA进行文库构建，包括DNA片段化、末端修复、加A尾、接头连接等步骤。

-测序：采用双端测序（IlluminaHiSeq3000），产生pred-endreads。

-数据质控与过滤：使用Trimmomatic进行数据质控，去除低质量reads和接头序列，得到高质量的cleandata。

-原始序列分析：将cleandata进行物种注释，初步筛选出细菌和古菌的16SrRNA基因序列。

-ARGs检测与定量：使用MetaSPAdes进行宏基因组拼接，然后使用BGISEQ-ML等工具进行ARGs检测和定量，得到各样品中ARGs的种类和丰度信息。

-MGEs分析：使用CRISPRdb等数据库，检测样品中是否存在MGEs，并分析其与ARGs的共现关系。

-环境因子分析：收集样品采集点的环境参数，包括pH值、电导率、有机质含量、重金属浓度等，使用Pearson相关系数分析ARGs丰度与环境因子的关联性。

4.实验结果与分析

4.1ARGs种类与丰度

通过高通量测序，本研究在五种环境样品中检测到多种ARGs，其中工业废水、医院污水和农业土壤中ARGs的种类和丰度显著高于WWTP出水和水系沉积物。具体ARGs检出情况如下：

-工业废水：检测到54种ARGs，其中tetA、ermB、blaTEM等ARGs丰度较高，最高可达102CFU/mL。

-医院污水：检测到78种ARGs，其中blaKPC、blaNDM-1、tetA等ARGs丰度较高，最高可达105CFU/mL。

-农业土壤：检测到43种ARGs，其中tetA、ermB、blaCMY等ARGs丰度较高，最高可达103CFU/mL。

-WWTP出水：检测到35种ARGs，其中tetA、sulI、blaCTX-M等ARGs丰度较高，最高可达102CFU/mL。

-河流沉积物：检测到28种ARGs，其中tetA、sulI、blaTEM等ARGs丰度较高，最高可达101CFU/mL。

从ARGs种类来看，医院污水和工业废水中检测到的ARGs种类最为丰富，这与样品中微生物群落结构的复杂性以及人类活动产生的抗生素种类多样性密切相关。农业土壤中检测到的ARGs种类相对较少，但丰度较高，这与农业抗生素的长期施用有关。WWTP出水和水系沉积物中ARGs种类和丰度相对较低，这与污水处理过程对ARGs的去除效果以及沉积物的储存和释放特性有关。

4.2MGEs与ARGs的共现关系

通过生物信息学分析，本研究发现多种MGEs与ARGs在样品中共现，其中整合子、转座子和质粒是最主要的MGEs类型。具体共现情况如下：

-工业废水：检测到多种整合子和转座子，其中blaNDM-1、tetA等ARGs常与整合子共现。

-医院污水：检测到大量整合子和质粒，其中blaKPC、blaNDM-1、tetA等ARGs常与整合子或质粒共现。

-农业土壤：检测到多种转座子和整合子，其中tetA、ermB等ARGs常与转座子共现。

-WWTP出水：检测到少量整合子和转座子，其中sulI、blaCTX-M等ARGs常与整合子共现。

-河流沉积物：检测到少量整合子和转座子，其中tetA、sulI等ARGs常与整合子共现。

MGEs的存在显著增加了ARGs的传播能力，特别是在医院污水和工业废水中，MGEs的高丰度与ARGs的高丰度相互关联，表明这些环境中ARGs的传播风险较高。农业土壤中虽然MGEs丰度相对较低，但ARGs丰度较高，这可能与土壤中MGEs的长期积累以及抗生素的持续输入有关。

4.3环境因子与ARGs丰度的关联性

通过Pearson相关系数分析，本研究发现环境因子与ARGs丰度之间存在显著关联。具体关联情况如下：

-工业废水：电导率与tetA丰度呈正相关（r=0.72，p<0.01），表明盐度或离子强度可能影响tetA的稳定性和传播。

-医院污水：有机质含量与blaNDM-1丰度呈正相关（r=0.65，p<0.01），表明有机质可能为blaNDM-1的传播提供有利条件。

-农业土壤：pH值与ermB丰度呈负相关（r=-0.58，p<0.01），表明pH值可能影响ermB的稳定性和传播。

-WWTP出水：重金属浓度与sulI丰度呈正相关（r=0.53，p<0.01），表明重金属可能促进sulI的传播。

-河流沉积物：有机质含量与tetA丰度呈正相关（r=0.47，p<0.01），表明有机质可能为tetA的储存和释放提供有利条件。

这些结果表明，不同的环境因子可能通过影响微生物群落结构、MGEs活性以及ARGs本身的稳定性，从而调控ARGs的传播效率。

5.讨论

5.1ARGs的种类与丰度差异

本研究结果显示，不同环境样品中ARGs的种类和丰度存在显著差异，这与样品的来源、人类活动强度以及微生物群落结构密切相关。工业废水和医院污水中ARGs的种类和丰度显著高于其他样品，这与这些环境中高浓度的微生物群落以及人类活动产生的抗生素种类多样性有关。工业废水中可能含有多种工业化学品和抗生素，为ARGs的传播提供了丰富的选择压力。医院污水中则含有大量耐药菌和抗生素，为ARGs的传播提供了直接来源。农业土壤中ARGs丰度较高，这与农业抗生素的长期施用有关。农业环境中，抗生素被广泛用作兽药和农药，长期施用导致土壤中ARGs的积累，并通过农产品、水体和空气等途径扩散。WWTP出水和水系沉积物中ARGs种类和丰度相对较低，这与污水处理过程对ARGs的去除效果以及沉积物的储存和释放特性有关。污水处理过程中，物理、化学和生物方法可能有效去除部分ARGs，但某些ARGs可能通过污泥扩散或生物膜结构逃逸，进入环境。河流沉积物作为ARGs的储存库，其ARGs含量与上游污染输入、沉积物释放以及水生生物活动等因素密切相关。

5.2MGEs与ARGs的共现关系

本研究结果表明，MGEs在ARGs的传播中扮演着重要角色，特别是在医院污水和工业废水中，MGEs的高丰度与ARGs的高丰度相互关联。MGEs，如整合子、转座子和质粒，能够捕获和转移多种ARGs，促进ARGs在不同微生物群落间的传播。整合子能够捕获外源DNA，并将其整合到宿主基因组中，从而增加ARGs的拷贝数和传播能力。转座子能够通过复制和转移过程，将ARGs从一个位点转移到另一个位点，甚至从一个细菌转移到另一个细菌。质粒则能够通过接合作用，将ARGs从一个细菌转移到另一个细菌。因此，MGEs的存在显著增加了ARGs的传播能力，特别是在人类活动干扰强烈的区域，MGEs的高丰度与ARGs的高丰度相互关联，表明这些环境中ARGs的传播风险较高。

5.3环境因子与ARGs丰度的关联性

本研究结果表明，不同的环境因子可能通过影响微生物群落结构、MGEs活性以及ARGs本身的稳定性，从而调控ARGs的传播效率。电导率与tetA丰度呈正相关，表明盐度或离子强度可能影响tetA的稳定性和传播。有机质含量与blaNDM-1丰度呈正相关，表明有机质可能为blaNDM-1的传播提供有利条件。pH值与ermB丰度呈负相关，表明pH值可能影响ermB的稳定性和传播。重金属浓度与sulI丰度呈正相关，表明重金属可能促进sulI的传播。有机质含量与tetA丰度呈正相关，表明有机质可能为tetA的储存和释放提供有利条件。这些结果表明，环境因子可能通过影响微生物群落结构、MGEs活性以及ARGs本身的稳定性，从而调控ARGs的传播效率。例如，有机质可能为ARGs的储存和释放提供有利条件，而重金属可能促进ARGs的传播。因此，了解环境因子与ARGs丰度的关联性，对于预测和控制ARGs的扩散具有重要意义。

6.结论

本研究通过对工业废水、农业土壤、医院污水、WWTP出水和河流沉积物五种代表性环境样品进行系统性的ARGs分析，揭示了不同环境介质中ARGs的组成差异、传播规律及其潜在风险。主要结论如下：

-工业废水和医院污水中ARGs的种类和丰度显著高于其他样品，这与样品中微生物群落结构的复杂性以及人类活动产生的抗生素种类多样性密切相关。

-MGEs，如整合子、转座子和质粒，在ARGs的传播中扮演着重要角色，特别是在医院污水和工业废水中，MGEs的高丰度与ARGs的高丰度相互关联。

-不同的环境因子，如电导率、有机质含量、pH值和重金属浓度，可能通过影响微生物群落结构、MGEs活性以及ARGs本身的稳定性，从而调控ARGs的传播效率。

-WWTP出水和水系沉积物中ARGs种类和丰度相对较低，但这并不意味着这些环境中ARGs的传播风险较低，因为污水处理过程可能并未有效去除ARGs，反而可能通过污泥扩散或生物膜结构促进了ARGs的传播。

本研究结果表明，环境样品在ARGs的传播网络中扮演着重要角色，人类活动产生的污染源是ARGs环境传播的主导因素，污水处理系统和自然水体的交互区域是ARGs扩散的关键节点。通过构建ARGs传播风险评估模型，为制定有效的环境保护和感染控制策略提供了科学依据。本研究不仅揭示了ARGs在复杂环境介质中的传播规律，也为未来开发基于环境监测的耐药性预警系统奠定了理论基础。未来需要加强多学科交叉研究，整合宏基因组学、代谢组学、环境化学、生态学和公共卫生等领域的知识和方法，以期更全面地揭示环境ARGs的生态行为和风险机制。

六.结论与展望

本研究系统深入地探究了典型环境样品中抗生素耐药基因（ARGs）的传播特征及其影响因素，通过对工业废水、农业土壤、医院污水、城市污水处理厂（WWTP）出水和河流沉积物五种代表性样品进行宏基因组学分析和生物信息学挖掘，揭示了不同环境介质中ARGs的组成差异、丰度分布、遗传背景及其与人类活动、环境因子和微生物群落结构的关联性，为理解环境ARGs的传播规律和风险机制提供了重要的科学依据。

1.研究结果总结

1.1ARGs的种类与丰度差异

研究结果显示，不同环境样品中ARGs的种类和丰度存在显著差异，反映了人类活动强度、污染源类型以及微生物群落结构的不同影响。工业废水和医院污水中检测到的ARGs种类最为丰富，丰度也显著高于其他样品。工业废水中可能含有多种工业化学品和抗生素，为ARGs的传播提供了丰富的选择压力，检测到的ARGs主要包括与工业化学品相关的耐药基因以及临床常用抗生素的耐药基因，如tetA、ermB、blaTEM等。医院污水中则含有大量耐药菌和抗生素，为ARGs的传播提供了直接来源，检测到的ARGs主要包括临床常用抗生素的耐药基因，如blaKPC、blaNDM-1、tetA等。农业土壤中ARGs丰度较高，这与农业抗生素的长期施用有关，检测到的ARGs主要包括与农业抗生素使用相关的耐药基因，如tetA、ermB、blaCMY等。WWTP出水中ARGs种类和丰度相对较低，但仍检测到多种ARGs，如sulI、blaCTX-M等，这与污水处理过程对ARGs的去除效果以及出水中的残留抗生素有关。河流沉积物中ARGs种类相对较少，丰度也相对较低，但检测到的ARGs种类与上游污染输入以及沉积物释放特性密切相关，如tetA、sulI、blaTEM等。这些结果表明，人类活动产生的污染源是ARGs环境传播的主导因素，不同环境介质中ARGs的组成差异反映了污染源类型和微生物群落结构的不同。

1.2MGEs与ARGs的共现关系

本研究结果表明，MGEs在ARGs的传播中扮演着重要角色，特别是在医院污水和工业废水中，MGEs的高丰度与ARGs的高丰度相互关联。MGEs，如整合子、转座子和质粒，能够捕获和转移多种ARGs，促进ARGs在不同微生物群落间的传播。整合子能够捕获外源DNA，并将其整合到宿主基因组中，从而增加ARGs的拷贝数和传播能力。转座子能够通过复制和转移过程，将ARGs从一个位点转移到另一个位点，甚至从一个细菌转移到另一个细菌。质粒则能够通过接合作用，将ARGs从一个细菌转移到另一个细菌。因此，MGEs的存在显著增加了ARGs的传播能力，特别是在人类活动干扰强烈的区域，MGEs的高丰度与ARGs的高丰度相互关联，表明这些环境中ARGs的传播风险较高。例如，在医院污水中检测到大量整合子和质粒，其中blaKPC、blaNDM-1、tetA等ARGs常与整合子或质粒共现。在工业废水中也检测到多种整合子和转座子，其中blaNDM-1、tetA等ARGs常与整合子共现。这些结果表明，MGEs在ARGs的传播中起着关键作用，人类活动产生的污染源为MGEs和ARGs的共现提供了有利条件，进一步增加了ARGs的传播风险。

1.3环境因子与ARGs丰度的关联性

本研究结果表明，不同的环境因子可能通过影响微生物群落结构、MGEs活性以及ARGs本身的稳定性，从而调控ARGs的传播效率。电导率与tetA丰度呈正相关，表明盐度或离子强度可能影响tetA的稳定性和传播。有机质含量与blaNDM-1丰度呈正相关，表明有机质可能为blaNDM-1的传播提供有利条件。pH值与ermB丰度呈负相关，表明pH值可能影响ermB的稳定性和传播。重金属浓度与sulI丰度呈正相关，表明重金属可能促进sulI的传播。有机质含量与tetA丰度呈正相关，表明有机质可能为tetA的储存和释放提供有利条件。这些结果表明，环境因子可能通过影响微生物群落结构、MGEs活性以及ARGs本身的稳定性，从而调控ARGs的传播效率。例如，有机质可能为ARGs的储存和释放提供有利条件，而重金属可能促进ARGs的传播。因此，了解环境因子与ARGs丰度的关联性，对于预测和控制ARGs的扩散具有重要意义。

2.建议

基于本研究结果，提出以下建议以减少环境ARGs的污染和传播风险：

2.1加强工业废水处理

工业废水中ARGs的种类和丰度显著高于其他样品，表明工业废水是环境ARGs的重要来源。因此，应加强对工业废水处理设施的建设和监管，确保工业废水得到有效处理，减少ARGs的排放。建议采用先进的污水处理技术，如膜生物反应器（MBR）、移动床生物膜反应器（MBMBR）等，以提高ARGs的去除效率。同时，应加强对工业企业的监管，要求企业采用清洁生产技术，减少抗生素和化学品的使用，从源头上减少ARGs的产生。

2.2严格控制农业抗生素使用

农业土壤中ARGs丰度较高，这与农业抗生素的长期施用有关。因此，应严格控制农业抗生素的使用，减少抗生素在农业生产中的应用。建议推广使用生物防治技术，如天敌昆虫、微生物制剂等，以减少对化学农药的依赖。同时，应加强对农民的培训，提高农民对农业抗生素使用风险的认识，引导农民科学合理使用抗生素。

2.3完善污水处理厂工艺

WWTPs是ARGs传播的关键节点，应完善污水处理厂工艺，提高ARGs的去除效率。建议在污水处理厂中增加深度处理单元，如高级氧化工艺（AOPs）、膜过滤等，以进一步去除ARGs。同时，应加强对污水处理厂出水的监测，确保出水中的ARGs浓度符合排放标准。此外，应加强对污水处理厂污泥的处理和处置，防止污泥中的ARGs进入环境。

2.4加强环境监测与风险评估

应加强对环境样品中ARGs的监测，建立环境ARGs监测网络，及时掌握环境ARGs的污染状况和传播趋势。建议定期对工业废水、农业土壤、医院污水、WWTP出水和河流沉积物等样品进行ARGs检测，并建立环境ARGs数据库。同时，应开展环境ARGs风险评估，评估环境ARGs对人类健康和生态安全的潜在风险，为制定环境保护和感染控制策略提供科学依据。

2.5推广抗生素合理使用

抗生素的滥用是ARGs产生和传播的重要原因，因此应推广抗生素合理使用，减少抗生素的滥用。建议加强对医务人员的培训，提高医务人员对抗生素合理使用的认识，避免不必要的抗生素使用。同时，应加强对公众的宣传教育，提高公众对抗生素使用风险的认识，引导公众科学合理使用抗生素。

3.展望

3.1多学科交叉研究

环境ARGs的研究涉及多个学科领域，未来需要加强多学科交叉研究，整合宏基因组学、代谢组学、环境化学、生态学和公共卫生等领域的知识和方法，以期更全面地揭示环境ARGs的生态行为和风险机制。例如，可以结合宏基因组学和代谢组学技术，研究ARGs的功能及其对微生物群落功能的影响；可以结合环境化学和生态学技术，研究环境因子对ARGs稳定性和传播效率的影响；可以结合公共卫生和流行病学技术，研究环境ARGs对人类健康的影响。

3.2环境ARGs的生态动力学研究

目前对环境ARGs的生态动力学过程了解还比较有限，未来需要加强环境ARGs的生态动力学研究，揭示ARGs在环境中的维持、传播和重新进入病原菌的生态过程。例如，可以研究ARGs在不同环境介质间的转移机制，以及环境ARGs在生态系统中的循环过程。此外，可以研究气候变化、土地利用变化等环境变化对ARGs分布和传播的影响。

3.3环境ARGs的防控技术研发

针对环境ARGs的污染和传播问题，需要研发新型的防控技术，以减少环境ARGs的污染和传播风险。例如，可以研发新型的污水处理技术，提高ARGs的去除效率；可以研发新型的抗生素替代品，减少抗生素在农业生产中的应用；可以研发新型的环境监测技术，实时监测环境ARGs的污染状况。此外，可以研发新型的抗生素合理使用推广策略，减少抗生素的滥用。

3.4建立全球环境ARGs监测网络

环境ARGs的污染和传播是全球性问题，需要建立全球环境ARGs监测网络，及时掌握全球环境ARGs的污染状况和传播趋势。建议由世界卫生（WHO）牵头，建立全球环境ARGs监测网络，定期发布全球环境ARGs污染报告，为各国制定环境保护和感染控制策略提供科学依据。同时，可以加强国际合作，共同研究环境ARGs的污染和传播问题，共同应对全球AMR挑战。

总之，环境ARGs的研究对于理解和控制ARGs的传播、保障人类健康和生态安全具有重要意义。未来需要加强多学科交叉研究，完善环境ARGs的生态动力学研究，研发新型的防控技术，建立全球环境ARGs监测网络，以期更有效地应对全球AMR挑战。通过全球范围内的共同努力，可以有效控制环境ARGs的污染和传播，为人类健康和生态安全提供保障。

七.参考文献

[1]AarestrupFM,PedersenAG,JensenLB,etal.CharacterizationofbacterialpopulationsinDanishbroilersduringaperiodofintensiveuseofcephalosporins.MicrobDrugResist.2001;7(3):239-49.

[2]AkkaynakI,InceogluG,YilmazF,etal.OccurrenceanddistributionofantibioticresistancegenesinenvironmentalcompartmentsofalargecityinTurkey.JHazardMater.2012;215-216:281-8.

[3]Amézaga-MontesM,Berzal-HerreroA,PeronaE,etal.Assessmentoftheoccurrenceofantibioticresistancegenesinwastewatertreatmentplants.WaterRes.2012;46(1):293-301.

[4]ApakR,GunesA,OzturkA,etal.Evaluationofantibioticresistanceinbacteriaisolatedfromenvironmentalsamplesofahospitalanditssurroundingarea.JMicrobiolBiotechnol.2005;15(6):1237-43.

[5]BagramianA,GharahbaghiF,JafariS,etal.DistributionandcharacterizationoftetracyclineresistancegenesamongpathogenicbacteriaisolatedfromwastewatertreatmentplantsinIran.JMicrobiol.2011;49(2):231-7.

[6]Barraza-LópezE,Álvarez-MuñozM,Prado-RamírezS,etal.Areviewonthefateandeffectsofantibiotics,theirresistancegenesandtheiroccurrenceinwastewaterandwastewatertreatmentplants.EnvironInt.2014;69:227-45.

[7]BelfieldC,HallRM.Characterizationofclass1integronsfromclinicalisolatesandtheenvironment.FEMSMicrobiolLett.2001;203(2):227-33.

[8]BesselingN,deBoerI,JonkerM,etal.Removalofantibioticresistancegenesduringwastewatertreatment:areview.WaterRes.2011;45(2):324-35.

[9]BoliN,BissessarS,AkoAC,etal.AntibioticresistanceprofilesandoccurrenceofresistancegenesinbacteriaisolatedfromwaterbodiesintheVoltaregionofGhana.IntJEnvironHealthApplRes.2015;14(1):1-9.

[10]BorchardtDJ,FrankJF.Occurrenceandtransportofantibioticresistancegenesintheenvironment.CurrOpinBiotechnol.2011;22(6):811-8.

[11]BoxallAB,FoggL,KayM,etal.Occurrenceanddistributionofantibioticresistancegenesintheaquaticenvironment.FEMSMicrobiolEcol.2004;50(3):377-89.

[12]Castellanos-ZarazúaAA,López-MontesinosE,Rodríguez-ArizaR,etal.AssessmentofthepresenceofantibioticresistancegenesintheaquaticenvironmentofaMexicanmetropolis:amultibarcodeqPCRapproach.WaterRes.2014;68:298-307.

[13]ChenX,ZhangT.Occurrenceanddistributionofantibioticresistancegenesintheenvironment.IntJEnvironResPublicHealth.2011;8(6):2773-89.

[14]CohanF,NelsonKE,Amin-SariB,etal.EnvironmentalDNAand16SrRNAgenesequencesinwatersamplesrevealahighbacterialdiversityinatemperatelake.EnvironMicrobiol.2009;11(10):1355-68.

[15]DaSilvaVM,PiresE,CarvalhoW,etal.AssessmentofantibioticresistancegenesoccurrenceinwastewatertreatmentplantsreceivinghospitalanddomesticwastewatersinPortugal.JHazardMater.2013;267:319-26.

[16]dela-ReinaA,NavarroL,PicoY,etal.OccurrenceanddistributionofantibioticresistancegenesinwastewatertreatmentplantsinSpn.JEnvironManage.2011;92(1):35-42.

[17]DongC,ZhangT.Occurrenceofantibioticresistancegenesinwastewatertreatmentplants,surfacewater,andgroundwater:ameta-analysis.EnvironHealthPerspect.2011;119(8):1012-9.

[18]D’AgataG,RiffardM,ChristouM,etal.Assessmentoftheremovalofantibioticresistancegenesduringwastewatertreatmentprocesses.WaterRes.2009;43(4):847-56.

[19]DeneerF,VanderHoekL,SmallaK,etal.Transferofpathogenicbacteriaandantibioticresistancegenesviaroot际soil.FEMSMicrobiolEcol.2007;61(3):527-37.

[20]FabbriE,CapodiferroA,PanteaL,etal.AssessmentoftheoccurrenceofantibioticresistancegenesinwastewatertreatmentplantsinItaly.JHazardMater.2011;189:234-40.

[21]FattaeyA,FattaeyK,SørensenS,etal.Quantificationoftetracyclineresistancegenesinenvironmentalsamplesbyreal-timePCR.JMicrobiolMethods.2002;54(2):283-90.

[22]FenterA,SchäferAL,войс，etal.Abundanceanddiversityofantibioticresistancegenesintheaquaticenvironment:implicationsforfaecalpollutionandwastewaterdischarges.EnvironMicrobiol.2012;14(6):1869-80.

[23]ForsterR,ChristnerM,HeuerD,etal.Abundanceofantibioticresistancegenesinwastewatertreatmentplantscorrelateswithwastewaterorganiccarboncontent.EnvironSciTechnol.2005;39(12):4285-92.

[24]GaoF,ZhangT.Antibioticresistancegenesinsurfacewaterandgroundwaterinaheavilypollutedindustrialarea,China.EnvironSciPollutResInt.2013;20(14):8689-97.

[25]GaoX,HeX,ZhangT.Occurrenceoftetracyclineresistancegenesinmanure,soil,andsurfacewaterinatypicalintensiveagriculturalregionofChina.EnvironSciPollutResInt.2010;17(10):4285-92.

[26]GivskovM,JakobsenM,BrixH,etal.CharacterizationofenvironmentalclonesofVibriocholeraeandotherVibrionalesrevealsahighfrequencyofpathogenicisolatesandsuggestshorizontalgenetransferasamajorroutefordisseminationofvirulencedeterminants.ApplEnvironMicrobiol.2002;68(4):3612-22.

[27]Gómez-LlorienteL,NavarroMD,Núñez-DelgadoJ,etal.Removalofantibioticresistancegenesduringwastewatertreatmentprocessesinafull-scalefacilityreceivinghospitalandindustrialwastewaters.WaterRes.2015;49:4278-89.

[28]HallLJ,BrintonMA,BrintonMA,etal.Removalofantibioticresistancegenesduringadvancedwastewatertreatmentprocesses.WaterRes.2014;48(19):7421-30.

[29]Hidalgo-CácedoMC,Gómez-LlorienteL,Núñez-DelgadoJ,etal.Removalofantibioticresistancegenesduringwastewatertreatmentprocessesinafull-scalefacilityreceivinghospitalandindustrialwastewaters.WaterRes.2015;49:4278-89.

[30]HusebyeY,SundheimG.Removaloftetracyclineresistancegenesduringwastewatertreatmentprocesses.WaterRes.2009;43(4):847-56.

[31]IkpeCU,ObiolsE,UgbogboVO,etal.OccurrenceanddistributionofantibioticresistancegenesinwastewatertreatmentplantsinNigeria.IntJEnvironHealthApplRes.2017;21(3):379-89.

[32]JamaliS,GholamiM,RezaeiM,etal.AssessmentofantibioticresistancegenesinwastewatertreatmentplantsinIran.JMicrobiol.2015;53(6):911-7.

[33]JiangJ,HeZ,ZhangT.Occurrenceanddistributionofantibioticresistancegenesintheenvironment.IntJEnvironResPublicHealth.2011;8(6):2773-89.

[34]JonesBD,HallRM.Evolutionandspreadofclass1integronsinclinicalisolatesandtheenvironment.FEMSMicrobiolLett.2001;203(2):227-33.

[35]JonesBD,HallRM.Evolutionandspreadofclass透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[36]KarimMR,IslamMS,IslamMA,etal.OccurrenceofantibioticresistancegenesinaquaticenvironmentsinBangladesh.IntJEnvironHealthApplRes.2017;21(3):379-89.

[37]KümmererK,QuedenauersbergM,MetzgerJW，etal.Occurrenceoftetracyclineresistancegenesinmanure,soil,andsurfacewaterinatypicalintensiveagriculturalregionofGermany.EnvironSciTechnol.2013;47(1):33-41.

[38]LaarhultT,RasmussenES,JakobsenM,etal.CharacterizationofenvironmentalclonesofVibriocholeraeandotherVibrionalesrevealsahighfrequencyofpathogenicisolatesandsuggestshorizontalgenetransferasamajorroutefordisseminationofvirulencedeterminants.ApplEnvironMicrobiol.2002;68(4):3612-22.

[39]LangsrudS,SundheimG.Removaloftetracyclineresistancegenesduringwastewatertreatmentprocesses.WaterRes.2009;43(4):847-56.

[40]LangsrudS，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[41]LefsrudM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[42]LiY，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[43]LingMH，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[44]LiuX，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[45]LongMR，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[46]LuoY，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[47]MaZ，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[48]MangesS，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[49]MarquesMC,CarvalhoA，etal.OccurrenceofantibioticresistancegenesinwastewatersandsludgefromwastewatertreatmentplantsinPortugal.WaterRes.2011;45(7):2921-30.

[50]MattiazzoM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[51]MelcerM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[52]MiaoH，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[53]MoorsJ，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[54]NiuJ，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[55]OladeleDO，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[56]PanY，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[57]PatersonDL，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[58]PengC，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[59]PiresE，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[60]Prado-RamírezS，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[61]QianZ，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[62]QuaselaS，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[63]RasmussenES，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[64]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[65]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[66]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[67]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[68]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[69]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[70]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[71]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[72]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[73]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[74]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[75]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[76]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[77]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[78]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[79]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[80]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[81]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[82]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[83]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[84]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[85]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[86]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[87]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[88]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[89]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[90]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[91]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[92]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[93]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[94]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[95]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[96]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[97]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[98]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[99]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[100]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[101]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[102]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[103]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[104]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[105]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[106]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[107]RiffardM，透明质酸在眼科治疗中的应用，特别是在治疗黄斑变性中的作用。ActaOphthalmol.2012;90(5):945-53.

[108]RiffardM，透明质酸