Coding Exon-structure Aware Realigner (CESAR) Utilizes Genome Alignments for Accurate Comparative Gene Annotation

Overview

Journal Nucleic Acids Res

Publisher Oxford University Press

Specialty Biochemistry

Date 2016 Mar 27

PMID 27016733

Citations 29

Authors

Virag Sharma

Anas Elghafari

Michael Hiller

Affiliations

Soon will be listed here.

Abstract

Identifying coding genes is an essential step in genome annotation. Here, we utilize existing whole genome alignments to detect conserved coding exons and then map gene annotations from one genome to many aligned genomes. We show that genome alignments contain thousands of spurious frameshifts and splice site mutations in exons that are truly conserved. To overcome these limitations, we have developed CESAR (Coding Exon-Structure Aware Realigner) that realigns coding exons, while considering reading frame and splice sites of each exon. CESAR effectively avoids spurious frameshifts in conserved genes and detects 91% of shifted splice sites. This results in the identification of thousands of additional conserved exons and 99% of the exons that lack inactivating mutations match real exons. Finally, to demonstrate the potential of using CESAR for comparative gene annotation, we applied it to 188 788 exons of 19 865 human genes to annotate human genes in 99 other vertebrates. These comparative gene annotations are available as a resource (http://bds.mpi-cbg.de/hillerlab/CESAR/). CESAR (https://github.com/hillerlab/CESAR/) can readily be applied to other alignments to accurately annotate coding genes in many other vertebrate and invertebrate genomes.

Citing Articles

Conservation assessment of human splice site annotation based on a 470-genome alignment.

Minkin I, Salzberg S, Salzberg S bioRxiv. 2023; .

PMID: 38076842 PMC: 10705407. DOI: 10.1101/2023.12.01.569581.

High-quality haploid genomes corroborate 29 chromosomes and highly conserved synteny of genes in Hyles hawkmoths (Lepidoptera: Sphingidae).

Hundsdoerfer A, Schell T, Patzold F, Wright C, Yoshido A, Marec F BMC Genomics. 2023; 24(1):443.

PMID: 37550607 PMC: 10405479. DOI: 10.1186/s12864-023-09506-y.

ncOrtho: efficient and reliable identification of miRNA orthologs.

Langschied F, Leisegang M, Brandes R, Ebersberger I Nucleic Acids Res. 2023; 51(13):e71.

PMID: 37260093 PMC: 10359484. DOI: 10.1093/nar/gkad467.

Integrating gene annotation with orthology inference at scale.

Kirilenko B, Munegowda C, Osipova E, Jebb D, Sharma V, Blumer M Science. 2023; 380(6643):eabn3107.

PMID: 37104600 PMC: 10193443. DOI: 10.1126/science.abn3107.

Building the Chordata Olfactory Receptor Database using more than 400,000 receptors annotated by Genome2OR.

Han W, Wu Y, Zeng L, Zhao S Sci China Life Sci. 2022; 65(12):2539-2551.

PMID: 35696018 DOI: 10.1007/s11427-021-2081-6.

References

Curwen V, Eyras E, Andrews T, Clarke L, Mongin E, Searle S . The Ensembl automatic gene annotation system. Genome Res. 2004; 14(5):942-50. PMC: 479124. DOI: 10.1101/gr.1858004. View

Lindblad-Toh K, Garber M, Zuk O, Lin M, Parker B, Washietl S . A high-resolution map of human evolutionary constraint using 29 mammals. Nature. 2011; 478(7370):476-82. PMC: 3207357. DOI: 10.1038/nature10530. View

. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012; 489(7414):57-74. PMC: 3439153. DOI: 10.1038/nature11247. View

Lu D, Brown R, Arumugam M, Brent M . Pairagon: a highly accurate, HMM-based cDNA-to-genome aligner. Bioinformatics. 2009; 25(13):1587-93. PMC: 2732315. DOI: 10.1093/bioinformatics/btp273. View

Speir M, Zweig A, Rosenbloom K, Raney B, Paten B, Nejad P . The UCSC Genome Browser database: 2016 update. Nucleic Acids Res. 2015; 44(D1):D717-25. PMC: 4702902. DOI: 10.1093/nar/gkv1275. View

Kuhn R, Karolchik D, Zweig A, Wang T, Smith K, Rosenbloom K . The UCSC Genome Browser Database: update 2009. Nucleic Acids Res. 2008; 37(Database issue):D755-61. PMC: 2686463. DOI: 10.1093/nar/gkn875. View

Brent M, Guigo R . Recent advances in gene structure prediction. Curr Opin Struct Biol. 2004; 14(3):264-72. DOI: 10.1016/j.sbi.2004.05.007. View

Hiller M, Huse K, Szafranski K, Jahn N, Hampe J, Schreiber S . Widespread occurrence of alternative splicing at NAGNAG acceptors contributes to proteome plasticity. Nat Genet. 2004; 36(12):1255-7. DOI: 10.1038/ng1469. View

Robinson G, Hackett K, Purcell-Miramontes M, Brown S, Evans J, Goldsmith M . Creating a buzz about insect genomes. Science. 2011; 331(6023):1386. DOI: 10.1126/science.331.6023.1386. View

10.

Schneider A, Cannarozzi G, Gonnet G . Empirical codon substitution matrix. BMC Bioinformatics. 2005; 6:134. PMC: 1173088. DOI: 10.1186/1471-2105-6-134. View

11.

Clarke S, VanderMeer J, Wenger A, Schaar B, Ahituv N, Bejerano G . Human developmental enhancers conserved between deuterostomes and protostomes. PLoS Genet. 2012; 8(8):e1002852. PMC: 3410860. DOI: 10.1371/journal.pgen.1002852. View

12.

Dewey C . Whole-genome alignment. Methods Mol Biol. 2012; 855:237-57. DOI: 10.1007/978-1-61779-582-4_8. View

13.

Alexandersson M, Cawley S, Pachter L . SLAM: cross-species gene finding and alignment with a generalized pair hidden Markov model. Genome Res. 2003; 13(3):496-502. PMC: 430255. DOI: 10.1101/gr.424203. View

14.

Gross S, Brent M . Using multiple alignments to improve gene prediction. J Comput Biol. 2006; 13(2):379-93. DOI: 10.1089/cmb.2006.13.379. View

15.

Stedman H, Kozyak B, Nelson A, Thesier D, Su L, Low D . Myosin gene mutation correlates with anatomical changes in the human lineage. Nature. 2004; 428(6981):415-8. DOI: 10.1038/nature02358. View

16.

Stanke M, Diekhans M, Baertsch R, Haussler D . Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008; 24(5):637-44. DOI: 10.1093/bioinformatics/btn013. View

17.

Hsu F, Kent W, Clawson H, Kuhn R, Diekhans M, Haussler D . The UCSC Known Genes. Bioinformatics. 2006; 22(9):1036-46. DOI: 10.1093/bioinformatics/btl048. View

18.

Hiller M, Agarwal S, Notwell J, Parikh R, Guturu H, Wenger A . Computational methods to detect conserved non-genic elements in phylogenetically isolated genomes: application to zebrafish. Nucleic Acids Res. 2013; 41(15):e151. PMC: 3753653. DOI: 10.1093/nar/gkt557. View

19.

Gross S, Do C, Sirota M, Batzoglou S . CONTRAST: a discriminative, phylogeny-free approach to multiple informant de novo gene prediction. Genome Biol. 2007; 8(12):R269. PMC: 2246271. DOI: 10.1186/gb-2007-8-12-r269. View

20.

Davydov E, Goode D, Sirota M, Cooper G, Sidow A, Batzoglou S . Identifying a high fraction of the human genome to be under selective constraint using GERP++. PLoS Comput Biol. 2010; 6(12):e1001025. PMC: 2996323. DOI: 10.1371/journal.pcbi.1001025. View