Evidence Directed Gene Construction for Eukaryotes
Don Gilbert, gilbertd at indiana edu, 2010 .. 2013
Informant evidence for gene models is transcript and protein
data. Transcripts from long (e.g. dbEST) and short reads (e.g.
Illumina), are mapped to the genome assembly with GMAP (long)
or GSNAP (short).
EST assemblies are constructed using PASA, from all transcript
data. RNA-seq assemblies are first constructed from aligned short
reads with Cufflinks, and combined with EST reads for a full
Protein genes from related sequenced genomes are BLASTX aligned to
repeat/transposon-soft-masked genome, then refined with Exonerate to
protein gene models.
Gene models are predicted with evidence-directed AUGUSTUS
predictor. Augustus is trained for gene parameters
using full cDNA genes derirved from the PASA EST assemblies.
Several prediction sets from different evidence sets and
parameters are combined, selecting highest evidence scored model
per locus, and avoiding artifactual gene splits and joins.
This consensus gene set has the best match to EST and protein
Predicted genes are UTR-extended and/or improved by PASA. Genes
are annotated with Uniprot descriptions, and classified by
evidence scores including transposable elements.
NCBI BLAST: blast.ncbi.nlm.nih.gov/