CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes

MOTIVATION The numbers of finished and ongoing genome projects are increasing at a rapid rate, and providing the catalog of genes for these new genomes is a key challenge. Obtaining a set of well-characterized genes is a basic requirement in the initial steps of any genome annotation process. An accurate set of genes is needed in order to learn about species-specific properties, to train gene-finding programs, and to validate automatic predictions. Unfortunately, many new genome projects lack comprehensive experimental data to derive a reliable initial set of genes. RESULTS In this study, we report a computational method, CEGMA (Core Eukaryotic Genes Mapping Approach), for building a highly reliable set of gene annotations in the absence of experimental data. We define a set of conserved protein families that occur in a wide range of eukaryotes, and present a mapping procedure that accurately identifies their exon-intron structures in a novel genomic sequence. CEGMA includes the use of profile-hidden Markov models to ensure the reliability of the gene structures. Our procedure allows one to build an initial set of reliable gene annotations in potentially any eukaryotic genome, even those in draft stages. AVAILABILITY Software and data sets are available online at http://korflab.ucdavis.edu/Datasets.

[1]  Sean R. Eddy,et al.  Profile hidden Markov models , 1998, Bioinform..

[2]  Michael Y. Galperin,et al.  Using the COG Database to Improve Gene Recognition in Complete Genomes , 2004, Genetica.

[3]  Michael R Brent,et al.  Genome annotation past, present, and future: how to define an ORF at each locus. , 2005, Genome research.

[4]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[5]  R. Durbin,et al.  GeneWise and Genomewise. , 2004, Genome research.

[6]  C. Burge,et al.  Computational inference of homologous gene structures in the human genome. , 2001, Genome research.

[7]  H. Akashi,et al.  Translational selection and molecular evolution. , 1998, Current opinion in genetics & development.

[8]  Darren A. Natale,et al.  The COG database: an updated version includes eukaryotes , 2003, BMC Bioinformatics.

[9]  Jeff Shrager,et al.  Chlamydomonas reinhardtii at the Crossroads of Genomics , 2003, Eukaryotic Cell.

[10]  M. Borodovsky,et al.  Recognition of genes in DNA sequence with ambiguities. , 1993, Bio Systems.

[11]  R. Guigó,et al.  GeneID in Drosophila. , 2000, Genome research.

[12]  M. Borodovsky,et al.  Gene identification in novel eukaryotic genomes by self-training algorithm , 2005, Nucleic acids research.

[13]  S. Karlin,et al.  Prediction of complete gene structures in human genomic DNA. , 1997, Journal of molecular biology.

[14]  D. Higgins,et al.  T-Coffee: A novel method for fast and accurate multiple sequence alignment. , 2000, Journal of molecular biology.

[15]  X. Gu,et al.  Intron gain and loss in segmentally duplicated genes in rice , 2006, Genome Biology.

[16]  J. Harrow,et al.  GENCODE: producing a reference annotation for ENCODE , 2006, Genome Biology.

[17]  Jian Wang,et al.  The Genome Sequence of the Malaria Mosquito Anopheles gambiae , 2002, Science.

[18]  Ian Korf,et al.  Gene finding in novel genomes , 2004, BMC Bioinformatics.

[19]  T. Andrews,et al.  The Ensembl automatic gene annotation system. , 2004, Genome research.

[20]  Paul Richardson,et al.  The Draft Genome of Ciona intestinalis: Insights into Chordate and Vertebrate Origins , 2002, Science.

[21]  P. Pevzner,et al.  Gene recognition via spliced sequence alignment. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[22]  R. Guigó,et al.  Evaluation of gene structure prediction programs. , 1996, Genomics.

[23]  Ian Korf,et al.  Integrating genomic homology into gene structure prediction , 2001, ISMB.

[24]  R. Guigó,et al.  Comparative gene prediction in human and mouse. , 2003, Genome research.

[25]  Li Li,et al.  ToxoDB: accessing the Toxoplasma gondii genome , 2003, Nucleic Acids Res..