Comparison of RefSeq protein-coding regions in human and vertebrate genomes.
BACKGROUND
Advances in high-throughput sequencing technology have yielded a large number of publicly available vertebrate genomes, many of which are selected for inclusion in NCBI's RefSeq project and subsequently processed by NCBI's eukaryotic annotation pipeline. Genome annotation results are affected by differences in available support evidence and may be impacted by annotation pipeline software changes over time. The RefSeq project has not previously assessed annotation trends across organisms or over time. To address this deficiency, we have developed a comparative protocol which integrates analysis of annotated protein-coding regions across a data set of vertebrate orthologs in genomic sequence coordinates, protein sequences, and protein features.


RESULTS
We assessed an ortholog dataset that includes 34 annotated vertebrate RefSeq genomes including human. We confirm that RefSeq protein-coding gene annotations in mammals exhibit considerable similarity. Over 50% of the orthologous protein-coding genes in 20 organisms are supported at the level of splicing conservation with at least three selected reference genomes. Approximately 7,500 ortholog sets include at least half of the analyzed organisms, show highly similar sequence and conserved splicing, and may serve as a minimal set of mammalian "core proteins" for initial assessment of new mammalian genomes. Additionally, 80% of the proteins analyzed pass a suite of tests to detect proteins that lack splicing conservation and have unusual sequence or domain annotation. We use these tests to define an annotation quality metric that is based directly on the annotated proteins thus operates independently of other quality metrics such as availability of transcripts or assembly quality measures. Results are available on the RefSeq FTP site [http://ftp.ncbi.nlm.nih.gov/refseq/supplemental/ProtCore/SM1.txt].


CONCLUSIONS
Our multi-factored analysis demonstrates a high level of consistency in RefSeq protein representation among vertebrates. We find that the majority of the RefSeq vertebrate proteins for which we have calculated orthology are good as measured by these metrics. The process flow described provides specific information on the scope and degree of conservation for the analyzed protein sequences and annotations and will be used to enrich the quality of RefSeq records by identifying targets for further improvement in the computational annotation pipeline, and by flagging specific genes for manual curation.
National Center for Biotechnology Information, U,S, National Library of Medicine, 8600 Rockville Pike, Bethesda, MD 20894, USA. pruitt@ncbi.nlm.nih.gov.
Jessica H
JH
Fong
Terence D
TD
Murphy
Kim D
KD
Pruitt
D000595:Amino Acid Sequence; D000818:Animals; D017124:Conserved Sequence / Q000235:genetics; D030541:Databases, Genetic; D030561:Databases, Nucleic Acid; D016678:Genome / Q000235:genetics*; D015894:Genome, Human / Q000235:genetics*; D006801:Humans; D058977:Molecular Sequence Annotation; D016366:Open Reading Frames / Q000235:genetics*; D017434:Protein Structure, Tertiary; D012326:RNA Splicing / Q000235:genetics; D012689:Sequence Homology, Nucleic Acid; D013045:Species Specificity; D014714:Vertebrates / Q000235:genetics*; D000511:alpha-Macroglobulins / Q000235:genetics
D000511:alpha-Macroglobulins