Where can I download HMM profile for different species?
0
0
Entering edit mode
8.4 years ago

I would like to ask is there any database allow me to download HMM for different species such as: 1. Pan troglodytes 2. Pongo abelii 3. Rhinopithecus rexellana

hmm profile species genome sequence • 3.3k views
ADD COMMENT
1
Entering edit mode

It seems you can download this from HMMER. If so, get more information from their paper HMMER web server: 2015 update.

ADD REPLY
0
Entering edit mode

I can't find any hmm link to download the above 3 species :/

ADD REPLY
1
Entering edit mode

Which kind of HMM profile are you looking for ? => I mean is it for coding genes ? For what purpose do you need it ? More precisely, which tool do you plan to use with ?

ADD REPLY
0
Entering edit mode

i would like to include the hmm in snap ab initio gene predictor since the package doesn't provide the other species

ADD REPLY
1
Entering edit mode

Each ab-initio tool has its own hmm profile format. And they don't necessarily use the same kind of information. So, I doubt you can use the profile coming from another tool as HMMER as suggested. I guess you are trying to annotate the gene coding. So the simplest is to create the profile yourself. Download the gff and the genome of your species (on Ensembl as example). Put the fasta sequence within the gff file as it should be (Introduced by the ##FASTA line at the end of the file). Then prepare the needed files to train snap by using the maker2zff.pl tool that is coming with the MAKER genome annotation tool package. You will get two output files:genome.ann and genome.dna Then to create the profile you will have to launch those commands:

fathom -categorize 1000 genome.ann genome.dna
    fathom -export 1000 -plus uni.ann uni.dna
    forge export.ann export.dna
    hmm-assembler.pl SPECIES . > SPECIES.hmm
ADD REPLY
0
Entering edit mode

It sounds like a reasonable approach!

ADD REPLY
0
Entering edit mode

Which HMM do you need?

ADD REPLY
0
Entering edit mode

the 3 species in my question

ADD REPLY
0
Entering edit mode

Yes I understand, but do you need the HMM for the nucleotides seqs of all genes, for the proteins? Also, I am not sure you have hmm specific for these species, maybe you can try to use Pfam domains and start with that for gene annotation.

ADD REPLY

Login before adding your answer.

Traffic: 1885 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6