Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Highly Accessed Methodology article

Improving the specificity of high-throughput ortholog prediction

Debra L Fulton12, Yvonne Y Li13, Matthew R Laird1, Benjamin GS Horsman1, Fiona M Roche1 and Fiona SL Brinkman1*

  • * Corresponding author: Fiona SL Brinkman brinkman@sfu.ca

  • † Equal contributors

Author Affiliations

1 Department of Molecular Biology and Biochemistry, Simon Fraser University, Burnaby, BC, Canada

2 Department of Medical Genetics, University of British Columbia, Vancouver, BC, Canada

3 Canada's Michael Smith Genome Sciences Centre, 570 W. 7th Avenue, Vancouver, BC, Canada

For all author emails, please log on.

BMC Bioinformatics 2006, 7:270  doi:10.1186/1471-2105-7-270

Published: 28 May 2006

Abstract

Background

Orthologs (genes that have diverged after a speciation event) tend to have similar function, and so their prediction has become an important component of comparative genomics and genome annotation. The gold standard phylogenetic analysis approach of comparing available organismal phylogeny to gene phylogeny is not easily automated for genome-wide analysis; therefore, ortholog prediction for large genome-scale datasets is typically performed using a reciprocal-best-BLAST-hits (RBH) approach. One problem with RBH is that it will incorrectly predict a paralog as an ortholog when incomplete genome sequences or gene loss is involved. In addition, there is an increasing interest in identifying orthologs most likely to have retained similar function.

Results

To address these issues, we present here a high-throughput computational method named Ortholuge that further evaluates previously predicted orthologs (including those predicted using an RBH-based approach) – identifying which orthologs most closely reflect species divergence and may more likely have similar function. Ortholuge analyzes phylogenetic distance ratios involving two comparison species and an outgroup species, noting cases where relative gene divergence is atypical. It also identifies some cases of gene duplication after species divergence. Through simulations of incomplete genome data/gene loss, we show that the vast majority of genes falsely predicted as orthologs by an RBH-based method can be identified. Ortholuge was then used to estimate the number of false-positives (predominantly paralogs) in selected RBH-predicted ortholog datasets, identifying approximately 10% paralogs in a eukaryotic data set (mouse-rat comparison) and 5% in a bacterial data set (Pseudomonas putida – Pseudomonas syringae species comparison). Higher quality (more precise) datasets of orthologs, which we term "ssd-orthologs" (

    s
upporting-
    s
pecies-
    d
ivergence-orthologs), were also constructed. These datasets, as well as Ortholuge software that may be used to characterize other species' datasets, are available at http://www.pathogenomics.ca/ortholuge/ webcite (software under GNU General Public License).

Conclusion

The Ortholuge method reported here appears to significantly improve the specificity (precision) of high-throughput ortholog prediction for both bacterial and eukaryotic species. This method, and its associated software, will aid those performing various comparative genomics-based analyses, such as the prediction of conserved regulatory elements upstream of orthologous genes.