前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >circRNA分析利器KNIFE试用笔记

circRNA分析利器KNIFE试用笔记

作者头像
生信技能树
发布2018-09-21 16:28:10
1.4K0
发布2018-09-21 16:28:10
举报
文章被收录于专栏:生信技能树

根据circRNA分析软件benchmark文献推荐,KNIFE准确度最高,circRNA分析优先推荐使用KNIFE工具。

  • 参考文献:A comprehensive overview and evaluation of circular RNA detection tools

1. 下载修新版KNIFE软件

从release中选择v1.4稳定版本下载,不建议使用git clone下载,开发版还有一些bug未解决。

代码语言:javascript
复制
cd ~/biosoft/
wget https://github.com/lindaszabo/KNIFE/archive/v1.4.tar.gz
tar zxvf v1.4.tar.gz && cd KNIFE-v1.4
cd circularRNApipeline_Standalone/analysis
# 添加执行权限,否则运行过程中会由于部分软件无执行权限报错。
sudo chmod a+x *

2. 确保所有R、python、和shell脚本可以执行

如果无法运行./completeRun.sh,需要将sh默认执行程序设置为非dash

代码语言:javascript
复制
sudo dpkg-reconfigure dash

3. 安装下面软件,并添加至环境路径

  • Bowtie2 2.2.1+ (tested with 2.2.1)
  • Bowtie 0.12.7+ (tested with 0.12.7)
  • perl
  • python 2.7.3+ (tested with 2.7.5) and the numpy and scipy libraries
  • R 3.0.2+ (tested with 3.0.2) and the data.table package
  • samtools 1.2+ if you want alignment files output as bam instead of sam files

4. 下载bowtie索引文件

a. 创建索引文件夹

创建以下两个文件夹,分别用于存放bowtie和bowtie2索引文件。 circularRNA/pipeline_index和 circularRNA/pipeline_denovo_scripts/index

b. 下载bowtie索引文件并解压到相应文件夹

Human (hg19), Mouse (mm10), Rat (rn5) and Drosophila (dm3)有可用的junction索引文件下载。

作者同时打包了所有的转录组、基因组和核糖体RNA索引、fasta文件和gtf文件。

需要将bowtie1索引文件解压至以下文件夹:

代码语言:javascript
复制
cd /path/to/circularRNApipeline_*/denovo_scripts/index
tar zxvf genomeId_BT1index.tar.gz

将bowtie2索引文件解压至以下文件夹:

代码语言:javascript
复制
cd /path/to/circularRNApipeline_*/index
tar zxvf genomeId_BT2index.tar.gz

gtf文件需要下载并解压缩至以下文件夹:

代码语言:javascript
复制
cd /path/to/circularRNApipeline_*/denovo_scripts
gunzip genomeId_genes.gtf.gz

KNIFE提供的hg19参考基因组和索引文件可以通过下面链接下载: https://mega.nz/#F!RtsCHCQb!fyxYNWjoCef5Ie361vUxiA

如果内存够大,也可以用自己的机器构建索引。

5. 按照testOutput/README指导进行测试

(1)下载release1.1中提供的两个fastq.gz文件
代码语言:javascript
复制
cd /path/to/KNIFE/testData
nohup wget -c https://github.com/lindaszabo/KNIFE/releases/download/v1.1/SRR1027187_1.fq.gz &
nohup wget -c https://github.com/lindaszabo/KNIFE/releases/download/v1.1/SRR1027187_2.fq.gz &
(2) 测试KNIFE是否安装好(单机运算)
代码语言:javascript
复制
sh completeRun.sh /home/jshi/biosoft/KNIFE/testData/ complete /home/jshi/biosoft/KNIFE/new testData 8 phred64 circReads 40 1 2>&1

可以将输出文件和github中给出的结果进行比较,看是否完成了环状RNA分析。

值得注意的是,全部分析流程都被包装到了一个shell脚本里。

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2018-08-26,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 生信技能树 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 1. 下载修新版KNIFE软件
  • 2. 确保所有R、python、和shell脚本可以执行
  • 3. 安装下面软件,并添加至环境路径
  • 4. 下载bowtie索引文件
    • a. 创建索引文件夹
      • b. 下载bowtie索引文件并解压到相应文件夹
      • 5. 按照testOutput/README指导进行测试
        • (1)下载release1.1中提供的两个fastq.gz文件
          • (2) 测试KNIFE是否安装好(单机运算)
          领券
          问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档