Email updates

Keep up to date with the latest news and content from Genome Biology and BioMed Central.

Open Access Research

Interkingdom gene fusions

Yuri I Wolf, Alexey S Kondrashov and Eugene V Koonin*

Author Affiliations

National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Rockville Pike, Bethesda, MD 20894, USA

For all author emails, please log on.

Genome Biology 2000, 1:research0013-research0013.13  doi:10.1186/gb-2000-1-6-research0013

The electronic version of this article is the complete one and can be found online at: http://genomebiology.com/2000/1/6/research/0013


Received:5 June 2000
Revisions received:11 September 2000
Accepted:6 November 2000
Published:4 December 2000

© 2000 GenomeBiology.com

Abstract

Background

Genome comparisons have revealed major lateral gene transfer between the three primary kingdoms of life - Bacteria, Archaea, and Eukarya. Another important evolutionary phenomenon involves the evolutionary mobility of protein domains that form versatile multidomain architectures. We were interested in investigating the possibility of a combination of these phenomena, with an invading gene merging with a pre-existing gene in the recipient genome.

Results

Complete genomes of fifteen bacteria, four archaea and one eukaryote were searched for interkingdom gene fusions (IKFs); that is, genes coding for proteins that apparently consist of domains originating from different primary kingdoms. Phylogenetic analysis supported 37 cases of IKF, each of which includes a 'native' domain and a horizontally acquired 'alien' domain. IKFs could have evolved via lateral transfer of a gene coding for the alien domain (or a larger protein containing this domain) followed by recombination with a native gene. For several IKFs, this scenario is supported by the presence of a gene coding for a second, stand-alone version of the alien domain in the recipient genome. Among the genomes investigated, the greatest number of IKFs has been detected in Mycobacterium tuberculosis, where they are almost always accompanied by a stand-alone alien domain. For most of the IKF cases detected in other genomes, the stand-alone counterpart is missing.

Conclusions

The results of comparative genome analysis show that IKF formation is a real, but relatively rare, evolutionary phenomenon. We hypothesize that IKFs are formed primarily via the proposed two-stage mechanism, but other than in the Actinomycetes, in which IKF generation seems to be an active, ongoing process, most of the stand-alone intermediates have been eliminated, perhaps because of functional redundancy.

Background

Comparative genome analysis has revealed major lateral gene transfer between the three primary kingdoms of life, Bacteria, Archaea, and Eukarya [1,2,3,4]. The best recognized form of lateral gene flux is the transfer of numerous genes from mitochondria and chloroplasts to eukaryotic nuclear genomes [5]. Far beyond that, however, the role of lateral gene exchange, along with lineage-specific gene loss, as one of the principal factors of evolution, at least among prokaryotes, is obvious from the fact that the great majority of conserved families of orthologous genes show a 'patchy' phyletic distribution [6,7]. In many cases, such families are shared by phylogenetically distant species (for example, bacteria and archaea), while they are missing in some of the more closely related species (for example, bacteria from the same lineage). Correlations have been noticed between the preferred routes of gene transfer and the lifestyles of the organisms involved. Thus, massive gene exchange seems to have occurred between archaeal and bacterial hyperthermophiles [8,9], whereas certain parasitic bacteria, for example, chlamydia and spirochetes, appear to have acquired significantly more eukaryotic genes than free-living bacteria [10,11,12].

Another evolutionary trend that is predominant in eukaryotes, but is important also in bacteria and archaea, involves the evolutionary mobility of protein domains that combine to form variable multidomain architectures [13,14,15,16,18]. Domain fusion is one of the foundations of most forms of regulation and signal transduction in the cell. Examples include prokaryotic transcriptional regulators, most of which consist of the DNA-binding helix-turn-helix domain fused to a variety of small-molecule-binding domains [19], the two-component signal transduction system that is based on fusions of histidine kinases with sensor domains and of receiver domains with DNA-binding domains [20], and the sugar phosphotransferase (PTS) systems that include complex fusions of several enzymes [21]. In the evolution of eukaryotes, domain fusion takes the form of domain accretion, whereby proteins from complex organisms (such as animals) that are involved in various forms of regulation and signal transduction tend to accrue multiple domains that facilitate the formation of complex networks of interactions [22].

We were interested in exploring the possibility of a meeting between these two major evolutionary phenomena - lateral gene exchange and gene fusion - which would result in the formation of multidomain proteins in which different domains display distinct evolutionary provenance. In particular, we sought to identify fusions between domains originating from different primary kingdoms - Bacteria, Archaea and Eukarya - which we term interkingdom gene (domain) fusions (IKFs), and obtain clues to the pathways of IKF origin through comparative genome analysis. We show that, although IKF in general is a rare phenomenon, one bacterial lineage, the Actinomycetes, displays a significantly increased frequency of such events; we also propose a probable mechanism for IKF formation.

Results and discussion

To identify IKFs, all protein sequences encoded in the analyzed genomes were compared to the non-redundant protein database, and those proteins in which distinct parts showed the greatest similarity to homologs from different primary kingdoms were identified (see the Materials and methods section). In most cases, the reported alignments were highly statistically significant, leaving no doubt that true homologs were detected (Table 1). On the few occasions when the database search statistics in themselves were not fully convincing (for example, the OB-fold nucleic acid-binding domain in the Bacillus subtilis protein YhcN and the methyltransferase domain in the YabN protein, also from B. subtilis), the homologous relationship was validated by detection of the salient sequence motifs known to be involved in the corresponding protein functions (data not shown). Such motif analysis was performed for all analyzed domains in order not only to validate homology, but also to distinguish between active and inactivated forms of enzymes. Figure 1 shows multiple alignments of two domains involved in an IKF, illustrating the conservation of the characteristic functional motifs and the specific similarity between each of the domains of the IKF protein (in this case from Aquifex aeolicus) and their archaeal and bacterial homologs, respectively.

thumbnailFigure 1. Multiple alignments of two domains comprising an interkingdom domian fusion. Alignments of (a) the PHP-hydrolase domain [4] and (b) the pyruvate formate lyase activating enzyme domain of the IKF protein aq_2060 from A. aeolicus. The sequences of the aq_2060 domains are placed with the most similar sequences of the corresponding stand-alone enzymes, bacterial ones in the case of PHP-hydrolase and archaeal ones in the case of the pyruvate formate lyase activating enzyme. The phylogenetic trees produced form these alignments are shown in Figure 2c. The numbers in parentheses show the lengths of regions between the aligned blocks that are not shown. The consensus includes amino acid residues and residue classes that are conserved in 75% of the aligned sequences; the residue classes are as follows: h, hydrophobic; l, aliphatic; a, aromatic; s, small; u, tiny; p, polar; b, big; t, residues with high turn-forming propensity. Asterisks show the predicted active site residues; note the replacements in some of the sequences that are predicted to be inactivated versions of the respective enzymes (see text). The alignments were colored using the BOXSHADE program [30]; individual residues conserved in at least 50% of the aligned sequences are in red; residues similar to the conserved ones and groups of conserved similar residues are in blue.

Table 1. Interkingdom domain fusions and their probable origins

In several cases, the chimeric origin of a gene was obvious at a qualitative level because no homolog of the 'alien' domain with comparable sequence similarity was detected in the recipient superkingdom (Table 1, Figure 2a,b). For the rest of the candidate IKFs, phylogenetic tree analysis was performed to corroborate the origin of the invading domain by horizontal transfer; statistically significant grouping of a candidate IKF domain with homologs from the donor superkingdom provides such evidence (Figure 2c,d). The overall number of confirmed IKFs is relatively small - 37 in 21 compared genomes (about 0.1% of the genes) - compared to the total number of likely interkingdom gene transfers. For completely sequenced bacterial genomes this has been conservatively estimated as 1-2% of the genes, with a greater fraction (2-10%) detected in archaea and hyperthermophilic bacteria ([23], and K.S. Makarova, L. Aravind and E.V.K., unpublished observations). Examination of the clusters of orthologous groups (COGs) of proteins from complete genomes [6], in which multidomain proteins are split into the constituent domains if the orthologs of the latter are present as stand-alone forms in some of the genomes, shows that IKFs constitute only a small fraction of all fusions of evolutionarily mobile domains (Figure 3). Generally, the small number of identified IKFs compared to the total number of inferred horizontal transfer events and the total number of domain fusions could be compatible with a random model of domain fusion subsequent to lateral gene transfer.

thumbnailFigure 2. Examples of phylogenetic trees supporting the contribution of interkingdom horizontal gene transfer to the emergence of interkingdom domain fusions. The names of proteins from different primary kingdoms are color-coded: black, bacterial; pink, archaeal; green, eukaryotic; the domains involved in the apparent IKF are shown in red. Red circles show nodes with bootstrap support >70%, and yellow circles show nodes with 50-70% support. The bar unit corresponds to 0.1 substitutions per site (10 PAM). (a) IKF: Rv1683 (gi| 7476858) from M. tuberculosis. Fusion of a bacterial poly(3-hydroxy-butyrate) (PHB) synthase and eukaryotic very long chain acyl-CoA synthetase. Note the absence of eukaryotic homologs in the PHB synthase tree and of bacterial homologs other than the two from M. leprae in the acyl-CoA synthetase tree. (b) IKF: yeast YOL066c (gi|6324506). Fusion of a eukaryotic pesudouridylate synthetase with a bacterial pyrimidine deaminase. Note the absence of eukaryotic homologs, other than that from S. pombe, in the pyrimidine deaminase tree. (c) IKF: aq_2060 (gi|2984285) from Aquifex aeolicus. This protein is a fusion of a PHP superfamily hydrolase of apparent bacterial origin and a pyruvate formate-lyase activating enzyme of archaeal origin. (d) IKF: yeast YOL055c (gi|1419865), YPL258c (gi|2132251) and YPR121w (gi|2132289) from S. cerevisiae. Fusion of a eukaryotic phosphomethylpyrimidine kinase and a bacterial transcriptional activator. Species abbreviations: Bac.meg., Bacillus megaterium; Chr.vin., Chromatium vinosum; Thi.vi., Thiocystis violacea; Am.med., Amycolatopsis mediterranei; Coch.het., Cochliobolus heterostrophus; Dme, Drosophila melanogaster; Cel, Caenorhabditis elegans; Mus, Mus musculus; Spo, Schizosaccharomyces pombe; Ath, Arabidopsis thaliana; Strep.co., Streptomyces coelicolor; The.nea., Thermotoga neapolitana; Bac. am, Bacillus amyloliquefaciens; Shi.fl., Shigella flexneri; Hsa, Homo sapiens.

thumbnailFigure 3. Overall numbers of domain fusions estimated using the COGs and interkingdom domain fusions encoded in completely sequenced genomes. The data for estimating the overall number of domain fusions were from the current COG release [6], which does not include several bacterial and archaeal species (for example, Aeropyrum pernix and Deinococcus radiodurans) that have been analyzed in the present work (Table 1). Accordingly, the data for these genomes are not shown in the figure. Species name abbreviations: Af, Archaeoglobus fulgidus; Mj, Methanococcus jannaschii; Mth, Methanobacterium thermoautotrophicum; Ph, Pyrococcus horikoshii; Sc, Saccharomyces cerevisiae; Aa, Aquifex aeolicus; Tm, Thermotoga maritima; Ssp, Synechocystis sp.; Ec, Escherichia coli; Bs, Bacillus subtilis; Mtu, Mycobacterium tuberculosis; Hi, Haemophilus influenzae; Hp, Helicobacter pylori; Mg; Mycoplasma genitalium; Mp, Mycoplasma pneumoniae; Bb, Borrelia burgdorferi; Tp, Treponema pallidum; Ct, Chlamydia trachomatis; Cp, Chlamydophila pneumoniae; Rp, Rickettsia prowazekii.

However, the distribution of IKFs among genomes is distinctly non-random, suggesting that such a simple model may be incorrect. Specifically, 12 IKFs were detected in Mycobacterium tuberculosis and 10 were found in the yeast Saccharomyces cerevisiae, but only a small number or none was identified in each of the other bacterial and archaeal genomes (Figure 2, Table 1). The excess of IKFs in Mycobacterium is particularly notable, given that the fraction of genes horizontally transferred from archaea and eukaryotes in the mycobacterial genome is only slightly greater than that in most of the other bacteria, and considerably lower than that in the hyperthermophilic bacteria Aquifex and Thermotoga (K.S. Makarova, L. Aravind and E.V.K., unpublished observations). Similarly, whereas the overall number of domain fusions in M. tuberculosis is greater than in most other bacteria, the difference is insufficient to account for the over-representation of IKFs; furthermore, the cyanobacterium Synechocystis sp. has an even greater overall number of fusions but does not have any detectable IKFs (Figure 3). At present, we cannot provide a defendable biological explanation for the comparatively high frequency of IKF in Mycobacterium. It is tempting to interpret this trend in terms of adaptation of this bacterium to its relatively recently occupied parasitic niche, but examination of the individual IKF cases does not offer immediate clues in mycobacterial biology. The yeast IKFs clearly represent relatively recent horizontal transfers distinct from the gene influx from the mitochondria following the establishment of endosymbiosis because, under the protocol of IKF detection used here, only those alien domains were identified that have no counterparts in other eukaryotes.

Most of the IKFs are unique, but B. subtilis, M. tuberculosis and yeast each also encode families of two to three paralogous IKFs, which apparently have evolved by duplication subsequent to the respective fusion events (Table 1). Strikingly, the same IKF, the three-domain uridine kinase, is shared by Treponema pallidum and Thermotoga maritima (Table 1). Given that these two bacteria are not specifically related and that Borrelia burgdorferi, the second spirochete whose genome has been sequenced, encodes a typical bacterial uridine kinase, the presence of a common IKF in Treponema and Thermotoga cannot be realistically attributed to vertical inheritance of this gene from a common ancestor. It thus probably reflects horizontal transfer of the gene encoding the three-domain protein subsequent to its emergence in either the spirochetes or the Thermotogales.

Two evolutionary issues pertaining to IKFs need to be addressed, namely the mechanism(s) of their origin and the selective forces responsible for their preservation. From general considerations, it seems likely that IKFs have evolved via a two-step process, which involves lateral transfer of the complete gene coding for the IKF's alien portion, followed by domain fusion. This scenario rests on the assumption that the acquired foreign gene is selectively advantageous, because otherwise it would have been inactivated by mutations before recombination could take place. Under this mechanism, the alien portion of an IKF is likely to be present in the recipient genome also as a stand-alone gene. A clear-cut case of such a duplication of a horizontally transferred domain has been noticed in Chlamydia, whose genomes encode the SWI domain, implicated in chromatin condensation, both as a stand-alone protein and as the carboxy-terminal portion of topoisomerase I [10]. Apart from this case, the IKFs fall into two readily discernible classes, namely those from Mycobacterium and all the rest. M. tuberculosis (the only complete genome of an actinomycete available) possesses considerably more IKFs than any other bacterial or archaeal species (see above), and typically, the alien portions of these proteins show high level of similarity to the homologs from the donor superkingdom (eukaryotes). Most significantly, there is also, with a single exception, a stand-alone counterpart in the mycobacterial genome; in some cases, such a counterpart is seen only in a closely related species, M. leprae, and in one case, it is found in Streptomyces, a distantly related actinomycete (Table 1). In the other genomes, the IKFs are generally less similar to the apparent donor and, with a few exceptions, stand-alone versions of the alien domains are missing (Table 1). The hypothesis that seems to be most compatible with these observations is that IKFs indeed evolve via a stand-alone, horizontally transferred intermediate, but in the case of ancient IKFs, these intermediates are typically eliminated during evolution, perhaps because their function becomes redundant with the formation of the IKF. The IKFs identified in actinomycetes appear to result from relatively recent gene fusion events so that the original, stand-alone transferred genes are still present in the genome.

The IKFs include a variety of protein functions. Only some of these are well understood such as, for example, those of the bifunctional nucleotide and coenzyme metabolism enzymes that are particularly abundant in yeast (Table 1). In other cases, the function of an IKF-encoded protein could be predicted only tentatively on the basis of the functions of its constituent domains (Table 1). The selective advantage of the formation of multidomain proteins, at least as far as enzymes are involved, lies in the possibility of effective coupling of the reactions catalyzed by the different domains [16]; this may be generalized also for functional coordination of non-enzymatic domains. Fusion may result in the addition of a regulatory function to an enzymatic one. For example, it appears most likely that the RNA-binding TGS domain [24] in the uridine kinases of Treponema pallidum and Thermotoga maritima is involved in autoregulation of translation. The unusual aspect of the IKFs appears to be the compatibility of evolutionarily distant domains.

Examination of the phyletic distribution of the multidomain architectures of IKFs may help in pinpointing the evolutionary stage at which the fusion (but not necessarily the preceding horizontal gene transfer) has occurred. For example, the fusion of the SWI domain with topoisomerase belongs after the radiation of Chlamydia from other bacterial lineages, but before the radiation of Chlamydia pneumoniae and Chlamydia trachomatis (Table 1). The majority of IKFs detected in the yeast S. cerevisiae are also present in Schizosaccharomyces pombe and/or other ascomycetes (Table 1, and data not shown), but not in any other eukaryotes, and accordingly, they should have evolved at a relatively early stage of fungal evolution, but not before the fungal clade diverged from the rest of the eukaryotic crown group.

Finally, it should be noted that formation of some of the IKFs might have required more complex rearrangements of the contributing proteins than simple domain fusion. Figure 4 shows the domain architectures of proteins that contribute domains to two IKFs. In each case, a simple fusion between genes encoding the respective individual domains is insufficient to explain the emergence of the IKF. For example, the uridine kinase example mentioned above (Figure 4a) should have involved isolation of the TGS-HxxxH domains of threonyl-tRNA synthetase before or concomitantly with their fusion with the uridine kinase. The specific molecular mechanism could have involved selective duplication of the upstream portion of the threonyl-tRNA synthetase gene. Similarly, the sialic acid synthase homologous domain, which is fused to hydroxymethylpyrimidine phosphate kinase in A. pernix and pyrococci, appears to have been derived from two-domain proteins that additionally contain a helix-turn-helix DNA-binding domain (Figure 4b). These hypotheses of a complex mechanism of gene fusion involved in the emergence of IKFs are based on a limited sample of sequenced genomes. An alternative possibility is that, before the postulated horizontal transfer event, the recipient domain(s) has been encoded by a stand-alone gene; such genes that do not contain the fused alien domain may yet be discovered in newly sequenced genomes. In fact, a stand-alone version of the sialic acid synthase homologous domain is seen in Methanobacterium, although it is considerably less similar to the IKF than the version fused to the HTH domain (Figure 4b).

thumbnailFigure 4. Multidomain architectures of interkingdom fusion proteins and their homologs (examples). (a) The three-domain uridine kinase; (b) the sialic acid synthase homologous domain fused to hydroxymethylpyrimidine phosphate kinase. Domain name abbreviations: TTRS, threonyl-tRNA synthetase; UDK, uridine kinase; TGS and H3H, amino-terminal domains of TTRS; HMP-PK, hydroxymethylpyrimidine phosphate kinase; SISH, sialic acid synthase homologous domain; HTH, helix-turn-helix DNA-binding domain. Different shades represent distinct sequence families of each domain. Species name abbreviations: Tp, Treponema pallidum; Tm, Thermotoga maritima; Mth, Methanobacterium thermoautotrophicum; Mj, Methanococcus jannaschii; Ap, Aeropyrum pernix; Ph, Pyrococcus horikoshii; Pa, Pyrococcus abyssii.

The identification of IKFs underscores the complexity of the evolutionary process as revealed by comparison of multiple genomes. In and by itself, this phenomenon may not have a unique biological significance, but it reveals the overlap between two major evolutionary trends, horizontal gene transfer and protein domain rearrangement, and shows that domains, rather then entire proteins (genes), should be considered fundamental units of genetic material exchange.

Materials and methods

Protein sequences encoded in 21 complete genomes of archaea, bacteria and the yeast Saccharomyces cerevisiae were extracted from the Genome division of the Entrez retrieval system [25]. Each protein encoded in these genomes was used as the query in a comparison against the non-redundant protein sequence database (National Center for Biotechnology Information, NIH, Bethesda, USA) using the BLASTP program [26]. For each query, the set of local similarities detected by BLASTP was automatically (using a Perl script written for this purpose) screened for putative IKFs, that is situations in which the query did not have full-size homologs outside its immediate taxonomic group (for example, the Proteobacteria for Escherichia coli) and in which different regions of the query showed the greatest similarity to proteins from different primary kingdoms. The pseudocode for the script follows:

The script itself is available as an 1. The candidate IKF cases were further examined to detect situations where one or more distinct regions of the query could be classified as 'native' or 'alien' either on the basis of the lack of close homologs from the respective primary kingdom or using phylogenetic analysis. Multiple sequence alignments were generated using the ClustalW program [27], and when necessary, manually corrected to ensure the proper alignment of conserved motifs typical of the respective domains. Phylogenetic trees were constructed using the PROTDIST and FITCH programs of the PHYLIP package [28]. Trees were made separately for each domain of a putative IKF, and its mixed ancestry was considered confirmed if the affinities of the domains with different primary kingdoms were supported by bootstrap values of at least 50%. Additional iterative database searches were performed using the PSI-BLAST program [26,29] in order to predict functions of the individual domains of the identified IKFs in cases when these were not immediately clear.

Additional data

The following additional data are included with the online version of this paper: the 1 used to screen local similarities for putative IKFs.

Additional data file 1.

Format: TXT Size: 9KB Download fileOpen Data

References

  1. Doolittle WF: Lateral genomics.

    Trends Cell Biol 1999, 9:M5-M8. PubMed Abstract | Publisher Full Text OpenURL

  2. Doolittle WF: Phylogenetic classification and the universal tree.

    Science 1999, 284:2124-2129. PubMed Abstract | Publisher Full Text OpenURL

  3. 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. PubMed Abstract | Publisher Full Text OpenURL

  4. Ochman H, Lawrence JG, Groisman EA: Lateral gene transfer and the nature of bacterial innovation.

    Nature 2000, 405:299-304. PubMed Abstract | Publisher Full Text OpenURL

  5. Gray MW: Evolution of organellar genomes.

    Curr Opin Genet Dev 1999, 9:678-687. PubMed Abstract | Publisher Full Text OpenURL

  6. Tatusov RL, Galperin MY, Natale DA, Koonin EV: The COG database: a tool for genome-scale analysis of protein functions and evolution.

    Nucleic Acids Res 2000, 28:33-36. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families.

    Science 1997, 278:631-637. PubMed Abstract | Publisher Full Text OpenURL

  8. Aravind LR, Tatusov L, Wolf YI, Walker DR, Koonin EV: Evidence for massive gene exchange between archaeal and bacterial hyperthermophiles.

    Trends Genet 1998, 14:442-444. PubMed Abstract | Publisher Full Text OpenURL

  9. Nelson KE, Clayton RA, Gill SR, Gwinn ML, Dodson RJ, Haft DH, Hickey EK, Peterson JD, Nelson WC, Ketchum KA, et al.: Evidence for lateral gene transfer between Archaea and bacteria from genome sequence of Thermotoga maritima.

    Nature 1999, 399:323-329. PubMed Abstract | Publisher Full Text OpenURL

  10. Stephens RS, Kalman S, Lammel C, Fan J, Marathe R, Aravind L, Mitchell W, Olinger L, Tatusov RL, Zhao Q, et al.: Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis.

    Science 1998, 282:754-759. PubMed Abstract | Publisher Full Text OpenURL

  11. Subramanian G, Koonin EV, Aravind L: Comparative genome analysis of pathogenic spirochetes - Borrelia burgdorferi and Treponema pallidum.

    Infect Immun 2000, 68:1633-1648. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Wolf YI, Aravind L, Koonin EV: Rickettsiae and Chlamydiae: evidence of horizontal gene transfer and gene exchange.

    Trends Genet 1999, 15:173-175. PubMed Abstract | Publisher Full Text OpenURL

  13. Doolittle RF: The multiplicity of domains in proteins.

    Annu Rev Biochem 1995, 64:287-314. PubMed Abstract | Publisher Full Text OpenURL

  14. Enright AJ, Ilipoulos I, Kyrpides NC, Ouzounis CA: Protein interaction maps for complete genomes based on gene fusion events.

    Nature 1999, 402:86-90. PubMed Abstract | Publisher Full Text OpenURL

  15. Galperin MY, Koonin EV: Who is your neighbor: new computational approaches in functional genomics.

    Nat Biotechnol 2000, 18:609-613. PubMed Abstract | Publisher Full Text OpenURL

  16. Marcotte EM, Pellegrini M, Ng HL, Rice DW, Yeates TO, Eisenberg D: Detecting protein function and protein-protein interactions from genome sequences.

    Science 1999, 285:751-753. PubMed Abstract | Publisher Full Text OpenURL

  17. Marcotte EM, Pellegrini M., Thompson MJ, Yeates TO, Eisenberg D: A combined algorithm for genome-wide prediction of protein function.

    Nature 1999, 402:83-86. PubMed Abstract | Publisher Full Text OpenURL

  18. Aravind L, Koonin EV: DNA-binding proteins and evolution of transcription regulation in the archaea.

    Nucleic Acids Res 1999, 27:4658-4670. PubMed Abstract | Publisher Full Text OpenURL

  19. Grebe TW, Stock JB: The histidine protein kinase superfamily.

    Adv Microb Physiol 1999, 41:139-227. PubMed Abstract | Publisher Full Text OpenURL

  20. Saier MH Jr, Reizer J: The bacterial phosphotransferase system: new frontiers 30 years later.

    Mol Microbiol 1994, 13:755-764. PubMed Abstract | Publisher Full Text OpenURL

  21. Koonin EV, Aravind L, Kondrashov AS: The impact of comparative genomics on our understanding of evolution.

    Cell 2000, 101:573-576. PubMed Abstract | Publisher Full Text OpenURL

  22. 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. PubMed Abstract | Publisher Full Text OpenURL

  23. Wolf YI, Aravind L, Grishin NV, Koonin EV: Evolution of aminoacyl-tRNA synthetases - analysis of unique domain architectures and phylogenetic trees reveals a complex history of horizontal gene transfer events.

    Genome Res 1999, 9:689-710. PubMed Abstract | Publisher Full Text OpenURL

  24. National Center for Biotechnology Information [http://www.ncbi.nlm.nih.gov/Entrez/Genome/org.html] webcite

  25. 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. PubMed Abstract | Publisher Full Text OpenURL

  26. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.

    Nucleic Acids Res 1994, 22:4673-4680. PubMed Abstract | Publisher Full Text OpenURL

  27. Felsenstein J: Inferring phylogenies from protein sequences by parsimony, distance, and likelihood methods.

    Methods Enzymol 1996, 266:418-427. PubMed Abstract | Publisher Full Text OpenURL

  28. Altschul SF, Koonin EV: PSI-BLAST - a tool for making discoveries in sequence databases.

    Trends Biochem Sci 1998, 23:444-447. PubMed Abstract | Publisher Full Text OpenURL

  29. BOXSHADE [http://www.ch.embnet.org/software/BOX_form.html] webcite