- Research article
- Open Access
The COG database: an updated version includes eukaryotes
- Roman L Tatusov1Email author,
- Natalie D Fedorova1,
- John D Jackson1,
- Aviva R Jacobs1,
- Boris Kiryutin1,
- Eugene V Koonin1,
- Dmitri M Krylov1,
- Raja Mazumder2,
- Sergei L Mekhedov1,
- Anastasia N Nikolskaya2,
- B Sridhar Rao1,
- Sergei Smirnov1,
- Alexander V Sverdlov1,
- Sona Vasudevan1,
- Yuri I Wolf1,
- Jodie J Yin1 and
- Darren A Natale2
© Tatusov et al; licensee BioMed Central Ltd. 2003
Received: 20 May 2003
Accepted: 11 September 2003
Published: 11 September 2003
The availability of multiple, essentially complete genome sequences of prokaryotes and eukaryotes spurred both the demand and the opportunity for the construction of an evolutionary classification of genes from these genomes. Such a classification system based on orthologous relationships between genes appears to be a natural framework for comparative genomics and should facilitate both functional annotation of genomes and large-scale evolutionary studies.
We describe here a major update of the previously developed system for delineation of Clusters of Orthologous Groups of proteins (COGs) from the sequenced genomes of prokaryotes and unicellular eukaryotes and the construction of clusters of predicted orthologs for 7 eukaryotic genomes, which we named KOGs after euk aryotic o rthologous g roups. The COG collection currently consists of 138,458 proteins, which form 4873 COGs and comprise 75% of the 185,505 (predicted) proteins encoded in 66 genomes of unicellular organisms. The euk aryotic o rthologous g roups (KOGs) include proteins from 7 eukaryotic genomes: three animals (the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster and Homo sapiens), one plant, Arabidopsis thaliana, two fungi (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and the intracellular microsporidian parasite Encephalitozoon cuniculi. The current KOG set consists of 4852 clusters of orthologs, which include 59,838 proteins, or ~54% of the analyzed eukaryotic 110,655 gene products. Compared to the coverage of the prokaryotic genomes with COGs, a considerably smaller fraction of eukaryotic genes could be included into the KOGs; addition of new eukaryotic genomes is expected to result in substantial increase in the coverage of eukaryotic genomes with KOGs. Examination of the phyletic patterns of KOGs reveals a conserved core represented in all analyzed species and consisting of ~20% of the KOG set. This conserved portion of the KOG set is much greater than the ubiquitous portion of the COG set (~1% of the COGs). In part, this difference is probably due to the small number of included eukaryotic genomes, but it could also reflect the relative compactness of eukaryotes as a clade and the greater evolutionary stability of eukaryotic genomes.
The updated collection of orthologous protein sets for prokaryotes and eukaryotes is expected to be a useful platform for functional annotation of newly sequenced genomes, including those of complex eukaryotes, and genome-wide evolutionary studies.
The rapid accumulation of genome sequences is a major challenge to researchers attempting to extract the maximum functional and evolutionary information from the new genomes. To avoid informational overflow from the constant influx of new genome sequences, a comprehensive evolutionary classification of the genes from all sequenced genomes is required. Such classifications are based on two fundamental notions from evolutionary biology: orthology and paralogy, which describe the two fundamentally different types of homologous relationships between genes [1–4]. Orthologs are homologous genes derived by vertical descent from a single ancestral gene in the last common ancestor of the compared species. Paralogs, in contrast, are homologous genes, which, at some stage of evolution of the respective gene family, have evolved by duplication of an ancestral gene. The notions of orthology and paralogy are intimately linked because, if a duplication (s) occurred after the speciation event that separated the compared species, orthology becomes a relationship between sets of paralogs (co-orthologs), rather than individual genes. A classic case of the interplay between orthologous and paralogous relationships is seen in the globin family: all animal globins, including myoglobin, are paralogs, but they are all co-orthologs of the plant leghemoglobin(s) .
Deciphering orthologous and paralogous relationships among genes is critical for both the functional and the evolutionary aspects of comparative genomics [4, 5]. Orthologs typically occupy the same functional niche in different species, whereas paralogs tend to evolve toward functional diversification. Therefore, robustness of genome annotation depends on accurate identification of orthologs. Similarly, knowing which homologous genes are orthologs and which are paralogs is required for constructing evolutionary scenarios involving, along with vertical inheritance, lineage-specific gene loss and horizontal gene transfer.
In principle, identification of orthologs requires phylogenetic analysis of entire families of homologous proteins, which is expected to isolate orthologous protein sets in distinct clades [6–8]. However, on the scale of complete genomes, such analysis is both extremely labor-intensive and error-prone due to the inherent artifacts of phylogenetic tree construction. Therefore shortcuts have been developed by introducing the notion of a genome-specific best hit (BeT). A BeT is the protein in a target genome, which is most similar to a given protein from the query genome [9, 10]. The underlying premise is that orthologs are more similar to each other than they are to any other protein from the respective genomes. In multiple-genome comparisons, pairs of potential orthologs identified via BeTs can be joined to form clusters of orthologs represented in all or a subset of the analyzed genomes [9, 11]. This approach to the identification of orthologous protein sets meets with two obvious complications. Firstly, many proteins belong to lineage-specific expansions, i.e., have evolved via duplication(s) after the divergence of the compared species [12–14]. In these cases, deciphering (co)orthologous relationships can be a hard task and clusters of orthologs that include such expansions should be treated with particular caution. The second complication is caused by the fact that many proteins exist in multidomain forms encoded by a single gene in some species and as products of two or more stand-alone genes in others. In protein clustering, multidomain proteins may connect distinct clusters of orthologs resulting in artifactual lumping.
The approach to the identification of orthologous protein sets based on clustering of consistent BeTs has been implemented in the collection of Clusters of Orthologous Groups (COGs) of proteins [9, 15]. The COG construction protocol included an automatic procedure for detecting candidate sets of orthologs, manual splitting of multidomain proteins into the component domains, and subsequent manual curation and annotation. The COGs started with 6 prokaryotic genomes and one genome of a unicellular eukaryote, yeast Saccharomyces cerevisiae . Subsequent updates increased the number of prokaryotic genomes in the COGs to 43 . The procedure for COG construction required that each COG included proteins from at least three sufficiently distant species. This conservative approach notwithstanding, ~60 to ~85% of the proteins encoded in prokaryotic genomes were included in the COGs.
The COG system, which includes the COGNITOR program for adding new members to COGs (RLT, unpublished results), has become a widely used tool for computational genomics. The most important applications of the COGs are functional annotation of newly sequenced genomes [16–20] and genome-wide evolutionary analyses [21–25].
Here, we present a major update to the COGs, with over 63 sequenced prokaryotic genomes and three genomes of unicellular prokaryotes now included. Furthermore, the COG system is extended to complex, multicellular eukaryotes by constructing clusters of probable orthologs, which we named KOGs (euk aryotic o rthologous g roups) for 7 sequenced genomes of animals, fungi, microsporidia, and plants.
Results and discussion
Update of the COGs
Coverage of unicellular organisms in COGs
Number of annotated proteins
Number (and percentage) of proteins in COGs
Number of COGs that include the given species
Escherichia coli K12
Escherichia coli O157:H7
Escherichia coli O157:H7 EDL933
Neisseria meningitidis MC58
Neisseria meningitides Z2491
Helicobacter pylori 26695
Helicobacter pylori J99
Low-GC Gram-positive bacteria
Mycobacterium tuberculosis H37Rv
Mycobacterium tuberculosis CDC1551
Construction of KOGs for 7 sequenced eukaryotic genomes
Eukaryotic KOGs were constructed from annotated proteins encoded in the genomes of three animals (Homo sapiens , the fruit fly Drosophila melanogaster , and the nematode Caenorhabditis elegans) , the green plant Arabidopsis thaliana (thale cress) , two fungi (budding yeast Saccharomyces cerevisiae  and fission yeast Schizosaccharomyces pombe , and the microsporidian Encephalitozoon cuniculi ). The basic procedure for KOG construction was the same as the procedure previously employed for prokaryotic genomes (Refs. [9, 15] and see Materials and Methods). Given the abundance of multidomain architectures among eukaryotic proteins and the fact that apparent orthologs often differ in domain composition [32, 39], the protocol based on the BeT analysis was amended with domain identification using the RPS-BLAST program . Proteins assigned to a KOG by the initial KOG construction procedure were kept in that KOG without splitting them into individual domains if they shared a common core of domains. In addition, proteins, which consisted solely of widespread, "promiscuous" domains (e.g., SH2, SH3, WD40 repeats or TPR repeats) and did not show clear-cut orthologous relationships, were assigned to Fuzzy Orthologous Groups (FOGs). In addition to KOGs and FOGs, we also identified provisional clusters of orthologs represented in two genomes (TWOGs) by detecting bi-directional BeTs between proteins not included in KOGs or FOGs and assigning additional members by examination of the BLAST search outputs. Finally, lineage-specific expansions (LSEs) of paralogs among the proteins from each genome not included in KOGs, FOGs, and TWOGs were detected by using the clustering procedure described previously  accompanied by a newly developed procedure for finding tight protein clusters (BK and RLT, unpublished results). The construction of TWOGs and LSEs involved more extensive case by case evaluation than the KOG construction due to the lack of well established procedures to generate these types of clusters; nevertheless, these clusters should be considered preliminary until further validation.
Representation of the 7 analyzed eukaryotic species in KOGs
Number of annotated proteins
Number of proteins in KOGs (%)
Compared to prokaryotes, a considerably smaller fraction of eukaryotic genes could be included into KOGs (Tables 1 and 2). Thus, the apparent difference in coverage with highly conserved clusters of orthologs (C/KOGs) between prokaryotes and eukaryotes, particularly complex ones, is probably due to the relatively small number of eukaryotic genomes included in this analysis and is expected to level off with the growth of the eukaryotic genome collection. This view is compatible with the observed dependence of the KOG coverage on the number of genes (Table 1), which suggests that the KOGs are still far from saturation.
The functions of human globins and globin homologs, primarily in oxygen delivery to different tissues, at different developmental stages have been studied in great detail . In contrast, the dramatic proliferation of globin-like proteins in the nematode C. elegans, while noticed, in part, in previous work , is not well understood. To our knowledge, KOG3378 is the most complete current representation of this lineage-specific expansion of globin-like paralogs; the experimental study of these genes is expected to reveal novel aspects of invertebrate physiology.
Another notable observation comes from the analysis of the yeast members of KOG3378. A BLAST search of the non-redundant protein sequence database (NCBI, NIH, Bethesda) and examination of the domain composition of the S. cerevisiae protein YGR234w shows that this protein (named flavohemoglobin) consists of the globin domain fused to a flavodoxin reductase domain and is highly similar to a variety of oxidoreductases from several bacterial species and some lower eukaryotes (e.g., slime molds and other protists), which have the same domain composition ( and data not shown). The S. pombe flavohemoglobin belongs to the same protein family but is not the closest relative of the S. cerevisiae flavohemoglobin (data not shown). These observations strongly suggest that the yeast flavohemoglobin genes have been acquired from bacteria via horizontal gene transfer and hence have an evolutionary history that is distinct and independent from those of plant and animal globins. Notably, the second member of this KOG from S. cerevisie YNL234w is not at all a close paralog of the flavohemoglobins. The only identifiable domain in this large protein is the globin domain, which is most similar to vertebrate neuroglobins. These observations illustrate an important general point to be kept in mind when perusing the KOGs: although a given set of proteins may have been legitimately brought together in the same KOG in the context of eukaryotic genome comparison, on some occasions, different KOG members have different evolutionary trajectories.
Prokaryotic and eukaryotic orthologous gene sets: evolutionary connections and functional differences
The two sets of orthologous genes overlap because the three species of unicellular eukaryotes were included in both sets; the proteins from these species obviously form connections between prokaryotic orthologous sets (COGs) and eukaryotic orthologous sets (KOGs). Such connections, suggestive of orthologous relationships, were established between 1253 COGs, each of which included at least one protein from a unicellular eukaryote (not counting COGs that consisted exclusively of eukaryotic proteins), and 2000 eukaryotic KOGs. The greater number of eukaryotic KOGs involved in this relationship is due to the fact that, on many occasions, several proteins from unicellular eukaryotes that are part of the same COG have their distinct orthologs in other eukaryotes and, accordingly, belong to several KOGs. Only relatively small fractions of the prokaryotic COGs (27% of the COGs that include at least one prokaryotic species) and eukaryotic KOGs (34% of the KOGs and TWOGs) comprised sets of putative orthologs represented in both prokaryotes and eukaryotes. This emphasizes the distinction between the repertoires of genes that are conserved in prokaryotes and in eukaryotes and the considerable amount of innovation in both groups of organisms. However, these numbers give the low bound of the shared clusters of orthologs because some of the KOGs are not represented in the relatively small genomes of unicellular eukaryotes, primarily due to gene loss in the latter, but have prokaryotic counterparts.
Using phyletic patterns to examine gene function and evolution
The second case in point that we consider here is a search for COGs, which are represented in the causative agent of plague, Yersinia pestis , but not in other Proteobacteria (the taxon to which Y. pestis belongs) or eukaryotes; this query retrieves 7 COGs (Fig. 5b). These genes probably have been acquired by Y. pestis via horizontal gene transfer. On a more practical note, some of these genes could be potential targets for highly selective anti-bacterial agents. It is noticeable that three of these genes are predicted to be involved in cell wall metabolism (COGs 2152, 2401, and 3867), whereas the functions of others remain uncharacterized.
The collection of COGs from prokaryotes and unicellular eukaryotes was substantially amended to include 66 species and eukaryotic orthologous groups (KOGs) for 7 species were constructed. The prokaryotic COG system already covers most of the globular proteins encoded in bacterial and archaeal genomes. Eukaryotic KOGs include a lower fraction of the encoded proteins but this difference is expected to level off with the growth of the eukaryotic genome collection. The eukaryotic KOG analysis revealed a substantial conserved core of eukaryotic genes as well as major lineage-specific variations. Lineage-specific expansion of paralogous families within the KOGs and expansion of families that do not have orthologs in other compared genomes make major contributions to the eukaryotic gene repertoire. Only a minority of eukaryotic KOGs have readily detectable prokaryotic counterparts and the same holds for prokaryotic COGs, emphasizing the extent of innovation in both the eukaryotic and prokaryotic divisions of life. The wide scatter of the phyletic patterns among the KOGs testifies to the importance of lineage-specific gene loss in the evolution of eukaryotic genomes.
The current collection of eukaryotic KOGs includes 7 genomes whose sequences had been available as of July 1, 2002. Manual correction and annotation of KOGs is a labor-intensive process, which precluded immediate inclusion of the genomes of the mouse , fugu fish , mosquito , the urochordate Ciona instestinalis , and the malarial parasite Plasmodium falciparum , which have become available since that date. However, once the basic system is established, it is expected that inclusion of these and other newly sequenced genomes in the KOG system proceeds at a greater pace.
The C/KOG system can be employed for functional annotation of genes from new genomes by using the COGNITOR program and for research into genome evolution. The utility of the system for both of these purposes should increase progressively with the inclusion of new genomes, particularly those of early-branching eukaryotes.
Protein sets for new genomes
The protein sets for all newly included bacterial and archaeal genomes, the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, the microsporidian Encaephalitozoon cuniculi, the thale cress Arabidopsis thaliana, and the fruit fly Drosophila melanogaster were extracted from the Genome division of the (NCBI, NIH, Bethesda). The protein sequences for the nematode Caenorhabditis elegans were from the WormPep67 database, the sequences for Homo sapiens were from the NCBI build 30.
Addition of new genomes to the COGs
The new genomes were added to the COGs by using the COGNITOR program, with the results validated manually, essentially as described previously [9, 15]. After the completion of the validation process, the remaining proteins were subject to the COG construction procedure, in order to detect new COGs that could not be formed without the added genomes; the validation and annotation steps were repeated with the newly detected COGs.
Sequence analysis, construction and annotation of KOGs
The construction of KOGs followed the previously outlined strategy based on sets of consistent BeTs [9, 15], but included additional steps that reflected specific features of eukaryotic proteins. Briefly, the procedure was as follows. 1. Detection and masking of widespread, typically repetitive domains, which was performed by using the RPS-BLAST program and the PSSMs for the respective domains from the CDD collection . These domains, namely, PPR (pfam01535), WD40 (pfam00400), IG (pfam00047), IGc1, Igv, IG_like, RRM (pfam00076), ANK (pfam00023), myosin tail (pfam01576), Fn3 (pfam00041), CA, (IG), ANK, kelch (pfam01344), OAD_kelch, SH3 (pfam00018), intermediate filaments (pfam00038), C2H2 finger (pfam00096), PDZ (pfam00595), POZ (pfam00651), PH (pfam00169), ZnF-C4 (pfam00105), spectrin (pfam00435), Sushi (pfam00084), TPR (pfam00017), BTB, LRR_CC, LY, ARM, SH2, and CH, were detected and masked prior to applying the COG construction procedure. Masking these domains was required to ensure the robust classification of the eukaryotic orthologous clusters with the KOG detection procedure because hits between these common, "promiscuous" domains resulted in spurious lumping of numerous non-orthologous proteins. 2. All-against-all comparison of protein sequences from the analyzed genomes by using the gapped BLAST program , with filtering for low sequence complexity regions performed using the SEG program . 3. Detection of triangles of mutually consistent, genome-specific best hits (BeTs). 4. Merging triangles with a common side to form crude, preliminary KOGs. 5. Case by case analysis of each candidate KOG. This analysis serves to eliminate the false-positives that are incorporated in the KOGs during the automatic steps and included, primarily, examination of the domain composition of KOG members, which was determined using the RPS-BLAST program and the CDD collection of position-specific scoring matrices (PSSMs) for individual domains . Generally, proteins were kept in the same KOG when they shared a conserved core domain architecture. However, in cases when KOGs were artificially bridged by multidomain proteins, the latter were split into individual domains (or arrays of domains) and steps (1)-(4) were repeated with these sequences; this results in the assignment of individual domains to KOGs in accordance with their distinct evolutionary affinities. 6. Assignment of proteins containing promiscuous domains. In cases when a sequence assigned to a KOG contained one or more masked promiscuous domains, these domains were restored and became part of the respective KOG. Proteins containing promiscuous domains but not assigned to any KOG were classified in Fuzzy Orthologous Groups (FOGs) named after the respective domains. 7. Examination of large KOGs, which included multiple members from all or several of the compared genomes by using phylogenetic trees, cluster analysis with the BLASTCLUST program ftp://ftp.ncbi.nih.gov/blast/, comparison of domain architectures, and visual inspection of alignments; as a result, some of these protein sets were split into two or more smaller ones that were included in the final set of KOGs.
The KOGs were annotated on the basis of the annotations available through GenBank and other public databases, which were critically assessed against the primary literature. For proteins that are currently annotated as "hypothetical" or "unknown", iterative sequence similarity searches with the PSI-BLAST program , the results of the RPS-BLAST searches, additional domain architecture analysis performed by using the SMART system , and comparison to the COG database by using the COGNITOR program (RLT, unpublished results) were employed to identify distant homologs with experimentally characterized functions and/or structures. The known and predicted functions of KOGs were classified into 23 categories (see legend to Fig. 4); these were modified from the functional classification previously employed for prokaryotic COGs  by including several specific eukaryotic categories.
Availability of the results
We thank L. Aravind, David Lipman, Kira Makarova and Wei Yang for useful discussions, and Igor Garkavtsev for his contributions at the initial stages of the KOG project.
- Fitch WM: Distinguishing homologous from analogous proteins. Systematic Zoology 1970, 19: 99–106.View ArticlePubMedGoogle Scholar
- Fitch WM: Homology a personal view on some of the problems. Trends Genet 2000, 16: 227–231. 10.1016/S0168-9525(00)02005-9View ArticlePubMedGoogle Scholar
- Henikoff S, Greene EA, Pietrokovski S, Bork P, Attwood TK, Hood L: Gene families: the taxonomy of protein paralogs and chimeras. Science 1997, 278: 609–614. 10.1126/science.278.5338.609View ArticlePubMedGoogle Scholar
- Sonnhammer EL, Koonin EV: Orthology, paralogy and proposed classification for paralog subtypes. Trends Genet 2002, 18: 619–620. 10.1016/S0168-9525(02)02793-2View ArticlePubMedGoogle Scholar
- Wilson CA, Kreychman J, Gerstein M: Assessing annotation transfer for genomics: quantifying the relations between protein sequence, structure and function through traditional and probabilistic scores. J Mol Biol 2000, 297: 233–249. 10.1006/jmbi.2000.3550View ArticlePubMedGoogle Scholar
- Sicheritz-Ponten T, Andersson SG: A phylogenomic approach to microbial evolution. Nucleic Acids Res 2001, 29: 545–552. 10.1093/nar/29.2.545PubMed CentralView ArticlePubMedGoogle Scholar
- Zmasek CM, Eddy SR: RIO: Analyzing proteomes by automated phylogenomics using resampled inference of orthologs. BMC Bioinformatics 2002, 3: 14. 10.1186/1471-2105-3-14PubMed CentralView ArticlePubMedGoogle Scholar
- Storm CE, Sonnhammer EL: Automated ortholog inference from phylogenetic trees and calculation of orthology reliability. Bioinformatics 2002, 18: 92–99. 10.1093/bioinformatics/18.1.92View ArticlePubMedGoogle Scholar
- Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families. Science 1997, 278: 631–637. 10.1126/science.278.5338.631View ArticlePubMedGoogle Scholar
- Huynen MA, Bork P: Measuring genome evolution. Proc Natl Acad Sci U S A 1998, 95: 5849–5856. 10.1073/pnas.95.11.5849PubMed CentralView ArticlePubMedGoogle Scholar
- Montague MG, Hutchison CA 3rd: Gene content phylogeny of herpesviruses. Proc Natl Acad Sci U S A 2000, 97: 5334–5339. 10.1073/pnas.97.10.5334PubMed CentralView ArticlePubMedGoogle Scholar
- Jordan IK, Makarova KS, Spouge JL, Wolf YI, Koonin EV: Lineage-specific gene expansions in bacterial and archaeal genomes. Genome Res 2001, 11: 555–565. 10.1101/gr.GR-1660RPubMed 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: 1041–1052. 10.1006/jmbi.2000.5197View ArticlePubMedGoogle Scholar
- Lespinet O, Wolf YI, Koonin EV, Aravind L: The role of lineage-specific gene family expansion in the evolution of eukaryotes. Genome Res 2002, 12: 1048–1059. 10.1101/gr.174302PubMed CentralView ArticlePubMedGoogle Scholar
- Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, Kiryutin B, Galperin MY, Fedorova ND, Koonin EV: The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res 2001, 29: 22–28. 10.1093/nar/29.1.22PubMed CentralView ArticlePubMedGoogle Scholar
- Slesarev AI, Mezhevaya KV, Makarova KS, Polushin NN, Shcherbinina OV, Shakhova VV, Belova GI, Aravind L, Natale DA, Rogozin IB, Tatusov RL, Wolf YI, Stetter KO, Malykh AG, Koonin EV, Kozyavkin SA: The complete genome of hyperthermophile Methanopyrus kandleri AV19 and monophyly of archaeal methanogens. Proc Natl Acad Sci U S A 2002, 99: 4644–4649. 10.1073/pnas.032671499PubMed CentralView ArticlePubMedGoogle Scholar
- Natale DA, Shankavaram UT, Galperin MY, Wolf YI, Aravind L, Koonin EV: Genome annotation using clusters of orthologous groups of proteins (COGs) – towards understanding the first genome of a Crenarchaeon. Genome Biology 2001, 5: RESEARCH0009.Google Scholar
- Nolling J, Breton G, Omelchenko MV, Makarova KS, Zeng Q, Gibson R, Lee HM, Dubois J, Qiu D, Hitti J, Wolf YI, Tatusov RL, Sabathe F, Doucette-Stamm L, Soucaille P, Daly MJ, Bennett GN, Koonin EV, Smith DR: Genome sequence and comparative analysis of the solvent-producing bacterium Clostridium acetobutylicum. J Bacteriol 2001, 183: 4823–4838. 10.1128/JB.183.16.4823-4838.2001PubMed CentralView ArticlePubMedGoogle Scholar
- McClelland M, Sanderson KE, Spieth J, Clifton SW, Latreille P, Courtney L, Porwollik S, Ali J, Dante M, Du F, Hou S, Layman D, Leonard S, Nguyen C, Scott K, Holmes A, Grewal N, Mulvaney E, Ryan E, Sun H, Florea L, Miller W, Stoneking T, Nhan M, Waterston R, Wilson RK: Complete genome sequence of Salmonella enterica serovar Typhimurium LT2. Nature 2001, 413: 852–856. 10.1038/35101614View ArticlePubMedGoogle Scholar
- Makarova KS, Aravind L, Wolf YI, Tatusov RL, Minton KW, Koonin EV, Daly MJ: Genome of the extremely radiation-resistant bacterium Deinococcus radiodurans viewed from the perspective of comparative genomics. Microbiol Mol Biol Rev 2001, 65: 44–79. 10.1128/MMBR.65.1.44-79.2001PubMed CentralView ArticlePubMedGoogle Scholar
- Jordan IK, Kondrashov FA, Rogozin IB, Tatusov RL, Wolf YI, Koonin EV: Constant relative rate of protein evolution and detection of functional diversification among bacterial, archaeal and eukaryotic proteins. Genome Biol 2001, 2: RESEARCH0053. 10.1186/gb-2001-2-12-research0053PubMed CentralView ArticlePubMedGoogle Scholar
- Yanai I, Derti A, DeLisi C: Genes linked by fusion events are generally of the same functional category: A systematic analysis of 30 microbial genomes. Proc Natl Acad Sci U S A 2001, 98: 7940–7945. 10.1073/pnas.141236298PubMed CentralView ArticlePubMedGoogle Scholar
- Lecompte O, Ripp R, Puzos-Barbe V, Duprat S, Heilig R, Dietrich J, Thierry JC, Poch O: Genome evolution at the genus level: comparison of three complete genomes of hyperthermophilic archaea. Genome Res 2001, 11: 981–993. 10.1101/gr.GR1653RPubMed CentralView ArticlePubMedGoogle Scholar
- Koonin EV, Makarova KS, Aravind L: Horizontal gene transfer in prokaryotes: quantification and classification. Annu Rev Microbiol 2001, 55: 709–742. 10.1146/annurev.micro.55.1.709View ArticlePubMedGoogle Scholar
- Jordan IK, Rogozin IB, Wolf YI, Koonin EV: Essential genes are more evolutionarily conserved than are nonessential genes in bacteria. Genome Res 2002, 12: 962–968. 10.1101/gr.87702. Article published online before print in May 2002PubMed CentralView ArticlePubMedGoogle Scholar
- Wernegreen JJ: Genome evolution in bacterial endosymbionts of insects. Nat Rev Genet 2002, 3: 850–861. 10.1038/nrg931View ArticlePubMedGoogle Scholar
- Casjens S: Borrelia genomes in the year 2000. J Mol Microbiol Biotechnol 2000, 2: 401–410.PubMedGoogle Scholar
- Koonin EV, Mushegian AR, Galperin MY, Walker DR: Comparison of archaeal and bacterial genomes: computer analysis of protein sequences predicts novel functions and suggests a chimeric origin for the archaea. Mol Microbiol 1997, 25: 619–637. 10.1046/j.1365-2958.1997.4821861.xView ArticlePubMedGoogle Scholar
- Snel B, Bork P, Huynen MA: Genomes in flux: the evolution of archaeal and proteobacterial gene content. Genome Res 2002, 12: 17–25. 10.1101/gr.176501View ArticlePubMedGoogle Scholar
- Wolf YI, Rogozin IB, Grishin NV, Tatusov RL, Koonin EV: Genome trees constructed using five different approaches suggest new major bacterial clades. BMC Evol Biol 2001, 1: 8. 10.1186/1471-2148-1-8PubMed CentralView ArticlePubMedGoogle Scholar
- Mirkin BG, Fenner TI, Galperin MY, Koonin EV: Algorithms for computing parsimonious evolutionary scenarios for genome evolution, the last universal common ancestor and dominance of horizontal gene transfer in the evolution of prokaryotes. BMC Evol Biol 2003, 3: 2. 10.1186/1471-2148-3-2PubMed CentralView ArticlePubMedGoogle Scholar
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, Funke R, Gage D, Harris K, Heaford A, Howland J, Kann L, Lehoczky J, LeVine R, McEwan P, McKernan K, Meldrim J, Mesirov JP, Miranda C, Morris W, Naylor J, Raymond C, Rosetti M, Santos R, Sheridan A, Sougnez C, Stange-Thomann N, Stojanovic N, Subramanian A, Wyman D, Rogers J, Sulston J, Ainscough R, Beck S, Bentley D, Burton J, Clee C, Carter N, Coulson A, Deadman R, Deloukas P, Dunham A, Dunham I, Durbin R, French L, Grafham D, Gregory S, Hubbard T, Humphray S, Hunt A, Jones M, Lloyd C, McMurray A, Matthews L, Mercer S, Milne S, Mullikin JC, Mungall A, Plumb R, Ross M, Shownkeen R, Sims S, Waterston RH, Wilson RK, Hillier LW, McPherson JD, Marra MA, Mardis ER, Fulton LA, Chinwalla AT, Pepin KH, Gish WR, Chissoe SL, Wendl MC, Delehaunty KD, Miner TL, Delehaunty A, Kramer JB, Cook LL, Fulton RS, Johnson DL, Minx PJ, Clifton SW, Hawkins T, Branscomb E, Predki P, Richardson P, Wenning S, Slezak T, Doggett N, Cheng JF, Olsen A, Lucas S, Elkin C, Uberbacher E, Frazier M, et al.: Initial sequencing and analysis of the human genome. Nature 2001, 409: 860–921. 10.1038/35057062View ArticlePubMedGoogle Scholar
- Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG, Scherer SE, Li PW, Hoskins RA, Galle RF, et al.: The genome sequence of Drosophila melanogaster. Science 2000, 287: 2185–2195. 10.1126/science.287.5461.2185View ArticlePubMedGoogle Scholar
- Consortium TCeS: Genome sequence of the nematode C. elegans: a platform for investigating biology. Science 1998, 282: 2012–2018. 10.1126/science.282.5396.2012View ArticleGoogle Scholar
- Initiative: TAG Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 2000, 408: 796–815. 10.1038/35048692View ArticleGoogle Scholar
- Goffeau A, Barrell BG, Bussey H, Davis RW, Dujon B, Feldmann H, Galibert F, Hoheisel JD, Jacq C, Johnston M, Louis EJ, Mewes HW, Murakami Y, Philippsen P, Tettelin H, Oliver SG: Life with 6000 genes. Science 1996, 274: 563–567. 10.1126/science.274.5287.546View ArticleGoogle Scholar
- Wood V, Gwilliam R, Rajandream MA, Lyne M, Lyne R, Stewart A, Sgouros J, Peat N, Hayles J, Baker S, Basham D, Bowman S, Brooks K, Brown D, Brown S, Chillingworth T, Churcher C, Collins M, Connor R, Cronin A, Davis P, Feltwell T, Fraser A, Gentles S, Goble A, Hamlin N, Harris D, Hidalgo J, Hodgson G, Holroyd S, Hornsby T, Howarth S, Huckle EJ, Hunt S, Jagels K, James K, Jones L, Jones M, Leather S, McDonald S, McLean J, Mooney P, Moule S, Mungall K, Murphy L, Niblett D, Odell C, Oliver K, O'Neil S, Pearson D, Quail MA, Rabbinowitsch E, Rutherford K, Rutter S, Saunders D, Seeger K, Sharp S, Skelton J, Simmonds M, Squares R, Squares S, Stevens K, Taylor K, Taylor RG, Tivey A, Walsh S, Warren T, Whitehead S, Woodward J, Volckaert G, Aert R, Robben J, Grymonprez B, Weltjens I, Vanstreels E, Rieger M, Schafer M, Muller-Auer S, Gabel C, Fuchs M, Dusterhoft A, Fritzc C, Holzer E, Moestl D, Hilbert H, Borzym K, Langer I, Beck A, Lehrach H, Reinhardt R, Pohl TM, Eger P, Zimmermann W, Wedler H, Wambutt R, Purnelle B, Goffeau A, Cadieu E, Dreano S, Gloux S, et al.: The genome sequence of Schizosaccharomyces pombe. Nature 2002, 415: 871–880. 10.1038/nature724View ArticlePubMedGoogle Scholar
- Katinka MD, Duprat S, Cornillot E, Metenier G, Thomarat F, Prensier G, Barbe V, Peyretaillade E, Brottier P, Wincker P, Delbac F, El Alaoui H, Peyret P, Saurin W, Gouy M, Weissenbach J, Vivares CP: Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi. Nature 2001, 414: 450–453. 10.1038/35106579View ArticlePubMedGoogle Scholar
- Koonin EV, Aravind L, Kondrashov AS: The impact of comparative genomics on our understanding of evolution. Cell 2000, 101: 573–576.View ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Anderson JB, DeWeese-Scott C, Fedorova ND, Geer LY, He S, Hurwitz DI, Jackson JD, Jacobs AR, Lanczycki CJ, Liebert CA, Liu C, Madej T, Marchler GH, Mazumder R, Nikolskaya AN, Panchenko AR, Rao BS, Shoemaker BA, Simonyan V, Song JS, Thiessen PA, Vasudevan S, Wang Y, Yamashita RA, Yin JJ, Bryant SH: CDD: a curated Entrez database of conserved domain alignments. Nucleic Acids Res 2003, 31: 383–387. 10.1093/nar/gkg087PubMed CentralView ArticlePubMedGoogle Scholar
- Makarova KS, Aravind L, Galperin MY, Grishin NV, Tatusov RL, Wolf YI, Koonin EV: Comparative genomics of the Archaea (Euryarchaeota): evolution of conserved protein families, the stable core, and the variable shell. Genome Res 1999, 9: 608–628.PubMedGoogle Scholar
- Hedges SB: The origin and evolution of model organisms. Nat Rev Genet 2002, 3: 838–849. 10.1038/nrg929View ArticlePubMedGoogle Scholar
- Pesce A, Bolognesi M, Bocedi A, Ascenzi P, Dewilde S, Moens L, Hankeln T, Burmester T: Neuroglobin and cytoglobin. Fresh blood for the vertebrate globin family. EMBO Rep 2002, 3: 1146–1151. 10.1093/embo-reports/kvf248PubMed CentralView ArticlePubMedGoogle Scholar
- Neuwald AF, Liu JS, Lipman DJ, Lawrence CE: Extracting protein alignment models from the sequence database. Nucleic Acids Res 1997, 25: 1665–1677. 10.1093/nar/25.9.1665PubMed CentralView ArticlePubMedGoogle Scholar
- Gardner PR, Gardner AM, Martin LA, Dou Y, Li T, Olson JS, Zhu H, Riggs AF: Nitric-oxide dioxygenase activity and function of flavohemoglobins. sensitivity to nitric oxide and carbon monoxide inhibition. J Biol Chem 2000, 275: 31581–31587. 10.1074/jbc.M004141200View ArticlePubMedGoogle Scholar
- Huynen MJ, Snel B: Gene and context: integrative approaches to genome analysis. Adv Prot Chem 2000, 54: 345–379.View ArticleGoogle Scholar
- Huynen M, Snel B, Lathe W, Bork P: Exploitation of gene context. Curr Opin Struct Biol 2000, 10: 366–370. 10.1016/S0959-440X(00)00098-1View ArticlePubMedGoogle Scholar
- Wolf YI, Rogozin IB, Kondrashov AS, Koonin EV: Genome alignment, evolution of prokaryotic genome organization and prediction of gene function using genomic context. Genome Res 2001, 11: 356–372. 10.1101/gr.GR-1619RView ArticlePubMedGoogle Scholar
- Olsen GJ, Woese CR: Archaeal genomics: an overview. Cell 1997, 89: 991–994.View ArticlePubMedGoogle Scholar
- Winkler HH, Neuhaus HE: Non-mitochondrial ATP transport. Trends Biochem Sci 1999, 24: 64–68. 10.1016/S0968-0004(98)01334-6View ArticlePubMedGoogle Scholar
- Wolf YI, Aravind L, Koonin EV: Rickettsiae and Chlamydiae: evidence of horizontal gene transfer and gene exchange. Trends Genet 1999, 15: 173–175. 10.1016/S0168-9525(99)01704-7View ArticlePubMedGoogle Scholar
- Parkhill J, Wren BW, Thomson NR, Titball RW, Holden MT, Prentice MB, Sebaihia M, James KD, Churcher C, Mungall KL, Baker S, Basham D, Bentley SD, Brooks K, Cerdeno-Tarraga AM, Chillingworth T, Cronin A, Davies RM, Davis P, Dougan G, Feltwell T, Hamlin N, Holroyd S, Jagels K, Karlyshev AV, Leather S, Moule S, Oyston PC, Quail M, Rutherford K, Simmonds M, Skelton J, Stevens K, Whitehead S, Barrell BG: Genome sequence of Yersinia pestis, the causative agent of plague. Nature 2001, 413: 523–527. 10.1038/35097083View ArticlePubMedGoogle Scholar
- Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, Antonarakis SE, Attwood J, Baertsch R, Bailey J, Barlow K, Beck S, Berry E, Birren B, Bloom T, Bork P, Botcherby M, Bray N, Brent MR, Brown DG, Brown SD, Bult C, Burton J, Butler J, Campbell RD, Carninci P, Cawley S, Chiaromonte F, Chinwalla AT, Church DM, Clamp M, Clee C, Collins FS, Cook LL, Copley RR, Coulson A, Couronne O, Cuff J, Curwen V, Cutts T, Daly M, David R, Davies J, Delehaunty KD, Deri J, Dermitzakis ET, Dewey C, Dickens NJ, Diekhans M, Dodge S, Dubchak I, Dunn DM, Eddy SR, Elnitski L, Emes RD, Eswara P, Eyras E, Felsenfeld A, Fewell GA, Flicek P, Foley K, Frankel WN, Fulton LA, Fulton RS, Furey TS, Gage D, Gibbs RA, Glusman G, Gnerre S, Goldman N, Goodstadt L, Grafham D, Graves TA, Green ED, Gregory S, Guigo R, Guyer M, Hardison RC, Haussler D, Hayashizaki Y, Hillier LW, Hinrichs A, Hlavina W, Holzer T, Hsu F, Hua A, Hubbard T, Hunt A, Jackson I, Jaffe DB, Johnson LS, Jones M, Jones TA, Joy A, Kamal M, Karlsson EK, et al.: Initial sequencing and comparative analysis of the mouse genome. Nature 2002, 420: 520–562. 10.1038/nature01262View ArticlePubMedGoogle Scholar
- Aparicio S, Chapman J, Stupka E, Putnam N, Chia JM, Dehal P, Christoffels A, Rash S, Hoon S, Smit A, Gelpke MD, Roach J, Oh T, Ho IY, Wong M, Detter C, Verhoef F, Predki P, Tay A, Lucas S, Richardson P, Smith SF, Clark MS, Edwards YJ, Doggett N, Zharkikh A, Tavtigian SV, Pruss D, Barnstead M, Evans C, Baden H, Powell J, Glusman G, Rowen L, Hood L, Tan YH, Elgar G, Hawkins T, Venkatesh B, Rokhsar D, Brenner S: Whole-genome shotgun assembly and analysis of the genome of Fugu rubripes. Science 2002, 297: 1301–1310. 10.1126/science.1072104View ArticlePubMedGoogle Scholar
- Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, Salzberg SL, Loftus B, Yandell M, Majoros WH, Rusch DB, Lai Z, Kraft CL, Abril JF, Anthouard V, Arensburger P, Atkinson PW, Baden H, de Berardinis V, Baldwin D, Benes V, Biedler J, Blass C, Bolanos R, Boscus D, Barnstead M, Cai S, Center A, Chatuverdi K, Christophides GK, Chrystal MA, Clamp M, Cravchik A, Curwen V, Dana A, Delcher A, Dew I, Evans CA, Flanigan M, Grundschober-Freimoser A, Friedli L, Gu Z, Guan P, Guigo R, Hillenmeyer ME, Hladun SL, Hogan JR, Hong YS, Hoover J, Jaillon O, Ke Z, Kodira C, Kokoza E, Koutsos A, Letunic I, Levitsky A, Liang Y, Lin JJ, Lobo NF, Lopez JR, Malek JA, McIntosh TC, Meister S, Miller J, Mobarry C, Mongin E, Murphy SD, O'Brochta DA, Pfannkoch C, Qi R, Regier MA, Remington K, Shao H, Sharakhova MV, Sitter CD, Shetty J, Smith TJ, Strong R, Sun J, Thomasova D, Ton LQ, Topalis P, Tu Z, Unger MF, Walenz B, Wang A, Wang J, Wang M, Wang X, Woodford KJ, Wortman JR, Wu M, Yao A, Zdobnov EM, Zhang H, Zhao Q, et al.: The genome sequence of the malaria mosquito Anopheles gambiae. Science 2002, 298: 129–149. 10.1126/science.1076181View ArticlePubMedGoogle Scholar
- Dehal P, Satou Y, Campbell RK, Chapman J, Degnan B, De Tomaso A, Davidson B, Di Gregorio A, Gelpke M, Goodstein DM, Harafuji N, Hastings KE, Ho I, Hotta K, Huang W, Kawashima T, Lemaire P, Martinez D, Meinertzhagen IA, Necula S, Nonaka M, Putnam N, Rash S, Saiga H, Satake M, Terry A, Yamada L, Wang HG, Awazu S, Azumi K, Boore J, Branno M, Chin-Bow S, DeSantis R, Doyle S, Francino P, Keys DN, Haga S, Hayashi H, Hino K, Imai KS, Inaba K, Kano S, Kobayashi K, Kobayashi M, Lee BI, Makabe KW, Manohar C, Matassi G, Medina M, Mochizuki Y, Mount S, Morishita T, Miura S, Nakayama A, Nishizaka S, Nomoto H, Ohta F, Oishi K, Rigoutsos I, Sano M, Sasaki A, Sasakura Y, Shoguchi E, Shin-i T, Spagnuolo A, Stainier D, Suzuki MM, Tassy O, Takatori N, Tokuoka M, Yagi K, Yoshizaki F, Wada S, Zhang C, Hyatt PD, Larimer F, Detter C, Doggett N, Glavina T, Hawkins T, Richardson P, Lucas S, Kohara Y, Levine M, Satoh N, Rokhsar DS: The draft genome of Ciona intestinalis: insights into chordate and vertebrate origins. Science 2002, 298: 2157–2167. 10.1126/science.1080049View ArticlePubMedGoogle Scholar
- Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, Carlton JM, Pain A, Nelson KE, Bowman S, Paulsen IT, James K, Eisen JA, Rutherford K, Salzberg SL, Craig A, Kyes S, Chan MS, Nene V, Shallom SJ, Suh B, Peterson J, Angiuoli S, Pertea M, Allen J, Selengut J, Haft D, Mather MW, Vaidya AB, Martin DM, Fairlamb AH, Fraunholz MJ, Roos DS, Ralph SA, McFadden GI, Cummings LM, Subramanian GM, Mungall C, Venter JC, Carucci DJ, Hoffman SL, Newbold C, Davis RW, Fraser CM, Barrell B: Genome sequence of the human malaria parasite Plasmodium falciparum. Nature 2002, 419: 498–511. 10.1038/nature01097View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25: 3389–3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMedGoogle Scholar
- Wootton JC, Federhen S: Analysis of compositionally biased regions in sequence databases. Methods Enzymol 1996, 266: 554–571.View ArticlePubMedGoogle Scholar
- Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci U S A 1998, 95: 5857–5864. 10.1073/pnas.95.11.5857PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.