OrthoFinder

GitHub上对此描述的也是很详细: https://github.com/davidemms/OrthoFinder

一、下载安装

1.conda安装

conda install -c bioconda orthofinder

(本人第一次利用conda安装成功了,第二次可能是因为镜源连接不稳定等原因,安装不成功,第二次尝试了源码下载安装;也可能存在orthofinder与conda中的Python的版本不对应的情况,所以,查看orthofinder的安装环境)

2.源码下载安装

参考https://davidemms.github.io/orthofinder_tutorials/alternative-ways-of-getting-OrthoFinder.html
在ReadMe文件中,存在安装的方式,是依赖于其他软件的(MCL, FastME, DIAMOND/MMseqs2/BLAST+)(但可能在github下载的时候已经将这些都包含里面了,不需要再自行下载这些软件)。
*还存在一些问题,在使用的过程中,不知道为什么需要添加 ./ *

wget https://github.com/davidemms/OrthoFinder/releases/download/2.5.4/OrthoFinder.tar.gz
tar xzf OrthoFinder.tar.gz
./OrthoFinder/orthofinder -h ##运行帮助文档,查看是否安装成功
OrthoFinder/orthofinder -f OrthoFinder/ExampleData ###用里面的例子数据运行尝试

若是需要下载依赖的文件,参考//www.greatytc.com/p/16e0bbb2ba19

二、基础使用

1.准备数据

准备没有isoform的蛋白数据(没有重复)(此处选择的是最长转录本数据),将数据放在一个新建的目录下data_pep,数据文件夹下至少需要两个文件:【若是存在query,则不要忘记添加query文件

2.基础运行

基础的快速的运行:

orthofinder -f data_pep##可以添加上-t 200
###Run full OrthoFinder analysis on FASTA format proteomes in <dir>
 ## orthofinder [options] -f <dir>

简单的参数选择设置:

 -t <int>        Number of parallel sequence search threads [Default = 104]
 -a <int>        Number of parallel analysis threads
 -d              Input is DNA sequences
 -M <txt>        Method for gene tree inference. Options 'dendroblast' & 'msa'
                 [Default = dendroblast]
 -S <txt>        Sequence search program [Default = diamond]
                 Options: blast, diamond, diamond_ultra_sens, blast_gz, mmseqs, blast_nucl
 -A <txt>        MSA program, requires '-M msa' [Default = mafft]
                 Options: mafft, muscle
 -T <txt>        Tree inference method, requires '-M msa' [Default = fasttree]
                 Options: fasttree, raxml, raxml-ng, iqtree

通常采用的运行命令:

orthofinder -f data_pep -S diamond -M msa -T raxml -t 200 -a 200

三、报错

第一次使用出现了报错(系统对于用户打开文件数量的限制):

ERROR: The system limits on the number of files a process can open is too low. 
For 240 species OrthoFinder needs to be able to open at least r=57700 files. 
Please increase the limit and restart OrthoFinder

查看用户对于打开文件的限制;若需要改动则需要root对用户进行设置:增加系统对该进程打开的文件

ulimit -Hn  #打开文件数量的硬限制
ulimit -Sn  #打开文件数量的软限制
ulimit -n  #允许当前用户进程打开的文件数限制
cat /proc/sys/fs/file-max  #查看能够打开的文件的最大限制

正常情况下,用户可以设置Sn的最大值为Hn,ulimit -Sn 4096
若需要继续增加数值,则root进一步设置文件能够打开的最大的限制:

vim /etc/security/limits.conf  # 在文件的最后添加一行
your_user_name hard nofile 57700
your_user_name soft nofile 57700

然后在编辑/etc/pam.d/login文件,在文件后边添加如下行:

session required /lib/security/pam_limits.so

四、结果分析

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

推荐阅读更多精彩内容