13518219792
这篇文章主要介绍“自动化运行OrthoMCL的方法”，在日常操作中，相信很多人在自动化运行OrthoMCL的方法问题上存在疑惑，小编查阅了各式资料，整理出简单好用的操作方法，希望对大家解答”自动化运行OrthoMCL的方法”的疑惑有所帮助！接下来，请跟着小编一起来学习吧！
成都创新互联是一家集网站建设,海拉尔企业网站建设,海拉尔品牌网站建设,网站定制,海拉尔网站建设报价,网络营销,网络优化,海拉尔网站推广为一体的创新建站企业，帮助传统企业提升企业形象加强企业竞争力。可充分满足这一群体相比中小企业更为丰富、高端、多元的互联网需求。同时我们时刻保持专业、时尚、前沿，时刻以成就客户成长自我，坚持不断学习、思考、沉淀、净化自己，让我们为更多的企业打造出实用型网站。
安装OrthoMCL Pipeline
直接去github下载zip安装包，或者直接使用git命令下载
git clone https://github.com/apetkau/orthomcl-pipeline.git
安装perl依赖包
使用cpan或者cpanm安装依赖包
$ cpanm BioPerl DBD::MySQL DBI Parallel::ForkManager YAML::Tiny Set::Scalar Text::Table Exception::Class Test::Most Test::Warn Test::Exception Test::Deep Moose SVG Algorithm::Combinatorics
安装blast，MCL及OrthoMCL
安装好2.2.26版本的BLAST (blastall, formatdb, 不是blast+, 更旧的版本可能不会正确运行)
下载编译安装MCL
下载编译安装OrthoMCL
确保blastall, formatdb, MCl和OrthoMCL都加入了环境变量
检查依赖的软件是否安装齐全
出现下面的提示则安装齐全了
$ perl scripts/orthomcl-pipeline-setup.pl
Checking for Software dependencies...
Checking for OthoMCL ... OK
Checking for formatdb ... OK
Checking for blastall ... OK
Checking for mcl ... OK
Wrote new configuration to orthomcl-pipeline/scripts/../etc/orthomcl-pipeline.conf
Wrote executable file to orthomcl-pipeline/scripts/../bin/orthomcl-pipeline
Please add directory orthomcl-pipeline/scripts/../bin to PATHOrthomcl-pipeline配置文件
配置文件在orthomcl-pipeline-installed-path/etc/orthomcl-pipeline.conf
blast: F: 'm S' b: 100000 e: 1e-5 v: 100000 filter: max_percent_stop: 20 min_length: 10 mcl: inflation: 1.5 path: blastall: '/usr/bin/blastall' formatdb: '/usr/bin/formatdb' mcl: '/usr/local/bin/mcl' orthomcl: '/home/aaron/software/orthomcl/bin' scheduler: fork split: 4
根据计算资源可以把split改为合适的数字，split类似于运行的线程数
配置数据库
需要安装好MySQL数据库
1. 如果之前已经创建过orthomcl数据库，可以用下面的perl脚本进行配置
$ perl scripts/orthomcl-setup-database.pl --user orthomcl --password orthomcl --host localhost --database orthomcl --outfile orthomcl.conf --no-create-database Connecting to database orthmcl on host orthodb with user orthomcl ...OK Config file **orthomcl.conf** created.
2. 如果之前没有配置过orthomcl数据库，需要先以root运行mysql，创建orthomcl用户
$ mysql -u root -p
Enter password:
mysql> GRANT SELECT, INSERT, UPDATE, DELETE, CREATE, CREATE VIEW, INDEX, DROP on *.* to orthomcl; #创建用户并授权
mysql> set password for orthomcl@localhost = password('orthomcl'); #设置用户密码
mysql> quit;当用户创建后，再用perl脚本配置
$ perl scripts/orthomcl-setup-database.pl --user orthomcl --password orthomcl --host localhost --database orthomcl --outfile orthomcl.conf Connecting to mysql and creating database **orthmcldb** on host orthodb with user orthomcl ...OK database orthmcl created ...OK Config file **orthomcl.conf** created.
脚本会生成orthomcl.conf文件，这个文件之后就不用坐任何修改了。
测试
$ perl t/test_pipeline.pl -m orthomcl.conf -s fork -t /tmp Test using scheduler fork TESTING NON-COMPLIANT INPUT TESTING FULL PIPELINE RUN 3 README: Tests case of one gene (in 1.fasta and 2.fasta) not present in other files. ok 1 - Expected matched returned groups file ...
出现上面的信息说明安装成功了！！！
运行
$ ./bin/orthomcl-pipeline Error: no input-dir definedUsage: orthomcl-pipeline -i [input dir] -o [output dir] -m [orthmcl config] [Options] ...
可以加上—nocompliant就可以不显示没必要的提示信息了
$ ./bin/orthomcl-pipeline -i ~/test_data/zbl -o ~/test_data/zbl-out -m orthomcl.conf --nocompliant
zbl文件夹里面放入fasta格式的蛋白序列文件，一个基因组一个文件，例如genome1.fasta, genome2.fasta, 文件后缀名必需是.fasta。zbl-out文件夹里面是所有的结果，包括聚类完成的groups文件，orthologs文件，inparalogs文件以及coorthologs文件。
到此，关于“自动化运行OrthoMCL的方法”的学习就结束了，希望能够解决大家的疑惑。理论与实践的搭配能更好的帮助大家学习，快去试试吧！若想继续学习更多相关知识，请继续关注创新互联网站，小编会继续努力为大家带来更多实用的文章！
