《Environmental Microbiology Reports》：Metagenomic Characterisation of Antibiotic Resistance in Anaerobic Digesters and Their Implications for Environmental Health
编辑推荐：
中文翻译：环境系统中的抗微生物耐药性（antimicrobial resistance, AMR）是一项新兴的“同一健康”（One Health）挑战，厌氧消化器可能作为耐药相关遗传特征的储库。本研究分析了从NCBI Sequence Read Archive
中文翻译：环境系统中的抗微生物耐药性（antimicrobial resistance, AMR）是一项新兴的“同一健康”（One Health）挑战，厌氧消化器可能作为耐药相关遗传特征的储库。本研究分析了从NCBI Sequence Read Archive检索获得的四个公开可用的厌氧消化器鸟枪法宏基因组数据集（M2、M3、M12和M17），以表征微生物功能、耐药相关注释和抗生素靶点基因座。这些数据集采用Illumina HiSeq 2000测序生成，并通过质量评估、组装、功能注释和耐药基因组（resistome）分析进行处理。共鉴定出1082个耐药相关注释，其中M2中丰度最高（348个），其次是M12（276个）、M17（259个）和M3（205个）。检测到的注释与多种抗菌药物类别相关，包括氨基糖苷类、喹诺酮类、β-内酰胺类、四环素类、糖肽类、大环内酯类和磺胺类。多个频繁检测到的基因座，包括rpoB、rpoC、gyrA、gyrB、EF-Tu、EF-G、Ddl和KasA，被解释为抗生素靶点或持家基因座，而非已确认的获得性耐药基因。探索性共丰度分析确定了功能相关基因座之间的强关联，包括rpoB-rpoC和gyrA-gyrB，这可能反映了共享的基因组存在或功能关系。这些发现强调了厌氧消化器的耐药组潜力，同时强调了需要宿主分辨的宏基因组学、移动遗传元件分析和突变水平表征来阐明环境AMR风险。
论文解读文章
一、研究背景与问题
抗生素耐药性已被公认为最紧迫的全球健康挑战之一，对人类健康、环境完整性和动物福利产生深远影响。耐药性在常用抗生素中的扩散导致治疗失败、住院时间延长和医疗成本上升，给全球卫生系统带来巨大压力。虽然抗微生物耐药性（antimicrobial resistance, AMR）的临床维度已得到充分认识，但其环境驱动因素，特别是在农业生态系统中，正受到越来越多的关注。农业系统已成为AMR的关键储库和传播途径，其来源包括污水处理厂（WWTPs）出水、牲畜粪便以及厌氧消化器。厌氧消化（anaerobic digestion, AD）被广泛用于处理农业残留物、食物垃圾、动物粪便和污水污泥等有机废弃物，兼具废弃物管理和可再生能源生产的双重效益。然而，AD系统也可能成为抗生素耐药基因（antibiotic resistance genes, ARGs）和移动遗传元件（mobile genetic elements, MGEs）的热点区域，特别是在消化产物被重新用作肥料时，可能促进耐药决定因子的环境传播。尽管对废弃物处理环境中AMR的关注日益增加，厌氧消化器对环境耐药组的贡献仍未得到充分解析。这一空白之所以重要，是因为厌氧消化系统的消化产物常被用作有机肥料，从而为耐药相关遗传物质进入农业土壤和更广泛的生态网络创造了可能途径。此外，还需要更清晰地区分获得性抗生素耐药基因与抗生素靶点或持家基因座，因为仅检测靶基因并不能确认耐药性，除非证明了赋予耐药性的突变或移动基因关联。
二、研究内容与意义
本研究利用来自NCBI的宏基因组测序数据，调查了地理上不同地点采集的厌氧消化产物样品中ARGs和MGEs的丰度与多样性，核心研究问题是：“厌氧消化器系统中抗生素耐药基因和移动遗传元件的多样性与丰度如何，其对环境健康和同一健康有何影响？”通过阐明厌氧消化器的微生物动态和耐药谱，本研究旨在表征公开可用的厌氧消化器鸟枪法宏基因组的功能谱、耐药相关注释和共丰度模式，同时区分获得性耐药决定因子与抗生素靶点或持家基因座，以对其环境和同一健康相关性提供审慎评估。研究结果表明，厌氧消化器系统含有多种功能基因、耐药相关注释和抗生素靶点基因座，支持了厌氧消化器可能对环境耐药组有所贡献的观点，特别是在消化产物在农业环境中再利用的情况下。然而，由于多个频繁检测到的基因座（包括rpoB、rpoC、gyrA、gyrB、EF-Tu和EF-G）是持家基因或抗生素靶点基因而非确定的获得性ARGs，研究结果应谨慎解读。这些基因座的存在本身并不能在没有突变水平、等位基因水平或移动遗传元件证据的情况下确认抗生素耐药性。
三、关键技术方法
研究人员从NCBI Sequence Read Archive（SRA）检索了四个公开可用的全基因组鸟枪法宏基因组数据集，对应BioSample登录号SAMN27412227、SAMN27412222、SAMN03742276和SAMN27412223，分别标记为M2、M3、M12和M17。其中M2、M3和M17来自美国伊萨卡，M12来自中国香港，均采用Illumina HiSeq 2000平台测序。原始序列数据经过质量评估后，使用BV-BRC基因组组装服务进行组装，并用Unicycler v0.4.8进行优化，Samtools v1.17处理比对，Pilon v1.2.3进行错误校正和间隙填充，QUAST v5.2.0评估组装质量。功能注释使用BV-BRC流程中的RASTtk工具包进行，并将注释的编码序列映射到KEGG数据库。耐药相关注释通过BV-BRC注释流程中的综合抗生素耐药数据库（CARD）鉴定，并进一步分为推定获得性ARGs和抗生素靶点或持家基因座两类。分类学分析使用Kraken v2.1.2和MiniKraken2数据库（2023年3月版本）进行。数据分析使用Python v3.11中的pandas、numpy、seaborn和matplotlib库，共丰度模式通过Spearman相关系数评估并在Gephi v0.9.2中可视化。
四、研究结果
3.1 基因组组装特征与注释统计
四个厌氧消化器宏基因组组装共产生138.36 Mb的基因组长度，包含12,000个重叠群和146,445个预测编码序列（CDSs）。组装大小在各样品间存在差异，从M3的22.64 Mb到M2的45.66 Mb不等。M2的CDS数量最高（48,058个），M3最低（23,851个）。组装质量指标也存在差异：重叠群N50值从M12的20,284 bp到M3的83,331 bp不等；GC含量从M17的45.30%到M12的53.88%。所有组装中共预测出72,120个功能蛋白和74,325个假定蛋白。四个宏基因组含有2,205个RNA基因和3,694个CRISPR相关基因。针对耐药数据库的注释在所有样品中鉴定出1,082个耐药相关基因，M2最高（348个），其次是M12（276个）、M17（259个）和M3（205个）。
3.2 厌氧消化器宏基因组中的功能基因分布
功能注释显示，代谢相关基因在所有宏基因组中构成最丰富的子系统类别，数量从M3的3,565个到M2的6,618个不等。蛋白质加工是第二大功能类别，各样品间有1,562至2,991个基因。能量代谢相关基因也高度丰富，特别是在M2（1,847个）和M12（1,822个）中。其他功能类别包括细胞过程、DNA加工、应激反应、防御和毒力、膜转运以及RNA加工。M2在大多数功能类别中通常表现出最高丰度，而M3的注释基因最少。细胞包膜功能、调控和细胞信号传导以及杂项功能在所有宏基因组中是最不丰富的子系统类别。
3.3 耐药相关注释分类为获得性ARGs和抗生素靶点基因座
检测到的耐药相关注释被分为两个解释性类别。第一类包括推定获得性ARGs，如与酶促失活、外排介导耐药、靶点保护和其他公认耐药机制相关的注释。第二类包括抗生素靶点或持家基因座，如rpoB、rpoC、gyrA、gyrB、EF-Tu、EF-G、Ddl、KasA和S12p。这些基因座在生物学上很重要，因为它们是已知的抗生素靶点，但仅检测到它们并不能证明耐药性，除非确认了耐药相关突变、等位基因水平变异或与移动遗传元件的关联。因此，这些基因座被解释为抗生素靶点或耐药相关基因座，而非确定的获得性ARGs。
3.4 耐药相关基因和抗生素靶点基因座的分布
耐药相关基因座的相对丰度在四个厌氧消化器宏基因组中存在差异。EF-Tu、EF-G、rpoB、rpoC、gyrA、gyrB、Ddl、KasA、S12p、Iso-tRNA和MirA在所有样品中均被检测到，但其比例代表存在差异。M2在多个基因座上通常表现出最高相对丰度，包括EF-Tu、EF-G和Iso-tRNA，而M12在rpoB和rpoC上表现出更高的代表性。M3和M17在各检测基因座上贡献了不同比例，反映了厌氧消化器微生物组中耐药组组成的差异。
3.5 耐药相关基因注释在抗生素类别间的分布
耐药相关基因注释在所有厌氧消化器宏基因组中与多种抗生素类别相关。氨基糖苷类、β-内酰胺类、四环素类、磺胺类、喹诺酮类、糖肽类、肽类抗生素和大环内酯类是最常代表的类别。抗生素类别相关注释的相对丰度在各样品间存在差异，M2和M12通常在多个抗生素类别上表现出更广泛的代表性，而M3在若干类别上的代表性相对较低。
3.6 耐药相关和抗生素靶点基因座间的探索性共丰度模式
探索性相关分析用于检验检测到的耐药相关注释和抗生素靶点基因座之间的共丰度模式。rpoB和rpoC之间以及gyrA和gyrB之间观察到强正相关。这些关系在生物学上是预期的，因为rpoB和rpoC编码细菌RNA聚合酶亚基，而gyrA和gyrB编码DNA旋转酶亚基。因此，这些相关性不应被解释为耐药传播、共选择或水平基因转移的直接证据，而更可能反映共享的基因组存在、功能关联或所分析宏基因组中细菌群落组成的相似性。
五、讨论与结论
讨论部分指出，厌氧消化系统不仅因其在可再生能源生产和废弃物增值中的作用而被认可，而且作为可能影响抗微生物耐药性环境动态的重要微生物生态系统也日益受到关注。研究结果揭示了各样品间组装质量、功能基因分布和耐药组组成的显著差异，反映了厌氧消化系统内微生物群落结构和环境条件的差异。测序策略的差异似乎影响组装结果，配对末端数据集比更新的Illumina HiSeq数据集产生更碎片化的组装。大量假定蛋白的存在表明厌氧消化器内很大一部分微生物功能潜力仍未表征。代谢功能在所有注释基因中占很大比例，支持厌氧消化的连续水解、发酵、产乙酸和产甲烷阶段。应激反应、防御和毒力相关基因的丰富存在可能反映了微生物对波动底物组成、抑制性代谢物积累、营养竞争和微生物拮抗等环境胁迫的适应。检测到的耐药相关注释与人类医学、兽医实践和农业中广泛使用的抗生素类别相关。然而，由于本研究未进行突变水平变异调用，这些基因座的检测被谨慎解释为耐药相关或抗生素靶点潜力，而非表型AMR的直接证据。共丰度分析中观察到的强正相关（如rpoB-rpoC和gyrA-gyrB）更可能反映共享基因组存在、共同分类分布或功能相互依赖，而非耐药传播的直接证据。研究结论指出，未来研究应整合获得性ARG重点分析、宿主分辨的宏基因组重建、质粒和整合子检测、重叠群水平ARG-MGE关联以及抗生素靶基因的变异水平分析。此类方法将有助于确定厌氧消化器中的耐药相关注释是否代表可转移的环境AMR风险，还是反映常驻微生物群落的背景基因组组成。建议在“同一健康”监测框架内对厌氧消化系统和消化产物进行常规监测，特别是在消化产物施用于农业土壤的情况下。
