叶绿体基因组组装工具的系统比较

A systematic comparison of chloroplast genome assembly tools

作者信息Jan A Freudenthal, Simon Pfaff, Niklas Terhoeven, Arthur Korte, Markus J Ankenbrand, Frank Förster
PMID32988404
发布时间2020-09-28
DOI10.1186/s13059-020-02153-6

实验完整度

中

研究包含模拟和真实数据的定量评估、一致性检验及新基因组组装,但未涉及生物学功能验证。

主要模型

拟南芥(Arabidopsis thaliana)模拟数据 369个真实测序数据集 105个无参考叶绿体的新数据集

重点核对

工具版本及提交哈希(表4) 输入数据(模拟或真实) 线程数(1、2、4、8) 内存限制(10 GiB) 时间限制(48小时)

摘要

背景:叶绿体是使植物能够进行光合作用的细胞内细胞器。它们通过原核细胞被真核宿主细胞共生整合而产生,并且仍然含有带有独特基因组信息的自身基因组。质体基因组包含必需基因,并经常用于生物技术或系统发育学。已经开发了能够评估质体基因组的不同组装器。这些组装器通常使用全基因组测序实验的数据,这些数据通常包含完整叶绿体基因组的读数。结果:不同组装工具的性能从未被系统比较过。在这里,我们提出了对七个叶绿体组装工具的基准测试,这些工具能够在超过60%的已知真实数据集中成功。我们的结果显示,测试的组装器在生成完整叶绿体基因组序列和计算需求方面存在显著差异。对来自未知质体基因组的物种的105个数据集的检查导致组装了20个新的叶绿体基因组。结论:我们为每个测试工具创建了docker镜像,这些镜像免费提供给科学界,并确保分析的可重复性。这些容器允许使用标准计算基础设施对叶绿体基因组数据集进行分析和筛选。因此,对基因组测序数据中的叶绿体进行大规模筛选是可行的。

实验结论

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

研究问题
现有开源叶绿体基因组组装工具在默认参数下,从全基因组测序数据中组装完整叶绿体基因组的性能和效率如何?
核心机制
不同组装工具在从WGS数据中提取叶绿体读段和组装环状基因组(包括IR区域)的策略上存在差异,导致性能差异。
主要证据
在模拟数据和真实数据集上,通过比对参考基因组计算得分,评估组装质量,并比较了运行时间、内存使用等。GetOrganelle在模拟和真实数据上均表现最佳,Fast-Plast次之。
研究意义
本研究提供了叶绿体组装工具的基准测试结果和推荐指南,并提供了Docker容器以提高可重复性,使得大规模叶绿体基因组筛选成为可能。

研究路径

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

1

模拟数据生成

创建无测序错误和生物变异的理想数据集,用于公平比较各组装工具的基本能力。

基于A. thaliana (TAIR10)叶绿体和核心基因组,使用seqkit滑窗方法生成不同基因组与叶绿体读数比例(0:1, 1:10, 1:100, 1:1000)和不同读长(150bp和250bp)的模拟数据,并包含2百万读段对版本。

2

真实数据收集

获取广泛绿色植物的真实测序数据,用于评估工具在真实数据上的性能。

从NCBI SRA中筛选符合条件(绿色植物、WGS策略、Illumina平台、双端、随机选择、公开)的数据集,每个具有CpBase参考叶绿体的物种选择一个数据集,共369个数据集。

3

组装工具执行

运行七个选定的组装工具,记录性能指标并生成组装结果。

使用统一的wrapper脚本定义输入文件(forward.fq, reverse.fq)、输出文件(output.fa)和线程数环境变量。所有工具封装在Docker镜像中,在HPC上通过Singularity运行。每个组装任务使用4线程、10 GiB RAM,时间限制48小时。

4

组装质量评估

量化各工具组装结果的质量,以比较其性能。

使用minimap2将组装结果与参考基因组比对,计算得分(包括完整性、正确性、重复分辨率、连续性)。对于新的数据,使用单重叠群、长度≥130kbp、IR≥17kbp作为成功标准。

5

一致性验证

评估各工具结果的稳定性。

对真实数据重复运行所有组装,并比较两次运行获得的得分。

研究方法

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

产品清单

实验环节名称品牌货号
GetOrganelle--587c1c51c34e270eb9178a42a77a5150157e6925
IOGA--c460ea9d9fe176fec2bd76d369b0cbb36793b2bf
NOVOPlasty--6af0894f8ea1d76a1b71df9cb762cf6e48dceac1
chloroExtractor--87364e48ec84a3f6ee91fc8d995b0bda5a0fa82d
Chloroplast assembly protocol--250d16ac02005d6a5939bf182b3d2995d0e88229
Fast-Plast--7e32b2e797fd1f49d32d6559e8345afefbaff803
ORG.Asm--830313acae3ca773b63f6bea9fc6d017e021bde5
minimap2--v2.16
seqkit----

关键环节

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

环节核对要点
组装执行
工具版本(提交哈希)
阅读提示:表4
组装执行
运行参数(线程数、内存、时间限制)
阅读提示:Methods: Our setup
数据准备
输入数据(模拟或真实)
阅读提示:Methods: Data
性能评估
评分标准
阅读提示:Methods: Quantitative