Inparanoid 寻找物种间直系同源基因

InParanoid 是一个用于寻找直系同源基因(Orthologs)的计算机程序,其目标是发现不同物种之间具有相似功能的基因。

软件下载安装:

git clone https://bitbucket.org/sonnhammergroup/inparanoid.git

或者

wget   https://bitbucket.org/sonnhammergroup/inparanoid/get/d0690ac396ac.zip

在运行的时候,系统有报错,报错内容为:

Can't locate Parallel/ForkManager.pm in @INC (you may need to install the Parallel::ForkManager module) (@INC contains: /build/Cellar/anaconda2/lib/site_perl/5.26.2/x86_64-linux-thread-multi /build/Cellar/anaconda2/lib/site_perl/5.26.2 /build/Cellar/anaconda2/lib/5.26.2/x86_64-linux-thread-multi /build/Cellar/anaconda2/lib/5.26.2 .) at inparanoid.pl line 116.BEGIN failed--compilation aborted at inparanoid.pl line 116.

缺少 Parallel::ForkManager 模块。

下载安装 perl 的 Parallel::ForkManager 模块:

wget   https://cpan.metacpan.org/authors/id/Y/YA/YANICK/Parallel-ForkManager-2.02.tar.gz

tar    zxvf   Parallel-ForkManager-2.02.tar.gz

perl   MANIFEST   INSTALL_BASE=/home/xx/

make

make   install

将路径添加到  ~/.bashrc  中的环境变量:

export   PERL5LIB=/home/xx/lib/perl5:$PERL5LIB

就 ok 了。

环境依赖

软件的运行需要提前安装好 DIAMOND 和 BLAST。

## 安装 DIAMOND 

wget http://github.com/bbuchfink/diamond/releases/download/v2.0.8/diamond-linux64.tar.gz

tar xzf diamond-linux64.tar.gz

## 高版本的软件需要 GLIBC_2.17,系统上没有该版本,低版本的软件没有这个问题。

## 安装 BLAST

wget ftp://ftp.ncbi.nlm.nih.gov/blast/executables/legacy.NOTSUPPORTED/2.2.18/blast-2.2.18-x64-linux.tar.gz 

tar xvfz blast-2.2.18-x64-linux.tar.gz

软件运行

## 运行示例文件,全部使用默认参数

perl   inparanoid.pl   -input-dir   ./testInput/

可选参数:

- f1                物种 A 的蛋白质序列的 FASTA 文件。

- f2                物种 B 的蛋白质序列的 FASTA 文件。

- outgroup     作为外类群的物种 C 的蛋白质序列的 FASTA 文件(可选,默认无).

- input-dir        包含多个物种蛋白质序列的 FASTA 文件的目录。将运行所有物种之间的比对。如果使用了此选项,则将 -f1 和 -f2 留空。

- out-dir        指定输出文件的目录.

- seq-tool        用于序列相似性比对的工具。可选项为 Diamond 或 Blast(默认为 Diamond)。

- 2pass        运行两遍比对策略。不适用于 Diamond,但推荐用于 Blast(默认为 False)。

- bootstrap        运行引导法以估计正交基因的置信度(默认为 False)。

- seedscore        计算 seedscore 以估计正交基因的置信度(默认为 False)。

- score-cutoff        设置位分数阈值。任何低于此分数的匹配将被忽略(默认为 40)。

- seq-cutoff        设置序列重叠阈值。匹配区域应至少覆盖较长序列的此部分。匹配区域是从第一个片段的开始到最后一个片段的结束(默认为 0.5)。

- seg-cutoff        设置片段覆盖率阈值。匹配片段必须覆盖较长序列的此部分(默认为 0.25)。

- conf-cutoff        置信度阈值。在这个置信度或更高的情况下,将包括内源同源基因。默认值为 0.05。

- grp-cutoff        组重叠阈值。如果一个组中的正交基因在另一个组中的置信度超过此值,那么这两个组将被合并。默认值为 0.5。

- sensitivity    Diamond 的敏感度设置。Diamond 提供了多种敏感度模式,可以根据具体需求选择。可选项包括 mid-sensitive(中等敏感度)、sensitive(敏感度)、more-sensitive(更高敏感度)、very-sensitive(非常敏感度)和 ultra-sensitive(超敏感度)。默认值为 very-sensitive。

结果文件

SQLtable 文件以制表符分割,每行是一个蛋白序列,每一列内容分别是:

Group-id      Score    Species   Confidence-score    Protein-name

序列分析原始文件:

ProteinA      ProteinB     Score     LengthA     LengthB      Length-MatchingRegionA      Length-MatchingRegionB      Total-Length-of-MatchA     Total-Length-of-MatchB     Location-of-Matches

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,222评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,455评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 157,720评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,568评论 1 284
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,696评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,879评论 1 290
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,028评论 3 409
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,773评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,220评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,550评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,697评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,360评论 4 332
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,002评论 3 315
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,782评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,010评论 1 266
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,433评论 2 360
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,587评论 2 350

推荐阅读更多精彩内容