《Applied and Environmental Microbiology》：Whole-genome surveillance supports hazard profiling of Escherichia coli lineages in recycled water treatment systems
编辑推荐：
在气候变化背景下，使用处理后的废水对于可持续水资源管理日益重要，然而基于大肠杆菌（Escherichia coli）计数的常规监测手段，在揭示菌株多样性和相关公共卫生危害方面提供的信息有限。本研究应用纵向全基因组测序（whole-genome sequenci
在气候变化背景下，使用处理后的废水对于可持续水资源管理日益重要，然而基于大肠杆菌（Escherichia coli）计数的常规监测手段，在揭示菌株多样性和相关公共卫生危害方面提供的信息有限。本研究应用纵向全基因组测序（whole-genome sequencing, WGS）技术，对一个再生水设施从进水到最终出水的整个处理流程中收集的180株大肠杆菌分离株进行了分析。基因组分析揭示了广泛的菌株水平异质性，包括跨越8个系统发育群的88种序列类型（sequence types, STs），且进水的多样性高于处理后的出水。与澳大利亚背景基因组的系统发育比较表明，这些菌株与来自伴侣动物、野生鸟类、人类和牲畜的菌株聚类，提示存在多个潜在来源库而非单一主导来源，尽管来源贡献并非确定性的。尽管总大肠杆菌载量减少了>90%，但从上游和下游阶段回收的分离株展现出大致可比的毒力因子（virulence factor, VF）和抗生素抗性基因（antimicrobial resistance gene, ARG）图谱，这表明在所培养的分离株集合中，丰度的降低超过了基因组组成的改变。为了评估其运行相关性，研究人员构建了一个原型化的基因组学知情危害框架，整合了毒力决定因子、ARGs、质粒相关的移动性以及特定回用场景的暴露背景。使用该框架，92.8%的分离株被归类为低危害，7.2%为中等危害，未发现符合高或严重危害分类标准的分离株。这些发现表明，对指示生物进行基因组分析可以揭示常规计数法无法捕捉到的种群结构和危害异质性，并为将基因组信息纳入再生水系统的危害知情监测提供了实践基础。
论文解读文章
研究背景与目的
在全球气候变化和人口增长的驱动下，水资源需求日益增加，推动了对废水回收利用技术的广泛应用。再生水在农业灌溉、工业及环境补给等领域扮演着越来越重要的角色。为确保其微生物安全性，传统监测方法主要依赖对大肠杆菌（Escherichia coli, E. coli）的培养计数，以此作为粪便污染和处理效能的监管指标。然而，大肠杆菌物种内部存在显著的基因组多样性，既包含环境持久性的共生菌株，也携带临床相关的毒力因子（virulence factor, VF）和抗生素抗性基因（antimicrobial resistance gene, ARG）。因此，相同的菌落计数可能对应截然不同的公共健康风险。现有的常规监测无法区分这些具有不同潜在危害的遗传谱系。为解决这一问题，研究人员提出将全基因组测序（whole-genome sequencing, WGS）技术整合到再生水系统的监测中，以期获得比传统二元合规性指标更高的分辨率，从而更准确地评估微生物危害。
研究方法概述
本研究对澳大利亚维多利亚州一座日处理量约600兆升（ML）、年产再生水约30,000 ML的污水厂进行了为期一年半（2024年1月至2025年6月）的纵向采样。研究人员从三个关键点位采集水样：进入再生水处理系统的进水（W1）、经过紫外线（UV）消毒或UV加氯消毒后的中间处理水（W2），以及最终储存在接触池中准备用于A类再生水回用的出水（W3）。通过对所有样品进行培养，共获得180株大肠杆菌分离株。关键技术方法包括：
- 1.
全基因组测序与组装：采用牛津纳米孔（Oxford Nanopore）平台进行长读长测序，并使用Flye和Medaka等工具进行de novo组装和质量评估。
- 2.
群体基因组学与系统发育分析：通过多位点序列分型（multi-locus sequence typing, MLST）和Clermont分型确定序列类型（sequence type, ST）和系统发育群。利用核心基因组单核苷酸多态性（single nucleotide polymorphism, SNP）构建最大似然系统发育树，并与来自澳大利亚的121个背景参考基因组（相似度≥99%）进行比较，进行潜在来源溯源。
- 3.
功能基因注释：使用ABRicate工具分别比对VFDB数据库和CARD数据库，鉴定毒力因子和抗生素抗性基因。利用MOB-suite软件预测质粒及其移动性类别。
- 4.
基因组危害评估框架构建：开发了一套综合评分体系，对毒力潜力、ARG临床危害等级、多重耐药广度（multidrug resistance, MDR）和质粒介导的水平转移潜力进行量化，并结合不同回用场景的暴露途径权重，形成最终的暴露调整危害指数（exposure-adjusted hazard index, EAHI）。
主要研究结果
空间-时间动态与环境驱动因素
- •
结论：沿处理工艺，大肠杆菌载量显著下降，进水（W1）浓度最高，处理后（W2, W3）浓度极低。季节变化对大肠杆菌载量的影响不显著。总体数据中，大肠杆菌浓度与pH呈正相关，与温度呈弱负相关，但这些关联在单个采样点内并不显著，表明处理过程是影响大肠杆菌丰度的首要因素。
基因组多样性与潜在来源归因
- •
结论：180株分离株显示出极高的基因组异质性，包含88种STs和8个系统发育群。其中，B1群占主导地位（66.7%）。进水（W1）的ST丰富度（56个）高于出水（W3，42个）。通过与澳大利亚背景菌株的比较，研究人员发现这些分离株最常与伴侣动物（尤其是犬类，37.8%）和野生鸟类（24.4%）来源的菌株聚类，其次是人类和家畜，表明再生水系统中的大肠杆菌可能源自多种宿主库，而非单一优势来源。出水（W3）中与野生鸟类相关的谱系比例增加，可能反映了处理后环境的生态选择或与周边野生动物的相互作用。
组成与种群结构
- •
结论：进水和出水之间的系统发育群组成和ST分布存在显著差异。进水（W1）的α多样性更高。ST162和ST155是两个最常见的ST。单核苷酸多态性（SNP）分析显示，ST162中存在跨位点和时间的近相同基因组（≤1个SNP），暗示短期存留或反复引入；而ST155则表现出高度的基因组分歧，缺乏紧密相关的跨位点配对，更符合反复引入的模式。
毒力决定簇与抗生素抗性谱
- •
结论：共检测到195种不同的VF和47种ARGs。尽管位点间的基因丰富度有差异，但每个基因组携带的VF和ARG数量中位数在不同处理阶段（W1和W3）之间无显著差异。这表明处理过程虽然大幅降低了细菌丰度，但并未显著改变存活菌群中每株菌所携带的毒力和抗性基因负担。几乎所有分离株都携带外排泵和β-内酰胺类抗性决定簇，而其他临床意义更大的抗性类别（如喹诺酮类、磺胺类）检出率较低且在处理前后相似。
基因组危害评估框架
- •
结论：应用综合基因组危害评估框架后，绝大多数分离株（92.8%）被归类为低危害，仅有7.2%（13株）被归类为中等危害，未发现高或严重危害菌株。值得注意的是，中等危害菌株在出水（W3）中的比例（14.3%）显著高于进水（W1，1.2%）。这一结果表明，尽管低危害菌株在所有点位均占主导，但处理过程可能选择性地富集了某些具有更高基因组危害潜力的谱系。
讨论与结论总结
讨论部分总结：研究表明，传统的大肠杆菌计数法无法揭示菌株水平的多样性、持久性以及与公共健康相关的基因组特征。基因组分析能够区分低风险的背景检出与携带较高危害潜力的菌株。处理工艺是控制大肠杆菌丰度的主要因素，但对菌群基因组组成的影响相对较小。源追踪结果显示菌群来源复杂，突出了再生水系统受多源输入影响的特点。本研究所开发的基因组危害框架是一个有效的筛选工具，能够识别出常规监测无法区分的、具有潜在更高危害的菌株，为优化处理性能和风险管理提供了新思路。研究也存在局限性，例如单一场址的研究设计、各采样点分离株数量的不均衡，以及该框架尚未经过流行病学数据的验证等。
研究结论部分翻译：总体而言，对再生水系统中大肠杆菌的基因组分析揭示了高的谱系多样性、有限但可检测的位点间结构分化，以及以低危害为主的种群，同时在各处理阶段检测到一小部分中等危害的分离株。这些发现表明，基因组解析分析提高了对大肠杆菌检出的解读能力，超越了传统的计数方法，并为将基因组数据纳入再生水系统的危害知情监测提供了实践基础。
