- Open Access
QuartetS-DB: a large-scale orthology database for prokaryotes and eukaryotes inferred by evolutionary evidence
© Yu et al.; licensee BioMed Central Ltd. 2012
- Received: 20 December 2011
- Accepted: 28 May 2012
- Published: 22 June 2012
The concept of orthology is key to decoding evolutionary relationships among genes across different species using comparative genomics. QuartetS is a recently reported algorithm for large-scale orthology detection. Based on the well-established evolutionary principle that gene duplication events discriminate paralogous from orthologous genes, QuartetS has been shown to improve orthology detection accuracy while maintaining computational efficiency.
QuartetS-DB is a new orthology database constructed using the QuartetS algorithm. The database provides orthology predictions among 1621 complete genomes (1365 bacterial, 92 archaeal, and 164 eukaryotic), covering more than seven million proteins and four million pairwise orthologs. It is a major source of orthologous groups, containing more than 300,000 groups of orthologous proteins and 236,000 corresponding gene trees. The database also provides over 500,000 groups of inparalogs. In addition to its size, a distinguishing feature of QuartetS-DB is the ability to allow users to select a cutoff value that modulates the balance between prediction accuracy and coverage of the retrieved pairwise orthologs. The database is accessible at https://applications.bioanalysis.org/quartetsdb.
QuartetS-DB is one of the largest orthology resources available to date. Because its orthology predictions are underpinned by evolutionary evidence obtained from sequenced genomes, we expect its accuracy to continue to increase in future releases as the genomes of additional species are sequenced.
- Orthology detection
- Orthology database
The ability to sequence complete genomes is transforming life-science research, especially as it becomes faster, less expensive, and more widely available. Access to this vast amount of genomic data allows us to perform large-scale comparative genomic analyses, so that experimental knowledge can be transferred from well-studied organisms to newly sequenced genomes. Critical to this transfer of knowledge from one organism to another is the concept of homology. Homologous genes evolve from a common ancestral gene, and can be characterized as orthologs or paralogs [1, 2]. Orthologous genes evolve through speciation events and, as such, are expected to have retained their molecular functions. Although the functional similarity of orthologs is still being debated [3, 4], orthology has been used to infer gene function , identify conserved regulatory and functional regions across multiple species , predict novel signaling pathways , identify common, broad-spectrum protein targets across pathogenic organisms , and select pathogen-specific drug targets without human orthologs to minimize potential detrimental effects in patients . Conversely, paralogous genes evolve through duplication events, which offer an opportunity for genes to escape selection pressures and undergo mutations leading to new functional roles [10, 11]. Even though, in general, paralogs are not expected to share similar molecular functions, a certain type (termed inparalogs), which evolved from recent duplications subsequent to a given speciation event , are also likely to share similar functions. Some orthology databases, such as InParanoid , provide within-species inparalog information.
With the recent growth of databases containing complete genome sequences, there has been a parallel growth in the development of new orthology databases, derived from algorithms employing various evolutionary assumptions and approximations. In the past five years, we have observed a steady growth in the number of accessible orthology databases, including many successive updates, where each release covers a larger number of complete genome sequences. For example, OrthoDB’s coverage of complete genomes increased from 57 in its original release to 191 in its last update [13, 14], PhylomeDB increased from 443 to 717 [15, 16], eggNOG increased from 373 to 630 [17, 18], and OMA’s coverage increased nearly 7-fold from 150 to 1000 genomes between 2005 and 2010 . (By the time of this writing, PhylomeDB, eggNOG, and OMA had increased their coverage to 1415, 1133, and 1109 genomes, respectively).
To address the accelerated growth in the number of available genome sequences, orthology detection methods must become more computationally efficient, even when extensive high-performance computational resources are available. One possibility, suggested by Altenhoff et al.  is to develop algorithms that bypass portions of the computationally intensive, all-against-all comparative procedure widely used by methods based on bi-directional best hits (BBH). While minimizing such costly pairwise comparisons would be computationally advantageous, care must be taken not to compromise detection accuracy. Hence, the challenge is to develop new methods that can handle large-scale applications (e.g., thousands of genomes) while balancing often diametrically opposed objectives: detection accuracy and computational efficiency.
Recently, we developed a novel orthology detection method, termed QuartetS, that attempts to balance detection accuracy and computational efficiency . QuartetS is a sequence similarity-based orthology prediction method grounded on the well-established evolutionary concept that gene duplication events distinguish orthologous from paralogous relationships . It provides accurate predictions by evaluating whether putative orthologous pairs identified by BLAST analysis as BBH across two genomes have originated from a duplication event. This is typically achieved by analyzing phylogenetic gene trees formed by the two genes of interest and two genes from a third genome, for all available genomes. However, because the construction of phylogenetic trees is computationally intensive, to achieve computational efficiency QuartetS approximates this process via an analytic expression that considers pairwise sequence similarities identified by BLAST. Since evolutionary evidence is extracted from all genomes, we expect the accuracy of QuartetS to continue to increase as additional organisms are sequenced.
In a systematic, large-scale comparison of 624 bacterial genomes using both function- and phylogeny-based metrics , we showed that QuartetS slightly, but consistently, outperformed the highly specific Orthologous Matrix (OMA) method . This is notable because in a recent comparative study , OMA outperformed a large number of orthology detection methods, and its database is one of the largest of its kind, including orthology predictions from 1000 complete genomes .
Phylum distribution of the 1621 species in QuartetS-DB
Number of species
The QuartetS method
provides the likelihood of a duplication event inferred by z1 and z2 along the evolution of the two target genes x and y. Thus, we infer that x and y are orthologs when α is smaller than a specified cutoff value Ω, with a smaller Ω leading to the identification of fewer, more accurate orthologs.
Identifying putative orthologous pairs x and y flagged by BLAST analysis as BBH across two species;
Computing α for a set of gene pairs z1 and z2 in a third species that form BBH with x and y, and computing all possible α for each and every species in a genome database; and
Predicting x and y as orthologs if α ≤ Ω, for all α.
This strategy has been shown to produce an acceptable tradeoff between detection accuracy and computational efficiency. The computational cost to fetch the sequence similarity values Si,j adds less than 0.5% to the cost of obtaining BBH pairs with BLAST, while the prediction accuracy is similar to that obtained with precisely constructed gene trees .
Acquisition of protein sequence data
We collected the protein sequence data of complete genomes from the NCBI RefSeq database release 43 . This included nearly five million proteins in 1457 prokaryotic species (1365 bacterial and 92 archaeal) and over two million proteins in 164 eukaryotic species. In QuartetS-DB we identified proteins using NCBI GenInfo Identifier (GI) numbers. However, because users may wish to query the database using naming conventions linked to records in other resources, we also downloaded additional protein identifiers, including the NCBI Gene ID, RefSeq accession, Locus tag, and UniProt accession, from the RefSeq database (ftp://ftp.ncbi.nih.gov/refseq/release/release-catalog/) and mapped them to the corresponding protein GI numbers. We obtained protein annotations, such as gene symbol, functional description, and Gene Ontology (GO) terms (GO description and GO accession numbers), from the same site. We named all 1621 species in the database according to RefSeq records and uniquely identified them by their NCBI taxonomy identifier. It should be noted that the RefSeq database contains multiple protein isoforms translated by alternative splicing of the same gene. In QuartetS-DB, we kept all of these isoforms (whose functions might be slightly different) and treated them as separate proteins. This allowed us to identify “orthologous” (functionally similar) protein isoforms in different species.
The core product of QuartetS-DB is the collection of pairwise orthologs predicted between each pair of genomes from the 1621 species used in developing the database. To obtain these pairwise orthologs, we first performed an all-against-all BLAST analysis for all seven million proteins from these 1621 genomes to generate BBH pairs and predict putative orthologs. Putative pairwise orthologs are BBH pairs that satisfied two conditions: 1) the alignment region must cover at least 50% of the length of each sequence, and 2) the bit-score of the pair alignment had to exceed a cutoff value of 50, which is equivalent to a 10-5 E-value cutoff in our database. Currently, QuartetS-DB contains 4.42 million putative pairwise orthologs. We then used the QuartetS method to compute the parameter α for each putative pairwise ortholog and stored the largest value in the database. This allows users to modulate the accuracy of the retrieved pairwise orthologs between two species by specifying a cutoff value, Ω, and only retrieving those orthologs with a corresponding α ≤ Ω. Users can retrieve fewer, but more accurate (i.e., a lower level of false positives) ortholog predictions by specifying a small Ω, or they can retrieve more, but less accurate predictions by specifying a large Ω.
Comparative studies of multiple species require the identification of orthologs across more than two species, which can be inferred from pairwise orthologs by constructing orthologous groups. We constructed more than 300,000 orthologous groups by post-processing the pairwise ortholog predictions (using Ω = 20) and clustering these using the Markov Cluster (MCL) program (version 08–213, downloaded from http://micans.org/mcl). MCL is an unsupervised clustering algorithm , previously used by OrthoMCL , which clusters pairwise orthologs into orthologous groups, so that each group contains proteins (considered to be orthologs among themselves) from multiple (≥ 3) distinct species.
Functional annotation for orthologous groups
We employed a simple consensus rule based on the annotations of the individual proteins forming each group to functionally annotate the orthologous groups. When an annotation taxonomy was assigned to more than 80% of the proteins in an orthologous group, that taxonomy, which includes gene symbol, functional description, GO description, and GO accession number, was used to characterize the function of the group. Otherwise, the orthologous group was not annotated with a particular taxonomy. Many orthologous groups were not annotated either because their proteins lacked an annotation in the NCBI RefSeq database or due to inconsistent protein annotations that did not satisfy the consensus rule.
Gene tree generation for each orthologous group
We constructed gene trees for each orthologous group that contained proteins belonging to four or more species (more than 236,000 gene trees in total). First, we performed multiple sequence alignments of all proteins in a group using muscle (version 3.8.31, accessible at http://www.drive5.com/muscle/) with the default settings. We then used the results to construct gene trees with the RAxML program (version 7.0.4, accessible at http://www.exelixis-lab.org/) using the JTT model. The gene trees are saved in Newick format, so that they can be viewed and exported using the PhyloWidget plug-in , and can be downloaded as images or as Newick format text files.
QuartetS-DB also provides inparalog relationships within each one of the 1621 species in the database. These inparalogs are defined as paralogs whose duplication events occurred after the most recent speciation event that separated each species from the other species in the database. Briefly, we constructed a group of inparalogs for a species by evaluating the results of an all-against-all BLAST analysis and identifying a set of proteins in the species such that the alignment and bit-score of any two proteins in the set satisfied the two conditions described in the “Pairwise orthologs” Section and the sequence similarity between any two proteins in the set was higher than the sequence similarity between any one protein in the set and any protein in any one of the other species in QuartetS-DB. We constructed a total of 515,055 groups of inparalogs for the 1621 species. In practice, this information extends the one-to-one pairwise orthology relationships to many-to-many co-orthology relationships because inparalogs of each protein in an orthologous pair are considered to be co-orthologs of the other protein (and its inparalogs) in the pair.
Summary of the content of QuartetS-DB
Number of species
Number of proteins
Number of putative pairwise orthologs1
Number of pairwise orthologs (Cutoff Ω = 20)
Number of orthologous groups (OG)
Number of OGs with gene symbols2
Number of OGs with functional descriptions2
Number of OGs with GO annotations2
Number of gene trees (for OG with size ≥ 4)
Median size of OGs
Number of inparalog groups
At its core, QuartetS-DB contains a list of pairwise orthologous proteins, orthologous groups (consisting of orthologous proteins present in more than two species), and inparalog groups. A Web-based interface provides the ability to query the database and to access and download results, which include brief functional descriptions of individual proteins and orthologous and inparalog groups. Additional detailed information for individual proteins is also available via links to external resources, such as NCBI and UniProt. The QuartetS-DB database primarily supports two types of studies: multiple genomes (where the user is interested in identifying a list of orthologs across two or more specific genomes) and individual proteins (where the user is interested in identifying a list of species that contains orthologs of a specific protein).
Study of multiple genomes
Users may want to identify all orthologs among multiple species to conduct genome-wide comparative studies. For example, when studying human gene functions through mouse models, it is often necessary to obtain human-mouse orthologs so that knowledge about the functions of the mouse genes can be transferred to the corresponding human orthologs. As another example, the identification of orthologs among multiple bacterial genomes is usually the first step in gene association studies involving different phenotypes, such as bacterial pathogenicity and virulence.
Study of individual proteins
Users can study individual proteins by using the query function in the “Orthologous Groups” page of the interface. For example, if users would like to infer the function of a specific protein, they can query the protein identifier, such as the protein GI, RefSeq accession, Gene ID, Locus tag, or UniProt accession, to obtain an orthologous group that contains this protein and its orthologs. The annotated orthologous proteins in the group, as well as the group’s functional description (if available), can help to infer the function of the query protein (Figure 2c). The group’s pre-computed gene tree (Figure 2d) may support studies of the evolution of the query protein and its orthologs. All gene trees are displayed using the PhyloWidget plug-in, and can be downloaded as images or as Newick format text files.
In addition to the ability to query a specific protein, users can directly query orthologous groups by their annotations, such as group symbol, functional description, and GO terms. This allows users to find proteins with particular functions and identify species that contain these proteins. For example, if users select two species, Homo sapiens and Mus musculus, and query groups that contain “telomerase” in the description, they will obtain group QTS_86077 that is described as “telomerase reverse transcriptase” and contains orthologs in the two species. But if they query “dnaC,” no orthologous group is returned, because dnaC is a bacterial protein that has no orthologs in either one of the two species. This feature of QuartetS-DB may be helpful in the selection of pathogen-specific drug targets that do not have human orthologs, so that potential detrimental effects in patients can be minimized.
By providing access to orthology predictions among 1621 complete genomes based on a recently developed, accurate, and computationally efficient method (QuartetS), QuartetS-DB constitutes one of the largest and reliable resources of its kind. Because its orthology predictions are underpinned by evolutionary evidence obtained from sequenced genomes, we expect its accuracy to continue to increase in future releases as the complete genomes from additional species are sequenced.
QuartetS-DB was primarily designed to support studies involving both the identification of orthologs between two specific species and among genomes of multiple species. For the former, it provides users with the ability to retrieve more than four million accurate pairwise orthologs and to select a cutoff value that modulates the balance between prediction accuracy and coverage of the retrieved pairwise orthologs, satisfying users with distinct needs. For the latter, it provides users with the ability to retrieve a list of all orthologs across multiple, user-specified genomes and to browse more than 236,000 precomputed gene trees of orthologous groups, both convenient features for comparative studies involving multiple genomes.
The database is freely accessible at https://applications.bioanalysis.org/quartetsdb. The Web interface operates with various browsers, including Safari, Internet Explorer 7 or later, Firefox 8 or later, and Google Chrome. We used the PhyloWidget Java Applet to display gene trees, which requires Sun Java 1.5 or greater. Detailed system requirements for this Java Applet can be found at the PhyloWidget Website (http://www.phylowidget.org). There are no other system requirements. We plan to update the database annually.
We acknowledge the anonymous reviewers for their constructive comments that helped improve our database.
This work was sponsored by the United States (U.S.) Department of Defense High Performance Computing Modernization Program, under the High Performance Computing Software Applications Institutes Initiative. Funding for open access charge is the same as the funding for the performed research.
- Koonin EV: Orthologs, paralogs, and evolutionary genomics. Annu Rev Genet 2005, 39: 309–338. 10.1146/annurev.genet.39.073003.114725View ArticlePubMedGoogle Scholar
- Chen F, Mackey AJ, Stoeckert CJ, Roos DS: OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res 2006, 34(Database issue):D363-D368.PubMed CentralView ArticlePubMedGoogle Scholar
- Studer RA, Robinson-Rechavi M: How confident can we be that orthologs are similar, but paralogs differ? Trends in genetics: TIG 2009, 25(5):210–216. 10.1016/j.tig.2009.03.004View ArticlePubMedGoogle Scholar
- Nehrt NL, Clark WT, Radivojac P, Hahn MW: Testing the ortholog conjecture with comparative functional genomic data from mammals. PLoS computational biology 2011, 7(6):e1002073. 10.1371/journal.pcbi.1002073PubMed CentralView ArticlePubMedGoogle Scholar
- Stein L: Genome annotation: from sequence to biology. Nat Rev Genet 2001, 2(7):493–503.View ArticlePubMedGoogle Scholar
- Taher L, McGaughey DM, Maragh S, Aneas I, Bessling SL, Miller W, Nobrega MA, McCallion AS, Ovcharenko I: Genome-wide identification of conserved regulatory function in diverged sequences. Genome Res 2011, 21(7):1139–1149. 10.1101/gr.119016.110PubMed CentralView ArticlePubMedGoogle Scholar
- Korcsmaros T, Szalay MS, Rovo P, Palotai R, Fazekas D, Lenti K, Farkas IJ, Csermely P, Vellai T: Signalogs: orthology-based identification of novel signaling pathway components in three metazoans. PLoS One 2011, 6(5):e19240. 10.1371/journal.pone.0019240PubMed CentralView ArticlePubMedGoogle Scholar
- Gerdes SY, Scholle MD, D’Souza M, Bernal A, Baev MV, Farrell M, Kurnasov OV, Daugherty MD, Mseeh F, Polanuyer BM, et al.: From genetic footprinting to antimicrobial drug targets: examples in cofactor biosynthetic pathways. J Bacteriol 2002, 184(16):4555–4572. 10.1128/JB.184.16.4555-4572.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Florez AF, Park D, Bhak J, Kim BC, Kuchinsky A, Morris JH, Espinosa J, Muskus C: Protein network prediction and topological analysis in Leishmania major as a tool for drug target selection. BMC Bioinformatics 2010, 11: 484. 10.1186/1471-2105-11-484PubMed CentralView ArticlePubMedGoogle Scholar
- Ohta T: Evolution by gene duplication revisited: differentiation of regulatory elements versus proteins. Genetica 2003, 118(2–3):209–216.View ArticlePubMedGoogle Scholar
- Serres MH, Kerr AR, McCormack TJ, Riley M: Evolution by leaps: gene duplication in bacteria. Biol Direct 2009, 4: 46. 10.1186/1745-6150-4-46PubMed CentralView ArticlePubMedGoogle Scholar
- Remm M, Storm CE, Sonnhammer EL: Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol 2001, 314(5):1041–1052. 10.1006/jmbi.2000.5197View ArticlePubMedGoogle Scholar
- Kriventseva EV, Rahman N, Espinosa O, Zdobnov EM: OrthoDB: the hierarchical catalog of eukaryotic orthologs. Nucleic Acids Res 2008, 36(Database issue):D271-D275.PubMed CentralPubMedGoogle Scholar
- Waterhouse RM, Zdobnov EM, Tegenfeldt F, Li J, Kriventseva EV: OrthoDB: the hierarchical catalog of eukaryotic orthologs in 2011. Nucleic Acids Res 2011, 39(Database issue):D283-D288.PubMed CentralView ArticlePubMedGoogle Scholar
- Huerta-Cepas J, Bueno A, Dopazo J, Gabaldon T: PhylomeDB: a database for genome-wide collections of gene phylogenies. Nucleic Acids Res 2008, 36(Database issue):D491-D496.PubMed CentralPubMedGoogle Scholar
- Huerta-Cepas J, Capella-Gutierrez S, Pryszcz LP, Denisov I, Kormes D, Marcet-Houben M, Gabaldon T: PhylomeDB v3.0: an expanding repository of genome-wide collections of trees, alignments and phylogeny-based orthology and paralogy predictions. Nucleic Acids Res 2011, 39(Database issue):D556-D560.PubMed CentralView ArticlePubMedGoogle Scholar
- Jensen LJ, Julien P, Kuhn M, von Mering C, Muller J, Doerks T, Bork P: eggNOG: automated construction and annotation of orthologous groups of genes. Nucleic Acids Res 2008, 36(Database issue):D250-D254.PubMed CentralPubMedGoogle Scholar
- Muller J, Szklarczyk D, Julien P, Letunic I, Roth A, Kuhn M, Powell S, von Mering C, Doerks T, Jensen LJ, et al.: eggNOG v2.0: extending the evolutionary genealogy of genes with enhanced non-supervised orthologous groups, species and functional annotations. Nucleic Acids Res 2010, 38(Database issue):D190-D195.PubMed CentralView ArticlePubMedGoogle Scholar
- Altenhoff AM, Schneider A, Gonnet GH, Dessimoz C: OMA 2011: orthology inference among 1000 complete genomes. Nucleic Acids Res 2011, 39(Database issue):D289-D294.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu C, Zavaljevski N, Desai V, Reifman J: QuartetS: a fast and accurate algorithm for large-scale orthology detection. Nucleic Acids Res 2011, 39(13):e88. 10.1093/nar/gkr308PubMed CentralView ArticlePubMedGoogle Scholar
- Altenhoff AM, Dessimoz C: Phylogenetic and functional assessment of orthologs inference projects and methods. PLoS computational biology 2009, 5(1):e1000262. 10.1371/journal.pcbi.1000262PubMed CentralView ArticlePubMedGoogle Scholar
- Pruitt KD, Tatusova T, Klimke W, Maglott DR: NCBI Reference Sequences: current status, policy and new initiatives. Nucleic Acids Res 2009, 37(Database issue):D32-D36.PubMed CentralView ArticlePubMedGoogle Scholar
- Enright AJ, Van Dongen S, Ouzounis CA: An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res 2002, 30(7):1575–1584. 10.1093/nar/30.7.1575PubMed CentralView ArticlePubMedGoogle Scholar
- Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res 2003, 13(9):2178–2189. 10.1101/gr.1224503PubMed CentralView ArticlePubMedGoogle Scholar
- Jordan GE, Piel WH: PhyloWidget: web-based visualizations for the tree of life. Bioinformatics 2008, 24(14):1641–1642. 10.1093/bioinformatics/btn235View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.