GeMoMa: Difference between revisions

From Jstacs
Jump to navigationJump to search
Line 105: Line 105:
<tr style="vertical-align:top">
<tr style="vertical-align:top">
<td><font color="green">a</font></td>
<td><font color="green">a</font></td>
<td>assignment (The assignment file, which combines parts of the CDS to transcripts)</td>
<td>assignment (The assignment file, which combines parts of the CDS to transcripts, OPTIONAL)</td>
<td>FILE</td>
<td>FILE</td>
</tr>
</tr>
Line 172: Line 172:
<td><font color="green">e</font></td>
<td><font color="green">e</font></td>
<td>e-value (The e-value for filtering blast results, default = 100.0)</td>
<td>e-value (The e-value for filtering blast results, default = 100.0)</td>
<td>DOUBLE</td>
</tr>
<tr style="vertical-align:top">
<td><font color="green">ct</font></td>
<td>contig threshold (The threshold for evaluating contigs, valid range = [0.0, 1.0], default = 0.9)</td>
<td>DOUBLE</td>
</tr>
<tr style="vertical-align:top">
<td><font color="green">r</font></td>
<td>region threshold (The threshold for evaluating regions, valid range = [0.0, 1.0], default = 0.9)</td>
<td>DOUBLE</td>
</tr>
<tr style="vertical-align:top">
<td><font color="green">h</font></td>
<td>hit threshold (The threshold for adding additional hits, valid range = [0.0, 1.0], default = 0.9)</td>
<td>DOUBLE</td>
<td>DOUBLE</td>
</tr>
</tr>
Line 198: Line 213:
<td>align (A flag which allows to output a tab-delimited file, which contains the results in a blast-like format (deprecated), default = false)</td>
<td>align (A flag which allows to output a tab-delimited file, which contains the results in a blast-like format (deprecated), default = false)</td>
<td>BOOLEAN</td>
<td>BOOLEAN</td>
</tr>
<tr style="vertical-align:top">
<td><font color="green">genomic</font></td>
<td>genomic (A flag which allows to output a fasta file containing the genomic regions of the predictions, default = false)</td>
<td>BOOLEAN</td>
</tr>
<tr style="vertical-align:top">
<td><font color="green">prefix</font></td>
<td>prefix (A prefix to be used for naming the predictions, default = )</td>
<td>STRING</td>
</tr>
</tr>
<tr style="vertical-align:top">
<tr style="vertical-align:top">

Revision as of 17:17, 22 December 2015

by Jens Keilwagen, Michael Wenk, Jessica L. Erickson, Martin H. Schattat, Jan Grau, and Frank Hartung

Gene Model Mapper (GeMoMa) is a homology-based gene prediction program that uses the annotation of protein-coding genes in a reference genome to infer annotation of protein-coding genes in a target genome.

Abstract

Annotation of protein-coding genes is very important in many fields of research and application. Thereby, homology-based gene prediction programs allow for transferring knowledge from an annotated organism to an organism of interest.

Here, we present a homology-based gene prediction program called GeMoMa. GeMoMa utilizes the intron position conservation of related genes in different organisms. We assess the performance of GeMoMa comparing it with state-of-the-art competitors on plant and animal genomes using an extended best reciprocal hit approach. We find that it often makes more precise predictions than its competitors yielding an improvement of up to 622% more correct transcripts. Subsequently, we use RNA-seq data to compare the predictions of homology-based gene prediction programs, and find again that GeMoMa performs well.

Hence, we conclude that exploiting intron position conservation improves homology-based gene prediction, and we make GeMoMa freely available as command-line tool and Galaxy integration.

Paper

We have submitted the paper for review.

Download

GeMoMa is implemented in Java using Jstacs. You can download a zip file containing a readme, the GeMoMa jar file and some tiny scripts for running GeMoMa. The jar file allows for

  • running the command line interface (CLI) version
  • creating the XML file needed for the Galaxy integration.

Running the command line application

For running the command line application, Java v1.6 or later is required.

For preparing the data, we provide the tool Extractor. You can run Extractor from the command line with
java -jar GeMoMa-1.0.jar CLI Extractor [<parameter>=<value> ...]
The parameters comprise:

name comment type

a annotation (Reference annotation file (GFF), which contains gene models annotated in the reference genome) FILE
g genome (Reference genome file (FASTA)) FILE
gc genetic code (whether to use the default or a user-specified genetic code, range={default, user-specified}, default = default)

No parameters for selection "default"
Parameters for selection "user-specified":

u user code (user-specified genetic code) FILE
p proteins (whether the complete proteins sequences should returned as output, default = false) BOOLEAN
t transcripts (whether the complete transcripts sequences should returned as output, default = false) BOOLEAN
v verbose (A flag which allows to output wealth of additional information, default = false) BOOLEAN
outdir The output directory, defaults to the current working directory (.) STRING

For predicting gene models, we provide the tool GeMoMa. You can run GeMoMa from the command line with
java -jar GeMoMa-1.0.jar CLI GeMoMa [<parameter>=<value> ...]
The parameters comprise:

name comment type

t tblastn results (The sorted tblastn results) FILE
q query cds parts (The query cds parts file (FASTA), i.e., the cds parts that have been blasted) FILE
a assignment (The assignment file, which combines parts of the CDS to transcripts, OPTIONAL) FILE
tg target genome (The target genome file (FASTA), i.e., the target sequences in the blast run) FILE
alignment alignment (for computing the optimal alignment score of the complete prediction vs. the query protein, range={no, yes}, default = no)

No parameters for selection "no"
Parameters for selection "yes":

qp query proteins (The path to the query protein file (FASTA), OPTIONAL) FILE
g genetic code (whether to use the default or a user-specified genetic code, range={default, user-specified}, default = default)

No parameters for selection "default"
Parameters for selection "user-specified":

c code (user-specified genetic code) FILE
s substitution matrix (the substitution matrix used in the alignment, range={default, user-specified}, default = default)

No parameters for selection "default"
Parameters for selection "user-specified":

m matrix (user-specified substitution matrix) FILE
go gap opening (The gap opening cost in the alignment, default = 11) INT
ge gap extension (The gap extension cost in the alignment, default = 1) INT
mil maximum intron length (The maximum length of an intron, default = 15000) INT
i intron-loss-gain-penalty (The penalty used for intron loss and gain, default = 25) INT
e e-value (The e-value for filtering blast results, default = 100.0) DOUBLE
ct contig threshold (The threshold for evaluating contigs, valid range = [0.0, 1.0], default = 0.9) DOUBLE
r region threshold (The threshold for evaluating regions, valid range = [0.0, 1.0], default = 0.9) DOUBLE
h hit threshold (The threshold for adding additional hits, valid range = [0.0, 1.0], default = 0.9) DOUBLE
p predictions (The (maximal) number of predictions per transcript, default = 1) INT
selected selected (The path to list file, which allows to make only a predictions for the contained transcript ids, OPTIONAL) FILE
as avoid stop (A flag which allows to avoid stop codons in a transcript (except the last AS), default = true) BOOLEAN
approx approx (whether an approximation is used to compute the score for intron gain, default = true) BOOLEAN
align align (A flag which allows to output a tab-delimited file, which contains the results in a blast-like format (deprecated), default = false) BOOLEAN
genomic genomic (A flag which allows to output a fasta file containing the genomic regions of the predictions, default = false) BOOLEAN
prefix prefix (A prefix to be used for naming the predictions, default = ) STRING
v verbose (A flag which allows to output wealth of additional information per transcript, default = false) BOOLEAN
outdir The output directory, defaults to the current working directory (.) STRING

GeMoMa returns the predicted annotation as gff file and the predicted proteins as fasta file.