Pittsburgh Supercomputing Center 

Advancing the state-of-the-art in high-performance computing,
communications and data analytics.

SNAP

 

 

SNAP is a general purpose gene finding program suitable for both eukaryotic and prokaryotic genomes. SNAP is an acroynm for Semi-HMM-based Nucleic Acid Parser.

Installed on blacklight and biou.

Other resources that may be helpful include:

Korf I. Gene finding in novel Genomes. BMC Bioinformatics 2004, 5:59
Website: http://korflab.ucdavis.edu/software.html 

Running SNAP

1) Make SNAP commands availiable for use
a) blacklight:
The SNAP programs will be made availiable for use through the module command. To load the SNAP module enter:

module load snap-hmm

b) biou:
The SNAP programs are availiable through the Galaxy instance on biou.

To make the SNAP programs availiable through the command line, csh users should enter the following command:

% source /packages/bin/SETUP_BIO_SOFTWARE

To make the SNAP programs availiable through the command line, bash users should enter the following command:

% source /packages/bin/SETUP_BIO_SOFTWARE

2) SNAP Command line usage:

snap [options] <HMM file> <FASTA file > [options]

options:

-help report useful information
-lcmask treat lowercase as N
-plus predict on plus strand only
-minus predict on minus strand only
-gff output annotation as GFF
-ace output annotation as ACED
-quiet do not send progress to STDERR
-aa <file> create FASTA file of proteins
-tx <file> create FASTA file of transcripts
-xdef <file> external definitions
-name <string> name for the gene [default snap]

3) HMM model files for SNAP

To use SNAP, you must either build your own HMM model file or use a precompiled HMM model file.

a) Using pre-compiled HMM model file.

A number of precompiled HMM model files are included in the SNAP release. These files include:

Acanium.hmm
A.gambiae.hmm
A.mellifera.hmm
A.thaliana.hmm
At.hmm
B.malayi.hmm
B.mori.hmm
Ce.hmm
C.elegans.hmm
C.intestinalis.hmm
D.melanogaster.hmm
Dm.hmm
ixodesA.hmm
ixodesB.hmm
mam39.hmm
mam39-ro.hmm
mam46.hmm
mam46-ro.hmm
mam54.hmm
mam54-ro.hmm
mamiso.hmm
minimal.hmm
Nasonia.hmm
nGASP.hmm
nGASPr.hmm
O.sativa.hmm
Os.hmm
worm1.hmm
brugia - same as B.malayi.hmm ciona - same as C.intestinalis.hmm fly - same as D.melanogaster.hmm mosquito - same as A.gambiae.hmm rice - same as O.sativa.hmm thale - same as At.hmm worm - same as Ce.hmm

If your genome is represented above (or is a close relative of a genome represented above), you may use the pre-compiled HMM model file, with the -d option followed by the directory containing the pre-compiled HMM model file. The precompiled HMM model files can be found in the directory:

$SNAPHMM_HOME/HMM

For example to use the precompiled set for the D.melanogaster genome on a set of sequences contained in the file fasta.file, you would use the following on the command line:

% snap $SNAPHMM_HOME/HMM/D.melanogaster.hmm fasta.file

b) compiling your own HMM model files

Note: The author of SNAP ( would like to be contacted should you wish to train SNAP for a new genome. A minimal parameter estimation procedure is outlined below. There are a number of options for forge and hmm-assembler.pl that are not described below.

1) Sequences must be in FASTA format. It's a good idea if you don't have genes that are too related to each other.

2) Gene structures must be in ZFF short format. In the short format, the sequence records are separated by a definition line, just like FASTA. there are 4 fields: Label, Begin, End, Group. The 4th field is optional. Label is a controlled vocabulary including:

    Esngl - single exon gene
    Einit - initial exon
    Eterm - terminal exon
    Exon - generic or internal exon
    (see $SNAPHMM_HOME/Zoe/zoeFeature.h for a complete list).

All exons of a gene must share the same unique group name. The strand of the feature is implied in the coordinates, so if Begin > End, the feature is on the minus strand. Here's an example with two sequences, each containing a single gene on the plus strand:

    >sequence-1
    Einit 201 325 Y73E7A.6
    Eterm 2175 2319 Y73E7A.6
    >sequence-2
    Einit 201 462 Y73E7A.7
    Exon 1803 2031 Y73E7A.7
    Exon 2929 3031 Y73E7A.7
    Exon 3467 3624 Y73E7A.7
    Exon 4185 4406 Y73E7A.7
    Eterm 5103 5280 Y73E7A.7

The most important part of parameter estimation is preparing a training set. There are many ways to go about this. At the end, you want these in the ZFF short format. Save the ZFF as genome.ann and the FASTA as genome.dna.

3) The first step is to look at some features of the genes:

fathom genome.ann genome.dna -gene-stats

4) Next, you want to verify that the genes have no obvious errors:

fathom genome.ann genome.dna -validate

You may find some errors and warnings. Check these out in some kind of genome browser and remove those that are real errors.

5) Next, break up the sequences into fragments with one gene per sequence with the following command:

fathom -genome.ann genome.dna -categorize 1000

There will be up to 1000 bp on either side of the genes. You will find several new files.

    alt.ann, alt.dna (genes with alternative splicing)
    err.ann, err.dna (genes that have errors)
    olp.ann, olp.dna (genes that overlap other genes)
    wrn.ann, wrn.dna (genes with warnings)
    uni.ann, uni.dna (single gene per sequence)

6) Convert the uni genes to plus stranded with the command:

fathom uni.ann uni.dna -export 1000 -plus

You will find 4 new files:

    export.aa proteins corresponding to each gene
    export.ann gene structure on the plus strand
    export.dna DNA of the plus strand
    export.tx transcripts for each gene

7) The parameter estimation program, forge, creates a lot of files. You probably want to create a directory to keep things tidy before you execute the program.

mkdir params
cd params
forge ../export.ann ../export.dna
cd ..

8) Last is to build an HMM.

hmm-assembler.pl my-genome params > my-genome.hmm

4) PBS example files are availiable for running snap on blacklight

Stay Connected

Stay Connected with PSC!

facebook 32 twitter 32 google-Plus-icon