基于污水流行病学的疾病暴发监测中方法学变异性的核算：以SARS-CoV-2为例

《ANALYTICAL AND BIOANALYTICAL CHEMISTRY》：Accounting for methodological variabilities in disease outbreak surveillance using wastewater-based epidemiology: a case study with SARS-CoV-2

　　

基于污水的SARS-CoV-2监测已成为在社区层面监测健康状况和疾病暴发的成熟方法。尽管RNA定量的分析方案已广泛实施并优化，但跨所有阶段（如储存、样品处理和PCR定量）引入的变异性表征仍然有限，包括其对流行病学解释的影响。本研究提出一种分析框架，以量化污水中SARS-CoV-2 RNA测量完整工作流程中的变异性，包括储存稳定性、样品处理和PCR定量。该框架基于一个大型SARS-CoV-2数据集开发，涵盖两个病毒基因靶标（N1、E-Sarbeco），并与人类相关CrAssphage一起定量，总计来自四个水处理回收中心（WRCs）24个月内的20,124个RT-qPCR数据点。分析揭示，样品冷藏比冷冻产生更丰富的流行病学数据，且鉴于CrAssphage变异性较高，其作为归一化生物标志物的效用极小。对SARS-CoV-2靶标分析变异性的量化促成了复合指标（Total Estimated SARS-CoV-2）的开发，以精度降低为代价提高了灵敏度并减少了未检出。这一新型多阶段框架能够改善WBE数据的灵敏度和解释以及WBE方法学的表征，对病原体监测具有广泛适用性。

研究背景与问题提出

SARS-CoV-2大流行凸显了加强疾病暴发早期预警系统的紧迫性，而疫情初期缺乏可扩展的社区级监测手段阻碍了有效追踪。污水流行病学（Wastewater-Based Epidemiology, WBE）通过检测市政污水中的核酸等标志物，能够以非侵入、近实时、低成本方式反映社区感染趋势，但检测流程涉及样品储存、处理、RT-qPCR定量等多个顺序阶段，各阶段均引入潜在变异性。既往研究多聚焦孤立步骤（如病毒回收率、检测限），缺乏对端到端累积变异性的系统量化，导致难以区分污水浓度变化是真实感染动态还是分析假象。该研究发表于《ANALYTICAL AND BIOANALYTICAL CHEMISTRY》，旨在建立多阶段分析框架，刻画RT-qPCR定量SARS-CoV-2 N1、E-Sarbeco及CrAssphage的完整方法学变异性，并评估储存条件对流行病学信息完整性的影响。

主要技术方法概述

研究人员于2020年4月至2022年2月从英格兰西南部四个水处理回收中心（Paulton、Radstock、Bath、Bristol对应 catchment A-D）采集797份进水污水每日样品（DS），采用24小时复合采样。每份样品分三份40 mL生物重复（BR）经60℃巴氏灭活、离心去除悬浮颗粒、PEG-8000沉淀浓缩、TRIzol^TM/氯仿提取核酸，以多重RT-qPCR（Agilent AriaMX）定量N1（US CDC 2019-nCoV_N1）、E-Sarbeco（Charité）及CrAssphage，每个BR的核酸作技术重复（TR）。以几何均值逐层汇总TR→BR→DS浓度，剔除未检出后计算；通过DS:BR、BR:TR、DS:TR比值分布量化处理变异性（PV）、分析变异性（AV）和方法学变异性（MV），以置信范围（CR=UCL/LCL）和线性调整%CV表示；储存稳定性考察+4℃（4周7时间点）、?20℃（1年10时间点）及最多4次冻融（FT）循环；流行病学完整性以UKHSA临床新冠病例数（每百万人病例数 cpm）对照污水靶标浓度，以250 cpm为切点算灵敏度/特异度；复合指标Total Estimated SARS-CoV-2按靶标%CV反比加权几何均值融合N1与E-Sarbeco。

研究结果

影响储存条件对基因靶标定量的影响

在+4℃下保存28天，N1、E-Sarbeco、CrAssphage均无总体降解趋势，标准差排序CrAssphage>E-Sarbeco>N1，未检出率低；?20℃保存1年各靶点出现大量未检出且变异性更高，E-Sarbeco与CrAssphage尤甚；FT循环使E-Sarbeco中位跌至未检出，N1相对耐受但仍具时点内变异。结论：冷藏（≤7天）优于冷冻，FT应尽量避免。

储存条件下的流行病学完整性评估

以UKHSA新冠cpm为参照，+4℃保存的污水N1与E-Sarbeco浓度和当日临床病例数呈线性（Pearson r N1≈0.84、E-Sarbeco≈0.68），95%特异度下N1灵敏度约86%、E-Sarbeco 74–77%；?20℃（FT=1）样品浓度聚左、无法区分病例负荷，N1阳性预测值从98%降至33%。结论：冷藏保留流行病学判别力，冷冻侵蚀早期预警价值。

变异性表征

分析变异性（AV，BR:TR）显示N1的CR最小、%CV最低，CrAssphage最大；处理变异性（PV，DS:BR）略小于AV但靶点排序一致，CrAssphage CR达8.32×致其浓度120–8320 gc/L波动，不适合作人群归一化标志物；方法学变异性（MV，DS:TR）为AV与PV叠加，低浓度DS未检出比例更高，证实未检出为“非随机缺失”。以MV_95%推得提取核酸LOQ_gc/μL=1×MV_95%/MV_50%，再乘浓度因子1250 μL/L得污水LOQ_gc/L。

用于疾病流行的复合指标

Total Estimated SARS-CoV-2以1/(%CV)²为权重合并N1与E-Sarbeco的TR几何均值，再几何汇总至BR与DS。该复合指标较单靶标提升低浓度检出、降未检出率，代价为精度下降（CR展宽）。

讨论与结论翻译

讨论指出，既往WBE研究依赖孤立步骤质控，该研究首次以原生靶标在真实污水矩阵中量化端到端变异性，证明（i）+4℃短期冷藏优于?20℃冷冻与FT，因后者引入RNA降解与回收损失；（ii）CrAssphage加工变异性过高，未经 catchment 内变异性评估不可用作归一化内标；（iii）以置信范围CR替代%CV更贴合RT-qPCR乘法型变异；（iv）复合指标以精度换灵敏度，适合低流行期暴发探测。研究结论如下：本研究开发的多阶段框架量化了污水SARS-CoV-2 RNA从储存、处理到RT-qPCR的全链条变异性；明确冷藏较冷冻保留更丰富流行病学信号，CrAssphage不宜作归一化 Biomarker；所建 Total Estimated SARS-CoV-2 复合指标在牺牲精度下提升灵敏度、压低未检出；该框架可推广至其他病原体WBE监测的方法学表征与数据解释。