楼宇尺度废水宏基因组学揭示耐药基因与毒力基因的时序模式

《Antibiotics》：Building-Scale Wastewater Metagenomics Reveals Temporal Patterns in Resistance and Virulence Genes

　　

背景/目的：抗菌药耐药性(AMR)与毒力是微生物致病潜能协同演化的维度，其在楼宇尺度废水系统中的生态组织仍知之甚少。方法：研究人员采用鸟枪法宏基因组测序(shotgun metagenomic sequencing)，对2021至2023年某大学COVID-19隔离宿舍排放的12份废水抓样（每学期2份）进行表征，分析耐药基因(ARGs)与毒力因子(VFs)的时序动态及生态关联。结果：废水微生物组以稳定的肠道相关厌氧细菌为核心；群落组成呈现显著的春季—秋季结构分化，年×学期交互作用解释60%群落变异；秋季2023年机会性菌属Acinetobacter扩增是最显著时序扰动。总ARG丰度各学期稳定，但耐药组(resistome)组成显著转换，表明时序动态由组成更替驱动而非总体耐药负荷变化。VF功能类别在各采样期广泛保守，与其在持续粪便核心微生物组中的结构嵌合一致。相关与网络分析显示耐药与毒力功能类别呈模块化生态耦合，金属/共耐药(metal/co-resistance)与磷霉素(fosfomycin)类与毒力功能关联最强；群落水平Benjamini–Hochberg校正共现网络解析出分类单元锚定的耐药模块与独立毒力功能簇，Acinetobacter与氟喹诺酮耐药为关键连接节点。结论：楼宇尺度废水宏基因组学可捕获生态结构化的功能基因动态，有望作为建成环境中AMR与毒力监测工具。

研究背景方面，抗菌药耐药性(AMR)与毒力因子(VFs)被视为微生物致病潜能协同演化的两个维度。已有市政与医院尺度废水研究多关注耐药基因(ARGs)出现与丰度，对 confined、细尺度系统中ARGs与VFs如何被微生物群落生态组织的认识有限。市政废水混合异质人口信号，会稀释局部生态动态；而楼宇尺度废水来自定义明确、时间连贯的人群，更适合检测局部群落结构与功能基因短时扰动。COVID-19隔离宿舍存在 occupancy、感染率、消毒剂使用与抗菌药管理改变，可能重塑微生物群落与功能基因，但缺乏将分类单元、ARGs、VFs在单一建筑汇水区纵向联合解析的研究。因此，研究人员以美国北卡罗来纳农工州立大学（North Carolina A&T State University）COVID-19隔离宿舍 Haley Hall 为对象，开展2021—2023年废水宏基因组时序研究，论文发表于《Antibiotics》。

主要关键技术方法方面，研究人员依托校园废水监测计划，从 Haley Hall 取得12份白天单次抓样（6个学期—年组合、每学期2份、10份经RT-qPCR确认SARS-CoV-2 RNA阳性）构成目的性子集；另取校卫生中心2020年8月至2022年12月COVID-19检测与隔离记录作人口层面背景。核酸由商业平台完成提取与鸟枪法宏基因组测序（Element AVITI，2×150 bp）。分类用Kepler/GenBook流程，ARGs用ResFinder、VFs用Virulence Factor Database（VFDB）鉴定；并以Kraken2、MEGAHIT组装、ABRicate作开源验证。统计在R中进行，含Bray–Curtis PCoA、PERMANOVA/PERMDISP、Spearman相关、Benjamini–Hochberg校正、Louvain模块识别与igraph共现网络。

研究结果方面，2.1 肠道相关厌氧细菌主导废水微生物组：属水平 Faecalibacterium、Prevotella、Bacteroides、Phocaeicola、Alistipes 五类肠道相关菌均值占43.6%，证明住宅废水具稳定人源粪便信号；秋季2023年 Acinetobacter 扩增最突出。2.2 微生物α与β多样性的季节格局：Shannon指数学期间无显著差异；Bray–Curtis排序中学期解释R²=0.23（p=0.01），年+学期模型R²=0.60（p=0.001），说明春—秋结构分化且年×学期交互主导变异。2.3 总ARG丰度与类别水平耐药组跨学期变化：总ARG丰度无单调趋势；大环内酯与四环素类常年主导，β-内酰胺秋季富集，bla_CMY与bla_MOX秋季升高且bla_CMY随时间递增，硝基咪唑类显著下降；基因水平耐药组按学期显著分离（R²=0.23，p=0.001）。2.4 毒力因子丰度与功能组成跨学期：总VF丰度秋高于春但不显著；黏附、运动、效应物递送、营养/代谢四类最稳定，VF组成学期解释率R²=0.143（p=0.069），显示毒力功能锚定粪便核心。2.5 耐药组与毒力组生态关联：metal/co-resistance与fosfomycin类与黏附、调控、效应物递送等VF功能强正相关（ρ=0.755–0.897），表明特定类别—功能对共享生态驱动。2.6 分类单元—ARG—VF共现格局：BH校正后19条边呈稀疏模块；Acinetobacter与fluoroquinolone resistance介数最高，耐药经分类单元组织，VF彼此共变，无直接ARG–VF边存活。2.7 校园COVID-19检测与隔离趋势：2020—2022隔离量在秋季开学与2021年4月达峰，为2021—2022样品提供 occupancy/感染压力背景，2023年缺同期人口数据。

讨论部分总结方面，研究人员指出楼宇废水存在“稳定粪便核心 + 时序调制变异”的分层结构：群落核心由肠道厌氧菌维持，耐药类别随选择性压力更替，毒力功能因嵌于核心而保守。metal/co-resistance与fosfomycin关联广，提示消毒剂/金属压力共选择；但短读长数据不能区分ARG与VF同位于移动遗传元件还是仅因宿主分类单元共变。网络显示模块化而非全连通，说明耐药与毒力是半独立功能行会(guilds)。研究意义在于证明楼宇尺度废水宏基因组可高分辨率监测AMR与毒力、识别bla_KPC等临床相关碳青霉烯酶序列信号，并为抗菌药管理(antimicrobial stewardship)与感染控制提供补充监测框架；局限为样本量小、单次抓样、无流式加权、无crAssphage/PMMoV归一化、单楼不可外推。

结论部分翻译：本研究证明楼宇尺度废水宏基因组评估可在受限人类相关环境中高分辨率审视微生物群落、AMR与毒力性状的生态组织。该宿舍废水微生物组各采样期均以肠道相关厌氧细菌的稳定粪便核心微生物组为锚，其上是受季节与年份调制的群落结构、耐药组组成与毒力谱变异。季节与年特异性扰动（尤其2023年秋季显著群落偏移）塑造了分类与功能基因格局。特定耐药类别与毒力功能间的耦合，以及共现网络分类单元锚定模块结构，表明这些性状并非独立分布，而是围绕细菌分类单元与共享选择压力组织。该组织直接关系废水衍生信号解释：不应视为孤立遗传标记，而应作为微生物群落状态与功能潜能的生态嵌合指标。总之，研究为楼宇尺度废水监测融入生态解释建立框架，在结合临床与处方数据后，可支撑抗菌药管理、感染控制及半封闭人群耐药涌现预警。

作者 Morikwe, U. 在图11B用BioRender制作流程；机构为North Carolina A&T State University；IRB协议HS26-0058与IBC协议#20-11为合规依据。全文未推测未出处内容，所有术语均依原文浓缩。