温馨提示×

温馨提示×

您好,登录后才能下订单哦!

密码登录×
登录注册×
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》

如何使用GetOrganelle软件组装叶绿体基因组

发布时间:2021-07-19 11:40:56 来源:亿速云 阅读:921 作者:chen 栏目:大数据

本篇内容主要讲解“如何使用GetOrganelle软件组装叶绿体基因组”,感兴趣的朋友不妨来看看。本文介绍的方法操作简单快捷,实用性强。下面就让小编来带大家学习“如何使用GetOrganelle软件组装叶绿体基因组”吧!


现在做植物的叶绿体基因组基本上都是直接以新鲜叶片做材料,提取总DNA测序,构建二代测序文库,然后利用现成的软件组装叶绿体基因组,省去了提取叶绿体的步骤

利用总DNA测序的数据组装叶绿体基因组的软件很多,有一篇综述介绍这些工具

https://genomebiology.biomedcentral.com/articles/10.1186/s13059-020-02153-6如何使用GetOrganelle软件组装叶绿体基因组

今天这篇推文我使用 GetOrganelle 这个软件,软件的Github链接 https://github.com/Kinggerm/GetOrganelle

对应的论文

https://genomebiology.biomedcentral.com/articles/10.1186/s13059-020-02154-5如何使用GetOrganelle软件组装叶绿体基因组

软件的github主页对软件的使用方法介绍的很详细,这款软件是昆明植物所的老师开发的,还开设了qq群进行答疑,qq群号在gitbub的主页末尾可以看到

安装直接使用conda,非常容易

conda install -c bioconda getorganelle

想要在linux系统下使用conda命令需要安装miniconda或者Anaconda3,我之前录制过一期视频介绍linux系统安装 Anaconda3,不熟悉的可以找来看看 https://www.bilibili.com/video/BV1Ft4y1e7w2

软件安装好以后还得下载参考基因组,运行如下命令如何使用GetOrganelle软件组装叶绿体基因组

组装用到的命令是get_organelle_from_reads.py,使用到的参数有

  • -1 -2 分别制定双端测序数据的路径

  • -o 制定输出文件的文件名

  • -R -k 具体是什么意思我还不知道,-R按照帮助文档直接设置15应该就可以,-k后面接的数字也可以按照他帮助文档的来设置,现在双端测序通常是150bp,这个-k参数直接设置105和121就可以,这个数字少一点,速度应该会快一点。

  • -F 指定参考,如果是叶绿体基因组后面直接跟 embplant_pt 就可以

因为叶绿体基因组的拷贝数高,基本上2G的数据量就够组装得到完整基因组用的了。所以直接用head命令取全基因组测序数据的前2千万行就够了

head -n 20000000 data_R1.fastq > R1.fastq
head -n 20000000 data_R2.fastq > R2.fastq

最后是组装

get_organelle_from_reads.py -1 R1.fastq -2 R2.fq -o plastome_output -R 15 -k 105,121 -F embplant_pt

如果组装成功最后会得到两条序列,这两条序列差别在小单拷贝区方向不同,我自己的处理方式是选一些近缘的参考叶绿体基因组来构建进化树,去掉枝长明显过长的那一个。

到此,相信大家对“如何使用GetOrganelle软件组装叶绿体基因组”有了更深的了解,不妨来实际操作一番吧!这里是亿速云网站,更多相关内容可以进入相关频道进行查询,关注我们,继续学习!

向AI问一下细节

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

AI