Directly to Contents

Navigation for:     BRAKER
BRAKER

Availability
Download BRAKER2 from http://bioinf.uni-greifswald.de/augustus/binaries/BRAKER2.tar.gz. The most recent release is version 2.1.0, March 13th 2018. Example data for testing the BRAKER2 pipeline is available at http://bioinf.uni-greifswald.de/augustus/binaries/BRAKER2examples.tar.gz (1.2 GB).

BRAKER is a tool for fully automated genome annotation with GeneMark-ET and AUGUSTUS. BRAKER is a joint project of Georgia Institute of Technology, USA and Institute for Mathematics and Computer Science, University of Greifswald, Germany.

In its initial version, BRAKER1 was able to process genome and RNA-Seq data, only. BioTechniques reported on BRAKER1 on 02/05/2016.

BRAKER1 required two input files:

First, GeneMark-ET performs RNA-Seq supported iterative training and generates initial gene structures. Second, AUGUSTUS uses predicted genes for training and then integrates RNA-Seq read information as extrinsic evidence into final gene predictions.

With the most recent release, BRAKER2, training of GeneMark-EP from alignment information of proteins from species with a longer evolutionary distance to the target genome is possible. You may thus achieve good gene prediction accuracy of GeneMark-EP and AUGUSTUS in the absence of RNA-Seq data. BRAKER1 functionality is fully maintained by BRAKER2. You find more information about BRAKER2 on our poster and in the slides from our talk at Plants and Animal Genomes conference 2018.

Accuracy

We compare prediction accuracy of BRAKER1 on four model species genomes to accuracy of MAKER2 and CodingQuarry (only applicable to fungi). The following table is an excerpt from our publication BRAKER1: Unsupervised RNA-Seq-Based Genome Annotation with GeneMark-ET and AUGUSTUS, Table1:

LevelArabidopsis thaliana Caenorhabditis elegans Drosophila melanogaster Schizosaccharomyces pombe
BRAKER1MAKER2BRAKER1MAKER2BRAKER1MAKER2BRAKER1MAKER2Coding Quarry
Gene Sens.64.451.355.041.064.955.277.442.879.7
Gene Spec.52.052.255.230.859.446.380.568.772.6
Exon Sens.82.976.180.269.475.066.483.250.179.6
Exon Spec.79.076.185.362.381.766.983.271.481.7

Usage

For straight forward gene prediction in a genome file with a RNA-Seq alignment bam file, run the following command:

perl braker.pl --genome=genome.fa --bam=RNAseq.bam

For usage with protein data, choose one aligner out of GenomeThreader (gth), Exonerate (exonerate) and Spaln (spaln):

perl braker.pl --genome=genome.fa --bam=RNAseq.bam --prot_seq=protein.fa --prg=(gth|exonerate|spaln) --ALIGNMENT_TOOL_PATH=/path/to/aligner

For usage with protein alignment data from remotely related species (result of the GeneMark-EP protein mapping pipeline):

perl braker.pl --genome=genome.fa --hints=protein.gff --epmode

BRAKER2 produces three important output files in the working directory:

Note:

We have to correct one important reference in the BRAKER1 publication. In computations of the gene prediction accuracy for the D. melanogaster genome we used the r6.07 version of the fly genome and annotation. However, the Supplementary materials to the paper (available at the "Bioinformatics" journal website) incorrectly cite the earlier r5.55 version of the D. melanogaster genome.

Publications

Please cite the following publications when using BRAKER for your project:

CONTACT
Institute for Mathematics und Computer Sciences
Walther-Rathenau-Straße 47
17487 Greifswald
Germany
Tel.: +49 (0)3834 86 - 46 24
Fax: +49 (0)3834 86 - 46 40

mario.stanke@uni-greifswald.de