Question: define PFAM dir in get_homologues perl script
1
gravatar for dago
4.1 years ago by
dago2.5k
Germany
dago2.5k wrote:

I think this is really a thrivial question but I am stuck on it.

I would like to use get_homolohues to perform a comaprative genomic analysis. It requires some dependencies, such as MCL, Blast, hmmer. I am using Biolinux 8 and most of them are installed already.

I follow the installation instructions. When checking the binaries needed I got the following:

 perl get_homologues.pl -v

Checking required binaries and data sources, all set in phyTools.pm :
        EXE_BLASTP : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/bin/ncbi-blast-2.2.27+/bin/blastp)
        EXE_BLASTN : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/bin/ncbi-blast-2.2.27+/bin/blastn)
      EXE_FORMATDB : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/bin/ncbi-blast-2.2.27+/bin/makeblastdb)
           EXE_MCL : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/mcl-02-063/shmcl/mcl)
      EXE_MAKEHASH : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/COGsoft/COGmakehash/COGmakehash )
     EXE_READBLAST : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/COGsoft/COGreadblast/COGreadblast )
        EXE_COGLSE : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/COGsoft/COGlse/COGlse )
        EXE_COGTRI : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/COGsoft/COGtriangles/COGtriangles )
       EXE_HMMPFAM : Pfam-A.hmm file needs to be installed (ftp://ftp.sanger.ac.uk/pub/databases/Pfam/)

        EXE_INPARA : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/_cluster_makeInparalog.pl)
         EXE_ORTHO : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/_cluster_makeOrtholog.pl)
         EXE_HOMOL : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/_cluster_makeHomolog.pl)
    EXE_SPLITBLAST : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/_split_blast.pl)
  EXE_SPLITHMMPFAM : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/_split_hmmscan.pl)

So I just needed to install PFAM. I dowloaded as follow: 

wget -N ftp://ftp.sanger.ac.uk/pub/databases/Pfam/current_release/Pfam-A.hmm.gz

gunzip Pfam-A.hmm.gz

 

I then changed the directory in the appropriate perl scipt to my actual dir where I have the Pfam-A.hmm, which is

/usr/local/lib/Pfam-A.hmm

However, when I check again the binaries needed I still get the same message as above, telling me that Pfam-A is not installed.

Anyone can help me with this?

uname -a                                                                                                                                           [ 3:20PM]
Linux bl8vbox 3.13.0-39-generic #66-Ubuntu SMP Tue Oct 28 13:30:27 UTC 2014 x86_64 x86_64 x86_64 GNU/Linux

Update

I downloaded again the Pfam-A.hmm.gz and this time it worked.

I run the  get_homologues.pl with the option -D for Pfam search and I do not get any error. However, when I open the file resulted from the analysis XXX.fasta.pfam I do not have any data, but this message:

Error: Unrecognized format, trying to open hmm file /usr/local/lib/get_homologues-x86_64-20140930/db/Pfam-A.hmm for reading.

Any idea?

 

hmmer pfam biolinux perl • 2.3k views
ADD COMMENTlink modified 4.1 years ago by b.contreras.moreira10 • written 4.1 years ago by dago2.5k

Did you follow all the instructions given here? Did you install HMMER3? Did you run the commands given below:

Generate binary files for Pfam-A.hmm using 

/your/path/hmmer-3.0/src/hmmpress Pfam-A.hmm

Edit file lib/phyTools.pm and modify the relevant line to:

if( ! defined($ENV{"PFAMDB"}) ){ $ENV{"PFAMDB"} = "db/Pfam-A.hmm"; }
ADD REPLYlink modified 4.1 years ago • written 4.1 years ago by komal.rathi3.3k

Thanks a lot.

Hmmer Is already installed in Biolinux. I followed the instruction as in the manual. I

ADD REPLYlink written 4.1 years ago by dago2.5k
1
gravatar for b.contreras.moreira
4.1 years ago by
Spain
b.contreras.moreira10 wrote:

Hi Dago,

get_homologues is bundled with all the required sources/binaries, so after running install you should be ready to go if all you get are OKs. The install.pl script will also take care of pressing the PfamA.hmm file located at the $ENV{"PFAMDB"} location, as mentioned above. I do recommend proceeding this way, as thus the software is using compatible hmmpress and then hmmscan binaries. If you choose to press the hmms yourself then you should check that the hmmpress binary in biolinux is compatible with the bundled hmmscan binary,

hope this helps,

Bruno

ADD COMMENTlink written 4.1 years ago by b.contreras.moreira10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1218 users visited in the last hour