《Microbiology Spectrum》：Optimizing a culture-enriched hybrid metagenomics pipeline to assess the AMR footprint of livestock manure in anaerobic digestate
编辑推荐：
来自畜牧生产系统的环境样品，包括粪便和厌氧消化物，作为抗菌药物耐药基因（antimicrobial resistance genes, ARGs）的储库可能被低估，因为传统宏基因组方法可能忽略低丰度ARGs，且通常缺乏将这些基因与其微生物宿主及共定位的移动遗传
来自畜牧生产系统的环境样品，包括粪便和厌氧消化物，作为抗菌药物耐药基因（antimicrobial resistance genes, ARGs）的储库可能被低估，因为传统宏基因组方法可能忽略低丰度ARGs，且通常缺乏将这些基因与其微生物宿主及共定位的移动遗传元件（mobile genetic elements, MGEs）相关联的分辨率。研究人员评估了在有或无抗生素选择的情况下，培养富集宏基因组学（culture-enriched metagenomics, CEMG）是否能增强厌氧消化物中ARGs的检测，并通过混合短读长和长读长宏基因组组装提高ARG-MGE-宿主关联的分辨率。CEMG提高了ARG回收率；平均ARG丰度从宏基因组新鲜消化物（fresh digestate, FD）中的15.4 counts per million (CPM)升至无抗生素CEMG中的124 CPM和抗生素选择性CEMG中的160 CPM。在FD中仅检测到9个独特的ARGs，而CEMG回收了112个，包括具有临床重要性的ARGs，如糖肽耐药、β-内酰胺酶基因，以及赋予多种抗生素类别交叉耐药的cfr 23S rRNA甲基转移酶。抗生素选择诱导了ARG谱的靶向、类别特异性变化，其中与四环素耐药相关的ARGs在各处理中持续富集。混合宏基因组组装解析了784个ARGs的基因组背景，其中59.3%与至少一类MGEs共定位，主要是质粒和整合性接合元件/整合性可移动元件。杀菌剂和金属耐药基因经常与ARGs在同一重叠群上共现。总之，这些发现表明抗生素选择性培养富集通过提高低丰度ARGs的检测增强了耐药组监测，而混合组装为评估其移动性和宿主关联提供了关键的基因组背景。
论文解读：培养富集混合宏基因组学优化厌氧消化物中抗生素耐药性足迹的评估
抗菌药物在畜禽养殖中的广泛使用，使得牲畜粪便及其厌氧消化产物成为抗菌药物耐药基因（antimicrobial resistance genes, ARGs）的重要环境储库。然而，直接进行宏基因组测序在复杂环境样品中常因测序深度不足而遗漏低丰度但具有临床意义的ARGs，并且难以将ARGs与其微生物宿主及共定位的移动遗传元件（mobile genetic elements, MGEs）关联起来，从而可能低估环境耐药组的多样性与传播风险。为突破这一局限，研究人员在厌氧消化物中系统评估了培养富集宏基因组学（culture-enriched metagenomics, CEMG）结合混合短读长与长读长组装策略的效能，旨在提高ARG检出率并解析其基因组背景。该研究发表于《Microbiology Spectrum》。
研究人员从实验室长期运行的厌氧消化器中采集新鲜消化物（fresh digestate, FD），分别进行有氧和厌氧培养富集，设置无抗生素对照组以及氨苄青霉素、头孢噻呋、红霉素、林可霉素、四环素、万古霉素六种抗生素处理，每种抗生素设三个浓度梯度（1×、10×、100×基于大小调整的最低选择浓度）。结合Illumina短读长测序和Oxford Nanopore长读长测序，采用读长比对和混合组装两种互补分析路径，全面比较直接宏基因组与CEMG对耐药组的表征能力，并解析ARG-MGE-宿主关联。
研究主要结果分为以下五个方面。
（1）培养富集对消化物微生物多样性的影响。基于属水平Bray-Curtis距离的主坐标分析和PERMANOVA显示，富集群落与FD显著分离，氧气条件是最强驱动因素（R2=0.240，P=0.001）。有氧富集群落以Alcaligenes为主（约50.2%），厌氧富集群落以Proteus（约25%）、Streptococcus（约13%）等为主。FD的Shannon多样性指数（3.02）显著高于CEMG无抗生素组（1.63）和CEMG抗生素组（1.68），而抗生素种类和浓度对整体多样性影响有限，仅万古霉素表现出显著的浓度依赖性属数下降。
（2）CEMG导致的耐药组组成和多样性变化。基于CARD数据库的读长比对显示，CEMG抗生素组的平均ARG丰度从FD的15.4 CPM升至160.0 CPM，提高约10.4倍；CEMG无抗生素组也达124 CPM。全部样品共鉴定出116个高置信度ARGs，其中FD仅检出9个，而CEMG检出112个（约96.6%）。富集样品独有检出糖肽耐药相关基因家族、CTX-M/OXA/TEM型β-内酰胺酶基因、多类耐药甲基转移酶基因cfr、lsa型ABC-F蛋白及多种氨基糖苷修饰酶基因。四环素耐药ARGs在所有处理中持续占优势，其丰度相对FD增加约14倍。
（3）CEMG增强不同预测耐药机制相关ARGs的检出。FD以抗生素灭活（46.15%）和靶标保护（39.56%）为主，主要由林可酰胺核苷酸转移酶和msr型ABC-F蛋白驱动。厌氧富集后，外排泵和靶标保护成为主要机制，MFS外排泵基因家族占比最高；有氧富集则仍以抗生素灭活为主，ant(4′)氨基糖苷核苷酸转移酶、erm和cfr 23S rRNA甲基转移酶等高丰度表达。
（4）ARG-MGE共定位及ARG与重金属耐药基因共选择。混合组装共获得784个ARG实例，分布于651条独特contigs上，其中59.3%与至少一类MGE共定位。质粒是最常见的MGE类型（占65.2%的MGE相关contigs），其次为插入序列（32.9%）、整合性接合元件/整合性可移动元件（ICEs/IMEs，21.8%）。同时，ARG与杀菌剂和金属耐药基因在同一contig上频繁共现，厌氧富集条件下锌、砷耐药基因占比较高，有氧富集则显著富集西曲溴铵和三氯生耐药基因。
（5）ARG-MGE-宿主关联。对ARG-bearing contigs进行Kraken2分类发现，大多数归属于Bacillota（53.4%）、Pseudomonadota（29.5%）等。四环素耐药基因的宿主范围最广，涵盖Streptococcus、Proteus、Enterococcus、Citrobacter等属，且常与ICE/IME、插入序列及质粒共定位；氨基糖苷类耐药基因主要分布于Acinetobacter、Enterococcus等属且以质粒关联为主；β-内酰胺类耐药基因几乎仅由质粒携带，见于Escherichia、Citrobacter等属。
讨论部分指出，培养富集不仅提高了ARG的可检测性，还通过选择特定微生物类群重塑了观察到的耐药组结构，使得在直接宏基因组中不可见的可培养低丰度耐药菌群得以显现。四环素耐药基因在多种非对应抗生素富集处理中持续增加，提示共选择机制可能维持了耐药决定簇的存续。质粒和ICE/IME上ARG的高频共定位意味着厌氧消化后仍可能存在水平转移潜力，而ARG与金属/杀菌剂耐药基因的共现进一步支持金属选择压力在维持环境耐药组中的作用。研究结论强调，仅基于整体群落宏基因组的评估可能低估厌氧消化物中ARGs的多样性和移动传播潜力；将富集策略与直接宏基因组分析相结合，能够为评估畜禽源厌氧消化物相关抗菌药物耐药性的环境传播风险提供更全面的框架。
