转录组cufflinks套装的使用

cufflinks套装有很多,我们主要使用的只有三个

Cufflinks是用来处理tophat的输出的bam文件然后输出gtf文件

cuffmerge把多个样本的gtf文件合并的,也没啥子用,主要是测多个样本可能会需要

cuffdiff算出分组的bam文件里面的差异基因。

一:下载安装该软件

是二进制版本,找到网址,然后用wget下载,解压即可使用

转录组cufflinks套装的使用-图片1

二:准备数据

cufflinks和cuffmerge和cuffdiff需要的文件各不相同

对cufflinks来说,需要的是tophat输出的bam文件

对cuffmerge来说,需要的是多次运行cufflinks输出的多个gtf文件

对cuffdiff来说,需要的多个样本的tophat的bam文件

三:运行命令

cufflinks对人类来说几乎没有用,因为对人类研究的太完全了,一大把的gtf文件,随便都可以去下载到,不需要从tophat的bam文件里面重新输出gtf文件,我简单给一个脚本

转录组cufflinks套装的使用-图片2

对cuffmerge来说,人类也是不需要的,主要是新物种还可以用一下,不需要脚本,就是简单的一句话而已,我都懒得讲了。

重点:对于cuffdiff来说,这个非常重要咯,我给出好几个版本的命令,希望你能总结出规律

第一个例子

转录组cufflinks套装的使用-图片3

有些是必须参数,跟你项目调整即可

第二个例子

cuffdiff

-o FN  #这是设置结果输出目录

-b /home/immune/refer_genome/hg19/hg19 #这是参考基因组文件

-p 18 #这里cpu数量的控制

-L case,control #这是标签

-u merged_gtf/merged.gtf  #这是参考gtf文件

./FN_A549_1/accepted_hits.bam,./LA_A549_2/accepted_hits.bam,./FN_A549_3/accepted_hits.bam    #这是case组的bam文件地址,共三个bam文件

./FN_con_A549_1/accepted_hits.bam,./FN_con_549_2/accepted_hits.bam,./FN_con_A549_3/accepted_hits.bam #这是control组的bam文件地址,也是三个bam文件

第三个例子

转录组cufflinks套装的使用-图片4

对了,参考基因组还需要有个fai文件,但是cuffdiff本身也会创造,没有也无所谓的

转录组cufflinks套装的使用-图片5

结束之后的log如下显示

转录组cufflinks套装的使用-图片6

4,输出数据解释,这里主要解释cuffdiff的结果文件

转录组cufflinks套装的使用-图片7

输出的文件太多了,我也没怎么具体细看。

反正那个gene_exp.diff里面可以看到差异基因的信息

主要研究了gene,isoform,cds,tss这四个方面的差异,也是我们可能会关心的四个方面的信息,同时包含着他们是如何具体算出了的tracking文件。

原文来自:http://www.bio-info-trainee.com/166.html

发表评论

匿名网友

拖动滑块以完成验证