CircRNA 预测软件之CIRCexplorer2使用流程（三）

今天的文章中我给大家介绍环状RNA预测的另一款重要的软件——CIRCexplorer2. 该软件是由中国科学院上海生命科学研究院生物化学与细胞生物学研究所陈玲玲研究组和计算生物学研究所杨力研究组共同开发的， 他们团队已经在《Cell》, 《Mol Cell》和《Genome Res》等国际顶尖杂志发表多篇环状RNA研究文章。

CIRCexplorer2继承了CIRCexplorer主要功能并增加了很多新的特性。它支持TopHat2/TopHat-Fusion, STAR, MapSplice, BWA and segemehl等多种RNA aligners, 并且能够精确地注释预测到环状RNA。更为重要的是，该工具可以侦测到多种环状RNA可变剪切事件，并且能够de nove组装环状RNA全长转录本。图1描绘了CIRCexplorer2的工作流程：

图1：The schematic diagram of CIRCexplorer2

下面给大家介绍CIRCexplorer2的工作流程和命令行参数

1、  安装

CIRCexplorer2是基于Python编写的，因此可以用pip install工具直接安装在本地服务器上。代码如下：

git clone https://github.com/YangLab/CIRCexplorer2.git

cd CIRCexplorer2

pip install -r requirements.txt

### install scipy according to http://www.scipy.org/install.html

python setup.py install

在CIRCexplorer2 documentation中还提供了几种可选择的安装方式，读者有兴趣的话可以参考：http://circexplorer2.readthedocs.io/en/latest/

Tophat-Fusion pair-end alignment

尽管CIRCexplorer2有自己的基因组比对模块（CIRCexplorer2 align），并且支持TopHat2/TopHat-Fusion，STAR, MapSplice等多种RNA aligner，但对于双端测序数据（pair-end sequencing reads）来说，只支持TopHat-Fusion。因此，我们只介绍基于TopHat-Fusion的比对数据。

$tophat --library-type fr-unstranded -o tophat_fusion_out -p 15 --fusion-search --keep-fasta-order --bowtie1 --no-coverage-search ${bowtie_index} read1.fq read2.fq

###TopHat参数介绍###

--library-type fr-unstranded ： 标准Illumina文库类型。

-o ：结果输出目录，包含accepted_hits.bam等主要输出文件

-p ：多线程

--fusion-search ：允许fusion mapping

--keep-fasta-order ： 保持与基因组fasta文件相同的顺序

--bowtie1 ： 采用bowtie1算法进行比对

--no-coverage-search ：不允许基于coverage的算法寻找junctions

2、  Parsing for Circular RNA Fusion Junction Reads

$CIRCexplorer2 parse --pe -t TopHat-Fusion tophat_fusion/accepted_hits.bam –o circ_out 2> CIRCexplorer2_parse.log

### CIRCexplorer2 parse参数介绍###

--pe ：Parse paired-end alignment file (only for TopHat-Fusion).

-t ：指定aligner

3、  Annotating for Circular RNAs

$ CIRCexplorer2 annotate -r hg19_ref_all.txt -g hg19.fa circ_out > CIRCexplorer2_annotate.log

### CIRCexplorer2 annotate参数介绍###

-r ：基因注释文件

-g ： 基因组fasta文件

经过以上4个步奏，CIRCexplorer2会生成一个annotate目录，里面包括circ_fusion.txt文件，即是环状RNA注释文件，其内容包括以下几列：