![[Feedback]](/icons/banner/feedbackACT.gif){kind=icon}
![[For Subscribers]](/icons/banner/subscriberACT.gif){kind=icon}
![[Archive]](/icons/banner/archiveACT.gif){kind=icon}
![[Search]](/icons/banner/searchACT.gif){kind=icon}
![[Contents]](/icons/banner/tocACT.gif){kind=icon}
|
|
|
|
|
||
|
Published online before print
June 17, 2005, 10.1101/gr.3666505 Genome Res. 15:954-959, 2005 ©2005 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/05 $5.00
Letter The net of life: Reconstructing the microbial phylogenetic networkComputational Genomics Group, The European Bioinformatics Institute, EMBL Cambridge Outstation, Cambridge CB10 1SD, United Kingdom
It has previously been suggested that the phylogeny of microbial species might be better described as a network containing vertical and horizontal gene transfer (HGT) events. Yet, all phylogenetic reconstructions so far have presented microbial trees rather than networks. Here, we present a first attempt to reconstruct such an evolutionary network, which we term the "net of life." We use available tree reconstruction methods to infer vertical inheritance, and use an ancestral state inference algorithm to map HGT events on the tree. We also describe a weighting scheme used to estimate the number of genes exchanged between pairs of organisms. We demonstrate that vertical inheritance constitutes the bulk of gene transfer on the tree of life. We term the bulk of horizontal gene flow between tree nodes as "vines," and demonstrate that multiple but mostly tiny vines interconnect the tree. Our results strongly suggest that the HGT network is a scale-free graph, a finding with important implications for genome evolution. We propose that genes might propagate extremely rapidly across microbial species through the HGT network, using certain organisms as hubs.
Following the legacy of Darwin's Origin of Species (Darwin 1859
All these tree-like representations of evolution have an inherent drawback, dealing solely with vertical inheritance (Bapteste et al. 2004
Attempts to deal with this issue include algorithmic solutions for network-like tree reconstruction, mostly addressing recombination (but not HGT) as a form of nonvertical inheritance (Wang et al. 2001 Here we present a first attempt to reconstruct the history of the microbial world, recording both horizontal and vertical gene transfer. For a scaffold depicting vertical gene transfer we use established tree reconstruction methods, on which we document the instances of horizontal transfer that intertwine the tree. We discuss the major properties of this complex phylogenetic network based on a multitude of genome comparisons; demonstrate its scale-free, small-world nature; and discuss the patterns of gene propagation through the network.
Data To ensure that our results are not affected solely by the orthology data (see Methods), we used two data sets: OFAM (see Methods) and groups of orthologs defined by STRING (von Mering et al. 2003  ). Similarly, to avoid possible bias from a single tree reconstruction method, we derived genomic trees with three independent methods: gene content, average ortholog similarity, and genome conservation (see Methods) for OFAM data and gene content for STRING data.
The summary of the evolutionary events reconstructed with each method is presented in Table 1. It is evident that although HGT is readily detectable, the bulk of the genes are still transferred by vertical gene transfer, which is the most prevailing mode of inheritance (Kunin and Ouzounis 2003a
HGT vine width distribution We define the HGT vine width as a summary of all horizontal transfer events between two nodes on the tree, subsequently fixated within the genome. The distribution of HGT vine widths, or number of genes transferred between any two nodes on the tree, is shown in Figure 1. All data sets and trees produce virtually identical frequency distribution (Fig. 1), following a power law (Table 2A), with the STRING data shifted by an order of magnitude, due to lower coverage of genomes (Table 1).
Connectivity of the network To investigate the properties of the HGT network, we removed the underlying (vertical inheritance) tree from the net of life. Since our inference of HGT vine widths is probabilistic (see Methods), we had to select a meaningful threshold to depict the inferred events. Thus, to investigate the connectivity of the HGT network, we experimented with several thresholds, namely, one (a single HGT), five, and 10. Irrespectively of the tree used and data set, the HGT network displays small-world behavior, with the diameter of the network fluctuating between five and six.
When higher thresholds are chosen for the analysis, the network also demonstrates power law distribution of connectivity of nodes (Table 2B), once again irrespectively of the data set or the tree used (Fig. 2). This power-law signal is obscured at the lowest thresholds, where many nodes appear to have high connectivity. We suggest that this deviation from the power law is a result of noise inevitable when a probability model is examined at low thresholds, namely, possibly containing more false-positive instances. Our usage of thresholds higher than one for evidence of HGT is indeed reinforced by biological observations that genes often travel between organisms as groups rather than singletons (Boucher et al. 2003
HGT champions
When examining 165 microbial genomes for the network hubs, certain species came out on the top of the connectivity list with a remarkable consistency between the results obtained from different trees and data sets (Table 3). We found Pirellula sp., Bradyrhizobium japonicum, and Erwinia carotovora always at the top of the list of (terminal) nodes with the largest number of HGT partners. Interestingly, the original genome report for Pirellula sp. provides certain hints for HGT events in this species (Glockner et al. 2003 ). Furthermore, there is evidence for HGT between B. japonicum and E. carotovora in the literature (Streit et al. 2004). In conclusion, these hubs can serve as bacterial "gene banks," providing a medium to acquire and redistribute genes in the microbial communities, caused either by specific genetic mechanisms or by virtue of their close proximity to and interaction with other species in their environmental niches.
We have also examined HGT vines that are reported to be wide and consistent across data sets and trees. One of the widest HGT vines is observed between the Bradyrhizobium genus (or sometimes the broader Rhizobiales group) of Alpha Proteobacteria and the Beta Proteobacterium Ralstonia solanacearum. Phylogenetically distant, both these species are soil bacteria, penetrating plant roots and forming—symbiotic in case of Bradyrhizobium (Kiers et al. 2003 ) and parasitic in case of Ralstonia (Alfano and Collmer 2004; Genin and Boucher 2004)—relationships with plants. Both cause tumor-like structures, and possess complex molecular mechanisms to interact with the host plants (Sawada et al. 2003). Both bacteria are reported to have acquired large number of genes horizontally (Kaneko et al. 2002; Salanoubat et al. 2002). Careful analysis of the genes that are transferred between the two bacteria can help to understand the mechanisms of pathogen–host interactions in these species, as well in other cases of HGT detected between species with similar life styles.
The strongest limitation of the types of the network reconstruction presented here is the inability of the ancestral state inference methods to precisely establish the donor organism for a HGT event. Often a HGT event is inferred across nodes of the tree that existed at different time periods. In this case, GeneTrace determines the donor group of organisms rather than a particular donor species, and the prediction should be read as "the donor is a progeny of the node." Although this effect might influence the character of the inferred network, the consistency between the results of medium- and high-confidence HGT vine width thresholds, as well any input data used in this study, indicates that the properties of the phylogenetic network reported here are genuine and realistic. A method to correctly infer HGT donors should greatly improve reconstruction of the network.
The GeneTrace method applied here uses phylogenetic distribution as a marker of HGT events. However, HGT more often occurs between related organisms, followed by homologous recombination (Vulic et al. 1997 Another limitation is our inability to determine the correct path across organisms when multiple HGT events happened. Although the probabilistic schema described in the Methods section was designed to reduce the impact of this phenomenon, identification of the exact order and direction of HGT events would drastically improve reconstruction of the network. The hubs of the HGT network presented here might partially result from the phylogenetic coverage of the sequenced species. When the coverage is low, multiple HGT events accumulate on long branches, and an artificial "hub" might appear. Thus, the reconstruction and understanding of the net of life will improve with better phylogenetic representation of sequenced organisms.
The currently acceptable representation of phylogenetic data is in the form of a tree-like structure in a two-dimensional space, often referred to as a "dendrogram" (meaning tree-graph in Greek). This presentation has the limitation of an inherent inability to depict HGT events. We propose to represent the phylogenetic data in the form of a three-dimensional tree, where beyond a tree drawn in the conventional two-dimensional space, HGT vines require a third dimension. When convergence of gene content is particularly high, participating nodes can be drawn closer in the third dimension. An example of such drawing is shown in Figure 3, with real data from this study. The full tree is available in VRML format, including all species identifiers (Janssen et al. 2003 Our results suggest that the connectivity of microbial HGT network has a power-law behavior; i.e., the connectivity distribution appears as a decreasing straight line on a log-log scale (Fig. 2). A network in which connectivity of nodes distributes as a power-law has also scale-free and small-world properties. Scale-free networks display identical properties when any random subset of the complete network is sampled, suggesting that our conclusions should not be strongly affected by an ever-increasing number of genomes.
In a small-world network, the average shortest path between any two of its nodes (termed "network diameter") involves traversing only relatively few nodes. This has a profound ecological meaning and strong implications for genome evolution. In the context of the HGT network, a small-world structure means that a substantially beneficial gene appearing in any organism can swing across species barriers and reach any other organism via a very small number of HGT events. In fact, this prediction of our hypothesis has an independent verification from the "experiment" of antibiotics-resistance genes that are known to spread extremely rapidly across species (Jacoby 1996
In order to reconstruct the phylogenetic network of microbial species, we required a data set of orthologs across all currently sequenced species. We used BLASTP (Altschul et al. 1997 ) to find best bidirectional hits across 165 microbial genomes in COGENT database release 184 (Janssen et al. 2003). To eliminate paralogy, we used only bidirectional best hits across genomes. We then clustered these hits by using Markov clustering algorithm (MCL) (Enright et al. 2002). The exhaustive nature of this schema ensures that all genes that had at least one bidirectional best hit in another organism are represented (L. Goldovsky, P. Jenssen, D. Ahrén, B. Audit, I. Cases, N. Darzentas, A.J. Enright, N. López-Bigas, J.M. Peregrin-Alvarez, M. Smith, et al., in prep.). We call the resulting protein families used for the analysis described herein as the "OFAM" data set. This data set is accessible at http://cgg.ebi.ac.uk/services/ortho-fam/.
To ensure that the results are not an artifact of the orthology definition, we used orthology information for 110 species from the STRING database (von Mering et al. 2003 ), from which we cross-linked 106 species to COGENT, resulting in 98 prokaryotic species, after excluding Eukaryotes. STRING adopts the definition of orthologs as groups of homologs built from at least one triplet of best-matching pairs of sequences, also known as clusters of orthologous genes (COGs) (Tatusov et al. 1997).
To reconstruct the microbial phylogenetic network, we required a phylogenetic tree. There are many methods for the reconstruction of phylogenetic trees (see Introduction); however, none guarantees 100% accuracy. To avoid biases generated by any single tree, we used three methods of genome-based phylogenetic reconstruction, i.e., gene content (Korbel et al. 2002 ), average gene similarity (Clarke et al. 2002), and genome conservation (Kunin et al. 2005). The first method derives phylogenetic distances from conservation of gene content, the second uses only sequence similarity between genomes, and the third combines the two measures to achieve maximum precision and contrast (Kunin et al. 2005). While being based on complete genomes, all these methods produce phylogenies that are remarkably similar to the classical 16S rRNA trees. All methods are implemented as it appears on the Genome Phylogeny Server (http://cgg.ebi.ac.uk/cgi-bin/gps/GPS.pl) and described elsewhere (Kunin et al. 2005). Only results consistent across different trees and with consistently high jacknife scores (Kunin and Ouzounis 2003b) are considered robust. For STRING data, we used a gene content tree constructed according to (Korbel et al. 2002).
Just as there are many methods to reconstruct phylogenetic trees, there are several available methods to identify HGT events. We could not use methods that are based on identification of biased GC content or codon usage, as these can only identify recently acquired genes and are not designed to reconstruct early events. We thus used GeneTrace—a method that identifies HGT from the phylogenetic distribution of protein families on the tree of life (Kunin and Ouzounis 2003b
A limitation of the GeneTrace approach to reconstructing HGT events is its inability to distinguish between the donor and the acceptor genomes (Kunin and Ouzounis 2003b Consider a situation when a protein family appears twice in distant sections of a tree. In this case, at least one HGT event may be necessary to explain the phylogenetic distribution of the family. Consider now a protein family that has three dispersed roots within a tree. Then, at least two horizontal transfer events are necessary to explain the distribution. However, simple linking of all nodes creates three possible edges for horizontal transfer. Assuming equal probability for all possible scenarios, we then assign the value of 2/3 as a probability for each possible event to be depicted correctly (and 1/3 for an incorrect detection). Thus, while the minimal number of edges required to connect all nodes (n) by HGT is n–1, the number of all possible connections is n(n–1)/2. This gives us the probability that each of the edges describes a valid HGT event as (n–1)/(n(n–1)/2), or 2/n. Thus, to each node that connects independent origins of a protein family, previously labeled by GeneTrace as arising from HGT, we assign a probability of 2/n. To describe the inferred sum of all HGT events between two nodes within an evolutionary net, we sum up all probabilities of transfer for each gene family transferred between the two nodes and term the resulting edge as "vine" and the weight of the edge as "vine width."
We thank members of the Computational Genomics Group for useful discussions. CAO acknowledges support from the UK Medical Research Council and IBM Research.
Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.3666505. Article published online before print in June 2005.
1 Present address: DOE Joint Genome Institute, Walnut Creek, California 94598, USA.
2 Corresponding author. [Supplemental material is available online at www.genome.org.]
Alfano, J.R. and Collmer, A. 2004. Type III secretion system effector proteins: Double agents in bacterial disease and plant defense. Annu. Rev. Phytopathol. 42: 385-414.[CrossRef][Medline]
Altschul, S.F., Madden, T.L., Schaffer, A.A., Zhang, J., Zhang, Z., Miller, W., and Lipman, D.J. 1997. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 25: 3389-3402. Bapteste, E., Boucher, Y., Leigh, J., and Doolittle, W.F. 2004. Phylogenetic reconstruction and lateral gene transfer. Trends Microbiol. 12: 406-411.[CrossRef][Medline] Boucher, Y., Douady, C.J., Papke, R.T., Walsh, D.A., Boudreau, M.E., Nesbo, C.L., Case, R.J., and Doolittle, W.F. 2003. Lateral gene transfer and the origins of prokaryotic groups. Annu. Rev. Genet. 37: 283-328.[CrossRef][Medline] Brown, J.R., Douady, C.J., Italia, M.J., Marshall, W.E., and Stanhope, M.J. 2001. Universal trees based on large combined protein sequence data sets. Nat. Genet. 28: 281-285.[CrossRef][Medline]
Clarke, G.D., Beiko, R.G., Ragan, M.A., and Charlebois, R.L. 2002. Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores. J. Bacteriol. 184: 2072-2080. Darwin, C. 1859. On the origin of species by means of natural selection. J. Murray, London.
Doolittle, R.F. 1981. Similar amino acid sequences: Chance or common ancestry? Science 214: 149-159.
———. 1999. Phylogenetic classification and the universal tree. Science 284: 2124-2129.
Enright, A.J., Van Dongen, S., and Ouzounis, C.A. 2002. An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 30: 1575-1584.
Fitz-Gibbon, S.T. and House, C.H. 1999. Whole genome-based phylogenetic analysis of free-living microorganisms. Nucleic Acids Res. 27: 4218-4222.
Fox, G.E., Stackebrandt, E., Hespell, R.B., Gibson, J., Maniloff, J., Dyer, T.A., Wolfe, R.S., Balch, W.E., Tanner, R.S., Magrum, L.J., et al. 1980. The phylogeny of prokaryotes. Science 209: 457-463. Genin, S. and Boucher, C. 2004. Lessons learned from the genome analysis of Ralstonia solanacearum. Annu. Rev. Phytopathol. 42: 107-134.[CrossRef][Medline]
Glockner, F.O., Kube, M., Bauer, M., Teeling, H., Lombardot, T., Ludwig, W., Gade, D., Beck, A., Borzym, K., Heitmann, K., et al. 2003. Complete genome sequence of the marine planctomycete Pirellula sp. strain 1. Proc. Natl. Acad. Sci. 100: 8298-8303. Gusfield, D., Eddhu, S., and Langley, C. 2004. Optimal, efficient reconstruction of phylogenetic networks with constrained recombination. J. Bioinform. Comput. Biol. 2: 173-213.[CrossRef][Medline] Jacoby, G.A. 1996. Antimicrobial-resistant pathogens in the 1990s. Annu. Rev. Med. 47: 169-179.[CrossRef][Medline]
Janssen, P., Enright, A.J., Audit, B., Cases, I., Goldovsky, L., Harte, N., Kunin, V., and Ouzounis, C.A. 2003. COmplete GENome Tracking (COGENT): A flexible data environment for computational genomics. Bioinformatics 19: 1451-1452. Kaneko, T., Nakamura, Y., Sato, S., Minamisawa, K., Uchiumi, T., Sasamoto, S., Watanabe, A., Idesawa, K., Iriguchi, M., Kawashima, K., et al. 2002. Complete genomic sequence of nitrogen-fixing symbiotic bacterium Bradyrhizobium japonicum USDA110. DNA Res. 9: 189-197.[Abstract] Kiers, E.T., Rousseau, R.A., West, S.A., and Denison, R.F. 2003. Host sanctions and the legume-rhizobium mutualism. Nature 425: 78-81.[CrossRef][Medline] Korbel, J.O., Snel, B., Huynen, M.A., and Bork, P. 2002. SHOT: A web server for the construction of genome phylogenies. Trends Genet. 18: 158-162.[CrossRef][Medline]
Kunin, V. and Ouzounis, C.A. 2003a. The balance of driving forces during genome evolution in prokaryotes. Genome Res. 13: 1589-1594.
———. 2003b. GeneTRACE: Reconstruction of gene content of ancestral species. Bioinformatics 19: 1412-1416.
Kunin, V., Ahren, D., Goldovsky, L., Janssen, P., and Ouzounis, C.A. 2005. Measuring genome conservation across taxa: Divided strains and united kingdoms. Nucleic Acids Res. 33: 616-621.
Lin, J. and Gerstein, M. 2000. Whole-genome trees based on the occurrence of folds and orthologs: Implications for comparing genomes on different levels. Genome Res. 10: 808-818. Makarenkov, V. and Legendre, P. 2004. From a phylogenetic tree to a reticulated network. J. Comput. Biol. 11: 195-212.[CrossRef][Medline] Martin, W. 1999. Mosaic bacterial chromosomes: A challenge en route to a tree of genomes. Bioessays 21: 99-104.[CrossRef][Medline] Nakamura, Y., Itoh, T., Matsuda, H., and Gojobori, T. 2004. Biased biological functions of horizontally transferred genes in prokaryotic genomes. Nat. Genet. 36: 760-766.[CrossRef][Medline]
Piel, W.H., Sanderson, M.J., and Donoghue, M.J. 2003. The small-world dynamics of tree networks and data mining in phyloinformatics. Bioinformatics 19: 1162-1168. Rokas, A., Williams, B.L., King, N., and Carroll, S.B. 2003. Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature 425: 798-804.[CrossRef][Medline] Salanoubat, M., Genin, S., Artiguenave, F., Gouzy, J., Mangenot, S., Arlat, M., Billault, A., Brottier, P., Camus, J.C., Cattolico, L., et al. 2002. Genome sequence of the plant pathogen Ralstonia solanacearum. Nature 415: 497-502.[CrossRef][Medline] Sawada, H., Kuykendall, L.D., and Young, J.M. 2003. Changing concepts in the systematics of bacterial nitrogen-fixing legume symbionts. J. Gen. Appl. Microbiol. 49: 155-179. Snel, B., Bork, P., and Huynen, M.A. 1999. Genome phylogeny based on gene content. Nat. Genet. 21: 108-110.[CrossRef][Medline]
———. 2002. Genomes in flux: The evolution of archaeal and proteobacterial gene content. Genome Res. 12: 17-25.
Streit, W.R., Schmitz, R.A., Perret, X., Staehelin, C., Deakin, W.J., Raasch, C., Liesegang, H., and Broughton, W.J. 2004. An evolutionary hot spot: The pNGR234b replicon of Rhizobium sp. strain NGR234. J. Bacteriol. 186: 535-542.
Tatusov, R.L., Koonin, E.V., and Lipman, D.J. 1997. A genomic perspective on protein families. Science 278: 631-637.
Tekaia, F., Lazcano, A., and Dujon, B. 1999. The genomic tree as revealed from whole proteome comparisons. Genome Res. 9: 550-557.
von Mering, C., Huynen, M., Jaeggi, D., Schmidt, S., Bork, P., and Snel, B. 2003. STRING: A database of predicted functional associations between proteins. Nucleic Acids Res. 31: 258-261.
Vulic, M., Dionisio, F., Taddei, F., and Radman, M. 1997. Molecular keys to speciation: DNA polymorphism and the control of genetic exchange in enterobacteria. Proc. Natl. Acad. Sci. 94: 9763-9767. Wang, L., Zhang, K., and Zhang, L. 2001. Perfect phylogenetic networks with recombination. J. Comput. Biol. 8: 69-78.[CrossRef][Medline]
http://cgg.ebi.ac.uk/services/ortho-fam/; OFAM data set. http://cgg.ebi.ac.uk/cgi-bin/gps/GPS.pl; Genome Phylogeny Server.
Received January 6, 2005; accepted in revised format May 2, 2005.  CiteULike  Connotea  Del.icio.us  Digg  Reddit  Technorati What's this?
This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Top
ABSTRACT
Results
