Harvesting evolutionary signals in a forest of prokaryotic gene trees

Phylogenomic studies produce increasingly large phylogenetic forests of trees with patchy taxonomical sampling. Typically, prokaryotic data generate thousands of gene trees of all sizes that are difficult, if not impossible, to root. Their topologies do not match the genealogy of lineages, as they are influenced not only by duplication, losses, and vertical descent but also by lateral gene transfer (LGT) and recombination. Because this complexity in part reflects the diversity of evolutionary processes, the study of phylogenetic forests is thus a great opportunity to improve our understanding of prokaryotic evolution. Here, we show how the rich evolutionary content of such novel phylogenetic objects can be exploited through the development of new approaches designed specifically for extracting the multiple evolutionary signals present in the forest of life, that is, by slicing up trees into remarkable bits and pieces: clans, slices, and clips. We harvested a forest of 6,901 unrooted gene trees comprising up to 100 prokaryotic genomes (41 archaea and 59 bacteria) to search for evolutionary events that a species tree would not account for. We identified 1) trees and partitions of trees that reflected the lifestyle of organisms rather than their taxonomy, 2) candidate lifestyle-specific genetic modules, used by distinct unrelated organisms to adapt to the same environment, 3) gene families, nonrandomly distributed in the functional space, that were frequently exchanged between archaea and bacteria, sometimes without major changes in their sequences. Finally, 4) we reconstructed polarized networks of genetic partnerships between archaea and bacteria to describe some of the rules affecting LGT between these two Domains.     

Phylogenomic networks

Phylogenomics is aimed at studying functional and evolutionary aspects of genome biology using phylogenetic analysis of whole genomes. Current approaches to genome phylogenies are commonly founded in terms of phylogenetic trees. However, several evolutionary processes are non tree-like in nature, including recombination and lateral gene transfer (LGT). Phylogenomic networks are a special type of phylogenetic network reconstructed from fully sequenced genomes. The network model, comprising genomes connected by pairwise evolutionary relations, enables the reconstruction of both vertical and LGT events. Modeling genome evolution in the form of a network enables the use of an extensive toolbox developed for network research. The structural properties of phylogenomic networks open up fundamentally new insights into genome evolution.     

Distant horizontal gene transfer is rare for multiple families of prokaryotic insertion sequences

Horizontal gene transfer in prokaryotes is rampant on short and intermediate evolutionary time scales. It poses a fundamental problem to our ability to reconstruct the evolutionary tree of life. Is it also frequent over long evolutionary distances? To address this question, we analyzed the evolution of 2,091 insertion sequences from all 20 major families in 438 completely sequenced prokaryotic genomes. Specifically, we mapped insertion sequence occurrence on a 16S rDNA tree of the genomes we analyzed, and we also constructed phylogenetic trees of the insertion sequence transposase coding sequences. We found only 30 cases of likely horizontal transfer among distantly related prokaryotic clades. Most of these horizontal transfer events are ancient. Only seven events are recent. Almost all of these transfer events occur between pairs of human pathogens or commensals. If true also for other, non-mobile DNA, the rarity of distant horizontal transfer increases the odds of reliable phylogenetic inference from sequence data.     

Does the 'Ring of Life' ring true?

In a recent stimulating paper, Rivera and Lake applied a new phylogenetic method to study the evolution of genomes, which challenges the classical representation of the Tree of Life. Acknowledging the evolutionary importance of lateral gene transfer, they used the conditioned genome approach to reconstruct the Tree of Life, and in the end proposed a Ring of Life. They explained that the Ring of Life structure is a result of a single fusion event between two prokaryotic genomes at the base of the eukaryotic tree, probably between the ancestors of a photosynthetic bacterium and an archaeon. Because this constitutes an important conclusion with regards to the evolutionary process and origin of the eukaryotic cell, their work deserves further attention before these conclusions can be accepted. Here we question the reconstruction and the meaning of the Ring of Life. In addition to general problems associated with gene-content-based phylogenetic analyses, we discuss some implicit premises and potential weaknesses of the conditioned genome method and conclude that, although Rivera and Lake's conclusions might be right, they have not been established by their current approach.     

Evolutionary patterns in prokaryotic genomes

Prokaryotic genomics is shifting towards comparative approaches to unravel how and why genomes change over time. Both phylogenetic and population genetics approaches are required to dissect the relative roles of selection and drift under these conditions. Lineages evolve adaptively by selection of changes in extant genomes and the way this occurs is being explored from a systemic and evolutionary perspective to understand how mutations relate with gene repertoire changes and how both are contextualized in cellular networks. Through an increased appreciation of genome dynamics in given ecological contexts, a more detailed picture of the genetic basis of prokaryotic evolution is emerging.     

Genomic and phylogenetic perspectives on the evolution of prokaryotes

Prokaryotes have been at the forefront of the genome sequencing revolution. Many genomes have been completely sequenced, revealing much about bacterial and archaeal genome content and organization. Yet, a meaningful evolutionary picture of prokaryotes still eludes us. Much of the problem lies in understanding the mode and tempo of genome evolution. Here phenylalanyl-tRNA synthetase is used as an example of the complex interplay among lateral gene transfer, operon recombination, and gene recruitment in the evolution of some prokaryotic genes. Promising new approaches to genomic analyses, which could add to our understanding prokaryotic evolution and help in their classification, are discussed.     

Transposable elements and viruses as factors in adaptation and evolution: an expansion and strengthening of the TE‐Thrust hypothesis

In addition to the strong divergent evolution and significant and episodic evolutionary transitions and speciation we previously attributed to TE‐Thrust, we have expanded the hypothesis to more fully account for the contribution of viruses to TE‐Thrust and evolution. The concept of symbiosis and holobiontic genomes is acknowledged, with particular emphasis placed on the creativity potential of the union of retroviral genomes with vertebrate genomes. Further expansions of the TE‐Thrust hypothesis are proposed regarding a fuller account of horizontal transfer of TEs, the life cycle of TEs, and also, in the case of a mammalian innovation, the contributions of retroviruses to the functions of the placenta. The possibility of drift by TE families within isolated demes or disjunct populations, is acknowledged, and in addition, we suggest the possibility of horizontal transposon transfer into such subpopulations. “Adaptive potential” and “evolutionary potential” are proposed as the extremes of a continuum of “intra‐genomic potential” due to TE‐Thrust. Specific data is given, indicating “adaptive potential” being realized with regard to insecticide resistance, and other insect adaptations. In this regard, there is agreement between TE‐Thrust and the concept of adaptation by a change in allele frequencies. Evidence on the realization of “evolutionary potential” is also presented, which is compatible with the known differential survivals, and radiations of lineages. Collectively, these data further suggest the possibility, or likelihood, of punctuated episodes of speciation events and evolutionary transitions, coinciding with, and heavily underpinned by, intermittent bursts of TE activity.     

Identification of Horizontal Gene Transfer from Phylogenetic Gene Trees

We suggest a new procedure to search for the genes with horizontal transfer events in their evolutionary history. The search is based on analysis of topology difference between the phylogenetic trees of gene (protein) groups and the corresponding phylogenetic species trees. Numeric values are introduced to measure the discrepancy between the trees. This approach was applied to analyze 40 prokaryotic genomes classified into 132 classes of orthologs. This resulted in a list of the candidate genes for which the hypothesis of horizontal transfer in evolution looks true.     

Identification of the horizontal gene transfer based on phylogenetic data

We suggest a new procedure to search for the genes with horizontal transfer events in their evolutionary history. The search is based on analysis of topology difference between the phylogenetic trees of gene (protein) groups and the corresponding phylogenetic species trees. Numeric values are introduced to measure the discrepancy between the trees. This approach was applied to analyze 40 prokaryotic genomes classified into 132 classes of orthologs. This resulted in a list of the candidate genes for which the hypothesis of horizontal transfer in evolution looks true.     

G+C3 structuring along the genome: a common feature in prokaryotes

The heterogeneity of gene nucleotide content in prokaryotic genomes is commonly interpreted as the result of three main phenomena: (1) genes undergo different selection pressures both during and after translation (affecting codon and amino acid choice); (2) genes undergo different mutational pressure whether they are on the leading or lagging strand; and (3) genes may have different phylogenetic origins as a result of lateral transfers. However, this view neglects the necessity of organizing genetic information on a chromosome that needs to be replicated and folded, which may add constraints to single gene evolution. As a consequence, genes are potentially subjected to different mutation and selection pressures, depending on their position in the genome. In this paper, we analyze the structuring of different codon usage measures along completely sequenced bacterial genomes. We show that most of them are highly structured, suggesting that genes have different base content, depending on their location on the chromosome. A peculiar pattern of genome structure, with a tendency toward an A+T-enrichment near the replication terminus, is found in most bacterial phyla and may reflect common chromosome constraints. Several species may have lost this pattern, probably because of genome rearrangements or integration of foreign DNA. We show that in several species, this enrichment is associated with an increase of evolutionary rate and we discuss the evolutionary implications of these results. We argue that structural constraints acting on the circular chromosome are not negligible and that this natural structuring of bacterial genomes may be a cause of overestimation in lateral gene transfer predictions using codon composition indices.     

Mammalian evolution and biomedicine: new views from phylogeny

Recent progress resolving the phylogenetic relationships of the major lineages of mammals has had a broad impact in evolutionary biology, comparative genomics and the biomedical sciences. Novel insights into the timing and historical biogeography of early mammalian diversification have resulted from a new molecular tree for placental mammals coupled with dating approaches that relax the assumption of the molecular clock. We highlight the numerous applications to come from a well-resolved phylogeny and genomic prospecting in multiple lineages of mammals, from identifying regulatory elements in mammalian genomes to assessing the functional consequences of mutations in human disease loci and those driving adaptive evolution.     

Homologs of eukaryotic Ras superfamily proteins in prokaryotes and their novel phylogenetic correlation with their eukaryotic analogs

Ras superfamily proteins are key regulators in a wide variety of cellular processes. Previously, they were considered to be specific to eukaryotes, and MglA, a group of obviously different prokaryotic proteins, were recognized as their only prokaryotic analogs or even ancestors. Here, taking advantage of quite a current accumulation of prokaryotic genomic databases, we have investigated the existence and taxonomic distribution of Ras superfamily protein homologs in a much wider prokaryotic range, and analyzed their phylogenetic correlation with their eukaryotic analogs. Thirteen unambiguous prokaryotic homologs, which possess the GDP/GTP-binding domain with all the five characteristic motifs of their eukaryotic analogs, were identified in 12 eubacteria and one archaebacterium, respectively. In some other archaebacteria, including four methanogenic archaebacteria and three Thermoplasmales, homologs were also found, but with the GDP/GTP-binding domains not containing all the five characteristic motifs. Many more MglA orthologs were identified than in previous studies mainly in delta-proteobacteria, and all were shown to have common unique features distinct from the Ras superfamily proteins. Our phylogenetic analysis indicated eukaryotic Rab, Ran, Ras, and Rho families have the closest phylogenetic correlation with the 13 unambiguous prokaryotic homologs, whereas the other three eukaryotic protein families (SRbeta, Sar1, and Arf) branch separately from them, but have a relatively close relationship with the methanogenic archaebacterial homologs and MglA. Although homologs were identified in a relative minority of prokaryotes with genomic databases, their presence in a relatively wide variety of lineages, their unique sequence characters distinct from those of eukaryotic analogs, and the topology of our phylogenetic tree altogether do not support their origin from eukaryotes as a result of lateral gene transfer. Therefore, we argue that Ras superfamily proteins might have already emerged at least in some prokaryotic lineages, and that the seven eukaryotic protein families of the Ras superfamily may have two independent prokaryotic origins, probably reflecting the 'fusion' evolutionary history of the eukaryotic cell.     

En route to a genome-based classification of Archaea and Bacteria?

Given the considerable promise whole-genome sequencing offers for phylogeny and classification, it is surprising that microbial systematics and genomics have not yet been reconciled. This might be due to the intrinsic difficulties in inferring reasonable phylogenies from genomic sequences, particularly in the light of the significant amount of lateral gene transfer in prokaryotic genomes. However, recent studies indicate that the species tree and the hierarchical classification based on it are still meaningful concepts, and that state-of-the-art phylogenetic inference methods are able to provide reliable estimates of the species tree to the benefit of taxonomy. Conversely, we suspect that the current lack of completely sequenced genomes for many of the major lineages of prokaryotes and for most type strains is a major obstacle in progress towards a genome-based classification of microorganisms. We conclude that phylogeny-driven microbial genome sequencing projects such as the Genomic Encyclopaedia of Archaea and Bacteria (GEBA) project are likely to rectify this situation.     

Lateral genomics

More than 20 complete prokaryotic genome sequences are now publicly available, each by itself an unparalleled resource for understanding organismal biology. Collectively, these data are even more powerful: they could force a dramatic reworking of the framework in which we understand biological evolution. It is possible that a single universal phylogenetic tree is not the best way to depict relationships between all living and extinct species. Instead a web- or net-like pattern, reflecting the importance of horizontal or lateral gene transfer between lineages of organisms, might provide a more appropriate visual metaphor. Here, I ask whether this way of thinking is really justified, and explore its implications.     

Lateral genomics

More than 20 complete prokaryotic genome sequences are now publicly available, each by itself an unparalleled resource for understanding organismal biology. Collectively, these data are even more powerful: they could force a dramatic reworking of the framework in which we understand biological evolution. It is possible that a single universal phylogenetic tree is not the best way to depict relationships between all living and extinct species. Instead a web- or net-like pattern, reflecting the importance of horizontal or lateral gene transfer between lineages of organisms, might provide a more appropriate visual metaphor. Here, I ask whether this way of thinking is really justified, and explore its implications.     

Lateral genomics

More than 20 complete prokaryotic genome sequences are now publicly available, each by itself an unparalleled resource for understanding organismal biology. Collectively, these data are even more powerful: they could force a dramatic reworking of the framework in which we understand biological evolution. It is possible that a single universal phylogenetic tree is not the best way to depict relationships between all living and extinct species. Instead a web- or net-like pattern, reflecting the importance of horizontal or lateral gene transfer between lineages of organisms, might provide a more appropriate visual metaphor. Here, I ask whether this way of thinking is really justified, and explore its implications.     

The slow road to the eukaryotic genome

The eukaryotic genome is a mosaic of eubacterial and archaeal genes in addition to those unique to itself. The mosaic may have arisen as the result of two prokaryotes merging their genomes, or from genes acquired from an endosymbiont of eubacterial origin. A third possibility is that the eukaryotic genome arose from successive events of lateral gene transfer over long periods of time. This theory does not exclude the endosymbiont, but questions whether it is necessary to explain the peculiar set of eukaryotic genes. We use phylogenetic studies and reconstructions of ancestral first appearances of genes on the prokaryotic phylogeny to assess evidence for the lateral gene transfer scenario. We find that phylogenies advanced to support fusion can also arise from a succession of lateral gene transfer events. Our reconstructions of ancestral first appearances of genes reveal that the various genes that make up the eukaryotic mosaic arose at different times and in diverse lineages on the prokaryotic tree, and were not available in a single lineage. Successive events of lateral gene transfer can explain the unusual mosaic structure of the eukaryotic genome, with its content linked to the immediate adaptive value of the genes its acquired. Progress in understanding eukaryotes may come from identifying ancestral features such as the eukaryotic splicesome that could explain why this lineage invaded, or created, the eukaryotic niche.     

Gene finding with a hidden Markov model of genome structure and evolution

Motivation: A growing number of genomes are sequenced. The differences in evolutionary pattern between functional regions can thus be observed genome-wide in a whole set of organisms. The diverse evolutionary pattern of different functional regions can be exploited in the process of genomic annotation. The modelling of evolution by the existing comparative gene finders leaves room for improvement. Results: A probabilistic model of both genome structure and evolution is designed. This type of model is called an Evolutionary Hidden Markov Model (EHMM), being composed of an HMM and a set of region-specific evolutionary models based on a phylogenetic tree. All parameters can be estimated by maximum likelihood, including the phylogenetic tree. It can handle any number of aligned genomes, using their phylogenetic tree to model the evolutionary correlations. The time complexity of all algorithms used for handling the model are linear in alignment length and genome number. The model is applied to the problem of gene finding. The benefit of modelling sequence evolution is demonstrated both in a range of simulations and on a set of orthologous human/mouse gene pairs. AVAILABILITY: Free availability over the Internet on www server: http://www.birc.dk/Software/evogene.     

Detection of Lateral Gene Transfer Events in the Prokaryotic tRNA Synthetases by the Ratios of Evolutionary Distances Method

The availability of large numbers of genomic sequences has demonstrated the importance of lateral gene transfer (LGT) in prokaryotic evolution. However, considerable uncertainty remains concerning the frequency of LGT compared to other evolutionary processes. To examine LGTs in ancient lineages of prokaryotes a method was developed that utilizes the ratios of evolutionary distances (RED) to distinguish between alternative evolutionary histories. The advantages of this approach are that the variability inherent in comparing protein sequences is transparent, the direction of LGT and the relative rates of evolution are readily identified, and it is possible to detect other types of evolutionary events. This method was standardized using 35 genes encoding ribosomal proteins that were believed to share a vertical evolution. Using RED-T, an original computer program designed to implement the RED method, the evolution of the genes encoding the 20 aminoacyl-tRNA synthetases was examined. Although LGTs were common in the evolution of the aminoacyl-tRNA synthetases, they were not sufficient to obscure the organismal phylogeny. Moreover, much of the apparent complexity of the gene tree was consistent with the formation of the paralogs in the ancestors to the modern lineages followed by more recent loss of one paralog or the other.     

Molecular evolutionary processes and conflicting gene trees: The hominoid case

Molecular evolutionary processes modify DNA over time, creating both newly derived substitutions shared by related descendant lineages (phylogenetic signal) and “false” similarities which confound phylogenetic reconstruction (homoplasy). However, some types of DNA regions, for example those containing tandem duplicate repeats, are preferentially subject to homoplasy-inducing processes such as sporadically occurring concerted evolution and DNA insertion/deletion. This added level of homoplasic “noise” can make DNA regions with repeats less reliable in phylogenetic reconstruction than those without repeats. Most molecular datasets which distinguish among African hominoids support a human-chimpanzee clade; the most notable exception is from the involucrin gene. However, phylogenetic resolution supporting a chimpanzee-gorilla clade is based entirely on involucrin DNA repeat regions. This is problematic because (1) involucrin repeats are difficult to align, and published alignments are contradictory; (2) involucrin repeats are subject to DNA insertion/deletion; (3) gorillas are polymorphic in that some do not have repeats reported to be synapomorphies linking chimpanzees and gorillas. Gene tree/species tree conflicts can occur due to the sorting of ancestrally polymorphic alleles during speciation. Because hominoid females transfer between groups, mitochondrial and nuclear gene flow occur to the same extent, and the probability of conflict between mitochondrial and nuclear gene trees is theoretically low. When hominoid intraspecific mitochondrial variability is taken into account [based on cytochrome oxidase subunit II (COII) gene sequences], humans and chimpanzees are most closely related, showing the same relative degree of separation from gorillas as when single individuals representing species are analyzed. Conflicting molecular phylogenies can be explained in terms of molecular evolutionary processes and sorting of ancient polymorphisms. This perspective can enhance our understanding of hominoid molecular phylogenies. © 1994 Wiley-Liss, Inc.