Nanopore 测序数据要不要质控?17 种预处理组合揭示变异检测与组装的取舍
先看结论
一项关于细菌基因组学的预印本基准研究比较了 17 种 Oxford Nanopore Technologies(ONT)测序数据预处理组合。研究显示,Nanopore 数据并不是“质控越多越好”:变异检测、从头组装和质粒恢复对读段处理的要求并不相同,流程选择应由下游任务决定。
ONT 测序能够产生较长的 DNA 读段,常用于细菌基因组的参考基因组变异检测和从头组装。这里的“质控”主要包括两类处理:一是去除接头、条形码等技术序列的 adapter trimming;二是根据读段长度或质量值筛选数据。前者可能去除残留的技术序列,后者可能提高输入读段的平均质量,但两者也都可能丢失对组装有价值的读段。
研究比较了哪些流程
研究围绕一个实际问题展开:随着 ONT basecalling(将测序仪采集的电信号转换为碱基序列)准确率提高,传统读段预处理是否仍然必要?作者比较了由三种接头处理工具 Dorado、Porechop_ABI 和 Barbell,以及三种质量处理工具 Filtlong、Chopper 和 SeqKit 构成的 17 种预处理组合,同时设置不去接头、不做质量处理的对照,并纳入同时执行两类处理的 Fastplong。
测试覆盖 13 个细菌分离株、100× 和 20× 两种测序深度,并分别使用 ONT 的 HAC 和 SUP 两种 basecalling 模型。变异检测以 F1 分数评价;F1 分数综合考虑精确率和召回率。组装评价则包括序列准确性、连续性、复制子恢复以及残留接头污染。这里的“复制子”可理解为染色体、质粒等能够独立复制的基因组单元。
组装的多个指标还被合并为加权算术平均值,权重来自微生物基因组学社区调查。这意味着“总体得分最高”不等于在每一个单项指标上都最好,实际选择仍取决于研究者更重视变异准确性、组装连续性,还是复制子和序列完整性的恢复。
研究观察到的主要结果
在 SUP basecalling 模型下,不做接头处理和质量处理的原始读段也能获得准确的 SNP 检测结果,并超过若干同时执行两类处理的流程。单独进行接头处理对 SNP 准确性没有明显影响。
按读段长度和质量联合筛选的 Filtlong 流程,在两种测序深度下通常取得最高的 SNP F1 分数;但在使用 Barbell 后再处理 20× 数据的组合中,Filtlong 的排名为倒数第二。这个例子说明,同一种质量过滤工具的效果也可能受到前一步接头处理方式和有效测序深度的影响。
组装结果为什么没有唯一赢家
对于从头组装,流程之间的权衡更加明显。每一种包含 Filtlong 的流程都产生了最低的错误率,但在 100× 和 20× 两种深度下均丢失了三个小质粒。对于需要完整恢复质粒的研究,这一结果意味着较低的碱基错误率并不自动代表更好的组装结果。
Barbell 接头处理提高了组装错误率,并在 20× 数据中降低了连续性。相反,跳过接头处理会使残留的接头和条形码序列进入组装结果。按照社区调查得到的权重汇总后,Porechop_ABI 加 Chopper 在两种测序深度下都取得最高总体得分,但它没有在任何单项指标上夺冠,因此不能据此断言它对所有项目都是最佳方案。
如何把结果转化为流程选择
如果目标是参考基因组 SNP 检测,SUP basecalling 下可以先把未处理读段作为基线:研究中,这类数据已经获得了准确的 SNP 结果,单独去接头也没有改变准确性。若更重视筛选后的读段质量,Filtlong 按长度和质量联合筛选在两种测序深度下通常取得最高的 SNP F1 分数,但需要同时检查它是否影响质粒或其他低丰度复制子的恢复。
如果目标是从头组装,则不应只看错误率或单一总体分数。应同时检查组装连续性、染色体和质粒等复制子的恢复情况,以及组装中是否残留接头和条形码序列。尤其在 20× 这类较低深度下,过滤掉短读段或低质量读段可能改善平均质量,却也可能进一步削弱小质粒和低覆盖区域的恢复。
实际分析时的操作顺序
建议保留原始 FASTQ 和 basecalling 模型信息,再分别建立“未处理”与候选质控流程。每条流程都应记录使用的接头处理工具、质量处理工具、输入读数、保留读数、总碱基数和读长分布;之后将同一批结果分别送入变异检测或组装,并独立检查 SNP F1、组装错误率、连续性、复制子恢复和残留技术序列。
需要复现这项比较时,不能只按工具名称照抄。所给研究摘要明确了 17 种流程的工具构成、13 个细菌分离株、100× 和 20× 深度,以及 HAC 和 SUP 两种模型,但没有列出各工具版本和逐项命令参数,也没有在所给资料中提供数据集或代码的公开入口。因此,摘要足以说明比较框架和主要结果,但不足以单独复现 17 种流程的全部细节。
结论与边界
这项研究的核心结论不是“应该使用某一个固定的质控工具”,而是质控策略必须服务于下游目标:变异检测可能受益于保留原始读段或采用质量过滤,组装则需要在错误率、连续性、接头污染和质粒恢复之间权衡。
研究比较的是 13 个细菌分离株、100× 与 20× 深度,以及 HAC 和 SUP 模型下的预处理组合。结果适合用来理解不同流程的取舍,不能直接视为对所有物种、文库类型、测序平台和数据批次的普遍保证。实际项目中,较稳妥的做法是保留未处理数据作为对照,并根据变异检测、完整组装或质粒恢复等目标选择和评估过滤策略。
原始来源: Evaluating Nanopore read trimming and quality filtering in bacterial genomics
摘要镜像: Life Science Network