The root of the universal tree and the origin of eukaryotes based on elongation factor phylogeny

The emerging view on the origin and early evolution of eukaryotic cells

The origin of the eukaryotic cell, with its compartmentalized nature and generally large size compared with bacterial and archaeal cells, represents a cornerstone event in the evolution of complex life on Earth. In a process referred to as eukaryogenesis, the eukaryotic cell is believed to have evolved between approximately 1.8 and 2.7 billion years ago from its archaeal ancestors, with a symbiosis with a bacterial (proto-mitochondrial) partner being a key event. In the tree of life, the branch separating the first from the last common ancestor of all eukaryotes is long and lacks evolutionary intermediates. As a result, the timing and driving forces of the emergence of complex eukaryotic features remain poorly understood. During the past decade, environmental and comparative genomic studies have revealed vital details about the identity and nature of the host cell and the proto-mitochondrial endosymbiont, enabling a critical reappraisal of hypotheses underlying the symbiotic origin of the eukaryotic cell. Here we outline our current understanding of the key players and events underlying the emergence of cellular complexity during the prokaryote-to-eukaryote transition and discuss potential avenues of future research that might provide new insights into the enigmatic origin of the eukaryotic cell.

Socialized mitochondria: mitonuclear crosstalk in stress

Traditionally, mitochondria are considered sites of energy production. However, recent studies have suggested that mitochondria are signaling organelles that are involved in intracellular interactions with other organelles. Remarkably, stressed mitochondria appear to induce a beneficial response that restores mitochondrial function and cellular homeostasis. These mitochondrial stress-centered signaling pathways have been rapidly elucidated in multiple organisms. In this review, we examine current perspectives on how mitochondria communicate with the rest of the cell, highlighting mitochondria-to-nucleus (mitonuclear) communication under various stresses. Our understanding of mitochondria as signaling organelles may provide new insights into disease susceptibility and lifespan extension.

Perspectives on Principles of Cellular Behavior from the Biophysics of Protists

Cells are the fundamental unit of biological organization. Although it may be easy to think of them as little more than the simple building blocks of complex organisms such as animals, single cells are capable of behaviors of remarkable apparent sophistication. This is abundantly clear when considering the diversity of form and function among the microbial eukaryotes, the protists. How might we navigate this diversity in the search for general principles of cellular behavior? Here, we review cases in which the intensive study of protists from the perspective of cellular biophysics has driven insight into broad biological questions of morphogenesis, navigation and motility, and decision making. We argue that applying such approaches to questions of evolutionary cell biology presents rich, emerging opportunities. Integrating and expanding biophysical studies across protist diversity, exploiting the unique characteristics of each organism, will enrich our understanding of general underlying principles.

A life off the beaten track in biomechanics: Imperfect elasticity, cytoskeletal glassiness, and epithelial unjamming

Textbook descriptions of elasticity, viscosity, and viscoelasticity fail to account for certain mechanical behaviors that typify soft living matter. Here, we consider three examples. First, strong empirical evidence suggests that within lung parenchymal tissues, the frictional stresses expressed at the microscale are fundamentally not of viscous origin. Second, the cytoskeleton (CSK) of the airway smooth muscle cell, as well as that of all eukaryotic cells, is more solid-like than fluid-like, yet its elastic modulus is softer than the softest of soft rubbers by a factor of 104-105. Moreover, the eukaryotic CSK expresses power law rheology, innate malleability, and fluidization when sheared. For these reasons, taken together, the CSK of the living eukaryotic cell is reminiscent of the class of materials called soft glasses, thus likening it to inert materials such as clays, pastes slurries, emulsions, and foams. Third, the cellular collective comprising a confluent epithelial layer can become solid-like and jammed, fluid-like and unjammed, or something in between. Esoteric though each may seem, these discoveries are consequential insofar as they impact our understanding of bronchospasm and wound healing as well as cancer cell invasion and embryonic development. Moreover, there are reasons to suspect that certain of these phenomena first arose in the early protist as a result of evolutionary pressures exerted by the primordial microenvironment. We have hypothesized, further, that each then became passed down virtually unchanged to the present day as a conserved core process. These topics are addressed here not only because they are interesting but also because they track the journey of one laboratory along a path less traveled by.

The universal Sua5/TsaC family evolved different mechanisms for the synthesis of a key tRNA modification

TsaC/Sua5 family of enzymes catalyzes the first step in the synthesis of N6-threonyl-carbamoyl adenosine (t⁶A) one of few truly ubiquitous tRNA modifications important for translation accuracy. TsaC is a single domain protein while Sua5 proteins contains a TsaC-like domain and an additional SUA5 domain of unknown function. The emergence of these two proteins and their respective mechanisms for t⁶A synthesis remain poorly understood. Here, we performed phylogenetic and comparative sequence and structure analysis of TsaC and Sua5 proteins. We confirm that this family is ubiquitous but the co-occurrence of both variants in the same organism is rare and unstable. We further find that obligate symbionts are the only organisms lacking sua5 or tsaC genes. The data suggest that Sua5 was the ancestral version of the enzyme while TsaC arose via loss of the SUA5 domain that occurred multiple times in course of evolution. Multiple losses of one of the two variants in combination with horizontal gene transfers along a large range of phylogenetic distances explains the present day patchy distribution of Sua5 and TsaC. The loss of the SUA5 domain triggered adaptive mutations affecting the substrate binding in TsaC proteins. Finally, we identified atypical Sua5 proteins in Archaeoglobi archaea that seem to be in the process of losing the SUA5 domain through progressive gene erosion. Together, our study uncovers the evolutionary path for emergence of these homologous isofunctional enzymes and lays the groundwork for future experimental studies on the function of TsaC/Sua5 proteins in maintaining faithful translation.

Phylogenomic Analyses of 2,786 Genes in 158 Lineages Support a Root of the Eukaryotic Tree of Life between Opisthokonts and All Other Lineages

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely.

Evidence for the Rapid and Divergent Evolution of Mycoplasmas: Structural and Phylogenetic Analysis of Enolases

Mycoplasmas are a group of prokaryotes without cell walls that have evolved through several rounds of degenerative evolution. With a low cell DNA G + C content and definitively long genetic lineages, mycoplasmas are thought to be in a state of rapid evolution. However, little associated evidence has been provided. Enolase is a key enzyme in glycolysis that is widely found in all species from the three domains, and it is evolutionarily conserved. In our previous studies, enolase acted as a virulence factor and participated in cell-surface adhesion in Mycoplasma hyopneumoniae. Furthermore, unique loop regions were first found in the crystal structure of Mhp Eno. Here, enolase structures from Mycoplasma pneumoniae and Mycoplasma bovis were determined. An extra helix 7 is specific and conservatively found in almost all mycoplasma enolases, as confirmed by crystal structures and sequence alignment. Particular motifs for helix 7, which is composed of F-K/G-K-L/F-K-X-A-I, have been proposed and could be regarded as molecular markers. To our surprise, the genetic distances between any two mycoplasma enolases were obviously longer than those between the two corresponding species themselves, indicating divergent evolution of mycoplasma enolases, whereas no horizontal gene transfer was detected in mycoplasma enolase genens. Furthermore, different evolutionary patterns were adopted by different loop regions of mycoplasma enolase. Enolases from different Mycoplasma species also showed different affinities for PLG and fibronectin. Our results indicate the rapid and divergent evolution of mycoplasma enolase and mycoplasmas. This study will also aid understanding the independent evolution of Mycoplasma species after separation from their common ancestor.

Rooting Species Trees Using Gene Tree-Species Tree Reconciliation

Interpreting phylogenetic trees requires a root, which provides the direction of evolution and polarizes ancestor-descendant relationships. But inferring the root using genetic data is difficult, particularly in cases where the closest available outgroup is only distantly related, which are common for microbes. In this chapter, we present a workflow for estimating rooted species trees and the evolutionary history of the gene families that evolve within them using probabilistic gene tree-species tree reconciliation. We illustrate the pipeline using a small dataset of prokaryotic genomes, for which the example scripts can be run using modest computer resources. We describe the rooting method used in this work in the context or other rooting strategies and discuss some of the limitations and opportunities presented by probabilistic gene tree-species tree reconciliation methods.

Immunogenic Proteins of Group B Streptococcus-Potential Antigens in Immunodiagnostic Assay for GBS Detection

Streptococcus agalactiae (Group B Streptococcus, GBS) is an opportunistic pathogen, which asymptomatically colonizes the gastrointestinal and genitourinary tract of up to one third of healthy adults. Nevertheless, GBS carriage in pregnant women may lead to several health issues in newborns causing life threatening infection, such as sepsis, pneumonia or meningitis. Recommended GBS screening in pregnant women significantly reduced morbidity and mortality in infants. Nevertheless, intrapartum antibiotic prophylaxis, recommended following the detection of carriage or in case of lack of a carriage test result for pregnant women who demonstrate certain risk factors, led to the expansion of the adverse phenomenon of bacterial resistance to antibiotics. In our paper, we reviewed some immunogenic GBS proteins, i.e., Alp family proteins, β protein, Lmb, Sip, BibA, FsbA, ScpB, enolase, elongation factor Tu, IMPDH, and GroEL, which possess features characteristic of good candidates for immunodiagnostic assays for GBS carriage detection, such as immunoreactivity and specificity. We assume that they can be used as an alternative diagnostic method to the presently recommended bacteriological cultivation and MALDI.

Inferring the deep past from molecular data

There is an expectation that analyses of molecular sequences might be able to distinguish between alternative hypotheses for ancient relationships, but the phylogenetic methods used and types of data analyzed are of critical importance in any attempt to recover historical signal. Here, we discuss some common issues that can influence the topology of trees obtained when using overly simple models to analyze molecular data that often display complicated patterns of sequence heterogeneity. To illustrate our discussion, we have used three examples of inferred relationships which have changed radically as models and methods of analysis have improved. In two of these examples, the sister-group relationship between thermophilic Thermus and mesophilic Deinococcus, and the position of long-branch Microsporidia among eukaryotes, we show that recovering what is now generally considered to be the correct tree is critically dependent on the fit between model and data. In the third example, the position of eukaryotes in the tree of life, the hypothesis that is currently supported by the best available methods is fundamentally different from the classical view of relationships between major cellular domains. Since heterogeneity appears to be pervasive and varied among all molecular sequence data, and even the best available models can still struggle to deal with some problems, the issues we discuss are generally relevant to phylogenetic analyses. It remains essential to maintain a critical attitude to all trees as hypotheses of relationship that may change with more data and better methods.

Multidomain ribosomal protein trees and the planctobacterial origin of neomura (eukaryotes, archaebacteria)

Palaeontologically, eubacteria are > 3× older than neomura (eukaryotes, archaebacteria). Cell biology contrasts ancestral eubacterial murein peptidoglycan walls and derived neomuran N-linked glycoprotein coats/walls. Misinterpreting long stems connecting clade neomura to eubacteria on ribosomal sequence trees (plus misinterpreted protein paralogue trees) obscured this historical pattern. Universal multiprotein ribosomal protein (RP) trees, more accurate than rRNA trees, are taxonomically undersampled. To reduce contradictions with genically richer eukaryote trees and improve eubacterial phylogeny, we constructed site-heterogeneous and maximum-likelihood universal three-domain, two-domain, and single-domain trees for 143 eukaryotes (branching now congruent with 187-protein trees), 60 archaebacteria, and 151 taxonomically representative eubacteria, using 51 and 26 RPs. Site-heterogeneous trees greatly improve eubacterial phylogeny and higher classification, e.g. showing gracilicute monophyly, that many 'rDNA-phyla' belong in Proteobacteria, and reveal robust new phyla Synthermota and Aquithermota. Monoderm Posibacteria and Mollicutes (two separate wall losses) are both polyphyletic: multiple outer membrane losses in Endobacteria occurred separately from Actinobacteria; neither phylum is related to Chloroflexi, the most divergent prokaryotes, which originated photosynthesis (new model proposed). RP trees support an eozoan root for eukaryotes and are consistent with archaebacteria being their sisters and rooted between Filarchaeota (=Proteoarchaeota, including 'Asgardia') and Euryarchaeota sensu-lato (including ultrasimplified 'DPANN' whose long branches often distort trees). Two-domain trees group eukaryotes within Planctobacteria, and archaebacteria with Planctobacteria/Sphingobacteria. Integrated molecular/palaeontological evidence favours negibacterial ancestors for neomura and all life. Unique presence of key pre-neomuran characters favours Planctobacteria only as ancestral to neomura, which apparently arose by coevolutionary repercussions (explained here in detail, including RP replacement) of simultaneous outer membrane and murein loss. Planctobacterial C-1 methanotrophic enzymes are likely ancestral to archaebacterial methanogenesis and β-propeller-α-solenoid proteins to eukaryotic vesicle coats, nuclear-pore-complexes, and intraciliary transport. Planctobacterial chaperone-independent 4/5-protofilament microtubules and MamK actin-ancestors prepared for eukaryote intracellular motility, mitosis, cytokinesis, and phagocytosis. We refute numerous wrong ideas about the universal tree.

Multidomain ribosomal protein trees and the planctobacterial origin of neomura (eukaryotes, archaebacteria)

Palaeontologically, eubacteria are > 3× older than neomura (eukaryotes, archaebacteria). Cell biology contrasts ancestral eubacterial murein peptidoglycan walls and derived neomuran N-linked glycoprotein coats/walls. Misinterpreting long stems connecting clade neomura to eubacteria on ribosomal sequence trees (plus misinterpreted protein paralogue trees) obscured this historical pattern. Universal multiprotein ribosomal protein (RP) trees, more accurate than rRNA trees, are taxonomically undersampled. To reduce contradictions with genically richer eukaryote trees and improve eubacterial phylogeny, we constructed site-heterogeneous and maximum-likelihood universal three-domain, two-domain, and single-domain trees for 143 eukaryotes (branching now congruent with 187-protein trees), 60 archaebacteria, and 151 taxonomically representative eubacteria, using 51 and 26 RPs. Site-heterogeneous trees greatly improve eubacterial phylogeny and higher classification, e.g. showing gracilicute monophyly, that many 'rDNA-phyla' belong in Proteobacteria, and reveal robust new phyla Synthermota and Aquithermota. Monoderm Posibacteria and Mollicutes (two separate wall losses) are both polyphyletic: multiple outer membrane losses in Endobacteria occurred separately from Actinobacteria; neither phylum is related to Chloroflexi, the most divergent prokaryotes, which originated photosynthesis (new model proposed). RP trees support an eozoan root for eukaryotes and are consistent with archaebacteria being their sisters and rooted between Filarchaeota (=Proteoarchaeota, including 'Asgardia') and Euryarchaeota sensu-lato (including ultrasimplified 'DPANN' whose long branches often distort trees). Two-domain trees group eukaryotes within Planctobacteria, and archaebacteria with Planctobacteria/Sphingobacteria. Integrated molecular/palaeontological evidence favours negibacterial ancestors for neomura and all life. Unique presence of key pre-neomuran characters favours Planctobacteria only as ancestral to neomura, which apparently arose by coevolutionary repercussions (explained here in detail, including RP replacement) of simultaneous outer membrane and murein loss. Planctobacterial C-1 methanotrophic enzymes are likely ancestral to archaebacterial methanogenesis and β-propeller-α-solenoid proteins to eukaryotic vesicle coats, nuclear-pore-complexes, and intraciliary transport. Planctobacterial chaperone-independent 4/5-protofilament microtubules and MamK actin-ancestors prepared for eukaryote intracellular motility, mitosis, cytokinesis, and phagocytosis. We refute numerous wrong ideas about the universal tree.

Evolution of P2A and P5A ATPases: ancient gene duplications and the red algal connection to green plants revisited

In a search for slowly evolving nuclear genes that may cast light on the deep evolution of plants, we carried out phylogenetic analyses of two well-characterized subfamilies of P-type pumps (P2A and P5A ATPases) from representative branches of the eukaryotic tree of life. Both P-type ATPase genes were duplicated very early in eukaryotic evolution and before the divergence of the present eukaryotic supergroups. Synapomorphies identified in the sequences provide evidence that green plants and red algae are more distantly related than are green plants and eukaryotic supergroups in which secondary or tertiary plastids are common, such as several groups belonging to the clade that includes Stramenopiles, Alveolata, Rhizaria, Cryptophyta and Haptophyta (SAR). We propose that red algae branched off soon after the first photosynthesizing eukaryote had acquired a primary plastid, while in another lineage that led to SAR, the primary plastid was lost but, in some cases, regained as a secondary or tertiary plastid.

The Diverse Functional Roles of Elongation Factor Tu (EF-Tu) in Microbial Pathogenesis

Elongation factor thermal unstable Tu (EF-Tu) is a G protein that catalyzes the binding of aminoacyl-tRNA to the A-site of the ribosome inside living cells. Structural and biochemical studies have described the complex interactions needed to effect canonical function. However, EF-Tu has evolved the capacity to execute diverse functions on the extracellular surface of both eukaryote and prokaryote cells. EF-Tu can traffic to, and is retained on, cell surfaces where can interact with membrane receptors and with extracellular matrix on the surface of plant and animal cells. Our structural studies indicate that short linear motifs (SLiMs) in surface exposed, non-conserved regions of the molecule may play a key role in the moonlighting functions ascribed to this ancient, highly abundant protein. Here we explore the diverse moonlighting functions relating to pathogenesis of EF-Tu in bacteria and examine putative SLiMs on surface-exposed regions of the molecule.

What is an archaeon and are the Archaea really unique?

The recognition of the group Archaea as a major branch of the tree of life (ToL) prompted a new view of the evolution of biodiversity. The genomic representation of archaeal biodiversity has since significantly increased. In addition, advances in phylogenetic modeling of multi-locus datasets have resolved many recalcitrant branches of the ToL. Despite the technical advances and an expanded taxonomic representation, two important aspects of the origins and evolution of the Archaea remain controversial, even as we celebrate the 40th anniversary of the monumental discovery. These issues concern (i) the uniqueness (monophyly) of the Archaea, and (ii) the evolutionary relationships of the Archaea to the Bacteria and the Eukarya; both of these are relevant to the deep structure of the ToL. To explore the causes for this persistent ambiguity, I examine multiple datasets and different phylogenetic approaches that support contradicting conclusions. I find that the uncertainty is primarily due to a scarcity of information in standard datasets-universal core-genes datasets-to reliably resolve the conflicts. These conflicts can be resolved efficiently by comparing patterns of variation in the distribution of functional genomic signatures, which are less diffused unlike patterns of primary sequence variation. Relatively lower heterogeneity in distribution patterns minimizes uncertainties and supports statistically robust phylogenetic inferences, especially of the earliest divergences of life. This case study further highlights the limitations of primary sequence data in resolving difficult phylogenetic problems, and raises questions about evolutionary inferences drawn from the analyses of sequence alignments of a small set of core genes. In particular, the findings of this study corroborate the growing consensus that reversible substitution mutations may not be optimal phylogenetic markers for resolving early divergences in the ToL, nor for determining the polarity of evolutionary transitions across the ToL.

Complex Evolutionary History of Translation Elongation Factor 2 and Diphthamide Biosynthesis in Archaea and Parabasalids

Diphthamide is a modified histidine residue which is uniquely present in archaeal and eukaryotic elongation factor 2 (EF-2), an essential GTPase responsible for catalyzing the coordinated translocation of tRNA and mRNA through the ribosome. In part due to the role of diphthamide in maintaining translational fidelity, it was previously assumed that diphthamide biosynthesis genes (dph) are conserved across all eukaryotes and archaea. Here, comparative analysis of new and existing genomes reveals that some archaea (i.e., members of the Asgard superphylum, Geoarchaea, and Korarchaeota) and eukaryotes (i.e., parabasalids) lack dph. In addition, while EF-2 was thought to exist as a single copy in archaea, many of these dph-lacking archaeal genomes encode a second EF-2 paralog missing key residues required for diphthamide modification and for normal translocase function, perhaps suggesting functional divergence linked to loss of diphthamide biosynthesis. Interestingly, some Heimdallarchaeota previously suggested to be most closely related to the eukaryotic ancestor maintain dph genes and a single gene encoding canonical EF-2. Our findings reveal that the ability to produce diphthamide, once thought to be a universal feature in archaea and eukaryotes, has been lost multiple times during evolution, and suggest that anticipated compensatory mechanisms evolved independently.

The Effect of Nonreversibility on Inferring Rooted Phylogenies

Most phylogenetic models assume that the evolutionary process is stationary and reversible. In addition to being biologically improbable, these assumptions also impair inference by generating models under which the likelihood does not depend on the position of the root. Consequently, the root of the tree cannot be inferred as part of the analysis. Yet identifying the root position is a key component of phylogenetic inference because it provides a point of reference for polarizing ancestor-descendant relationships and therefore interpreting the tree. In this paper, we investigate the effect of relaxing the unrealistic reversibility assumption and allowing the position of the root to be another unknown. We propose two hierarchical models that are centered on a reversible model but perturbed to allow nonreversibility. The models differ in the degree of structure imposed on the perturbations. The analysis is performed in the Bayesian framework using Markov chain Monte Carlo methods for which software is provided. We illustrate the performance of the two nonreversible models in analyses of simulated data using two types of topological priors. We then apply the models to a real biological data set, the radiation of polyploid yeasts, for which there is robust biological opinion about the root position. Finally, we apply the models to a second biological alignment for which the rooted tree is controversial: the ribosomal tree of life. We compare the two nonreversible models and conclude that both are useful in inferring the position of the root from real biological data.

Elongation factor Tu is a multifunctional and processed moonlighting protein

Many bacterial moonlighting proteins were originally described in medically, agriculturally, and commercially important members of the low G + C Firmicutes. We show Elongation factor Tu (Ef-Tu) moonlights on the surface of the human pathogens Staphylococcus aureus (Sa_Ef-Tu) and Mycoplasma pneumoniae (Mpn_Ef-Tu), and the porcine pathogen Mycoplasma hyopneumoniae (Mhp_Ef-Tu). Ef-Tu is also a target of multiple processing events on the cell surface and these were characterised using an N-terminomics pipeline. Recombinant Mpn_Ef-Tu bound strongly to a diverse range of host molecules, and when bound to plasminogen, was able to convert plasminogen to plasmin in the presence of plasminogen activators. Fragments of Ef-Tu retain binding capabilities to host proteins. Bioinformatics and structural modelling studies indicate that the accumulation of positively charged amino acids in short linear motifs (SLiMs), and protein processing promote multifunctional behaviour. Codon bias engendered by an A + T rich genome may influence how positively-charged residues accumulate in SLiMs.

Functional Constraints on Replacing an Essential Gene with Its Ancient and Modern Homologs

Genes encoding proteins that carry out essential informational tasks in the cell, in particular where multiple interaction partners are involved, are less likely to be transferable to a foreign organism. Here, we investigated the constraints on transfer of a gene encoding a highly conserved informational protein, translation elongation factor Tu (EF-Tu), by systematically replacing the endogenous tufA gene in the Escherichia coli genome with its extant and ancestral homologs. The extant homologs represented tuf variants from both near and distant homologous organisms. The ancestral homologs represented phylogenetically resurrected tuf sequences dating from 0.7 to 3.6 billion years ago (bya). Our results demonstrate that all of the foreign tuf genes are transferable to the E. coli genome, provided that an additional copy of the EF-Tu gene, tufB, remains present in the E. coli genome. However, when the tufB gene was removed, only the variants obtained from the gammaproteobacterial family (extant and ancestral) supported growth which demonstrates the limited functional interchangeability of E. coli tuf with its homologs. Relative bacterial fitness correlated with the evolutionary distance of the extant tuf homologs inserted into the E. coli genome. This reduced fitness was associated with reduced levels of EF-Tu and reduced rates of protein synthesis. Increasing the expression of tuf partially ameliorated these fitness costs. In summary, our analysis suggests that the functional conservation of protein activity, the amount of protein expressed, and its network connectivity act to constrain the successful transfer of this essential gene into foreign bacteria.IMPORTANCE Horizontal gene transfer (HGT) is a fundamental driving force in bacterial evolution. However, whether essential genes can be acquired by HGT and whether they can be acquired from distant organisms are very poorly understood. By systematically replacing tuf with ancestral homologs and homologs from distantly related organisms, we investigated the constraints on HGT of a highly conserved gene with multiple interaction partners. The ancestral homologs represented phylogenetically resurrected tuf sequences dating from 0.7 to 3.6 bya. Only variants obtained from the gammaproteobacterial family (extant and ancestral) supported growth, demonstrating the limited functional interchangeability of E. coli tuf with its homologs. Our analysis suggests that the functional conservation of protein activity, the amount of protein expressed, and its network connectivity act to constrain the successful transfer of this essential gene into foreign bacteria.

Robustness of predictions of extremely thermally stable proteins in ancient organisms

A number of studies have addressed the environmental temperatures experienced by ancient life. Computational studies using a nonhomogeneous evolution model have estimated ancestral G + C contents of ribosomal RNAs and the amino acid compositions of ancestral proteins, generating hypotheses regarding the mesophilic last universal common ancestor. In contrast, our previous study computationally reconstructed ancestral amino acid sequences of nucleoside diphosphate kinases using a homogeneous model and then empirically resurrected the ancestral proteins. The thermal stabilities of these ancestral proteins were equivalent to or greater than those of extant homologous thermophilic proteins, supporting the thermophilic universal ancestor theory. In this study, we reinferred ancestral sequences using a dataset from which hyperthermophilic sequences were excluded. We also reinferred ancestral sequences using a nonhomogeneous evolution model. The newly reconstructed ancestral proteins are still thermally stable, further supporting the hypothesis that the ancient organisms contained thermally stable proteins and therefore that they were thermophilic.

Probabilistic models of eukaryotic evolution: time for integration

In spite of substantial work and recent progress, a global and fully resolved picture of the macroevolutionary history of eukaryotes is still under construction. This concerns not only the phylogenetic relations among major groups, but also the general characteristics of the underlying macroevolutionary processes, including the patterns of gene family evolution associated with endosymbioses, as well as their impact on the sequence evolutionary process. All these questions raise formidable methodological challenges, calling for a more powerful statistical paradigm. In this direction, model-based probabilistic approaches have played an increasingly important role. In particular, improved models of sequence evolution accounting for heterogeneities across sites and across lineages have led to significant, although insufficient, improvement in phylogenetic accuracy. More recently, one main trend has been to move away from simple parametric models and stepwise approaches, towards integrative models explicitly considering the intricate interplay between multiple levels of macroevolutionary processes. Such integrative models are in their infancy, and their application to the phylogeny of eukaryotes still requires substantial improvement of the underlying models, as well as additional computational developments.

The two-domain tree of life is linked to a new root for the Archaea

One of the most fundamental questions in evolutionary biology is the origin of the lineage leading to eukaryotes. Recent phylogenomic analyses have indicated an emergence of eukaryotes from within the radiation of modern Archaea and specifically from a group comprising Thaumarchaeota/"Aigarchaeota" (candidate phylum)/Crenarchaeota/Korarchaeota (TACK). Despite their major implications, these studies were all based on the reconstruction of universal trees and left the exact placement of eukaryotes with respect to the TACK lineage unclear. Here we have applied an original two-step approach that involves the separate analysis of markers shared between Archaea and eukaryotes and between Archaea and Bacteria. This strategy allowed us to use a larger number of markers and greater taxonomic coverage, obtain high-quality alignments, and alleviate tree reconstruction artifacts potentially introduced when analyzing the three domains simultaneously. Our results robustly indicate a sister relationship of eukaryotes with the TACK superphylum that is strongly associated with a distinct root of the Archaea that lies within the Euryarchaeota, challenging the traditional topology of the archaeal tree. Therefore, if we are to embrace an archaeal origin for eukaryotes, our view of the evolution of the third domain of life will have to be profoundly reconsidered, as will many areas of investigation aimed at inferring ancestral characteristics of early life and Earth.

Genomewide comparison and novel ncRNAs of Aquificales

Background

The Aquificales are a diverse group of thermophilic bacteria that thrive in terrestrial and marine hydrothermal environments. They can be divided into the families Aquificaceae, Desulfurobacteriaceae and Hydrogenothermaceae. Although eleven fully sequenced and assembled genomes are available, only little is known about this taxonomic order in terms of RNA metabolism.

Results

In this work, we compare the available genomes, extend their protein annotation, identify regulatory sequences, annotate non-coding RNAs (ncRNAs) of known function, predict novel ncRNA candidates, show idiosyncrasies of the genetic decoding machinery, present two different types of transfer-messenger RNAs and variations of the CRISPR systems. Furthermore, we performed a phylogenetic analysis of the Aquificales based on entire genome sequences, and extended this by a classification among all bacteria using 16S rRNA sequences and a set of orthologous proteins.

Combining several in silico features (e.g. conserved and stable secondary structures, GC-content, comparison based on multiple genome alignments) with an in vivo dRNA-seq transcriptome analysis of Aquifex aeolicus, we predict roughly 100 novel ncRNA candidates in this bacterium.

Conclusions

We have here re-analyzed the Aquificales, a group of bacteria thriving in extreme environments, sharing the feature of a small, compact genome with a reduced number of protein and ncRNA genes. We present several classical ncRNAs and riboswitch candidates. By combining in silico analysis with dRNA-seq data of A. aeolicus we predict nearly 100 novel ncRNA candidates.

Archaea: the first domain of diversified life

The study of the origin of diversified life has been plagued by technical and conceptual difficulties, controversy, and apriorism. It is now popularly accepted that the universal tree of life is rooted in the akaryotes and that Archaea and Eukarya are sister groups to each other. However, evolutionary studies have overwhelmingly focused on nucleic acid and protein sequences, which partially fulfill only two of the three main steps of phylogenetic analysis, formulation of realistic evolutionary models, and optimization of tree reconstruction. In the absence of character polarization, that is, the ability to identify ancestral and derived character states, any statement about the rooting of the tree of life should be considered suspect. Here we show that macromolecular structure and a new phylogenetic framework of analysis that focuses on the parts of biological systems instead of the whole provide both deep and reliable phylogenetic signal and enable us to put forth hypotheses of origin. We review over a decade of phylogenomic studies, which mine information in a genomic census of millions of encoded proteins and RNAs. We show how the use of process models of molecular accumulation that comply with Weston's generality criterion supports a consistent phylogenomic scenario in which the origin of diversified life can be traced back to the early history of Archaea.

Close encounters of the third domain: the emerging genomic view of archaeal diversity and evolution

The Archaea represent the so-called Third Domain of life, which has evolved in parallel with the Bacteria and which is implicated to have played a pivotal role in the emergence of the eukaryotic domain of life. Recent progress in genomic sequencing technologies and cultivation-independent methods has started to unearth a plethora of data of novel, uncultivated archaeal lineages. Here, we review how the availability of such genomic data has revealed several important insights into the diversity, ecological relevance, metabolic capacity, and the origin and evolution of the archaeal domain of life.

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.

Rooted phylogeny of the three superkingdoms

The traditional bacterial rooting of the three superkingdoms in sequence-based gene trees is inconsistent with new phylogenetic reconstructions based on genome content of compact protein domains. We find that protein domains at the level of the SCOP superfamily (SF) from sequenced genomes implement with maximum parsimony fully resolved rooted trees. Such genome content trees identify archaea and bacteria (akaryotes) as sister clades that diverge from an akaryote common ancestor, LACA. Several eukaryote sister clades diverge from a eukaryote common ancestor, LECA. LACA and LECA descend in parallel from the most recent universal common ancestor (MRUCA), which is not a bacterium. Rather, MRUCA presents 75% of the unique SFs encoded by extant genomes of the three superkingdoms, each encoding a proteome that partially overlaps all others. This alone implies that the common ancestor to the superkingdoms was very complex. Such ancestral complexity is confirmed by phylogenetic reconstructions. In addition, the divergence of proteomes from the complex ancestor in each superkingdom is both reductive in numbers of unique SFs as well as cumulative in the abundance of surviving SFs. These data suggest that the common ancestor was not the first cell lineage and that modern global phylogeny is the crown of a "recently" re-rooted tree. We suggest that a bottlenecked survivor of an environmental collapse, which preceded the flourishing of the modern crown, seeded the current phylogenetic tree.

Being Aquifex aeolicus: Untangling a hyperthermophile's checkered past

Lateral gene transfer (LGT) is an important factor contributing to the evolution of prokaryotic genomes. The Aquificae are a hyperthermophilic bacterial group whose genes show affiliations to many other lineages, including the hyperthermophilic Thermotogae, the Proteobacteria, and the Archaea. Previous phylogenomic analyses focused on Aquifex aeolicus identified Thermotogae and Aquificae either as successive early branches or sisters in a rooted bacterial phylogeny, but many phylogenies and cellular traits have suggested a stronger affiliation with the Epsilonproteobacteria. Different scenarios for the evolution of the Aquificae yield different phylogenetic predictions. Here, we outline these scenarios and consider the fit of the available data, including three sequenced Aquificae genomes, to different sets of predictions. Evidence from phylogenetic profiles and trees suggests that the Epsilonproteobacteria have the strongest affinities with the three Aquificae analyzed. However, this pattern is shown by only a minority of encoded proteins, and the Archaea, many lineages of thermophilic bacteria, and members of genus Clostridium and class Deltaproteobacteria also show strong connections to the Aquificae. The phylogenetic affiliations of different functional subsystems showed strong biases: Most but not all genes implicated in the core translational apparatus tended to group Aquificae with Thermotogae, whereas a wide range of metabolic and cellular processes strongly supported the link between Aquificae and Epsilonproteobacteria. Depending on which sets of genes are privileged, either Thermotogae or Epsilonproteobacteria is the most plausible adjacent lineage to the Aquificae. Both scenarios require massive sharing of genes to explain the history of this enigmatic group, whose history is further complicated by specific affinities of different members of Aquificae to different partner lineages.

A congruent phylogenomic signal places eukaryotes within the Archaea

Determining the relationships among the major groups of cellular life is important for understanding the evolution of biological diversity, but is difficult given the enormous time spans involved. In the textbook ‘three domains’ tree based on informational genes, eukaryotes and Archaea share a common ancestor to the exclusion of Bacteria. However, some phylogenetic analyses of the same data have placed eukaryotes within the Archaea, as the nearest relatives of different archaeal lineages. We compared the support for these competing hypotheses using sophisticated phylogenetic methods and an improved sampling of archaeal biodiversity. We also employed both new and existing tests of phylogenetic congruence to explore the level of uncertainty and conflict in the data. Our analyses suggested that much of the observed incongruence is weakly supported or associated with poorly fitting evolutionary models. All of our phylogenetic analyses, whether on small subunit and large subunit ribosomal RNA or concatenated protein-coding genes, recovered a monophyletic group containing eukaryotes and the TACK archaeal superphylum comprising the Thaumarchaeota, Aigarchaeota, Crenarchaeota and Korarchaeota. Hence, while our results provide no support for the iconic three-domain tree of life, they are consistent with an extended eocyte hypothesis whereby vital components of the eukaryotic nuclear lineage originated from within the archaeal radiation.

Breaking through a phylogenetic impasse: a pair of associated archaea might have played host in the endosymbiotic origin of eukaryotes

For over a century, the origin of eukaryotes has been a topic of intense debate among scientists. Although it has become widely accepted that organelles such as the mitochondria and chloroplasts arose via endosymbiosis, the origin of the eukaryotic nucleus remains enigmatic. Numerous models for the origin of the nucleus have been proposed over the years, many of which use endosymbiosis to explain its existence. Proposals of microbes whose ancestors may have served as either a host or a guest in various endosymbiotic scenarios abound, none of which have been able to sufficiently incorporate the cell biological as well as phylogenetic data which links these organisms to the nucleus. While it is generally agreed that eukaryotic nuclei share more features in common with archaea rather than with bacteria, different studies have identified either one or the other of the two major groups of archaea as potential ancestors, leading to somewhat of a stalemate. This paper seeks to resolve this impasse by presenting evidence that not just one, but a pair of archaea might have served as host to the bacterial ancestor of the mitochondria. This pair may have consisted of ancestors of both Ignicoccus hospitalis as well as its ectosymbiont/ectoparasite ‘Nanoarchaeum equitans’.

Rooting the eukaryotic tree with mitochondrial and bacterial proteins

By exploiting the large body of genome data and the considerable progress in phylogenetic methodology, recent phylogenomic studies have provided new insights into the relationships among major eukaryotic groups. However, confident placement of the eukaryotic root remains a major challenge. This is due to the large evolutionary distance separating eukaryotes from their closest relatives, the Archaea, implying a weak phylogenetic signal and strong long-branch attraction artifacts. Here, we apply a new approach to the rooting of the eukaryotic tree by using a subset of genomic information with more recent evolutionary origin-mitochondrial sequences, whose closest relatives are α-Proteobacteria. For this, we identified and assembled a data set of 42 mitochondrial proteins (mainly encoded by the nuclear genome) and performed Bayesian and maximum likelihood analyses. Taxon sampling includes the recently sequenced Thecamonas trahens, a member of the phylogenetically elusive Apusozoa. This data set confirms the relationships of several eukaryotic supergroups seen before and places the eukaryotic root between the monophyletic "unikonts" and "bikonts." We further show that T. trahens branches sister to Opisthokonta with significant statistical support and question the bikont/excavate affiliation of Malawimonas species. The mitochondrial data set developed here (to be expanded in the future) constitutes a unique alternative means in resolving deep eukaryotic relationships.

The archaeal 'TACK' superphylum and the origin of eukaryotes

Although most hypotheses to explain the emergence of the eukaryotic lineage are conflicting, some consensus exists concerning the requirement of a genomic fusion between archaeal and bacterial components. Recent phylogenomic studies have provided support for eocyte-like scenarios in which the alleged 'archaeal parent' of the eukaryotic cell emerged from the Crenarchaeota/Thaumarchaeota. Here, we provide evidence for a scenario in which this archaeal parent emerged from within the 'TACK' superphylum that comprises the Thaumarchaeota, Crenarchaeota and Korarchaeota, as well as the recently proposed phylum 'Aigarchaeota'. In support of this view, functional and comparative genomics studies have unearthed an increasing number of features that are uniquely shared by the TACK superphylum and eukaryotes, including proteins involved in cytokinesis, membrane remodeling, cell shape determination and protein recycling.

How stands the Tree of Life a century and a half after The Origin?

We examine the Tree of Life (TOL) as an evolutionary hypothesis and a heuristic. The original TOL hypothesis has failed but a new "statistical TOL hypothesis" is promising. The TOL heuristic usefully organizes data without positing fundamental evolutionary truth.

A late methanogen origin for molybdenum-dependent nitrogenase

Mounting evidence indicates the presence of a near complete biological nitrogen cycle in redox-stratified oceans during the late Archean to early Proterozoic (c. 2.5-2.0 Ga). It has been suggested that the iron (Fe)- or vanadium (V)-dependent nitrogenase rather than molybdenum (Mo)-dependent form was responsible for dinitrogen fixation during this time because oceans were depleted in Mo and rich in Fe. We evaluated this hypothesis by examining the phylogenetic relationships of proteins that are required for the biosynthesis of the active site cofactor of Mo-nitrogenase in relation to structural proteins required for Fe-, V- and Mo-nitrogenase. The results are highly suggestive that among extant nitrogen-fixing organisms for which genomic information exists, Mo-nitrogenase is unlikely to have been associated with the Last Universal Common Ancestor. Rather, the origin of Mo-nitrogenase can be traced to an ancestor of the anaerobic and hydrogenotrophic methanogens with acquisition in the bacterial domain via lateral gene transfer involving an anaerobic member of the Firmicutes. A comparison of substitution rates estimated for proteins required for the biosynthesis of the nitrogenase active site cofactor and for a set of paralogous proteins required for the biosynthesis of bacteriochlorophyll suggests that Nif emerged from a nitrogenase-like ancestor approximately 1.5-2.2 Ga. An origin and ensuing proliferation of Mo-nitrogenase under anoxic conditions would likely have occurred in an environment where anaerobic methanogens and Firmicutes coexisted and where Mo was at least episodically available, such as in a redox-stratified Proterozoic ocean basin.

Archaeal phylogenomics provides evidence in support of a methanogenic origin of the Archaea and a thaumarchaeal origin for the eukaryotes

We have developed a machine-learning approach to identify 3537 discrete orthologue protein sequence groups distributed across all available archaeal genomes. We show that treating these orthologue groups as binary detection/non-detection data is sufficient to capture the majority of archaeal phylogeny. We subsequently use the sequence data from these groups to infer a method and substitution-model-independent phylogeny. By holding this phylogeny constrained and interrogating the intersection of this large dataset with both the Eukarya and the Bacteria using Bayesian and maximum-likelihood approaches, we propose and provide evidence for a methanogenic origin of the Archaea. By the same criteria, we also provide evidence in support of an origin for Eukarya either within or as sisters to the Thaumarchaea.

An actin-based cytoskeleton in archaea

In eukaryotic and bacterial cells, spatial organization is dependent upon cytoskeletal filaments. Actin is a main eukaryotic cytoskeletal element, involved in key processes such as cell shape determination, mechanical force generation and cytokinesis. We describe an archaeal cytoskeleton which forms helical structures within Pyrobaculum calidifontis cells, as shown by in situ immunostaining. The core components include an archaeal actin homologue, Crenactin, closely related to the eukaryotic counterpart. The crenactin gene belongs to a conserved gene cluster denoted Arcade (actin-related cytoskeleton in Archaea involved in shape determination). The phylogenetic distribution of arcade genes is restricted to the crenarchaeal Thermoproteales lineage, and to Korarchaeota, and correlates with rod-shaped and filamentous cell morphologies. Whereas Arcadin-1, -3 and -4 form helical structures, suggesting cytoskeleton-associated functions, Arcadin-2 was found to be localized between segregated nucleoids in a cell subpopulation, in agreement with possible involvement in cytokinesis. The results support a crenarchaeal origin of the eukaryotic actin cytoskeleton and, as such, have implications for theories concerning the origin of the eukaryotic cell.

Eukaryotic evolution: the importance of being archaebacterial

Approximately half of all eukaryotic genes show signs of prokaryotic origin. Genes derived from eubacteria are more abundant than those from archaebacteria, but the latter are functionally more important. This supports archaebacteria as founding ancestors of the eukaryotic nucleus.

Evolution of the Kdo2-lipid A biosynthesis in bacteria

Background

Lipid A is the highly immunoreactive endotoxic center of lipopolysaccharide (LPS). It anchors the LPS into the outer membrane of most Gram-negative bacteria. Lipid A can be recognized by animal cells, triggers defense-related responses, and causes Gram-negative sepsis. The biosynthesis of Kdo₂-lipid A, the LPS substructure, involves with nine enzymatic steps.

Results

In order to elucidate the evolutionary pathway of Kdo₂-lipid A biosynthesis, we examined the distribution of genes encoding the nine enzymes across bacteria. We found that not all Gram-negative bacteria have all nine enzymes. Some Gram-negative bacteria have no genes encoding these enzymes and others have genes only for the first four enzymes (LpxA, LpxC, LpxD, and LpxB). Among the nine enzymes, five appeared to have arisen from three independent gene duplication events. Two of such events happened within the Proteobacteria lineage, followed by functional specialization of the duplicated genes and pathway optimization in these bacteria.

Conclusions

The nine-enzyme pathway, which was established based on the studies mainly in Escherichia coli K12, appears to be the most derived and optimized form. It is found only in E. coli and related Proteobacteria. Simpler and probably less efficient pathways are found in other bacterial groups, with Kdo₂-lipid A variants as the likely end products. The Kdo₂-lipid A biosynthetic pathway exemplifies extremely plastic evolution of bacterial genomes, especially those of Proteobacteria, and how these mainly pathogenic bacteria have adapted to their environment.

Evolution of elongation factor G and the origins of mitochondrial and chloroplast forms

Protein synthesis elongation factor G (EF-G) is an essential protein with central roles in both the elongation and ribosome recycling phases of protein synthesis. Although EF-G evolution is predicted to be conservative, recent reports suggest otherwise. We have characterized EF-G in terms of its molecular phylogeny, genomic context, and patterns of amino acid substitution. We find that most bacteria carry a single "canonical" EF-G, which is phylogenetically conservative and encoded in an str operon. However, we also find a number of EF-G paralogs. These include a pair of EF-Gs that are mostly found together and in an eclectic subset of bacteria, specifically δ-proteobacteria, spirochaetes, and planctomycetes (the "spd" bacteria). These spdEFGs have also given rise to the mitochondrial factors mtEFG1 and mtEFG2, which probably arrived in eukaryotes before the eukaryotic last common ancestor. Meanwhile, chloroplasts apparently use an α-proteobacterial-derived EF-G rather than the expected cyanobacterial form. The long-term comaintenance of the spd/mtEFGs may be related to their subfunctionalization for translocation and ribosome recycling. Consistent with this, patterns of sequence conservation and site-specific evolutionary rate shifts suggest that the faster evolving spd/mtEFG2 has lost translocation function, but surprisingly, the protein also shows little conservation of sites related to recycling activity. On the other hand, spd/mtEFG1, although more slowly evolving, shows signs of substantial remodeling. This is particularly extensive in the GTPase domain, including a highly conserved three amino acid insertion in switch I. We suggest that subfunctionalization of the spd/mtEFGs is not a simple case of specialization for subsets of original activities. Rather, the duplication allows the release of one paralog from the selective constraints imposed by dual functionality, thus allowing it to become more highly specialized. Thus, the potential for fine tuning afforded by subfunctionalization may explain the maintenance of EF-G paralogs.

Evolution of the Archaea: emerging views on origins and phylogeny

Of the three domains of life, the Archaea are the most recently discovered and, from the perspective of systematics, perhaps the least understood. More than three decades after their discovery, there is still no overwhelming consensus as to their phylogenetic status, with diverse evidence supporting in varying degrees their monophyly, paraphyly, or even polyphyly. As a further complication, their evolutionary history is inextricably linked to the origin of Eukarya, one of the most challenging problems in evolutionary biology. This exclusive relationship between the eukaryal nucleocytoplasm and the Archaea is further supported by a new methodology for rooting the ribosomal Tree of Life based on amino acid composition. Novel approaches such as utilizing horizontal gene transfers as synchronizing events and branch length analysis of deep paralogs will help to clarify temporal relationships between these lineages, and may prove useful in evaluating the numerous conflicting hypotheses related to the evolution of the Archaea and Eukarya.