Harvesting evolutionary signals in a forest of prokaryotic gene trees

Genomes of Vibrio metoecus co-isolated with Vibrio cholerae extend our understanding of differences between these closely related species

BACKGROUND

Vibrio cholerae, the causative agent of cholera, is a well-studied species, whereas Vibrio metoecus is a recently described close relative that is also associated with human infections. The availability of V. metoecus genomes provides further insight into its genetic differences from V. cholerae. Additionally, both species have been co-isolated from a cholera-free brackish coastal pond and have been suggested to interact with each other by horizontal gene transfer (HGT).

RESULTS

The genomes of 17 strains from each species were sequenced. All strains share a large core genome (2675 gene families) and very few genes are unique to each species (< 3% of the pan-genome of both species). This led to the identification of potential molecular markers-for nitrite reduction, as well as peptidase and rhodanese activities-to further distinguish V. metoecus from V. cholerae. Interspecies HGT events were inferred in 21% of the core genes and 45% of the accessory genes. A directional bias in gene transfer events was found in the core genome, where V. metoecus was a recipient of three times (75%) more genes from V. cholerae than it was a donor (25%).

CONCLUSION

V. metoecus was misclassified as an atypical variant of V. cholerae due to their resemblance in a majority of biochemical characteristics. More distinguishing phenotypic assays can be developed based on the discovery of potential gene markers to avoid any future misclassifications. Furthermore, differences in relative abundance or seasonality were observed between the species and could contribute to the bias in directionality of HGT.

Horizontal gene transfer provides insights into the deep evolutionary history and biology of Trichinella

Evolution involves temporal changes in the characteristics of a species that are subsequently propagated or rejected through natural selection. In the case of parasites, host switching also plays a prominent role in the evolutionary process. These changes are rooted in genetic variation and gene flow where genes may be deleted, mutated (sequence), duplicated, rearranged and/or translocated and then transmitted through vertical gene transfer. However, the introduction of new genes is not driven only by Mendelian inheritance and mutation but also by the introduction of DNA from outside a lineage in the form of horizontal gene transfer between donor and recipient organisms. Once introduced and integrated into the biology of the recipient, vertical inheritance then perpetuates the newly acquired genetic factor, where further functionality may involve co-option of what has become a pre-existing physiological capacity. Upon sequencing the Trichinella spiralis (Clade I) genome, a cyanate hydratase (cyanase) gene was identified that is common among bacteria, fungi, and plants, but rarely observed among other eukaryotes. The sequence of the Trichinella cyanase gene clusters with those derived from the Kingdom Plantae in contrast to the genes found in some Clade III and IV nematodes that cluster with cyanases of bacterial origin. Phylogenetic analyses suggest that the Trichinella cyanase was acquired during the Devonian period and independently from those of other nematodes. These data may help inform us of the deep evolutionary history and ecological connectivity of early ancestors within the lineage of contemporary Trichinella. Further, in many extant organisms, cyanate detoxification has been largely superseded by energy requirements for metabolism. Thus, deciphering the function of Trichinella cyanase may provide new avenues for treatment and control.

Population Analysis of Vibrio cholerae in Aquatic Reservoirs Reveals a Novel Sister Species (Vibrio paracholerae sp. nov.) with a History of Association with Humans

Most efforts to understand the biology of Vibrio cholerae have focused on a single group, the pandemic-generating lineage harboring the strains responsible for all known cholera pandemics. Consequently, little is known about the diversity of this species in its native aquatic environment. To understand the differences in the V. cholerae populations inhabiting regions with a history of cholera cases and those lacking such a history, a comparative analysis of population composition was performed. Little overlap was found in lineage compositions between those in Dhaka, Bangladesh (where cholera is endemic), located in the Ganges Delta, and those in Falmouth, MA (no known history of cholera), a small coastal town on the United States east coast. The most striking difference was the presence of a group of related lineages at high abundance in Dhaka, which was completely absent from Falmouth. Phylogenomic analysis revealed that these lineages form a cluster at the base of the phylogeny for the V. cholerae species and were sufficiently differentiated genetically and phenotypically to form a novel species. A retrospective search revealed that strains from this species have been anecdotally found from around the world and were isolated as early as 1916 from a British soldier in Egypt suffering from choleraic diarrhea. In 1935, Gardner and Venkatraman unofficially referred to a member of this group as Vibrio paracholerae. In recognition of this earlier designation, we propose the name Vibrio paracholerae sp. nov. for this bacterium. Genomic analysis suggests a link with human populations for this novel species and substantial interaction with its better-known sister species. IMPORTANCE Cholera continues to remain a major public health threat around the globe. Understanding the ecology, evolution, and environmental adaptation of the causative agent (Vibrio cholerae) and tracking the emergence of novel lineages with pathogenic potential are essential to combat the problem. In this study, we investigated the population dynamics of Vibrio cholerae in an inland locality, which is known as endemic for cholera, and compared them with those of a cholera-free coastal location. We found the consistent presence of the pandemic-generating lineage of V. cholerae in Dhaka, where cholera is endemic, and an exclusive presence of a lineage phylogenetically distinct from other V. cholerae lineages. Our study suggests that this lineage represents a novel species that has pathogenic potential and a human link to its environmental abundance. The possible association with human populations and coexistence and interaction with toxigenic V. cholerae in the natural environment make this potential human pathogen an important subject for future studies.

Bipartite Network Analysis of Gene Sharings in the Microbial World

Extensive microbial gene flows affect how we understand virology, microbiology, medical sciences, genetic modification, and evolutionary biology. Phylogenies only provide a narrow view of these gene flows: plasmids and viruses, lacking core genes, cannot be attached to cellular life on phylogenetic trees. Yet viruses and plasmids have a major impact on cellular evolution, affecting both the gene content and the dynamics of microbial communities. Using bipartite graphs that connect up to 149,000 clusters of homologous genes with 8,217 related and unrelated genomes, we can in particular show patterns of gene sharing that do not map neatly with the organismal phylogeny. Homologous genes are recycled by lateral gene transfer, and multiple copies of homologous genes are carried by otherwise completely unrelated (and possibly nested) genomes, that is, viruses, plasmids and prokaryotes. When a homologous gene is present on at least one plasmid or virus and at least one chromosome, a process of "gene externalization," affected by a postprocessed selected functional bias, takes place, especially in Bacteria. Bipartite graphs give us a view of vertical and horizontal gene flow beyond classic taxonomy on a single very large, analytically tractable, graph that goes beyond the cellular Web of Life.

Genome-Wide Comparative Analysis of Phylogenetic Trees: The Prokaryotic Forest of Life

Genome-wide comparison of phylogenetic trees is becoming an increasingly common approach in evolutionary genomics, and a variety of approaches for such comparison have been developed. In this article we present several methods for comparative analysis of large numbers of phylogenetic trees. To compare phylogenetic trees taking into account the bootstrap support for each internal branch, the boot-split distance (BSD) method is introduced as an extension of the previously developed split distance (SD) method for tree comparison. The BSD method implements the straightforward idea that comparison of phylogenetic trees can be made more robust by treating tree splits differentially depending on the bootstrap support. Approaches are also introduced for detecting treelike and netlike evolutionary trends in the phylogenetic Forest of Life (FOL), i.e., the entirety of the phylogenetic trees for conserved genes of prokaryotes. The principal method employed for this purpose includes mapping quartets of species onto trees to calculate the support of each quartet topology and so to quantify the tree and net contributions to the distances between species. We describe the applications methods used to analyze the FOL and the results obtained with these methods. These results support the concept of the Tree of Life (TOL) as a central evolutionary trend in the FOL as opposed to the traditional view of the TOL as a "species tree."

Effect of the environment on horizontal gene transfer between bacteria and archaea

Background

Horizontal gene transfer, the transfer and incorporation of genetic material between different species of organisms, has an important but poorly quantified role in the adaptation of microbes to their environment. Previous work has shown that genome size and the number of horizontally transferred genes are strongly correlated. Here we consider how genome size confuses the quantification of horizontal gene transfer because the number of genes an organism accumulates over time depends on its evolutionary history and ecological context (e.g., the nutrient regime for which it is adapted).

Results

We investigated horizontal gene transfer between archaea and bacteria by first counting reciprocal BLAST hits among 448 bacterial and 57 archaeal genomes to find shared genes. Then we used the DarkHorse algorithm, a probability-based, lineage-weighted method (Podell & Gaasterland, 2007), to identify potential horizontally transferred genes among these shared genes. By removing the effect of genome size in the bacteria, we have identified bacteria with unusually large numbers of shared genes with archaea for their genome size. Interestingly, archaea and bacteria that live in anaerobic and/or high temperature conditions are more likely to share unusually large numbers of genes. However, high salt was not found to significantly affect the numbers of shared genes. Numbers of shared (genome size-corrected, reciprocal BLAST hits) and transferred genes (identified by DarkHorse) were strongly correlated. Thus archaea and bacteria that live in anaerobic and/or high temperature conditions are more likely to share horizontally transferred genes. These horizontally transferred genes are over-represented by genes involved in energy conversion as well as the transport and metabolism of inorganic ions and amino acids.

Conclusions

Anaerobic and thermophilic bacteria share unusually large numbers of genes with archaea. This is mainly due to horizontal gene transfer of genes from the archaea to the bacteria. In general, these transfers are from archaea that live in similar oxygen and temperature conditions as the bacteria that receive the genes. Potential hotspots of horizontal gene transfer between archaea and bacteria include hot springs, marine sediments, and oil wells. Cold spots for horizontal transfer included dilute, aerobic, mesophilic environments such as marine and freshwater surface waters.

The Dynamics of Genetic Interactions between Vibrio metoecus and Vibrio cholerae, Two Close Relatives Co-Occurring in the Environment

Vibrio metoecus is the closest relative of Vibrio cholerae, the causative agent of the potent diarrheal disease cholera. Although the pathogenic potential of this new species is yet to be studied in depth, it has been co-isolated with V. cholerae in coastal waters and found in clinical specimens in the United States. We used these two organisms to investigate the genetic interaction between closely related species in their natural environment. The genomes of 20 V. cholerae and 4 V. metoecus strains isolated from a brackish coastal pond on the US east coast, as well as 4 clinical V. metoecus strains were sequenced and compared with reference strains. Whole genome comparison shows 86-87% average nucleotide identity (ANI) in their core genes between the two species. On the other hand, the chromosomal integron, which occupies approximately 3% of their genomes, shows higher conservation in ANI between species than any other region of their genomes. The ANI of 93-94% observed in this region is not significantly greater within than between species, meaning that it does not follow species boundaries. Vibrio metoecus does not encode toxigenic V. cholerae major virulence factors, the cholera toxin and toxin-coregulated pilus. However, some of the pathogenicity islands found in pandemic V. cholerae were either present in the common ancestor it shares with V. metoecus, or acquired by clinical and environmental V. metoecus in partial fragments. The virulence factors of V. cholerae are therefore both more ancient and more widespread than previously believed. There is high interspecies recombination in the core genome, which has been detected in 24% of the single-copy core genes, including genes involved in pathogenicity. Vibrio metoecus was six times more often the recipient of DNA from V. cholerae as it was the donor, indicating a strong bias in the direction of gene transfer in the environment.

Production of poly(3-hydroxybutyrate) by Halomonas boliviensis in an air-lift reactor

Background

Microbial polyesters, also known as polyhydroxyalkanoates (PHAs), closely resemble physical and mechanical features of petroleum derived plastics. Recombinant Escherichia coli strains are being used in industrial production of PHAs in Stirred Tank Bioreactors (STRs). However, use of Air-Lift Reactors (ALRs) has been known to offer numerous technical operating options over STRs, and as such has been successfully implemented in many bioprocesses. Halomonas boliviensis is a halophilic bacterium that is known to assimilate various carbohydrates and convert them into a particular type of PHA known as poly(3-hydroxybutyrate) (PHB). Owing to this capability, it has been used to synthesize the polyester using hydrolysates of starch or wheat bran in stirred tank bioreactors.

Results

This research article firstly describes the production of PHB in shake flasks by H. boliviensis using different combinations of carbohydrates and partially hydrolyzed starch as carbon sources. The highest PHB yields, between 56 and 61 % (wt.), were achieved when either starch hydrolysate or a mixture of glucose and xylose were used as carbon sources. The starch hydrolysate obtained in this study was then used as carbon source in an ALR. The largest amount of PHB, 41 % (wt.), was attained after 24 hrs of cultivation during which maltose in the hydrolysate was assimilated more rapidly than glucose during active cell growth; however, the rate of assimilation of both the carbohydrates was found to be similar during synthesis of PHB. An incomplete pentose phosphate pathway, which lacks 6-phosphogluconate dehydrogenase, was deduced from the genome sequence of this bacterium and may result in the characteristic assimilation of glucose and maltose by the cells.

Conclusions

This study showed that the production of PHB by H. boliviensis using cheap substrates such as starch hydrolysate in a simple production system involving an ALR is feasible. Both maltose and glucose in the hydrolysate induce cell growth and PHB synthesis; most likely the cells balance adequately CoA and NAD(P)H during the assimilation of these carbohydrates. The combination of cheap substrates, simple production systems and the use of non-strict sterile conditions by the halophile H. boliviensis are desirable traits for large scale production of PHB, and should lead to a competitive bioprocess.

Electronic supplementary material

The online version of this article (doi:10.1186/s40709-015-0031-6) contains supplementary material, which is available to authorized users.

Horizontal gene transfer and genome evolution in Methanosarcina

Background

Genomes of Methanosarcina spp. are among the largest archaeal genomes. One suggested reason for that is massive horizontal gene transfer (HGT) from bacteria. Genes of bacterial origin may be involved in the central metabolism and solute transport, in particular sugar synthesis, sulfur metabolism, phosphate metabolism, DNA repair, transport of small molecules etc. Horizontally transferred (HT) genes are considered to play the key role in the ability of Methanosarcina spp. to inhabit diverse environments. At the moment, genomes of three Methanosarcina spp. have been sequenced, and while these genomes vary in length and number of protein-coding genes, they all have been shown to accumulate HT genes. However, previous estimates had been made when fewer archaeal genomes were known. Moreover, several Methanosarcinaceae genomes from other genera have been sequenced recently. Here, we revise the census of genes of bacterial origin in Methanosarcinaceae.

Results

About 5 % of Methanosarcina genes have been shown to be horizontally transferred from various bacterial groups to the last common ancestor either of Methanosarcinaceae, or Methanosarcina, or later in the evolution. Simulation of the composition of the NCBI protein non-redundant database for different years demonstrates that the estimates of the HGT rate have decreased drastically since 2002, the year of publication of the first Methanosarcina genome.

The phylogenetic distribution of HT gene donors is non-uniform. Most HT genes were transferred from Firmicutes and Proteobacteria, while no HGT events from Actinobacteria to the common ancestor of Methanosarcinaceae were found. About 50 % of HT genes are involved in metabolism. Horizontal transfer of transcription factors is not common, while 46 % of horizontally transferred genes have demonstrated differential expression in a variety of conditions. HGT of complete operons is relatively infrequent and half of HT genes do not belong to operons.

Conclusions

While genes of bacterial origin are still more frequent in Methanosarcinaceae than in other Archaea, most HGT events described earlier as Methanosarcina-specific seem to have occurred before the divergence of Methanosarcinaceae. Genes horizontally transferred from bacteria to archaea neither tend to be transferred with their regulators, nor in long operons.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0393-2) contains supplementary material, which is available to authorized users.

Modularized evolution in archaeal methanogens phylogenetic forest

Methanogens are methane-producing archaea that plays a key role in the global carbon cycle. To date, the evolutionary history of methanogens and closely related nonmethanogen species remains unresolved among studies conducted upon different genetic markers, attributing to horizontal gene transfers (HGTs). With an effort to decipher both congruent and conflicting evolutionary events, reconstruction of coevolved gene clusters and hierarchical structure in the archaeal methanogen phylogenetic forest, comprehensive evolution, and network analyses were performed upon 3,694 gene families from 41 methanogens and 33 closely related archaea. Our results show that 1) greater than 50% of genes are in topological dissonance with others; 2) the prevalent interorder HGTs, even for core genes, in methanogen genomes led to their scrambled phylogenetic relationships; 3) most methanogenesis-related genes have experienced at least one HGT; 4) greater than 20% of the genes in methanogen genomes were transferred horizontally from other archaea, with genes involved in cell-wall synthesis and defense system having been transferred most frequently; 5) the coevolution network contains seven statistically robust modules, wherein the central module has the highest average node strength and comprises a majority of the core genes; 6) different coevolutionary module genes boomed in different time and evolutionary lineage, constructing diversified pan-genome structures; 7) the modularized evolution is also closely related to the vertical evolution signals and the HGT rate of the genes. Overall, this study presented a modularized phylogenetic forest that describes a combination of complicated vertical and nonvertical evolutionary processes for methanogenic archaeal species.

Microbial communities evolve faster in extreme environments

Evolutionary analysis of microbes at the community level represents a new research avenue linking ecological patterns to evolutionary processes, but remains insufficiently studied. Here we report a relative evolutionary rates (rERs) analysis of microbial communities from six diverse natural environments based on 40 metagenomic samples. We show that the rERs of microbial communities are mainly shaped by environmental conditions, and the microbes inhabiting extreme habitats (acid mine drainage, saline lake and hot spring) evolve faster than those populating benign environments (surface ocean, fresh water and soil). These findings were supported by the observation of more relaxed purifying selection and potentially frequent horizontal gene transfers in communities from extreme habitats. The mechanism of high rERs was proposed as high mutation rates imposed by stressful conditions during the evolutionary processes. This study brings us one stage closer to an understanding of the evolutionary mechanisms underlying the adaptation of microbes to extreme environments.

To tree or not to tree? Genome-wide quantification of recombination and reticulate evolution during the diversification of strict intracellular bacteria

It is well known that horizontal gene transfer (HGT) is a major force in the evolution of prokaryotes. During the adaptation of a bacterial population to a new ecological niche, and particularly for intracellular bacteria, selective pressures are shifted and ecological niches reduced, resulting in a lower rate of genetic connectivity. HGT and positive selection are therefore two important evolutionary forces in microbial pathogens that drive adaptation to new hosts. In this study, we use genomic distance analyses, phylogenomic networks, tree topology comparisons, and Bayesian inference methods to investigate to what extent HGT has occurred during the evolution of the genus Rickettsia, the effect of the use of different genomic regions in estimating reticulate evolution and HGT events, and the link of these to host range. We show that ecological specialization restricts recombination occurrence in Rickettsia, but other evolutionary processes and genome architecture are also important for the occurrence of HGT. We found that recombination, genomic rearrangements, and genome conservation all show evidence of network-like evolution at whole-genome scale. We show that reticulation occurred mainly, but not only, during the early Rickettsia radiation, and that core proteome genes of every major functional category have experienced reticulated evolution and possibly HGT. Overall, the evolution of Rickettsia bacteria has been tree-like, with evidence of HGT and reticulated evolution for around 10–25% of the core Rickettsia genome. We present evidence of extensive recombination/incomplete lineage sorting (ILS) during the radiation of the genus, probably linked with the emergence of intracellularity in a wide range of hosts.

Romance of the three domains: how cladistics transformed the classification of cellular organisms

Cladistics is a biological philosophy that uses genealogical relationship among species and an inferred sequence of divergence as the basis of classification. This review critically surveys the chronological development of biological classification from Aristotle through our postgenomic era with a central focus on cladistics. In 1957, Julian Huxley coined cladogenesis to denote splitting from subspeciation. In 1960, the English translation of Willi Hennig’s 1950 work, Systematic Phylogenetics, was published, which received strong opposition from pheneticists, such as numerical taxonomists Peter Sneath and Robert Sokal, and evolutionary taxonomist, Ernst Mayr, and sparked acrimonious debates in 1960–1980. In 1977–1990, Carl Woese pioneered in using small subunit rRNA gene sequences to delimitate the three domains of cellular life and established major prokaryotic phyla. Cladistics has since dominated taxonomy. Despite being compatible with modern microbiological observations, i.e. organisms with unusual phenotypes, restricted expression of characteristics and occasionally being uncultivable, increasing recognition of pervasiveness and abundance of horizontal gene transfer has challenged relevance and validity of cladistics. The mosaic nature of eukaryotic and prokaryotic genomes was also gradually discovered. In the mid-2000s, high-throughput and whole-genome sequencing became routine and complex geneologies of organisms have led to the proposal of a reticulated web of life. While genomics only indirectly leads to understanding of functional adaptations to ecological niches, computational modeling of entire organisms is underway and the gap between genomics and phenetics may soon be bridged. Controversies are not expected to settle as taxonomic classifications shall remain subjective to serve the human scientist, not the classified.

How microbiology helps define the rhizome of life

In contrast to the tree of life (TOF) theory, species are mosaics of gene sequences with different origins. Observations of the extensive lateral sequence transfers in all organisms have demonstrated that the genomes of all life forms are collections of genes with different evolutionary histories that cannot be represented by a single TOF. Moreover, genes themselves commonly have several origins due to recombination. The human genome is not free from recombination events, so it is a mosaic like other organisms' genomes. Recent studies have demonstrated evidence for the integration of parasitic DNA into the human genome. Lateral transfer events have been accepted as major contributors of genome evolution in free-living bacteria. Furthermore, the accumulation of genomic sequence data provides evidence for extended genetic exchanges in intracellular bacteria and suggests that such events constitute an agent that promotes and maintains all bacterial species. Archaea and viruses also form chimeras containing primarily bacterial but also eukaryotic sequences. In addition to lateral transfers, orphan genes are indicative of the fact that gene creation is a permanent and unsettled phenomenon. Currently, a rhizome may more adequately represent the multiplicity and de novo creation of a genome. We wanted to confirm that the term “rhizome” in evolutionary biology applies to the entire cellular life history. This view of evolution should resemble a clump of roots representing the multiple origins of the repertoires of the genes of each species.

Evolutionary patterns of carbohydrate transport and metabolism in Halomonas boliviensis as derived from its genome sequence: influences on polyester production

BACKGROUND

Halomonas boliviensis is a halophilic bacterium that is included in the γ-Proteobacteria sub-group, and is able to assimilate different types of carbohydrates. H. boliviensis is also able to produce poly(3-hydroxybutyrate) (PHB) in high yields using glucose as the carbon precursor. Accumulation of PHB by microorganisms is induced by excess of intracellular NADH.The genome sequences and organization in microorganisms should be the result of evolution and adaptation influenced by mutation, gene duplication, horizontal gen transfer (HGT) and recombination. Furthermore, the nearly neutral theory of evolution sustains that genetic modification of DNA could be neutral or selected, albeit most mutations should be at the border between neutrality and selection, i.e. slightly deleterious base substitutions in DNA are followed by a slightly advantageous substitutions.

RESULTS

This article reports the genome sequence of H. boliviensis. The chromosome size of H. boliviensis was 4 119 979 bp, and contained 3 863 genes. A total of 160 genes of H. boliviensis were related to carbohydrate transport and metabolism, and were organized as: 70 genes for metabolism of carbohydrates; 47 genes for ABC transport systems and 43 genes for TRAP-type C4-dicarboxylate transport systems. Protein sequences of H. boliviensis related to carbohydrate transport and metabolism were selected from clusters of orthologous proteins (COGs). Similar proteins derived from the genome sequences of other 41 archaea and 59 bacteria were used as reference. We found that most of the 160 genes in H. boliviensis, c.a. 44%, were obtained from other bacteria by horizontal gene transfer, while 13% of the genes were acquired from haloarchaea and thermophilic archaea, only 34% of the genes evolved among Proteobacteria and the remaining genes encoded proteins that did not cluster with any of the proteins obtained from the reference strains. Furthermore, the diversity of the enzymes derived from these genes led to polymorphism in glycolysis and gluconeogenesis. We found further that an optimum ratio of glucose and sucrose in the culture medium of H. boliviensis favored cell growth and PHB production.

CONCLUSIONS

Results obtained in this article depict that most genetic modifications and enzyme polymorphism in the genome of H. boliviensis were mainly influenced by HGT rather than nearly neutral mutations. Molecular adaptation and evolution experienced by H. boliviensis were also a response to environmental conditions such as the type and amount of carbohydrates in its ecological niche. Consequently, the genome evolution of H. boliviensis showed to be strongly influenced by the type of microorganisms, genetic interaction among microbial species and its environment. Such trend should also be experienced by other prokaryotes. A system for PHB production by H. boliviensis that takes into account the evolutionary adaptation of this bacterium to the assimilation of combinations of carbohydrates suggests the feasibility of a bioprocess economically viable and environmentally friendly.

The rhizome of life: what about metazoa?

The increase in huge number of genomic sequences in recent years has contributed to various genetic events such as horizontal gene transfer (HGT), gene duplication and hybridization of species. Among them HGT has played an important role in the genome evolution and was believed to occur only in Bacterial and Archaeal genomes. As a result, genomes were found to be chimeric and the evolution of life was represented in different forms such as forests, networks and species evolution was described more like a rhizome, rather than a tree. However, in the last few years, HGT has also been evidenced in other group such as metazoa (for example in root-knot nematodes, bdelloid rotifers and mammals). In addition to HGT, other genetic events such as transfer by retrotransposons and hybridization between more closely related lineages are also well established. Therefore, in the light of such genetic events, whether the evolution of metazoa exists in the form of a tree, network or rhizome is highly questionable and needs to be determined. In the current review, we will focus on the role of HGT, retrotransposons and hybridization in the metazoan evolution.

Analysis of Porphyra membrane transporters demonstrates gene transfer among photosynthetic eukaryotes and numerous sodium-coupled transport systems

Membrane transporters play a central role in many cellular processes that rely on the movement of ions and organic molecules between the environment and the cell, and between cellular compartments. Transporters have been well characterized in plants and green algae, but little is known about transporters or their evolutionary histories in the red algae. Here we examined 482 expressed sequence tag contigs that encode putative membrane transporters in the economically important red seaweed Porphyra (Bangiophyceae, Rhodophyta). These contigs are part of a comprehensive transcriptome dataset from Porphyra umbilicalis and Porphyra purpurea. Using phylogenomics, we identified 30 trees that support the expected monophyly of red and green algae/plants (i.e. the Plantae hypothesis) and 19 expressed sequence tag contigs that show evidence of endosymbiotic/horizontal gene transfer involving stramenopiles. The majority (77%) of analyzed contigs encode transporters with unresolved phylogenies, demonstrating the difficulty in resolving the evolutionary history of genes. We observed molecular features of many sodium-coupled transport systems in marine algae, and the potential for coregulation of Porphyra transporter genes that are associated with fatty acid biosynthesis and intracellular lipid trafficking. Although both the tissue-specific and subcellular locations of the encoded proteins require further investigation, our study provides red algal gene candidates associated with transport functions and novel insights into the biology and evolution of these transporters.

E value cutoff and eukaryotic genome content phylogenetics

Genome content analysis has been used as a source of phylogenetic information in large prokaryotic tree of life studies. Recently the sequencing of many eukaryotic genomes has allowed for the similar use of genome content analysis for these organisms too. In this communication we examine the utility of genome content analysis for recovering phylogenetic patterns in several eukaryotic groups. By constructing multiple matrices using different e value cutoffs we examine the dynamics of altering the e value cutoff on five eukaryotic genome data sets. Our analysis indicates that the e value cutoff that is used as a criterion in the construction of the genome content matrix is a critical factor in both the accuracy and information content of the analysis. Strikingly, genome content by itself is not a reliable or accurate source of characters for phylogenetic analysis of the taxa in the five data sets we analyzed. We discuss two problems--small genome attraction and genome duplications as being involved in the rather poor performance of genome content data in recovering eukaryotic phylogeny.

Philosophy and evolution: minding the gap between evolutionary patterns and tree-like patterns

Ever since Darwin, the familiar genealogical pattern known as the Tree of Life (TOL) has been prominent in evolutionary thinking and has dominated not only systematics, but also the analysis of the units of evolution. However, recent findings indicate that the evolution of DNA, especially in prokaryotes and such DNA vehicles as viruses and plasmids, does not follow a unique tree-like pattern. Because evolutionary patterns track a greater range of processes than those captured in genealogies, genealogical patterns are in fact only a subset of a broader set of evolutionary patterns. This fact suggests that evolutionists who focus exclusively on genealogical patterns are blocked from providing a significant range of genuine evolutionary explanations. Consequently, we highlight challenges to tree-based approaches, and point the way toward more appropriate methods to study evolution (although we do not present them in technical detail). We argue that there is significant benefit in adopting wider range of models, evolutionary representations, and evolutionary explanations, based on an analysis of the full range of evolutionary processes. We introduce an ecosystem orientation into evolutionary thinking that highlights the importance of "type 1 coalitions" (functionally related units with genetic exchanges, aka "friends with genetic benefits"), "type 2 coalitions" (functionally related units without genetic exchanges), "communal interactions," and "emergent evolutionary properties." On this basis, we seek to promote the study of (especially prokaryotic) evolution with dynamic evolutionary networks, which are less constrained than the TOL, and to provide new ways to analyze an expanded range of evolutionary units (genetic modules, recombined genes, plasmids, phages and prokaryotic genomes, pangenomes, microbial communities) and evolutionary processes. Finally, we discuss some of the conceptual and practical questions raised by such network-based representation.

Genome-wide comparative analysis of phylogenetic trees: the prokaryotic forest of life

Genome-wide comparison of phylogenetic trees is becoming an increasingly common approach in evolutionary genomics, and a variety of approaches for such comparison have been developed. In this article, we present several methods for comparative analysis of large numbers of phylogenetic trees. To compare phylogenetic trees taking into account the bootstrap support for each internal branch, the Boot-Split Distance (BSD) method is introduced as an extension of the previously developed Split Distance method for tree comparison. The BSD method implements the straightforward idea that comparison of phylogenetic trees can be made more robust by treating tree splits differentially depending on the bootstrap support. Approaches are also introduced for detecting tree-like and net-like evolutionary trends in the phylogenetic Forest of Life (FOL), i.e., the entirety of the phylogenetic trees for conserved genes of prokaryotes. The principal method employed for this purpose includes mapping quartets of species onto trees to calculate the support of each quartet topology and so to quantify the tree and net contributions to the distances between species. We describe the application of these methods to analyze the FOL and the results obtained with these methods. These results support the concept of the Tree of Life (TOL) as a central evolutionary trend in the FOL as opposed to the traditional view of the TOL as a "species tree."

Ecological and temporal constraints in the evolution of bacterial genomes

Studies on the experimental evolution of microorganisms, on their in vivo evolution (mainly in the case of bacteria producing chronic infections), as well as the availability of multiple full genomic sequences, are placing bacteria in the playground of evolutionary studies. In the present article we review the differential contribution to the evolution of bacterial genomes that processes such as gene modification, gene acquisition and gene loss may have when bacteria colonize different habitats that present characteristic ecological features. In particular, we review how the different processes contribute to evolution in microbial communities, in free-living bacteria or in bacteria living in isolation. In addition, we discuss the temporal constraints in the evolution of bacterial genomes, considering bacterial evolution from the perspective of processes of short-sighted evolution and punctual acquisition of evolutionary novelties followed by long stasis periods.

Combined phylogenetic and genomic approaches for the high-throughput study of microbial habitat adaptation

High-throughput sequencing technologies provide new opportunities to address longstanding questions about habitat adaptation in microbial organisms. How have microbes managed to adapt to such a wide range of environments, and what genomic features allow for such adaptation? We review recent large-scale studies of habitat adaptation, with emphasis on those that utilize phylogenetic techniques. On the basis of current trends, we summarize methodological challenges faced by investigators, and the tools, techniques and analytical approaches available to overcome them. Phylogenetic approaches and detailed information about each environmental sample will be crucial as the ability to collect genome sequences continues to expand.

A rooted net of life

Abstract

Phylogenetic reconstruction using DNA and protein sequences has allowed the reconstruction of evolutionary histories encompassing all life. We present and discuss a means to incorporate much of this rich narrative into a single model that acknowledges the discrete evolutionary units that constitute the organism. Briefly, this Rooted Net of Life genome phylogeny is constructed around an initial, well resolved and rooted tree scaffold inferred from a supermatrix of combined ribosomal genes. Extant sampled ribosomes form the leaves of the tree scaffold. These leaves, but not necessarily the deeper parts of the scaffold, can be considered to represent a genome or pan-genome, and to be associated with members of other gene families within that sequenced (pan)genome. Unrooted phylogenies of gene families containing four or more members are reconstructed and superimposed over the scaffold. Initially, reticulations are formed where incongruities between topologies exist. Given sufficient evidence, edges may then be differentiated as those representing vertical lines of inheritance within lineages and those representing horizontal genetic transfers or endosymbioses between lineages.

Reviewers

W. Ford Doolittle, Eric Bapteste and Robert Beiko.

The Public Goods Hypothesis for the evolution of life on Earth

It is becoming increasingly difficult to reconcile the observed extent of horizontal gene transfers with the central metaphor of a great tree uniting all evolving entities on the planet. In this manuscript we describe the Public Goods Hypothesis and show that it is appropriate in order to describe biological evolution on the planet. According to this hypothesis, nucleotide sequences (genes, promoters, exons, etc.) are simply seen as goods, passed from organism to organism through both vertical and horizontal transfer. Public goods sequences are defined by having the properties of being largely non-excludable (no organism can be effectively prevented from accessing these sequences) and non-rival (while such a sequence is being used by one organism it is also available for use by another organism). The universal nature of genetic systems ensures that such non-excludable sequences exist and non-excludability explains why we see a myriad of genes in different combinations in sequenced genomes. There are three features of the public goods hypothesis. Firstly, segments of DNA are seen as public goods, available for all organisms to integrate into their genomes. Secondly, we expect the evolution of mechanisms for DNA sharing and of defense mechanisms against DNA intrusion in genomes. Thirdly, we expect that we do not see a global tree-like pattern. Instead, we expect local tree-like patterns to emerge from the combination of a commonage of genes and vertical inheritance of genomes by cell division. Indeed, while genes are theoretically public goods, in reality, some genes are excludable, particularly, though not only, when they have variant genetic codes or behave as coalition or club goods, available for all organisms of a coalition to integrate into their genomes, and non-rival within the club. We view the Tree of Life hypothesis as a regionalized instance of the Public Goods hypothesis, just like classical mechanics and euclidean geometry are seen as regionalized instances of quantum mechanics and Riemannian geometry respectively. We argue for this change using an axiomatic approach that shows that the Public Goods hypothesis is a better accommodation of the observed data than the Tree of Life hypothesis.

Of woods and webs: possible alternatives to the tree of life for studying genomic fluidity in E. coli

Background

We introduce several forest-based and network-based methods for exploring microbial evolution, and apply them to the study of thousands of genes from 30 strains of E. coli. This case study illustrates how additional analyses could offer fast heuristic alternatives to standard tree of life (TOL) approaches.

Results

We use gene networks to identify genes with atypical modes of evolution, and genome networks to characterize the evolution of genetic partnerships between E. coli and mobile genetic elements. We develop a novel polychromatic quartet method to capture patterns of recombination within E. coli, to update the clanistic toolkit, and to search for the impact of lateral gene transfer and of pathogenicity on gene evolution in two large forests of trees bearing E. coli. We unravel high rates of lateral gene transfer involving E. coli (about 40% of the trees under study), and show that both core genes and shell genes of E. coli are affected by non-tree-like evolutionary processes. We show that pathogenic lifestyle impacted the structure of 30% of the gene trees, and that pathogenic strains are more likely to transfer genes with one another than with non-pathogenic strains. In addition, we propose five groups of genes as candidate mobile modules of pathogenicity. We also present strong evidence for recent lateral gene transfer between E. coli and mobile genetic elements.

Conclusions

Depending on which evolutionary questions biologists want to address (i.e. the identification of modules, genetic partnerships, recombination, lateral gene transfer, or genes with atypical evolutionary modes, etc.), forest-based and network-based methods are preferable to the reconstruction of a single tree, because they provide insights and produce hypotheses about the dynamics of genome evolution, rather than the relative branching order of species and lineages. Such a methodological pluralism - the use of woods and webs - is to be encouraged to analyse the evolutionary processes at play in microbial evolution.

This manuscript was reviewed by: Ford Doolittle, Tal Pupko, Richard Burian, James McInerney, Didier Raoult, and Yan Boucher

Computational methods for Gene Orthology inference

Accurate inference of orthologous genes is a pre-requisite for most comparative genomics studies, and is also important for functional annotation of new genomes. Identification of orthologous gene sets typically involves phylogenetic tree analysis, heuristic algorithms based on sequence conservation, synteny analysis, or some combination of these approaches. The most direct tree-based methods typically rely on the comparison of an individual gene tree with a species tree. Once the two trees are accurately constructed, orthologs are straightforwardly identified by the definition of orthology as those homologs that are related by speciation, rather than gene duplication, at their most recent point of origin. Although ideal for the purpose of orthology identification in principle, phylogenetic trees are computationally expensive to construct for large numbers of genes and genomes, and they often contain errors, especially at large evolutionary distances. Moreover, in many organisms, in particular prokaryotes and viruses, evolution does not appear to have followed a simple 'tree-like' mode, which makes conventional tree reconciliation inapplicable. Other, heuristic methods identify probable orthologs as the closest homologous pairs or groups of genes in a set of organisms. These approaches are faster and easier to automate than tree-based methods, with efficient implementations provided by graph-theoretical algorithms enabling comparisons of thousands of genomes. Comparisons of these two approaches show that, despite conceptual differences, they produce similar sets of orthologs, especially at short evolutionary distances. Synteny also can aid in identification of orthologs. Often, tree-based, sequence similarity- and synteny-based approaches can be combined into flexible hybrid methods.