模拟数据生成
创建无测序错误和生物变异的理想数据集,用于公平比较各组装工具的基本能力。
基于A. thaliana (TAIR10)叶绿体和核心基因组,使用seqkit滑窗方法生成不同基因组与叶绿体读数比例(0:1, 1:10, 1:100, 1:1000)和不同读长(150bp和250bp)的模拟数据,并包含2百万读段对版本。
A systematic comparison of chloroplast genome assembly tools
研究包含模拟和真实数据的定量评估、一致性检验及新基因组组装,但未涉及生物学功能验证。
拟南芥(Arabidopsis thaliana)模拟数据 369个真实测序数据集 105个无参考叶绿体的新数据集
工具版本及提交哈希(表4) 输入数据(模拟或真实) 线程数(1、2、4、8) 内存限制(10 GiB) 时间限制(48小时)
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
创建无测序错误和生物变异的理想数据集,用于公平比较各组装工具的基本能力。
基于A. thaliana (TAIR10)叶绿体和核心基因组,使用seqkit滑窗方法生成不同基因组与叶绿体读数比例(0:1, 1:10, 1:100, 1:1000)和不同读长(150bp和250bp)的模拟数据,并包含2百万读段对版本。
获取广泛绿色植物的真实测序数据,用于评估工具在真实数据上的性能。
从NCBI SRA中筛选符合条件(绿色植物、WGS策略、Illumina平台、双端、随机选择、公开)的数据集,每个具有CpBase参考叶绿体的物种选择一个数据集,共369个数据集。
运行七个选定的组装工具,记录性能指标并生成组装结果。
使用统一的wrapper脚本定义输入文件(forward.fq, reverse.fq)、输出文件(output.fa)和线程数环境变量。所有工具封装在Docker镜像中,在HPC上通过Singularity运行。每个组装任务使用4线程、10 GiB RAM,时间限制48小时。
量化各工具组装结果的质量,以比较其性能。
使用minimap2将组装结果与参考基因组比对,计算得分(包括完整性、正确性、重复分辨率、连续性)。对于新的数据,使用单重叠群、长度≥130kbp、IR≥17kbp作为成功标准。
评估各工具结果的稳定性。
对真实数据重复运行所有组装,并比较两次运行获得的得分。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 组装工具 | 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 |