孟加拉国达卡医院废水中抗微生物耐药基因的宏基因组图谱与季节动态：基于鸟枪法宏基因组测序和纵向qPCR的跨季节研究

《Microbiology Spectrum》：Metagenomic profiling and temporal dynamics of antimicrobial resistance genes across hospital wastewaters in Dhaka, Bangladesh

摘要：孟加拉国医院废水（HWW）常未经充分处理即排入市政污水和水体，这引起了对医院可能作为临床重要抗微生物耐药基因（ARGs）储库和扩增点的担忧。然而，关于孟加拉国医院废水中这些基因组成、移动性和季节性行为的非培养依赖数据仍然有限。在此，研究人员将单季风时间点基线鸟枪法宏基因组测序与纵向定量实时PCR（qPCR）相结合，对达卡三家三级医院及附近一条社区污水线的废水在三个季节采样点进行了表征。基于基线宏基因组分析，各站点推定耐药组主要由氨基糖苷类、磺胺类、大环内酯-林可酰胺-链阳菌素（MLS）、四环素类和β-内酰胺类耐药性构成，主要与放线菌门（Actinomycetota）和假单胞菌门（Pseudomonadota）相关。医院出水含有推定临床重要基因，包括vanA簇、bla_VIM、bla_OXA-232、bla_IMP和bla_KPC。宏基因组分析将许多推定ARGs与质粒预测contigs、转座或整合模块相关联；毒力谱分析显示，与黏附、运动、铁获取和分泌相关的推定基因富集，尤其在MHWW-2中。来自基线宏基因组的探索性网络共现主要提示粪便来源的肠道共生菌以及病原体，作为ARGs输入废水系统的储库和推定载体。纵向qPCR显示，所选ARGs的绝对载量在初冬采样点较高，且1类整合子基因（intI1LC）在季节转换中追踪总ARG负荷。这些发现鉴定达卡医院废水主要为临床重要且潜在可移动耐药决定子的动态储库，并支持将其作为基于废水的抗微生物耐药性（AMR）监测靶点。
重要性：废水监测可提供AMR早期预警，但在孟加拉国，其应用很少超出培养或单基因PCR，限制了对新出现、可移动和设施特异性威胁的检测。本研究将宏基因组学与纵向qPCR整合，以连接存在哪些耐药基因、由谁携带以及它们在多所医院间如何季节性变化。研究人员表明，医院出水含有共享的背景耐药组，却在包括临床重要谱系的类群中不成比例地富集最后防线耐药性，且信号在初冬前上升。通过将高风险推定ARGs置于推定质粒、整合子和转座酶背景中，并与基线宏基因组分析所得的毒力库并列，研究人员的结果识别出可能促进其向医院外传播的合理基因组机制。这些发现初步将达卡医院出水定位为临床重要、可移动AMR的动态储库，并支持基于废水的流行病学用于持续监测和风险缓解。

论文主体解读

研究背景与问题：抗微生物耐药性（AMR）是重大公共卫生威胁，孟加拉国作为低收入和中等收入国家（LMIC），住院患者和普通人群中耐药感染均较普遍。2016至2023年，孟加拉国AMR增加11%，常规抗生素有效性下降82%。达卡人口超过1600万，政府医院压力大，医院获得性感染负担上升。医院废水（HWW）接收抗生素残留、生物杀灭剂和患者排泄物，可驱动耐药病原体选择和传播。既往孟加拉国研究多依赖培养、针对特定病原或单基因PCR，缺乏非培养依赖的全面耐药组和毒力组数据，也缺乏高风险耐药基因长期持续性和季节稳定性的数据。基于废水的流行病学（WBE）可非侵入性监测人群耐药水平。因此，研究人员结合鸟枪法宏基因组测序和靶向qPCR，调查达卡医院废水中AMR组成和季节动态。

研究开展、结论和意义：研究人员于2025年6月至11月，在三个季节采样点，对达卡南部三家政府三级医院及Azimpur社区污水线进行采样。基于单季风时间点基线宏基因组和纵向qPCR，研究人员发现医院出水较社区污水携带更广、更临床相关的推定耐药组，并含有与可移动遗传结构相关的ARGs；靶向ARGs绝对载量在初冬升高，1类整合子基因intI1LC与总ARG负荷同步变化。该研究将达卡医院废水初步定位为临床重要、潜在可移动AMR的动态储库，支持WBE用于持续监测和风险缓解。

关键技术方法（不超过250字）：样本队列来源为达卡南部三家政府三级医院，即Bangladesh Medical University Hospital（HWW-1）、Dhaka Medical College Hospital（HWW-2）、Mitford Medical Hospital（HWW-3），以及Azimpur社区污水（ACWW）；2025年6月、9月、11月分别代表季风、后季风和初冬。每个医院沿污水管线三个位置采集并混合。关键技术包括：单季风时间点基线鸟枪法宏基因组测序；Kraken2/Bracken进行分类；ShortBRED和RGI/CARD鉴定ARGs；PlasFlow/PlasClass及PlasmidFinder预测质粒；mobileOG-db鉴定移动遗传元件（MGEs）；VFDB鉴定毒力因子（VFs）；宏基因组组装基因组（MAGs）重建；纵向qPCR定量19个ARGs、intI1LC、intI2、intI3、转座酶和16S rRNA；常规PCR表征1类整合子CS区及基因盒；Spearman相关和Cytoscape网络分析。

研究结果：
样本采集（Sample collection）：2025年6月至11月，从孟加拉国达卡南部三家政府三级医院和Azimpur社区污水线采集HWW。Dhaka Medical College Hospital和Mitford Medical Hospital位于人口密集的旧达卡，患者量高。每个医院沿污水管线三个位置采集并混合成复合样本。
分类结构（单季风时间点基线鸟枪法宏基因组测序）：检测到37门、80纲、1311属、4849种。微生物组以Actinomycetota为主，占37.8%–76.4%，Pseudomonadota第二，占17.5%–25.8%。MHWW-2富集Bacillota和Bacteroidota。MHWW-2中ESKAPEE病原最高，包括Escherichia coli 2.64%、Klebsiella pneumoniae 1.97%、Enterococcus faecium 1.46%、Acinetobacter baumannii 0.72%。MHWW-1 Shannon多样性最低。
推定ARGs在废水中的流行（单季风时间点基线鸟枪法宏基因组测序）：ShortBRED鉴定218个非冗余推定AMR基因。三个HWW位点总ARG负荷高于ACWW，MHWW-2最高。氨基糖苷类、磺胺类、MLS、四环素类和β-内酰胺类约占总ARG丰度82%。sul1和qacEΔ1高丰度。检出vanA簇、bla_VIM、bla_OXA-232、bla_IMP、bla_KPC；mcr-5在MHWW-2和MHWW-3检出；tet(X)在所有位点低丰度。RGI鉴定132个推定ARGs，覆盖23类。MHWW-2β-内酰胺负荷最高，达14.13%，主要贡献者包括bla_OXA-205、bla_KPC-2、bla_NPS-1、bla_NDM-1、bla_OXA-1011、bla_CTX-M-15；bla_GES-5、bla_OXA-181、bla_IMP-26、bla_IMP-62仅在该位点检出。bla_NDM-1位于质粒预测contig并邻近转座/整合模块，但未进行直接水平转移实验。ACWW核心耐药组包括sul1、mrx、mphA、tet(A)、tet(C)、tet(O)及少量β-内酰胺酶。
推定ARGs与推定MGEs的基因组共现（单季风时间点基线鸟枪法宏基因组测序）：约三分之一推定ARGs位于质粒预测contig。鉴定49个多ARG contig，MHWW-2约占一半。β-内酰胺耐药与氨基糖苷或大环内酯耐药常见共现。临床重要推定ARGs常嵌入MGEs，如bla_KPC-2、bla_NDM-1、bla_VIM-4、bla_IMP-2c邻近tnp、tnpA、tnpR、int及ISKpn6、iSKpnG。MHWW-2和MHWW-3含完整vanA操纵子模块。qac基因常与ARGs共现，提示共选择。共鉴定29个质粒复制子标记，MHWW-2质粒谱最复杂。
推定VFs分布与丰度（单季风时间点基线鸟枪法宏基因组测序）：鉴定超过200个推定VFs，主要与运动和表面附着相关，包括IV型菌毛、pilG、pilI、pilT、pilQ、鞭毛和趋化蛋白。MHWW-2富集黏附、生物膜、铁获取和分泌相关模块，包括1型菌毛、curli、mrkB/mrkC、enterobactin系统entC/entD/entS及T2SS、T4SS成分。MHWW-3含east1毒素相关信号。
推定ARGs的分类宿主（单季风时间点基线鸟枪法宏基因组测序）：ARGs主要与Actinomycetota和Pseudomonadota相关。ACWW和MHWW-1以Actinomycetota相关ARGs为主，分别约45%和62%；MHWW-2转向Pseudomonadota和Bacillota，分别34%和31%。MHWW-2中ARGs集中于Escherichia、Enterococcus、Streptococcus、Acinetobacter和Bacteroidales。bla_NDM-1常仅解析到Gammaproteobacteria；bla_CTX-M-15、bla_OXA-181、bla_IMP-26解析到Enterobacterales；bla_IMP-62归属Pseudomonas；vanA簇主要与Enterococcus相关。探索性网络基于n=4基线样本，不能支持正式统计推断。
宏基因组组装基因组（MAG）（单季风时间点基线鸟枪法宏基因组测序）：重建并保留92个MAGs，覆盖53属。MHWW-2含7个推定ARG阳性MAGs和30个耐药基因，最高。ACWW和MHWW-1较少，MHWW-3仅1个四环素耐药相关MAG。vanA阳性MAG暂归Bacteroidales；MHWW-2中一MAG含bla_CTX-M-15和qnrS；Streptococcus相关MAG含多个β-内酰胺和大环内酯耐药基因；Escherichia相关MAG含多外排耐药决定子。ARG-VF共携带仅见于MHWW-2部分MAGs。
临床重要ARGs的时间动态（纵向qPCR监测）：19个ARGs定量显示，总ARG负荷从季风到初冬绝对增加，并伴随细菌生物量上升。HWW-1从4.66×10²增至4.37×10⁵，HWW-2从1.24×10⁴增至6.22×10⁶，HWW-3从3.0×10¹增至4.55×10⁵。HWW-1后季风相对负荷峰，HWW-2和HWW-3初冬相对富集最高。bla_TEM、ermB、floR、lnuB、tetM在所有位点和时间点持续检出。bla_KPC和bla_NDM初冬升高。HWW-2初冬dfrA12约5.8×10⁶，占该位点总ARG负荷94%。HWW-3从季风有限谱转为初冬广泛高负荷，bla_KPC峰值约1.6×10⁵。
MGEs的时间动态与表征（纵向qPCR监测）：intI1LC绝对丰度从季风到初冬增加数个数量级。HWW-1和HWW-2分别增加373倍和53倍；HWW-3从不可检测升至1.5×10⁵。intI2和intI3在季风低，后季风和初冬可检出并增加。基于三个医院三个季节点的qPCR数据，intI1LC与靶向ARG负荷显著相关，ρ=0.93，P<0.001；bla_NDM和bla_KPC分别与intI1LC相关，ρ=0.95和0.92。常规PCR显示后季风和初冬所有样本及HWW-2季风样本class 1 integron CS区阳性，条带700–1500 bp；巢式PCR提示基因盒含sul1、bla_OXA、dfrA1。转座酶Tnp4、Tnp5、Tnp7、Tps14所有样本检出，Tnp3仅后季风和初冬检出。

讨论总结与结论翻译：讨论指出，孟加拉国HWW处理不足，常直接排入水体并与社区污水混合，增加ARG传播可能。本研究显示医院出水较社区污水有更广、更临床相关的推定耐药组，并含可移动结构相关ARGs。微生物组以Actinomycetota为主，MHWW-2存在粪便来源严格厌氧菌，MHWW-3存在Methanothrix，提示停滞和沉积。碳青霉烯酶和ESBLs提示最后防线耐药储库。多ARG共定位提示共选择，MGEs可能作为稳定平台。vanA常与Enterococcus相关，但MAG分析在Bacteroidales也发现，原文提出可能为更广宿主、移动或技术伪影，但表明存在高风险潜在可移动vanA模块。ARG负荷与医院规模和患者量相关，MHWW-2最高。各医院耐药谱不同，HWW非同质暴露源，WBE可比较设施特异性热点。纵向qPCR显示初冬靶向ARG升高，可能由临床负荷增加和干燥少雨稀释减弱所致。整合子追踪总ARG负荷。局限包括：宏基因组仅季风基线，后续季节分析依赖靶向qPCR；短读长限制MGE分辨率，质粒分配为推定；16S rRNA归一化未考虑类群拷贝数变异，标准曲线采用代表性策略；仅一个社区参考；未直接评估季节稀释和环境混合。
研究结论翻译：基于这些发现，研究人员为LMIC提出两层监测框架。该框架包括通过核心标志物靶向qPCR进行高频哨点监测，核心标志物包括人为污染与移动组配对以及高后果ARGs；辅以周期性宏基因组测序以更新qPCR panel并解析新出现变体的基因组背景。研究人员的数据进一步表明，在先进废水处理尚未普及的地区，优先在高信号哨点节点（如HWW-2）进行现场缓解，对下游环境保护具有不成比例收益。WBE在发达国家广泛用于追踪病原暴发和耐药组特征，但在大多数发展中国家仍有限。据研究人员所知，本研究首次在孟加拉国纵向表征HWW中靶向ARG的季节变化。通过连接高分辨率基因组数据与实用监测，本研究为实施可持续AMR监测提供了实践框架。