1. KNIFE
根據(jù)circRNA分析軟件benchmark文獻(xiàn)推薦菱属，KNIFE準(zhǔn)確度最高召噩，circRNA分析優(yōu)先推薦使用KNIFE工具。
1. 下載修新版KNIFE軟件
從release中選擇v1.4穩(wěn)定版本下載栅螟，不建議使用git clone下載执俩，開發(fā)版還有一些bug未解決徐钠。
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
# 添加執(zhí)行權(quán)限，否則運(yùn)行過程中會由于部分軟件無執(zhí)行權(quán)限報(bào)錯(cuò)役首。
sudo chmod a+x *
2. 確保所有R尝丐、python、和shell腳本可以執(zhí)行
如果無法運(yùn)行./completeRun.sh衡奥，需要將sh默認(rèn)執(zhí)行程序設(shè)置為非dash
sudo dpkg-reconfigure dash
3. 安裝下面軟件爹袁，并添加至環(huán)境路徑
- 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. 創(chuàng)建索引文件夾
創(chuàng)建以下兩個(gè)文件夾，分別用于存放bowtie和bowtie2索引文件矮固。
circularRNApipeline_***/index和 circularRNApipeline_***/denovo_scripts/index
b. 下載bowtie索引文件并解壓到相應(yīng)文件夾
Human (hg19), Mouse (mm10), Rat (rn5) and Drosophila (dm3)有可用的junction索引文件下載失息。
作者同時(shí)打包了所有的轉(zhuǎn)錄組、基因組和核糖體RNA索引、fasta文件和gtf文件盹兢。
需要將bowtie1索引文件解壓至以下文件夾：
cd /path/to/circularRNApipeline_*/denovo_scripts/index
tar zxvf genomeId_BT1index.tar.gz
將bowtie2索引文件解壓至以下文件夾：
cd /path/to/circularRNApipeline_*/index
tar zxvf genomeId_BT2index.tar.gz
gtf文件需要下載并解壓縮至以下文件夾：
cd /path/to/circularRNApipeline_*/denovo_scripts
gunzip genomeId_genes.gtf.gz
KNIFE提供的hg19參考基因組和索引文件可以通過下面鏈接下載： https://mega.nz/#F!RtsCHCQb!fyxYNWjoCef5Ie361vUxiA
如果內(nèi)存夠大邻梆，也可以用機(jī)器構(gòu)建索引。
5. 按照testOutput/README指導(dǎo)進(jìn)行測試
（1）下載release1.1中提供的兩個(gè)fastq.gz文件
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是否安裝好(單機(jī)運(yùn)算)
sh completeRun.sh /home/jshi/biosoft/KNIFE/testData/ complete /home/jshi/biosoft/KNIFE/new testData 8 phred64 circReads 40 1 2>&1
可以將輸出文件和github中給出的結(jié)果進(jìn)行比較蛤迎，看是否完成了環(huán)狀RNA分析确虱。
