基因组大数据计算: CPU 和 GPU 加速方案深度评测
随着 NGS 测序通量的大幅提高,搭配高效 NGS 二级分析技术的精准解决方案快速融进基因组学的各个应用领域:遗传进化、临床诊断、分子育种、医药开发等。以下我们通过对基于 CPU 和 GPU 不同硬件平台的 NGS 二级分析方案进行详细评测,以期为基因组学研究领域的用户提供参考。
Sentieon 软件是面向 CPU 平台开发的,在不需要专用的编程语言,不依赖任何专用硬件的情况下进行快速基因变异检测分析,大幅降低了软件的配置、部署和维护成本。同时,Sentieon 针对第三代和第四代 Intel Xeon Scalable 处理器(以前代号为 Sapphire Rapids)进行了改进及优化,更适用于在多核系统上扩展,可以在 Xeon 平台上实现更快速地处理急诊科及 ICU 的紧急样本,也可在实验室环境下以更高的通量处理大量样本。
本次评测将比较面向 CPU 平台的 Sentieon 软件及面向 GPU 平台的 NVIDIA Clara Parabricks 的运算性能和分析准确性,以评估两者在基因组二级分析中的性价比及成本效益。
为什么没选择 GATK
GATK 是变异检测流程的行业金标准,但它是用 Java 写的,所以不是性能的金标准。伊利诺伊大学和梅奥诊所发表的文章已经证实,Sentieon 的运算性能明显优于 GATK,而且没有损失精度(https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6710408)。
同时基于最新来自 AWS 的官方评测,NVIDIA Clara Parabricks 相对 GATK 也有明显的性能优势(https://aws.amazon.com/cn/blogs/hpc/benchmarking-the-nvidia-clara-parabricks-germline-pipeline-on-aws/)
综上原因,我们没有将 GATK 纳入评测范围。我们的目标是比较 Sentieon 软件(用 C++编写,并针对 CPU 优化)与 Parabricks(用 CUDA 编写并针对 NVIDIA GPU 优化)。
测试项目
为了尽可能精确地比较 Sentieon 与 Parabricks 流程每个计算阶段的运行性能,我们按照两个流程中匹配的计算步骤(如下表),将 Parabricks 中 haplotypecaller、预处理和 fq2bam 各个阶段对应 Sentieon 流程的步骤标记出来,并统计两个流程对应计算阶段的运算时间。测试数据为 30X HG001 WGS 标准数据集。
性能评测
以下我们将展示 Sentieon 和 Parabricks 在不同硬件配置环境下的性能表现。其中 Intel Xeon 可扩展处理器提供了具有竞争力的性能,第四代 Intel Xeon 可扩展处理器(以前的代号为 Sapphire Rapids)提供最佳的整体性能。除了最重要的性能评测,我们详细对比了每个基因组的计算成本和功耗表现。
3rd Gen 硬件环境为 2* Intel 2.4GHz Intel Xeon Platinum 8368 CPU(152 核心, 超线程开启),256GB DDR4-3200 内存,1TB Intel 660p 和 2TB DC P4510 SSD。
4th Gen 硬件环境为 2* 4th Gen Intel Xeon Scalable CPU(原代号为 Sapphire Rapids,>40 核,超线程开启),Intel Pre-production BIOS,256GB DDR 内存(16(1DPC)/16 GB/4800 MT/s),1TB Intel D3-S4610 SSD。
3rd Gen 和 4th Gen 的操作系统均为. Ubuntu Linux* 20.04。
通过下表中各计算实例上 Sentieon vs. Parabricks 的性能对比,可以看出,3rd Intel Xeon 平台可在 40 分钟左右完成 30x WGS 的数据分析,与 GPU 平台速度相当。而 4th Intel Xeon 平台则跑出了 26.8 分钟的最佳成绩。
其中 NVIDIA Clara Parabricks 评测数据来源于 AWS 官方博客:https://aws.amazon.com/cn/blogs/hpc/benchmarking-the-nvidia-clara-parabricks-germline-pipeline-on-aws/
根据不同配置的硬件定价统计表可以发现,与 NVIDIA A100 Tensor Core 处理器(4.59 美元)相比,Intel Xeon Scalable CPU(1.54 美元)的每个基因组分析所需成本要低得多。如果第 4 代英特尔至强可扩展处理器具有类似的 AWS EC2 定价,则每个基因组的计算成本将降至不到 1 美元(2.1635 美元/小时 * 26.8 分钟 = 0.97 美元)。在此补充一下,这次基准测试中使用的 4th Gen Xeon Scalable CPU 是预发布硬件,因此最终产品的性能可能会有所提高。
在功耗方面,c6i.metal 实例中的两颗 Intel Xeon Platinum 8352M 处理器的功率是 370W,而 p4d.24xlarge 实例中的八个 NVIDIA A100 Tensor Core 的功率则高达 3,200W。Parabricks 要达到最佳性能,相较于 Intel 纯 CPU 硬件环境需要 8.6 倍的功率和 3.0 倍的成本,但分析效率相较于 Intel 3rd Gen Xeon Platinum 8352M CPU 来说,仅有 1.5 倍的性能提升。
精度评测
在分析精度方面,Parabricks 官方博客此前报道了其变异检测结果的准确性与 GATK 相当(F1 scores)。而 Sentieon 不但提供与 GATK 一致结果(一致性达 99%以上),更是 PrecisionFDA Truth Challenge 的连续冠军。在第二届 PrecisionFDA Truth Challenge 比赛中更是赢得了 4 个分项赛冠军。
关于对 Sentieon 软件的运行效率和准确度的行业认可,其遍布全球的用户和快速增长的数据处理量提供了广泛的例证。
结论
Sentieon 软件是通过改进算法模型实现性能加速(纯 CPU 环境,支持 X86/ARM),不依赖于昂贵高功耗的专用硬件配置(GPU/FPGA),不依赖专有编程语言;同时 Sentieon 软件针对几乎所有的短读长和长读测序平台进行了优化,是 FDA 多次公开挑战赛的连续赢家。本次评测展现了 Sentieon 软件在 Intel Xeon 平台上的卓越性能,是基因组二级分析的最佳解决方案。
评论