tRNAscan-SE安装使用教程

TIPs:

最近帮一个老师用SPORTS软件分析miRNA(也可以分析其他小RNA)序列的结构特征,可能会用到以下数据库,Few main database source:
mirbase database
rRNA database
GtRNAdb database
mitotRNAdb database
piRBase database
piRNABank
ensembl ncRNA database
rfam database
其中有些数据库中没有我需要的物种数据信息,需要借助其他软件分析提取,例如GtRNAdb database中没有鸽子tRNA信息,因此就需要下面的方法来预测tRNA出来。
推荐两个博客资源:1)https://blog.csdn.net/herokoking/article/details/77853662;2)https://blog.csdn.net/weixin_43569478/article/details/108079199

软件安装和测试

$ wget -c http://eddylab.org/infernal/infernal-1.1.2-linux-intel-gcc.tar.gz
$ tar xzvf infernal-1.1.2-linux-intel-gcc.tar.gz
$ cd infernal-1.1.2-linux-intel-gcc/
$ ./configure --prefix=/user/software/infernal-1.1.2-linux-intel-gcc
$ make && make install
$ cp binaries/*  bin

$ wget -c http://lowelab.ucsc.edu/software/tRNAscan-SE.tar.gz
$ tar -xvf tRNAscan-SE.tar.gz
$ mkdir tRNAscan-SE-2.0_installed # 跳过此操作安装在解压目录下也可以
$ ./configure --prefix=/user/software/tRNAscan-SE-2.0_installed
$ make && make install
# 添加bashrc【建议】或者terminal运行
    export PATH=/user/software/tRNAscan-SE-2.0_installed/bin:$PATH
    export PERL5LIB=/user/software/tRNAscan-SE-2.0_installed/lib:$PERL5LIB

# 需要编辑/user/software/tRNAscan-SE-2.0_installed/bin/tRNAscan-SE.conf文件,每个项目单独拷贝到相应目录下编辑保存
    bin_dir:/user/software/tRNAscan-SE-2.0_installed/bin
    lib_dir: /user/software/tRNAscan-SE-2.0_installed/lib
    infernal_dir: /user/software/infernal-1.1.2-linux-intel-gcc/bin

# 查看帮助说明和具体参数
$ tRNAscan-SE --help
tRNAscan-SE /user/database/GeZi_ncbi_1.0/final.fa -o final.out -f final.ss -s final.iso  -m final.stats -b final.bed -c /user/software/tRNAscan-SE-2.0_installed/bin/tRNAscan-SE.conf
# 该命令产生的结果与Demo文件夹中的结果文件类型一致,其中参数
  -o 结果文件
  -f 二级结构(tRNA)文件
  -s (Isotype specific)结果
  -m 统计文件(statistics)
  -b 输出bed格式结果,也可不设置

结果文件展示

final.out

Sequence Name: 输入的用于预测tRNA的序列名称
tRNA: 预测到的tRNA的个数(第n个)
tRNA Begin: tRNA基因的起始位置
Bounds End: tRNA 基因的终止位置
tRNA Type: tRNA 的类型,该tRNA转运的氨基酸的类型
Anticodon: 反义密码子
Intro Begin: 内含子的起始位置
Bounds End: 内含子的终止位置
Inf Score: cutoff score for reporting tRNAs (default=20)这里输出阈值以上的预测结果


final.out

final.ss

Str 代表的就是二级结构,其中每个< 和每个 > 是配对的,代表在二级结构中这两个碱基是连在一起的。
final.ss

final.iso

tRNAscan-SE 会根据反义密码子和 Isotype Model 两种方法预测tRNA的类型,对于每个tRNA, 分别给出反义密码子和 Isotype Model 预测的结果,以及两种结果预测的一致性,其中isotype model 会对每一种氨基酸进行打分,选得分最高的作为最终的结果
final.iso

final.stats

运行环境、参数等的设置在这里,另外还统计了基因组处理的过程、指标的统计。
final.stats

final.stats

bed输出我没跑,放上demo中的bed,估计我们文件产生的格式与之一致,至于bed文件现在看来作为常识了,自己百度吧。
bed文件

总结

  1. http://lowelab.ucsc.edu/tRNAscan-SE/index.html 可以网页版运行,很方便,但是注意网页中内置的参数等问题:
    网页运行
  2. 注意tRNA鉴定(search)的模式,或者如果不知道某些参数(score)的涵义等等,可以grep一下。
    search
    score
  3. 最后,聊以记录我学习该软件的过程,欢迎交流学习!
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 218,122评论 6 505
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,070评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,491评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,636评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,676评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,541评论 1 305
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,292评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,211评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,655评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,846评论 3 336
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,965评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,684评论 5 347
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,295评论 3 329
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,894评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,012评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,126评论 3 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,914评论 2 355

推荐阅读更多精彩内容