Horizontal Gene Transfer and the History of Life

Microbes acquire DNA from a variety of sources. The last decades, which have seen the development of genome sequencing, have revealed that horizontal gene transfer has been a major evolutionary force that has constantly reshaped genomes throughout evolution. However, because the history of life must ultimately be deduced from gene phylogenies, the lack of methods to account for horizontal gene transfer has thrown into confusion the very concept of the tree of life. As a result, many questions remain open, but emerging methodological developments promise to use information conveyed by horizontal gene transfer that remains unexploited today.The discovery of the existence of prokaryotic microbes dates back more than 300 years. Since then, our picture of our distant microscopic relatives has undergone several revolutions: from being the living “proofs” of the existence of spontaneous generation, they became later the “archaic” representatives of our distant ancestors, to finally be legitimately recognized as exceptionally diverse organisms, keystone to any ecosystem, including the most familiar and the most hostile environments on Earth. Similarly, although they were first seen as elementary and unbreakable bricks of life, they are now seen as genetically composite bodies, heavyweight champions of “gene robbery.” The most recent of these revolutions has indeed been the realization of their unparalleled ability to integrate genetic material coming from more or less evolutionarily distant organisms. This mechanism is called “horizontal gene transfer” as opposed to vertical transmission from mother to daughter cell.     

Evolution of Pentameric Ligand-Gated Ion Channels: Pro-Loop Receptors

Pentameric ligand-gated ion channels (pLGICs) are ubiquitous neurotransmitter receptors in Bilateria, with a small number of known prokaryotic homologues. Here we describe a new inventory and phylogenetic analysis of pLGIC genes across all kingdoms of life. Our main finding is a set of pLGIC genes in unicellular eukaryotes, some of which are metazoan-like Cys-loop receptors, and others devoid of Cys-loop cysteines, like their prokaryotic relatives. A number of such “Cys-less” receptors also appears in invertebrate metazoans. Together, those findings draw a new distribution of pLGICs in eukaryotes. A broader distribution of prokaryotic channels also emerges, including a major new archaeal taxon, Thaumarchaeota. More generally, pLGICs now appear nearly ubiquitous in major taxonomic groups except multicellular plants and fungi. However, pLGICs are sparsely present in unicellular taxa, suggesting a high rate of gene loss and a non-essential character, contrasting with their essential role as synaptic receptors of the bilaterian nervous system. Multiple alignments of these highly divergent sequences reveal a small number of conserved residues clustered at the interface between the extracellular and transmembrane domains. Only the “Cys-loop” proline is absolutely conserved, suggesting the more fitting name “Pro loop” for that motif, and “Pro-loop receptors” for the superfamily. The infered molecular phylogeny shows a Cys-loop and a Cys-less clade in eukaryotes, both containing metazoans and unicellular members. This suggests new hypotheses on the evolutionary history of the superfamily, such as a possible origin of the Cys-loop cysteines in an ancient unicellular eukaryote. Deeper phylogenetic relationships remain uncertain, particularly around the split between bacteria, archaea, and eukaryotes.     

Computational and phylogenetic validation of nematode horizontal gene transfer

Sequencing of expressed genes has shown that nematodes, particularly the plant-parasitic nematodes, have genes purportedly acquired from other kingdoms by horizontal gene transfer. The prevailing orthodoxy is that such transfer has been a driving force in the evolution of niche specificity, and a recent paper in BMC Evolutionary Biology that presents a detailed phylogenetic analysis of cellulase genes in the free-living nematode Pristionchus pacificus at the species, genus and family levels substantiates this hypothesis.     

Recurrent horizontal transfer of bacterial toxin genes to eukaryotes

In this work, we report likely recurrent horizontal (lateral) gene transfer events of genes encoding pore-forming toxins of the aerolysin family between species belonging to different kingdoms of life. Clustering based on pairwise similarity and phylogenetic analysis revealed several distinct aerolysin sequence groups, each containing proteins from multiple kingdoms of life. These results strongly support at least six independent transfer events between distantly related phyla in the evolutionary history of one protein family and discount selective retention of ancestral genes as a plausible explanation for this patchy phylogenetic distribution. We discuss the possible roles of these proteins and show evidence for a convergent new function in two extant species. We hypothesize that certain gene families are more likely to be maintained following horizontal gene transfer from commensal or pathogenic organism to its host if they 1) can function alone; and 2) are immediately beneficial for the ecology of the organism, as in the case of pore-forming toxins which can be utilized in multicellular organisms for defense and predation.     

Structural Evolution of the Protein Kinase–Like Superfamily

The protein kinase family is large and important, but it is only one family in a larger superfamily of homologous kinases that phosphorylate a variety of substrates and play important roles in all three superkingdoms of life. We used a carefully constructed structural alignment of selected kinases as the basis for a study of the structural evolution of the protein kinase–like superfamily. The comparison of structures revealed a “universal core” domain consisting only of regions required for ATP binding and the phosphotransfer reaction. Remarkably, even within the universal core some kinase structures display notable changes, while still retaining essential activity. Hence, the protein kinase–like superfamily has undergone substantial structural and sequence revision over long evolutionary timescales. We constructed a phylogenetic tree for the superfamily using a novel approach that allowed for the combination of sequence and structure information into a unified quantitative analysis. When considered against the backdrop of species distribution and other metrics, our tree provides a compelling scenario for the development of the various kinase families from a shared common ancestor. We propose that most of the so-called “atypical kinases” are not intermittently derived from protein kinases, but rather diverged early in evolution to form a distinct phyletic group. Within the atypical kinases, the aminoglycoside and choline kinase families appear to share the closest relationship. These two families in turn appear to be the most closely related to the protein kinase family. In addition, our analysis suggests that the actin-fragmin kinase, an atypical protein kinase, is more closely related to the phosphoinositide-3 kinase family than to the protein kinase family. The two most divergent families, α-kinases and phosphatidylinositol phosphate kinases (PIPKs), appear to have distinct evolutionary histories. While the PIPKs probably have an evolutionary relationship with the rest of the kinase superfamily, the relationship appears to be very distant (and perhaps indirect). Conversely, the α-kinases appear to be an exception to the scenario of early divergence for the atypical kinases: they apparently arose relatively recently in eukaryotes. We present possible scenarios for the derivation of the α-kinases from an extant kinase fold.     

Genetic exchange between kingdoms.

Bacterial conjugation with two evolutionarily divergent yeasts has been observed in the laboratory. Whether such trans-kingdom conjugation events, other than the well known Agrobacterium-plant cell interaction, actually occur in nature is not known. However, a few putative events have recently been uncovered by gene (or protein) sequence analysis, suggesting that horizontal gene transfer between phylogenetic kingdoms may be a real phenomenon.     

Serpins in prokaryotes

Members of the serpin (serine proteinase inhibitor) superfamily have been identified in higher multicellular eukaryotes (plants and animals) and viruses but not in bacteria, archaea, or fungi. Thus, the ancestral serpin and the origin of the serpin inhibitory mechanism remain obscure. In this study we characterize 12 serpin-like sequences in the genomes of prokaryotic organisms, extending this protein family to all major branches of life. Notably, these organisms live in dramatically different environments and some are evolutionarily distantly related. A sequence-based analysis suggests that all 12 serpins are inhibitory. Despite considerable sequence divergence between the proteins, in four of the 12 sequences the region of the serpin that determines proteinase specificity is highly conserved, indicating that these inhibitors are likely to share a common target. Inhibitory serpins are typically prone to polymerization upon heating; thus, the existence of serpins in the moderate thermophilic bacterium Thermobifida fusca, the thermophilic bacterium Thermoanaerobacter tengcongensis, and the hyperthermophilic archaeon Pyrobaculum aerophilum is of particular interest. Using molecular modeling, we predict the means by which heat stability in the latter protein may be achieved without compromising inhibitory activity.     

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ‐stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNA^Tyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences “U‐U‐x2‐C” and “U‐x‐G‐x2‐T,” thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.     

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the “scale” of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.     

Bryn Bridges and mutagenesis: exploring the intellectual space

The products of the SOS-regulated umuDC genes are required for most UV and chemical mutagenesis in Escherichia coli. Recently it has been recognized that UmuC is the founding member of a superfamily of novel DNA polymerases found in all three kingdoms of life. Key findings leading to these insights are reviewed, placing a particular emphasis on contributions made by Bryn Bridges and on his interest in the importance of interactions between the umuDC gene products and the replicative DNA polymerase.     

Bayesian Analysis of Congruence of Core Genes in Prochlorococcus and Synechococcus and Implications on Horizontal Gene Transfer

It is often suggested that horizontal gene transfer is so ubiquitous in microbes that the concept of a phylogenetic tree representing the pattern of vertical inheritance is oversimplified or even positively misleading. “Universal proteins” have been used to infer the organismal phylogeny, but have been criticized as being only the “tree of one percent.” Currently, few options exist for those wishing to rigorously assess how well a universal protein phylogeny, based on a relative handful of well-conserved genes, represents the phylogenetic histories of hundreds of genes. Here, we address this problem by proposing a visualization method and a statistical test within a Bayesian framework. We use the genomes of marine cyanobacteria, a group thought to exhibit substantial amounts of HGT, as a test case. We take 379 orthologous gene families from 28 cyanobacteria genomes and estimate the Bayesian posterior distributions of trees – a “treecloud” – for each, as well as for a concatenated dataset based on putative “universal proteins.” We then calculate the average distance between trees within and between all treeclouds on various metrics and visualize this high-dimensional space with non-metric multidimensional scaling (NMMDS). We show that the tree space is strongly clustered and that the universal protein treecloud is statistically significantly closer to the center of this tree space than any individual gene treecloud. We apply several commonly-used tests for incongruence/HGT and show that they agree HGT is rare in this dataset, but make different choices about which genes were subject to HGT. Our results show that the question of the representativeness of the “tree of one percent” is a quantitative empirical question, and that the phylogenetic central tendency is a meaningful observation even if many individual genes disagree due to the various sources of incongruence.     

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome‐scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence‐level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.     

Ancient gene duplications and the root(s) of the tree of life

Summary. Tracing organismal histories on the timescale of the tree of life remains one of the challenging tasks in evolutionary biology. The hotly debated questions include the evolutionary relationship between the three domains of life (e.g., which of the three domains are sister domains, are the domains para-, poly-, or monophyletic) and the location of the root within the universal tree of life. For the latter, many different points of view have been considered but so far no consensus has been reached. The only widely accepted rationale to root the universal tree of life is to use anciently duplicated paralogous genes that are present in all three domains of life. To date only few anciently duplicated gene families useful for phylogenetic reconstruction have been identified. Here we present results from a systematic search for ancient gene duplications using twelve representative, completely sequenced, archaeal and bacterial genomes. Phylogenetic analyses of identified cases show that the majority of datasets support a root between Archaea and Bacteria; however, some datasets support alternative hypotheses, and all of them suffer from a lack of strong phylogenetic signal. The results are discussed with respect to the impact of horizontal gene transfer on the ability to reconstruct organismal evolution. The exchange of genetic information between divergent organisms gives rise to mosaic genomes, where different genes in a genome have different histories. Simulations show that even low rates of horizontal gene transfer dramatically complicate the reconstruction of organismal evolution, and that the different most recent common molecular ancestors likely existed at different times and in different lineages. Correspondence and reprints: Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT 06269-3125, U.S.A. Present address: Genome Atlantic, Department of Biochemistry and Molecular Biology, Dalhousie University, Halifax, Nova Scotia, Canada.     

Uniform binding and negative catalysis at the origin of enzymes

Enzymes are well known for their catalytic abilities, some even reaching “catalytic perfection” in the sense that the reaction they catalyze has reached the physical bound of the diffusion rate. However, our growing understanding of enzyme superfamilies has revealed that only some share a catalytic chemistry while others share a substrate‐handle binding motif, for example, for a particular phosphate group. This suggests that some families emerged through a “substrate‐handle‐binding‐first” mechanism (“binding‐first” for brevity) instead of “chemistry‐first” and we are, therefore, left to wonder what the role of non‐catalytic binders might have been during enzyme evolution. In the last of their eight seminal, back‐to‐back articles from 1976, John Albery and Jeremy Knowles addressed the question of enzyme evolution by arguing that the simplest mode of enzyme evolution is what they defined as “uniform binding” (parallel stabilization of all enzyme‐bound states to the same degree). Indeed, we show that a uniform‐binding proto‐catalyst can accelerate a reaction, but only when catalysis is already present, that is, when the transition state is already stabilized to some degree. Thus, we sought an alternative explanation for the cases where substrate‐handle‐binding preceded any involvement of a catalyst. We find that evolutionary starting points that exhibit negative catalysis can redirect the reaction''s course to a preferred product without need for rate acceleration or product release; that is, if they do not stabilize, or even destabilize, the transition state corresponding to an undesired product. Such a mechanism might explain the emergence of “binding‐first” enzyme families like the aldolase superfamily.     

Conservation of Functionally Important Global Motions in an Enzyme Superfamily across Varying Quaternary Structures

The α‐d‐phosphohexomutase superfamily comprises enzymes involved in carbohydrate metabolism that are found in all kingdoms of life. Recent biophysical studies have shown for the first time that several of these enzymes exist as dimers in solution, prompting an examination of the oligomeric state of all proteins of known structure in the superfamily (11 different proteins; 31 crystal structures) via computational and experimental analyses. We find that these proteins range in quaternary structure from monomers to tetramers, with 6 of the 11 known structures being likely oligomers. The oligomeric state of these proteins not only is associated in some cases with enzyme subgroup (i.e., substrate specificity) but also appears to depend on domain of life, with the two archaeal proteins existing as higher‐order oligomers. Within the oligomers, three distinct interfaces are observed, one of which is found in both archaeal and bacterial proteins. Normal mode analysis shows that the topological arrangement of the oligomers permits domain 4 of each protomer to move independently as required for catalysis. Our analysis suggests that the advantages associated with protein flexibility in this enzyme family are of sufficient importance to be maintained during the evolution of multiple independent oligomers. This study is one of the first showing that global motions may be conserved not only within protein families but also across members of a superfamily with varying oligomeric structures.     

New Aspects on Lanosterol 14α-Demethylase and Cytochrome P450 Evolution: Lanosterol/Cycloartenol Diversification and Lateral Transfer

Sterol 14-demethylase (CYP51) is a member of the cytochrome P450 superfamily, widely found in animals, fungi, and plants but present in few prokaryotic groups. CYP51 is currently believed to be the ancestral cytochrome P450 that has been transferred from prokaryotes to eukaryotic kingdoms. We propose an alternate view of CYP51 evolution that has an impact on understanding the evolution of the entire CYP superfamily. Two hundred forty-nine bacterial and four archaeal CYP sequences have been aligned and a bacterial CYP tree designed, showing a separation of two branches. Prokaryotic CYP51s cluster to the minor branch, together with other eukaryote-like CYPs. Mycobacterial and methylococcal CYP51s cluster together (100% bootstrap probability), while Streptomyces CYP51 remains on a distant branch. A CYP51 phylogenetic tree has been constructed from 44 sequences resulting in a ((plant, bacteria),(animal, fungi)) topology (100% bootstrap probability). This is in accordance with the lanosterol/cycloartenol diversification of sterol biosynthesis. The lanosterol branch (nonphotosynthetic lineage) follows the previously proposed topology of animal and fungal orthologues (100% bootstrap probability), while plant and D. discoideum CYP51s belong to the cycloartenol branch (photosynthetic lineage), all in accordance with biochemical data. Bacterial CYP51s cluster within the cycloartenol branch (69% bootstrap probability), which is indicative of a lateral gene transfer of a plant CYP51 to the methylococcal/mycobacterial progenitor, suggesting further that bacterial CYP51s are not the oldest CYP genes. Lateral gene transfer is likely far more important than hitherto thought in the development of the diversified CYP superfamily. Consequently, bacterial CYPs may represent a mixture of genes with prokaryotic and eukaryotic origin.     

Recent emergence of the modern genetic code: a proposal

This article proposes that the genetic code was not fully formed before the divergence of life into three kingdoms. Rather, at least arginine and tryptophan evolved after the diversification of archaea, bacteria and eukaryotes, and were spread by horizontal gene transfer. Evidence for this hypothesis is based on data suggesting that enzymes for biosynthesis of arginine and tryptophan, and for arginine tRNA ligase, have shorter divergence times than the underlying lineages. Also, many of these genes display "star" phylogenies. This proposal is an extension of the idea that the genetic code was unified because of the evolutionary pressure from horizontal gene transfer. These considerations further undermine the need to postulate the existence of a "last common ancestor"; a simpler model would be that multiple lineages gave rise to life today.     

Serpins: structure,function and molecular evolution

The superfamily of serine proteinase inhibitors (serpins) are involved in a number of fundamental biological processes such as blood coagulation, complement activation, fibrinolysis, angiogenesis, inflammation and tumor suppression and are expressed in a cell-specific manner. The average protein size of a serpin family member is 350-400 amino acids, but gene structure varies in terms of number and size of exons and introns. Previous studies of all known serpins identified 16 clades and 10 orphan sequences. Vertebrate serpins can be conveniently classified into six sub-groups. We provide additional data that updates the phylogenetic analysis in the context of structural and functional properties of the proteins. From these, we can conclude that the functional classification of serpins relies on their protein structure and not on sequence similarity.     

Insights into the Phylogeny or Arylamine N-Acetyltransferases in Fungi

Previous studies have shown that Eumycetes fungi can acylate arylamine thanks to arylamine N-acetyltransferases, xenobiotic-metabolizing enzymes also found in animals and bacteria. In this article, we present the results of mining 96 available fungal genome sequences for arylamine N-acetyltransferase genes and propose their phylogeny. The filamentous Pezizomycotina are shown to possess many putative N-acetyltransferases, whilst these are often lacking in other fungal groups. The evolution of the N-acetyltransferases is best explained by the presence of at least one gene in the opisthokont ancestor of the fungi and animal kingdoms, followed by recurrent gene losses and gene duplications. A possible horizontal gene transfer event may have occurred from bacteria to the basidiomycetous yeast Malassezia globosa.     

Horizontal Transfer of a Subtilisin Gene from Plants into an Ancestor of the Plant Pathogenic Fungal Genus Colletotrichum

The genus Colletotrichum contains a large number of phytopathogenic fungi that produce enormous economic losses around the world. The effect of horizontal gene transfer (HGT) has not been studied yet in these organisms. Inter-Kingdom HGT into fungal genomes has been reported in the past but knowledge about the HGT between plants and fungi is particularly limited. We describe a gene in the genome of several species of the genus Colletotrichum with a strong resemblance to subtilisins typically found in plant genomes. Subtilisins are an important group of serine proteases, widely distributed in all of the kingdoms of life. Our hypothesis is that the gene was acquired by Colletotrichum spp. through (HGT) from plants to a Colletotrichum ancestor. We provide evidence to support this hypothesis in the form of phylogenetic analyses as well as a characterization of the similarity of the subtilisin at the primary, secondary and tertiary structural levels. The remarkable level of structural conservation of Colletotrichum plant-like subtilisin (CPLS) with plant subtilisins and the differences with the rest of Colletotrichum subtilisins suggests the possibility of molecular mimicry. Our phylogenetic analysis indicates that the HGT event would have occurred approximately 150–155 million years ago, after the divergence of the Colletotrichum lineage from other fungi. Gene expression analysis shows that the gene is modulated during the infection of maize by C. graminicola suggesting that it has a role in plant disease. Furthermore, the upregulation of the CPLS coincides with the downregulation of several plant genes encoding subtilisins. Based on the known roles of subtilisins in plant pathogenic fungi and the gene expression pattern that we observed, we postulate that the CPLSs have an important role in plant infection.