《Microorganisms》：Mycobacteriophages: An Integrated Platform for Rapid Diagnosis and Targeted Therapy of Tuberculosis Infection
编辑推荐：
传统污水宏基因组学能够实现人群层面的抗生素耐药性（antimicrobial resistance, AMR）监测，但对耐药基因的细菌宿主、基因组背景及移动性却只能提供有限的认识。本研究采用高通量染色体构象捕获（high-throughput chromoso
传统污水宏基因组学能够实现人群层面的抗生素耐药性（antimicrobial resistance, AMR）监测，但对耐药基因的细菌宿主、基因组背景及移动性却只能提供有限的认识。本研究采用高通量染色体构象捕获（high-throughput chromosome conformation capture, Hi-C）宏基因组学，这是一种邻近连接技术，用于解析宿主相关的AMR。研究人员在四个月的采样期内，从两座城市污水处理厂采集了六份进水污水样品。综合运用鸟枪法宏基因组测序和Hi-C测序来重建宏基因组组装基因组（metagenome-assembled genomes, MAGs）、质粒分箱和病毒分箱，并通过Hi-C接触将AMR、应激反应和毒力相关基因与其细菌宿主关联起来。从12877个基因组分箱中，有454个分箱达到质量阈值（139个MAGs、130个质粒分箱、185个病毒分箱），涵盖70个细菌和古菌科。AMRFinderPlus共识别出771个耐药及适应性相关基因命中（517个AMR基因、232个应激反应基因、22个毒力相关基因），以大环内酯-林可酰胺-链阳菌素B（macrolide–lincosamide–streptogramin B, MLSB）类（n = 117）、β-内酰胺类（n = 114）和四环素类（n = 96）耐药决定簇为主。Hi-C将490个基因命中与宿主科关联起来；其中肠杆菌科（Enterobacteriaceae）的宿主归属检出数最高（n = 131），而气单胞菌科（Aeromonadaceae）携带的AMR基因谱系最广（n = 102；72%为AMR基因）。宿主-质粒网络包含137个AMR阳性质粒分箱，通过166个Hi-C支持的关联连接至106个MAGs。有五个质粒重叠群（包括一个携带blaCMY-4的质粒）与来自两个不同细菌纲的MAGs相关联，提示存在跨谱系的质粒携带。病毒性AMR归属有限，仅在一个病毒重叠群上鉴定出erm(B)这一耐药决定簇。Hi-C宏基因组学将污水耐药组的分析从单纯的基因清单拓展到细菌宿主、质粒关联和基因组背景的解析，从而增强了基于污水的AMR监测在生态学解释方面的能力。
论文解读
一、研究背景与科学问题
抗生素耐药性（antimicrobial resistance, AMR）已成为全球重大公共卫生威胁，耐药决定簇在人类、动物和环境之间相互传播，构成了典型的“One Health”挑战。环境，尤其是城市污水系统，被视为耐药基因储存、进化与交换的重要场所。传统污水宏基因组学通过鸟枪法测序能够无培养地检测复杂微生物群落中的耐药基因并评估其相对丰度，但其核心局限在于：只能提供基因清单，无法确定特定耐药基因的宿主菌、其定位于染色体还是可移动遗传元件、以及在分类学边界间的传播潜力。该“基因-宿主-移动性”关联信息缺失，严重限制了环境耐药组数据向临床AMR风险评估的转化。高通量染色体构象捕获（Hi-C）宏基因组学利用细胞内DNA邻近连接原理，可将质粒、噬菌体和耐药基因归属至其微生物宿主，从而弥补这一缺口。然而，Hi-C此前在污水耐药组宿主解析方面的应用仍十分有限，尤其缺乏阿拉伯海湾地区的数据。为此，研究人员开展本研究，旨在联合Hi-C与鸟枪法宏基因组测序，解析污水中耐药基因的细菌宿主与基因组背景，探索耐药组的生态组织结构。
二、研究设计与技术方法
研究人员于2024年10月至2025年1月期间，从阿联酋迪拜两座污水处理厂（WWTP1和WWTP2）每月采集24小时复合进水样品，共获得八份样品，其中两份因Hi-C测序失败被排除，最终六份样品纳入分析。每份样品离心富集微生物生物量后提取总DNA，分别构建鸟枪法宏基因组文库和Hi-C文库，前者在Element AVITI平台采用双端150 bp测序，后者在Illumina NovaSeq平台测序。鸟枪法读段经fastp质控后用MEGAHIT进行de novo组装，Hi-C读段经BWA-MEM比对和SAMtools过滤后，利用ProxiMeta平台基于Hi-C邻近连接信号将重叠群聚类为宏基因组组装基因组（MAGs）、质粒分箱和病毒分箱。MAG质量用CheckM评估，分类学注释用GTDB-Tk（GTDB release v220）；病毒分箱用CheckV评估；质粒分箱通过与PLSDB数据库比对鉴定。微生物MAG质量阈值为完整度≥70%且污染度≤10%，质粒和病毒分箱阈值为完整度≥50%且污染度≤20%。宿主归属使用ProxiPhage工具，基于Hi-C长距离连接将质粒和病毒分箱与宿主MAG关联，并施加严格的连接数、连接比和拷贝数阈值过滤。耐药、应激反应和毒力相关基因采用AMRFinderPlus v3.10.5鉴定，并与Hi-C宿主归属数据整合。随后构建宿主-质粒关联网络和基因中心AMR互作网络，利用Python工具包进行可视化和统计分析。
三、主要研究结果
3.1 Hi-C宏基因组解卷积回收污水进水中的微生物、质粒和病毒基因组分箱
六份可分析样品共回收12877个基因组分箱（微生物1087个、质粒8841个、病毒2943个），总装箱序列长度4.32 Gb。经质量过滤后保留454个高质量分箱，包括139个MAGs、130个质粒分箱和185个病毒分箱。分箱回收量随测序深度变化，WWTP-1的两份样品分箱数远高于WWTP-2的单个样品。病毒分箱中有42个达到CheckV Complete（闭合基因组）级别，17个达高质量级别。值得注意的是，未过滤的微生物分箱质量跨度很大，中位CheckM完整度仅10.3%，表明Hi-C解卷积能有效富集高质量基因组。
3.2 微生物MAGs、质粒和病毒分箱的分类学组成
139个高质量MAGs分布于70个细菌和古菌科，其中伯克霍尔德菌科（Burkholderiaceae）和巴氏杆菌科（Pasteurellaceae）各9个最为丰富，其次为气单胞菌科和双歧杆菌科（各7个）。130个质粒分箱中91个（70.0%）可归属细菌宿主科，肠杆菌科占绝对优势（55个分箱，42.3%），在除WWTP-1-Jan外的所有样品中均检出；双歧杆菌科次之（11个）。185个病毒分箱中仅11个（5.9%）获得宿主归属。尽管宿主归属率差异显著，MAGs、质粒相关宿主和病毒相关宿主间的分类学模式一致，肠杆菌科和气单胞菌科同时出现在MAGs和质粒宿主中，提示它们是污水中可移动遗传元件的主要储库。
3.3 耐药组组成及AMR、应激反应和毒力决定簇的宿主分布
AMRFinderPlus共鉴定771个基因命中（517个AMR基因、232个应激反应基因、22个毒力相关基因），归并为31个基因类别。耐药组以MLSB类（117个）和β-内酰胺类（114个）为主，四环素类（96个）、氨基糖苷类（65个）次之。应激反应基因以重金属抗性为主，汞抗性（73个）、铜抗性（36个）、砷抗性（28个）和铜/银抗性（26个）共占约70%。毒力基因仅22个，全部来自WWTP-1-Dec样品。基因组分布方面，应激反应基因的质粒定位比例最高（58.6%），而AMR基因（39.5%）和毒力基因（40.9%）相对较低。Hi-C宿主归属将490个基因命中（AMR 323个、应激反应155个、毒力12个）连接至细菌宿主科。肠杆菌科拥有最多的宿主归属基因命中（131个），但以应激反应基因为主（66%）；气单胞菌科含102个基因命中，其中AMR基因占72%，覆盖β-内酰胺、氨基糖苷、四环素、喹诺酮、甲氧苄啶和磺胺类，并携带29个汞抗性基因；伯克霍尔德菌科以汞抗性为主的应激反应基因富集。
3.4 AMR阳性质粒的Hi-C宿主-质粒关联网络
重建的宿主-质粒网络包含137个AMR阳性质粒分箱，通过166个Hi-C支持的关联连接至106个MAGs，形成79个连通组件。绝大多数质粒分箱宿主连接度低（平均连接度1.21），而微生物MAGs连接度更高（平均1.57，最大10），呈现不对称网络结构。气单胞菌科贡献最多的宿主-质粒关联（30个关联、涉及28个质粒分箱），肠杆菌科次之（19个关联、17个质粒分箱）。最高连接度的MAGs包括一个肺炎克雷伯菌（Klebsiella）MAG和一个Tolumonas MAG，各连接10个质粒分箱。五个质粒重叠群与同一份样品中两个不同细菌纲的MAGs关联，其中携带blaCMY-4的质粒（k99_1341003，AmpC型β-内酰胺酶）同时关联梭菌纲（Clostridia）和互养菌纲（Synergistia）的MAGs；另外两个质粒分别携带氨基糖苷抗性基因aph(3″)-Ib和aph(6)-Id及四环素抗性基因tet(C)，均与γ-变形菌纲和α-变形菌纲MAGs关联。此外，在WWTP-2-Nov中，一个病毒重叠群携带erm(B)，归属至链球菌（Streptococcus）相关MAG。
3.5 基因中心宿主-AMR关联网络
基于30个最丰富的AMR基因构建的基因中心网络包含144个宿主关联检出，其中86个（59.7%）为质粒关联，57个（39.6%）为染色体关联，1个为病毒关联。四环素类和MLSB类各37个检出，质粒关联比例分别为70%和59%；氨基糖苷类16个检出，69%质粒关联。磺胺类、喹诺酮类、甲氧苄啶类、mph(A)和blaCMY-4绝大部分为质粒关联，而blaOXA（86%染色体）、vanR-O（100%）和emrD（100%）则几乎只出现在染色体上。宿主纲层面，45.8%的检出归属γ-变形菌纲，其中67%为质粒关联。
四、讨论与结论
讨论部分强调，Hi-C宏基因组学将污水耐药组分析从基因清单推进至宿主和基因组背景层面，揭示了传统丰度谱无法反映的生态组织差异。宿主归属显示，粪污微生物组中优势菌门存在显著功能分区：肠杆菌科虽检出最多基因命中，但以重金属抗性等应激反应基因为主；气单胞菌科则携带跨多药物类别的广泛AMR基因谱，这一结果与此前Stalder等利用Hi-C在污水中的开创性研究一致，强化了气单胞菌科作为污水耐药组重要环境储库的证据。宿主-质粒网络呈不对称结构，提示移动性耐药主要由携带多个AMR质粒的细菌宿主组织，而非单个质粒在无关菌群间的广泛传播；五个跨纲质粒关联事件提示潜在共享质粒携带，其中blaCMY-4由于赋予第三代头孢菌素耐药且不受临床β-内酰胺酶抑制剂抑制，具有流行病学意义。病毒相关AMR证据有限，仅检出一个erm(B)，可能与环境病毒宏基因组的固有挑战有关。应激反应基因与质粒的高度关联支持重金属共选择维持适应性可移动遗传元件的生态模型。研究局限性包括样本量小、两例样品技术性排除、MAG回收依赖阈值且未用正交方法验证Hi-C关联。结论指出，Hi-C宏基因组学可将耐药决定簇与细菌宿主和基因组背景关联，为污水AMR监测提供生态学视角，并支持进一步评估宿主分辨宏基因组学作为污水AMR监测补充方法的潜力。
