Oatk:利用HiFi Read进行细胞器组装

Oatk是我见过的一个非常实用、高效的植物细胞器基于PacBio HiFi测序的基因组组装工具。这里有两个限定语,一个是植物,一个是PacBio HiFi。如果不满足这两个需求,那就无法使用这个工具。

Oatk包括三个组件,syncasm用于从HiFi read中组装细胞器基因组。hmm_annotation会利用OatkDB的HMM profile对组装结果进行组装。最后的pathfinder 会根据注释信息和组装结果,得到环状的DNA序列。

Oatk软件编译本身仅有一个依赖,zlib,下载源码,就可以编译安装。

git clone https://github.com/c-zhou/oatk.git
cd oatk
make -j

在实际运行过程中,由于需要用到nhmmscan ,因此还需要额外安装了HMMER,可以使用conda安装。

另外还需要一个OatkDB,里面是作者预处理过的HMM profile,

git clone https://github.com/c-zhou/OatkDB.git

你可以设置一个环境变量,hmm_db_dir,用于后续分析。

hmm_db_dir=/path/to/OatkDB/v20230210/

Oatk的运行只需要一个命令,如下

input_hifi=/path/to/your/hifi_reads.fq

oatk -k 1001 -c 150 -t 8  \
    -m $hmm_db_dir/angiosperm_mito.fam \
    -p $hmm_db_dir/angiosperm_pltd.fam \
    -o oatk $input_hifi

其中-k和-c用于设置syncmer大小和深度。所谓的syncmer,类似Minimizers ,用于从序列中高效的选取K-mers(也就是长度为K的序列)。因此-k 1001 表示,kmer长度最小为1001 bp, -c 150 表示该类型的kmer至少要出现150次。

在后续自己分析设置参数时,-k可以保持1001,-c必须根据实际的测序深度来,作者推荐该值在平均测序深度的5-10倍,用于有效的过滤核基因组序列,确保只筛选出细胞器的序列。

-t 时线程数,根据实际的情况调整。

-m和-p对应的你实际物种的线粒体和叶绿体的HMM profile。

-o则是输出文件名的前缀。

最终的输出结果为

  • oatk.mito.fasta: 线粒体的fasta
  • oatk.pltd.fasta:叶绿体的fasta

整体来说,oatk和作者之前开发的工具一样,非常的简洁,高效,并且实用。正如他所用的编程语言C一样。如果你用Hi-C用来搭基因组的话,强烈推荐用他的yahs。

# 基因组 

评论

Your browser is out-of-date!

Update your browser to view this website correctly. Update my browser now

×