GeMoMa: Difference between revisions
No edit summary |
|||
Line 28: | Line 28: | ||
For preparing the data, we provide the tool Extractor. You can run Extractor from the command line with<br/> | For preparing the data, we provide the tool Extractor. You can run Extractor from the command line with<br/> | ||
<code>java -jar GeMoMa-1.1. | <code>java -jar GeMoMa-1.1.2.jar CLI Extractor [<parameter>=<value> ...]</code><br/> | ||
The parameters comprise: | The parameters comprise: | ||
Line 69: | Line 69: | ||
<td>transcripts (whether the complete transcripts sequences should returned as output, default = false)</td> | <td>transcripts (whether the complete transcripts sequences should returned as output, default = false)</td> | ||
<td>BOOLEAN</td> | <td>BOOLEAN</td> | ||
</tr> | |||
<tr style="vertical-align:top"> | |||
<td><font color="green">s</font></td> | |||
<td>selected (The path to list file, which allows to make only a predictions for the contained transcript ids, OPTIONAL)</td> | |||
<td>FILE</td> | |||
</tr> | </tr> | ||
<tr style="vertical-align:top"> | <tr style="vertical-align:top"> | ||
Line 83: | Line 88: | ||
For predicting gene models, we provide the tool GeMoMa. You can run GeMoMa from the command line with<br/> | For predicting gene models, we provide the tool GeMoMa. You can run GeMoMa from the command line with<br/> | ||
<code>java -jar GeMoMa-1.1. | <code>java -jar GeMoMa-1.1.2.jar CLI GeMoMa [<parameter>=<value> ...]</code><br/> | ||
The parameters comprise: | The parameters comprise: | ||
Line 237: | Line 242: | ||
GeMoMa returns the predicted annotation as gff file and the predicted proteins as fasta file. | GeMoMa returns the predicted annotation as gff file and the predicted proteins as fasta file. | ||
== Version history == | |||
[http://www.jstacs.de/download.php?which=GeMoMa GeMoMa 1.1.2] (05.02.2016) | |||
* GeMoMa bugfix (upstream, downstream sequence for splice site detection) | |||
* Extractor: new parameter s for selecting transcripts | |||
* improved Galaxy integration | |||
[http://www.jstacs.de/downloads/GeMoMa-1.1.1.zip GeMoMa 1.1.1] (01.02.2016) | |||
* initial release for paper |
Revision as of 15:01, 9 February 2016
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.1.2.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"
| ||||
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 | |||
s | selected (The path to list file, which allows to make only a predictions for the contained transcript ids, OPTIONAL) | FILE | |||
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.1.2.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"
| ||||
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"
| ||||
s | substitution matrix (the substitution matrix used in the alignment, range={default, user-specified}, default = default) No parameters for selection "default"
| ||||
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.
Version history
GeMoMa 1.1.2 (05.02.2016)
- GeMoMa bugfix (upstream, downstream sequence for splice site detection)
- Extractor: new parameter s for selecting transcripts
- improved Galaxy integration
GeMoMa 1.1.1 (01.02.2016)
- initial release for paper