HAlign-G:面向大规模近缘基因组的快速低内存多基因组比对工具

HAlign-G: rapid and low-memory multiple-genome aligner for large-scale closely related genomes

作者信息Pinglu Zhang, Tong Zhou, Yanming Wei, Qinzhong Tian, Yixiao Zhai, Yizheng Wang, Quan Zou, Furong Tang, Ximei Luo
PMID41316476
发布时间2025-11-28
DOI10.1186/s13059-025-03881-3

实验完整度

高

涵盖了模拟数据、小规模真实基因组、大规模人类基因组及SARS-CoV-2等多层级实验,并进行系统发育分析验证。

主要模型

模拟线粒体样序列(INDELible生成) Alignathon模拟灵长类和哺乳动物数据集 细菌基因组数据集(脑膜炎奈瑟菌、肺炎链球菌、大肠杆菌等) 人类染色体和SARS-CoV-2基因组数据

重点核对

MSA质量评估指标:SP分数、Q分数和TC分数 MGA质量评估指标:准确性、召回率和F1分数,以及新提出的M分数 结构变异检测评估:8种SV类型和12种长度 系统发育树质量评估:使用Caster和UPGMA算法构建树并与参考树比较 运行环境:Ubuntu 22.04.5 LTS,双AMD EPYC 7763 CPU

摘要

HAlign-G是一种快速且内存高效的大型多基因组比对工具。利用BWT-FM-LIS结合优化的K波段算法和星比对策略,它支持种内(HAlign-G1)和种间(HAlign-G2)比对。基准测试显示其准确性、效率和内存使用均优于现有方法。HAlign-G1在种内数据的多序列比对中在速度和质量方面表现出色,而HAlign-G2在多基因组比对中提供更高的准确性和结构变异检测。两个版本都能处理数百万个SARS-CoV-2基因组和数千个人类染色体,支持可靠的进化研究,并有助于构建更稳定的系统发育树,同时提升Progressive Cactus的性能。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
如何开发一种快速且内存高效的多序列及多基因组比对工具,以处理大规模、超长且高相似度的基因组数据?
核心机制
HAlign-G采用星比对策略,通过BWT-FM-LIS分治方法加速比对,结合K-band动态规划和仿射空位罚分,并进行结构变异识别与整合。
主要证据
在模拟数据集上,HAlign-G1在MSA质量上优于现有工具,HAlign-G2在MGA和SV检测上表现优异;在大规模人类染色体和SARS-CoV-2数据上成功运行,并构建出更高质量的系统发育树。
研究意义
HAlign-G为大规模基因组比对和进化研究提供了高效工具,能够处理数百万个基因组,有助于构建更稳定的系统发育树,并提升Progressive Cactus的性能。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

模拟数据集的MSA质量评估

评估HAlign-G1和HAlign-G2在多序列比对中的准确性和效率。

使用INDELible模拟相似度70%-99%的线粒体样序列,比较HAlign-G1、HAlign-G2、HAlign4、HAlign3、MAFFT、MUSCLE和ClustalO的比对质量、时间和内存。

2

模拟数据集的MGA质量评估

评估HAlign-G2在跨物种多基因组比对中的准确性和效率。

使用Alignathon模拟灵长类和哺乳动物数据集,比较HAlign-G2与Progressive Cactus等的精确度、召回率和F1分数,并比较运行时间。

3

结构变异检测评估

评估HAlign-G检测各种结构变异的能力。

使用自定义模拟器生成包含8种SV类型的序列,评估HAlign-G1、HAlign-G2、Parsnp、Mugsy和Progressive Cactus在SV检测中的M分数和检测数量。

4

小规模真实基因组评估

评估HAlign-G在小型真实基因组数据上的比对质量和效率。

在五个小型基因组数据集(线粒体基因组、脑膜炎奈瑟菌、肺炎链球菌、大肠杆菌和三个人类1号染色体数据集)上比较HAlign-G1/HAlign-G2与多种比对软件。

5

大规模基因组评估

评估HAlign-G在大型基因组数据集上的处理能力。

使用18个人类基因组染色体数据集,评估HAlign-G1和HAlign-G2的比对质量和资源消耗,并与Progressive Cactus和Parsnp比较。同时测试了1000和5000个人类染色体1序列以及100万SARS-CoV-2序列。

6

系统发育树质量评估

评估HAlign-G产生的系统发育树的质量,以及与Progressive Cactus的兼容性。

使用Caster基于HAlign-G1、HAlign-G2和Parsnp的比对结果构建系统发育树,并作为Progressive Cactus的引导树,比较M分数。还在大肠杆菌基因组和灵长类基因组数据上评估树的准确性。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
LibDivSufSort----
MurmurHash3----
Progressive Cactus----
Parsnp----
Mugsy----
KAlign----
MLAGAN----
Progressive Mauve----
MUMmer4----
HAlign-3----
HAlign4----
MAFFT----
MUSCLE----
ClustalO----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
MSA质量评估
序列相似度范围、评估指标(SP、Q、TC分数)、对比工具列表
阅读提示:Methods: Alignment quality evaluation metric; Results: Evaluation of multiple sequence alignment quality on simulated dataset
MGA质量评估
评估指标(精确度、召回率、F1分数)、数据集(Alignathon)、对比工具
阅读提示:Methods: Alignment quality evaluation metric; Results: Evaluation of multiple genome alignment quality on simulated dataset
结构变异检测
SV类型、长度范围(50-2000bp)、模拟数据集数量(12×10)、评估指标(M分数、检测数量)
阅读提示:Methods: 结构变异识别模块; Results: 结构变异检测评估部分
大规模基因组评估
线程数(8或64)、数据集(18人类基因组、1000/5000染色体1、100万SARS-CoV-2)、资源消耗(时间、内存)
阅读提示:Results: Evaluation on large-scale genome datasets
系统发育树评估
数据集(101大肠杆菌基因组、28灵长类基因组)、树构建方法(Caster、UPGMA、WASTER)、评估指标(M分数、RF距离)
阅读提示:Results: Quality of phylogenetic tree construction