title: Integration of mapped RNA-Seq reads into automatic training of eukaryotic gene finding algorithm.
abstract: We present a new approach to automatic training of a eukaryotic ab initio gene finding algorithm. With the advent of Next-Generation Sequencing, automatic training has become paramount, allowing genome annotation pipelines to keep pace with the speed of genome sequencing. Earlier we developed GeneMark-ES, currently the only gene finding algorithm for eukaryotic genomes that performs automatic training in unsupervised ab initio mode. The new algorithm, GeneMark-ET augments GeneMark-ES with a novel method that integrates RNA-Seq read alignments into the self-training procedure. Use of 'assembled' RNA-Seq transcripts is far from trivial; significant error rate of assembly was revealed in recent assessments. We demonstrated in computational experiments that the proposed method of incorporation of 'unassembled' RNA-Seq reads improves the accuracy of gene prediction; particularly, for the 1.3 GB genome of Aedes aegypti the mean value of prediction Sensitivity and Specificity at the gene level increased over GeneMark-ES by 24.5%. In the current surge of genomic data when the need for accurate sequence annotation is higher than ever, GeneMark-ET will be a valuable addition to the narrow arsenal of automatic gene prediction tools.
authors.affiliation: Joint Georgia Tech and Emory Wallace H. Coulter Department of Biomedical Engineering, Atlanta, GA, USA 30332.
authors.forename: Alexandre
authors.initials: A
authors.lastname: Lomsadze
authors.affiliation: Joint Georgia Tech and Emory Wallace H. Coulter Department of Biomedical Engineering, Atlanta, GA, USA 30332.
authors.forename: Paul D
authors.initials: PD
authors.lastname: Burns
authors.affiliation: Joint Georgia Tech and Emory Wallace H. Coulter Department of Biomedical Engineering, Atlanta, GA, USA 30332 School of Computational Science and Engineering, Georgia Tech, Atlanta, GA, USA 30332 Department of Bioinformatics, Moscow Institute of Physics and Technology, Moscow, Russia 141700 borodovsky@gatech.edu.
authors.forename: Mark
authors.initials: M
authors.lastname: Borodovsky
mesh_terms: D000465:Algorithms; D000818:Animals; D009033:Culicidae / Q000235:genetics; D004331:Drosophila melanogaster / Q000235:genetics; D020869:Gene Expression Profiling; D005796:Genes; D017344:Genes, Insect; D059014:High-Throughput Nucleotide Sequencing / Q000379:methods*; D016415:Sequence Alignment / Q000379:methods*; D017423:Sequence Analysis, RNA / Q000379:methods*