The universal tree of life: an update

Structures of Saccharolobus solfataricus initiation complexes with leaderless mRNAs highlight archaeal features and eukaryotic proximity

The archaeal ribosome is of the eukaryotic type. TACK and Asgard superphyla, the closest relatives of eukaryotes, have ribosomes containing eukaryotic ribosomal proteins not found in other archaea, eS25, eS26 and eS30. Here, we investigate the case of Saccharolobus solfataricus, a TACK crenarchaeon, using mainly leaderless mRNAs. We characterize the small ribosomal subunit of S. solfataricus bound to SD-leadered or leaderless mRNAs. Cryo-EM structures show eS25, eS26 and eS30 bound to the small subunit. We identify two ribosomal proteins, aS33 and aS34, and an additional domain of eS6. Leaderless mRNAs are bound to the small subunit with contribution of their 5'-triphosphate group. Archaeal eS26 binds to the mRNA exit channel wrapped around the 3' end of rRNA, as in eukaryotes. Its position is not compatible with an SD:antiSD duplex. Our results suggest a positive role of eS26 in leaderless mRNAs translation and possible evolutionary routes from archaeal to eukaryotic translation.

The Expanding Universe of Extensions and Tails: Ribosomal Proteins and Histones in RNA and DNA Complex Signaling and Dynamics

This short review bridges two biological fields: ribosomes and nucleosomes-two nucleoprotein assemblies that, along with many viruses, share proteins featuring long filamentous segments at their N- or C-termini. A central hypothesis is that these extensions and tails perform analogous functions in both systems. The evolution of these structures appears closely tied to the emergence of regulatory networks and signaling pathways, facilitating increasingly complex roles for ribosomes and nucleosome alike. This review begins by summarizing the structures and functions of ribosomes and nucleosomes, followed by a detailed comparison highlighting their similarities and differences, particularly in light of recent findings on the roles of ribosomal proteins in signaling and ribosome dynamics. The analysis seeks to uncover whether these systems operate based on shared principles and mechanisms. The nucleosome-ribosome analogy may offer valuable insights into unresolved questions in both fields. For instance, new structural insights from ribosomes might shed light on potential motifs formed by histone tails. From an evolutionary perspective, this study revisits the origins of signaling and regulation in ancient nucleoprotein assemblies, suggesting that tails and extensions may represent remnants of the earliest network systems governing signaling and dynamic control.

Stem Life: A Framework for Understanding the Prebiotic-Biotic Transition

Abiogenesis is frequently envisioned as a linear, ladder-like progression of increasingly complex chemical systems, eventually leading to the ancestors of extant cellular life. This "pre-cladistics" view is in stark contrast to the well-accepted principles of organismal evolutionary biology, as informed by paleontology and phylogenetics. Applying this perspective to origins, I explore the paradigm of "Stem Life," which embeds abiogenesis within a broader continuity of diversification and extinction of both hereditary lineages and chemical systems. In this new paradigm, extant life's ancestral lineage emerged alongside and was dependent upon many other complex prebiotic chemical systems, as part of a diverse and fecund prebiosphere. Drawing from several natural history analogies, I show how this shift in perspective enriches our understanding of Origins and directly informs debates on defining Life, the emergence of the Last Universal Common Ancestor (LUCA), and the implications of prebiotic chemical experiments.

The Last Universal Common Ancestor of Ribosome-Encoding Organisms: Portrait of LUCA

The existence of LUCA in the distant past is the logical consequence of the binary mechanism of cell division. The biosphere in which LUCA and contemporaries were living was the product of a long cellular evolution from the origin of life to the second age of the RNA world. A parsimonious scenario suggests that the molecular fabric of LUCA was much simpler than those of modern organisms, explaining why the evolutionary tempo was faster at the time of LUCA than it was during the diversification of the three domains. Although LUCA was possibly equipped with a RNA genome and most likely lacked an ATP synthase, it was already able to perform basic metabolic functions and to produce efficient proteins. However, the proteome of LUCA and its inferred metabolism remains to be correctly explored by in-depth phylogenomic analyses and updated datasets. LUCA was probably a mesophile or a moderate thermophile since phylogenetic analyses indicate that it lacked reverse gyrase, an enzyme systematically present in all hyperthermophiles. The debate about the position of Eukarya in the tree of life, either sister group to Archaea or descendants of Archaea, has important implications to draw the portrait of LUCA. In the second alternative, one can a priori exclude the presence of specific eukaryotic features in LUCA. In contrast, if Archaea and Eukarya are sister group, some eukaryotic features, such as the spliceosome, might have been present in LUCA and later lost in Archaea and Bacteria. The nature of the LUCA virome is another matter of debate. I suggest here that DNA viruses only originated during the diversification of the three domains from an RNA-based LUCA to explain the odd distribution pattern of DNA viruses in the tree of life.

Protein Fold Usages in Ribosomes: Another Glance to the Past

The analysis of protein fold usage, similar to codon usage, offers profound insights into the evolution of biological systems and the origins of modern proteomes. While previous studies have examined fold distribution in modern genomes, our study focuses on the comparative distribution and usage of protein folds in ribosomes across bacteria, archaea, and eukaryotes. We identify the prevalence of certain 'super-ribosome folds,' such as the OB fold in bacteria and the SH3 domain in archaea and eukaryotes. The observed protein fold distribution in the ribosomes announces the future power-law distribution where only a few folds are highly prevalent, and most are rare. Additionally, we highlight the presence of three copies of proto-Rossmann folds in ribosomes across all kingdoms, showing its ancient and fundamental role in ribosomal structure and function. Our study also explores early mechanisms of molecular convergence, where different protein folds bind equivalent ribosomal RNA structures in ribosomes across different kingdoms. This comparative analysis enhances our understanding of ribosomal evolution, particularly the distinct evolutionary paths of the large and small subunits, and underscores the complex interplay between RNA and protein components in the transition from the RNA world to modern cellular life. Transcending the concept of folds also makes it possible to group a large number of ribosomal proteins into five categories of urfolds or metafolds, which could attest to their ancestral character and common origins. This work also demonstrates that the gradual acquisition of extensions by simple but ordered folds constitutes an inexorable evolutionary mechanism. This observation supports the idea that simple but structured ribosomal proteins preceded the development of their disordered extensions.

Globin phylogeny, evolution and function, the newest update

Our globin census update allows us to refine our vision of globin origin, evolution, and structure to function relationship in the context of the currently accepted tree of life. The modern globin domain originates as a single domain, three-over-three α-helical folded structure before the diversification of the kingdoms of life (Bacteria, Archaea, Eukarya). Together with the diversification of prokaryotes, three monophyletic globin families (M, S, and T) emerged, most likely in Proteobacteria and Actinobacteria, displaying specific sequence and structural features, and spread by vertical and horizontal gene transfer, most probably already present in the last universal common ancestor (LUCA). Non-globin domains were added, and eventually lost again, creating multi-domain structures in key branches of M- (FHb and Adgb) and the vast majority of S globins, which with their coevolved multi-domain architectures, have predominantly "sensor" functions. Single domain T-family globins diverged into four major groups and most likely display functions related to reactive nitrogen and oxygen species (RNOS) chemistry, as well as oxygen storage/transport which drives the evolution of its major branches with their characteristic key distal residues (B10, E11, E7, and G8). M-family evolution also lead to distinctive major types (FHb and Fgb, Ngb, Adgb, GbX vertebrate Gbs), and shows the shift from high oxygen affinity controlled by TyrB10-Gln/AsnE11 likely related to RNOS chemistry in microorganisms, to a moderate oxygen affinity storage/transport function controlled by hydrophobic B10/E11-HisE7 in multicellular animals.

Chapter 4: A Geological and Chemical Context for the Origins of Life on Early Earth

Within the first billion years of Earth's history, the planet transformed from a hot, barren, and inhospitable landscape to an environment conducive to the emergence and persistence of life. This chapter will review the state of knowledge concerning early Earth's (Hadean/Eoarchean) geochemical environment, including the origin and composition of the planet's moon, crust, oceans, atmosphere, and organic content. It will also discuss abiotic geochemical cycling of the CHONPS elements and how these species could have been converted to biologically relevant building blocks, polymers, and chemical networks. Proposed environments for abiogenesis events are also described and evaluated. An understanding of the geochemical processes under which life may have emerged can better inform our assessment of the habitability of other worlds, the potential complexity that abiotic chemistry can achieve (which has implications for putative biosignatures), and the possibility for biochemistries that are vastly different from those on Earth.

Chapter 2: What Is Life?

The question "What is life?" has existed since the beginning of recorded history. However, the scientific and philosophical contexts of this question have changed and been refined as advancements in technology have revealed both fine details and broad connections in the network of life on Earth. Understanding the framework of the question "What is life?" is central to formulating other questions such as "Where else could life be?" and "How do we search for life elsewhere?" While many of these questions are addressed throughout the Astrobiology Primer 3.0, this chapter gives historical context for defining life, highlights conceptual characteristics shared by all life on Earth as well as key features used to describe it, discusses why it matters for astrobiology, and explores both challenges and opportunities for finding an informative operational definition.

Sua5 catalyzing universal t6A tRNA modification is responsible for multifaceted functions of the KEOPS complex in Cryptococcus neoformans

The N6-threonylcarbamoyl adenosine (t6A) tRNA modification is critical for ensuring translation fidelity across three domains of life. Our prior work highlighted the KEOPS complex, organized in a Pcc1-Kae1-Bud32-Cgi121 linear arrangement, not only serves an evolutionarily conserved role in t6A tRNA modification but also exerts diverse functional impacts on pathobiological characteristics in Cryptococcus neoformans, a leading cause of fungal meningitis worldwide. However, the extent to which the pleiotropic functions of the KEOPS complex are specifically tied to tRNA modification remains uncertain. To address this, we undertook a functional characterization of Sua5, responsible for generating the precursor threonylcarbamoyl-adenylate (TC-AMP) for t6A tRNA modification, using a reverse genetics approach. Comparative phenotypic analyses with KEOPS mutants revealed that Sua5 plays a vital role in multiple cellular processes, such as t6A tRNA modification, growth, sexual development, stress response, and virulence factor production, thus reflecting the multifaceted functions of the KEOPS complex. In support of this, sua5Δ bud32Δ double mutants showed phenotypes comparable to those of the corresponding single mutants. Intriguingly, a SUA5 allele lacking a mitochondria targeting sequence (SUA5MTSΔ) was sufficient to restore the wild-type phenotypes in the sua5Δ mutant, suggesting that Sua5's primary functional locus may be cytosolic, akin to the KEOPS complex. Further supporting this, the deletion of Qri7, a mitochondrial paralog of Kae1, had no discernible phenotypic impact on C. neoformans. We concluded that cytosolic t6A tRNA modifications, orchestrated by Sua5 and the KEOPS complex, are central to the regulation of diverse pathobiological functions in C. neoformans.IMPORTANCEUnderstanding cellular functions at the molecular level is crucial for advancing disease treatments. Our research reveals a critical connection between the KEOPS complex and Sua5 in Cryptococcus neoformans, a significant cause of fungal meningitis. While the KEOPS complex is known for its versatile roles in cellular processes, Sua5 is specialized in t6A tRNA modification. Our key finding is that the diverse roles of the KEOPS complex, ranging from cell growth and stress response to virulence, are fundamentally linked to its function in t6A tRNA modification. This conclusion is supported by the remarkable similarities between the impacts of Sua5 and KEOPS on these processes, despite their roles in different steps of the t6A modification pathway. This newfound understanding deepens our insight into fungal biology and opens new avenues for developing potential therapies against dangerous fungal diseases.

How Did Thylakoids Emerge in Cyanobacteria, and How Were the Primary Chloroplast and Chromatophore Acquired?

The emergence of thylakoid membranes in cyanobacteria is a key event in the evolution of all oxygenic photosynthetic cells, from prokaryotes to eukaryotes. Recent analyses show that they could originate from a unique lipid phase transition rather than from a supposed vesicular budding mechanism. Emergence of thylakoids coincided with the great oxygenation event, more than two billion years ago. The acquisition of semi-autonomous organelles, such as the mitochondrion, the chloroplast, and, more recently, the chromatophore, is a critical step in the evolution of eukaryotes. They resulted from primary endosymbiotic events that seem to share general features, i.e., an acquisition of a bacterium/cyanobacteria likely via a phagocytic membrane, a genome reduction coinciding with an escape of genes from the organelle to the nucleus, and, finally, the appearance of an active system translocating nuclear-encoded proteins back to the organelles. An intense mobilization of foreign genes of bacterial origin, via horizontal gene transfers, plays a critical role. Some third partners, like Chlamydia, might have facilitated the transition from cyanobacteria to the early chloroplast. This chapter further details our current understanding of primary endosymbiosis, focusing on primary chloroplasts, thought to have appeared over a billion years ago, and the chromatophore, which appeared around a hundred years ago.

Investigation of tRNA-based relatedness within the Planctomycetes-Verrucomicrobia-Chlamydiae (PVC) superphylum: a comparative analysis

The PVC superphylum is a diverse group of prokaryotes that require stringent growth conditions. RNA is a fascinating molecule to find evolutionary relatedness according to the RNA World Hypothesis. We conducted tRNA gene analysis to find evolutionary relationships in the PVC phyla. The analysis of genomic data (P = 9, V = 4, C = 8) revealed that the number of tRNA genes varied from 28 to 90 in Planctomycetes and Chlamydia, respectively. Verrucomicrobia has whole genomes and the longest scaffold (3 + 1), with tRNA genes ranging from 49 to 53 in whole genomes and 4 in the longest scaffold. Most tRNAs in the E. coli genome clustered with homologs, but approximately 43% clustered with tRNAs encoding different amino acids. Planctomyces, Akkermansia, Isosphaera, and Chlamydia were similar to E. coli tRNAs. In a phylum, tRNAs coding for different amino acids clustered at a range of 8 to 10%. Further analysis of these tRNAs showed sequence similarity with Cyanobacteria, Proteobacteria, Viridiplantae, Ascomycota and Basidiomycota (Eukaryota). This indicates the possibility of horizontal gene transfer or, otherwise, a different origin of tRNA in PVC bacteria. Hence, this work proves its importance for determining evolutionary relatedness and potentially identifying bacteria using tRNA. Thus, the analysis of these tRNAs indicates that primitive RNA may have served as the genetic material of LUCA before being replaced by DNA. A quantitative analysis is required to test these possibilities that relate the evolutionary significance of tRNA to the origin of life.

Main Factors Shaping Amino Acid Usage Across Evolution

The standard genetic code determines that in most species, including viruses, there are 20 amino acids that are coded by 61 codons, while the other three codons are stop triplets. Considering the whole proteome each species features its own amino acid frequencies, given the slow rate of change, closely related species display similar GC content and amino acids usage. In contrast, distantly related species display different amino acid frequencies. Furthermore, within certain multicellular species, as mammals, intragenomic differences in the usage of amino acids are evident. In this communication, we shall summarize some of the most prominent and well-established factors that determine the differences found in the amino acid usage, both across evolution and intragenomically.

The methanogen core and pangenome: conservation and variability across biology's growth temperature extremes

Temperature is a key variable in biological processes. However, a complete understanding of biological temperature adaptation is lacking, in part because of the unique constraints among different evolutionary lineages and physiological groups. Here we compared the genomes of cultivated psychrotolerant and thermotolerant methanogens, which are physiologically related and span growth temperatures from -2.5°C to 122°C. Despite being phylogenetically distributed amongst three phyla in the archaea, the genomic core of cultivated methanogens comprises about one-third of a given genome, while the genome fraction shared by any two organisms decreases with increasing phylogenetic distance between them. Increased methanogenic growth temperature is associated with reduced genome size, and thermotolerant organisms-which are distributed across the archaeal tree-have larger core genome fractions, suggesting that genome size is governed by temperature rather than phylogeny. Thermotolerant methanogens are enriched in metal and other transporters, and psychrotolerant methanogens are enriched in proteins related to structure and motility. Observed amino acid compositional differences between temperature groups include proteome charge, polarity and unfolding entropy. Our results suggest that in the methanogens, shared physiology maintains a large, conserved genomic core even across large phylogenetic distances and biology's temperature extremes.

Examining the diversity of structural motifs in fungal glycome

In this paper, we present the results of a systematic statistical analysis of the fungal glycome in comparison with the prokaryotic and protistal glycomes as described in the scientific literature and presented in the Carbohydrate Structure Database (CSDB). The monomeric and dimeric compositions of glycans, their non-carbohydrate modifications, glycosidic linkages, sizes of structures, branching degree and net charge are assessed. The obtained information can help elucidating carbohydrate molecular markers for various fungal classes which, in its turn, can be demanded for the development of diagnostic tools and carbohydrate-based vaccines against pathogenic fungi. It can also be useful for revealing specific glycosyltransferases active in a particular fungal species.

Calpains in cyanobacteria and the origin of calpains

Calpains are cysteine proteases involved in many cellular processes. They are an ancient and large superfamily of enzymes responsible for the cleavage and irreversible modification of a large variety of substrates. They have been intensively studied in humans and other mammals, but information about calpains in bacteria is scarce. Calpains have not been found among Archaea to date. In this study, we have investigated the presence of calpains in selected cyanobacterial species using in silico analyses. We show that calpains defined by possessing CysPC core domain are present in cyanobacterial genera Anabaena, Aphanizomenon, Calothrix, Chamaesiphon, Fischerella, Microcystis, Scytonema and Trichormus. Based on in silico protein interaction analysis, we have predicted putative interaction partners for identified cyanobacterial calpains. The phylogenetic analysis including cyanobacterial, other bacterial and eukaryotic calpains divided bacterial and eukaryotic calpains into two separate monophyletic clusters. We propose two possible evolutionary scenarios to explain this tree topology: (1) the eukaryotic ancestor or an archaeal ancestor of eukaryotes obtained calpain gene from an unknown bacterial donor, or alternatively (2) calpain gene had been already present in the last common universal ancestor and subsequently lost by the ancestor of Archaea, but retained by the ancestor of Bacteria and by the ancestor of Eukarya. Both scenarios would require multiple independent losses of calpain genes in various bacteria and eukaryotes.

Looking through the Lens of the Ribosome Biogenesis Evolutionary History: Possible Implications for Archaeal Phylogeny and Eukaryogenesis

Our understanding of microbial diversity and its evolutionary relationships has increased substantially over the last decade. Such an understanding has been greatly fueled by culture-independent metagenomics analyses. However, the outcome of some of these studies and their biological and evolutionary implications, such as the origin of the eukaryotic lineage from the recently discovered archaeal Asgard superphylum, is debated. The sequences of the ribosomal constituents are amongst the most used phylogenetic markers. However, the functional consequences underlying the analysed sequence diversity and their putative evolutionary implications are essentially not taken into consideration. Here, we propose to exploit additional functional hallmarks of ribosome biogenesis to help disentangle competing evolutionary hypotheses. Using selected examples, such as the multiple origins of halophily in archaea or the evolutionary relationship between the Asgard archaea and Eukaryotes, we illustrate and discuss how function-aware phylogenetic framework can contribute to refining our understanding of archaeal phylogeny and the origin of eukaryotic cells.

Accurate prediction of metagenome-assembled genome completeness by MAGISTA, a random forest model built on alignment-free intra-bin statistics

BACKGROUND

Although the total number of microbial taxa on Earth is under debate, it is clear that only a small fraction of these has been cultivated and validly named. Evidently, the inability to culture most bacteria outside of very specific conditions severely limits their characterization and further studies. In the last decade, a major part of the solution to this problem has been the use of metagenome sequencing, whereby the DNA of an entire microbial community is sequenced, followed by the in silico reconstruction of genomes of its novel component species. The large discrepancy between the number of sequenced type strain genomes (around 12,000) and total microbial diversity (10⁶-10¹² species) directs these efforts to de novo assembly and binning. Unfortunately, these steps are error-prone and as such, the results have to be intensely scrutinized to avoid publishing incomplete and low-quality genomes.

RESULTS

We developed MAGISTA (metagenome-assembled genome intra-bin statistics assessment), a novel approach to assess metagenome-assembled genome quality that tackles some of the often-neglected drawbacks of current reference gene-based methods. MAGISTA is based on alignment-free distance distributions between contig fragments within metagenomic bins, rather than a set of reference genes. For proper training, a highly complex genomic DNA mock community was needed and constructed by pooling genomic DNA of 227 bacterial strains, specifically selected to obtain a wide variety representing the major phylogenetic lineages of cultivable bacteria.

CONCLUSIONS

MAGISTA achieved a 20% reduction in root-mean-square error in comparison to the marker gene approach when tested on publicly available mock metagenomes. Furthermore, our highly complex genomic DNA mock community is a very valuable tool for benchmarking (new) metagenome analysis methods.

The expanding Asgard archaea and their elusive relationships with Eukarya

The discovery of Asgard archaea and the exploration of their diversity over the last 6 years have deeply impacted the scientific community working on eukaryogenesis, rejuvenating an intense debate on the topology of the universal tree of life (uTol). Here, we discuss how this debate is impacted by two recent publications that expand the number of Asgard lineages and eukaryotic signature proteins (ESPs). We discuss some of the main difficulties that can impair the phylogenetic reconstructions of the uTol and suggest that the debate about its topology is not settled. We notably hypothesize the existence of horizontal gene transfers between ancestral Asgards and proto-eukaryotes that could result in the observed abnormal behaviors of some Asgard ESPs and universal marker proteins. This hypothesis is relevant regardless of the scenario considered regarding eukaryogenesis. It implies that the Asgards were already diversified before the last eukaryotic common ancestor and shared the same biotopes with proto-eukaryotes. We suggest that some Asgards might be still living in symbiosis today with modern Eukarya.

Archaea: A Goldmine for Molecular Biologists and Evolutionists

The rebuttal of the prokaryote-eukaryote dichotomy and the elaboration of the three domains concept by Carl Woese and colleagues has been a breakthrough in biology. With the methodologies available at this time, they have shown that a single molecule, the 16S ribosomal RNA, could reveal the global organization of the living world. Later on, mining archaeal genomes led to major discoveries in archaeal molecular biology, providing a third model for comparative molecular biology. These analyses revealed the strong eukaryal flavor of the basic molecular fabric of Archaea and support rooting the universal tree between Bacteria and Arcarya (the clade grouping Archaea and Eukarya). However, in contradiction with this conclusion, it remains to understand why the archaeal and bacterial mobilomes are so similar and so different from the eukaryal one. These last years, the number of recognized archaea lineages (phyla?) has exploded. The archaeal nomenclature is now in turmoil and debates about the nature of the last universal common ancestor, the last archaeal common ancestor, and the topology of the tree of life are still going on. Interestingly, the expansion of the archaeal eukaryome, especially in the Asgard archaea, has provided new opportunities to study eukaryogenesis. In recent years, the application to Archaea of the new methodologies described in the various chapters of this book have opened exciting avenues to study the molecular biology and the physiology of these fascinating microorganisms.

Towards the Idea of Molecular Brains

How can single cells without nervous systems perform complex behaviours such as habituation, associative learning and decision making, which are considered the hallmark of animals with a brain? Are there molecular systems that underlie cognitive properties equivalent to those of the brain? This review follows the development of the idea of molecular brains from Darwin’s “root brain hypothesis”, through bacterial chemotaxis, to the recent discovery of neuron-like r-protein networks in the ribosome. By combining a structural biology view with a Bayesian brain approach, this review explores the evolutionary labyrinth of information processing systems across scales. Ribosomal protein networks open a window into what were probably the earliest signalling systems to emerge before the radiation of the three kingdoms. While ribosomal networks are characterised by long-lasting interactions between their protein nodes, cell signalling networks are essentially based on transient interactions. As a corollary, while signals propagated in persistent networks may be ephemeral, networks whose interactions are transient constrain signals diffusing into the cytoplasm to be durable in time, such as post-translational modifications of proteins or second messenger synthesis. The duration and nature of the signals, in turn, implies different mechanisms for the integration of multiple signals and decision making. Evolution then reinvented networks with persistent interactions with the development of nervous systems in metazoans. Ribosomal protein networks and simple nervous systems display architectural and functional analogies whose comparison could suggest scale invariance in information processing. At the molecular level, the significant complexification of eukaryotic ribosomal protein networks is associated with a burst in the acquisition of new conserved aromatic amino acids. Knowing that aromatic residues play a critical role in allosteric receptors and channels, this observation suggests a general role of π systems and their interactions with charged amino acids in multiple signal integration and information processing. We think that these findings may provide the molecular basis for designing future computers with organic processors.

Kinetics of the ancestral carbon metabolism pathways in deep-branching bacteria and archaea

The origin of life is believed to be chemoautotrophic, deriving all biomass components from carbon dioxide, and all energy from inorganic redox couples in the environment. The reductive tricarboxylic acid cycle (rTCA) and the Wood-Ljungdahl pathway (WL) have been recognized as the most ancient carbon fixation pathways. The rTCA of the chemolithotrophic Thermosulfidibacter takaii, which was recently demonstrated to take place via an unexpected reverse reaction of citrate synthase, was reproduced using a kinetic network model, and a competition between reductive and oxidative fluxes on rTCA due to an acetyl coenzyme A (ACOA) influx upon acetate uptake was revealed. Avoiding ACOA direct influx into rTCA from WL is, therefore, raised as a kinetically necessary condition to maintain a complete rTCA. This hypothesis was confirmed for deep-branching bacteria and archaea, and explains the kinetic factors governing elementary processes in carbon metabolism evolution from the last universal common ancestor.

Photosynthetic Systems Suggest an Evolutionary Pathway to Diderms

Bacteria are divided primarily into monoderms (with one cell membrane, and usually Gram-positive, due to a thick peptidoglycan layer) and diderms (with two cell membranes, and mostly Gram-negative, due to a thin peptidoglycan layer sandwiched between the two membranes). Photosynthetic species are spread among the taxonomic groups, some having type I reaction centers (RCI in monoderm phylum Firmicutes; and diderm phyla Acidobacteria and Chlorobi), others with type II reaction centers (RCII in monoderm phylum Chloroflexi; and diderm taxa Gemmatimonadetes, and alpha-, beta-, and gamma-Proteobacteria), and some containing both (RCI and RCII, only in diderm phylum Cyanobacteria). In most bacterial phylograms, photosystem types and diderm taxa are polyphyletic. A more parsimonious arrangement, which is supported by photosystem evolution, as well as additional sets of molecular characters, suggests that endosymbiotic events resulted in the formation of the diderms. In the model presented, monoderms readily form a monophyletic group, while diderms are produced by at least two endosymbiotic events, followed by additional evolutionary changes.

Evolution and origin of sliding clamp in bacteria, archaea and eukarya

The replication of DNA is an essential process in all domains of life. A protein often involved in replication is the sliding clamp. The sliding clamp encircles the DNA and helps replicative polymerase stay attached to the replication machinery increasing the processivity of the polymerase. In eukaryotes and archaea, the sliding clamp is called the Proliferating Cell Nuclear Antigen (PCNA) and consists of two domains. This PCNA forms a trimer encircling the DNA as a hexamer. In bacteria, the structure of the sliding clamp is highly conserved, but the protein itself, called beta clamp, contains three domains, which dimerize to form a hexamer. The bulk of literature touts a conservation of the structure of the sliding clamp, but fails to recognize the conservation of protein sequence among sliding clamps. In this paper, we have used PSI blast to the second iteration in NCBI to show a statistically significant sequence homology between Pyrococcus furiosus PCNA and Kallipyga gabonensis beta clamp. The last two domains of beta clamp align with the two domains of PCNA. This homology data demonstrates that PCNA and beta clamp arose from a common ancestor. In this paper, we have further used beta clamp and PCNA sequences from diverse bacteria, archaea and eukarya to build maximum likelihood phylogenetic tree. Most, but not all, species in different domains of life harbor one sliding clamp from vertical inheritance. Some of these species that have two or more sliding clamps have acquired them from gene duplication or horizontal gene transfer events.

Introduction-First Billion Years: Habitability

The physical processes active during the first billion years (FBY) of Earth's history, such as accretion, differentiation, and impact cratering, provide constraints on the initial conditions that were conducive to the formation and establishment of life on Earth. This motivated the Lunar and Planetary Institute's FBY topical initiative, which was a four-part conference series intended to look at each of these physical processes to study the basic structure and composition of our Solar System that was set during the FBY. The FBY Habitability conference, held in September 2019, was the last in this series and was intended to synthesize the initiative; specifically, to further our understanding of the origins of life, planetary and environmental habitability, and the search for life beyond Earth. The conference included discussions of planetary habitability and the potential emergence of life on bodies within our Solar System, as well as extrasolar systems by applying our knowledge of the Solar System's FBY, and in particular Earth's early history. To introduce this Special Collection, which resulted from work discussed at the conference, we provide a review of the main themes and a synopsis of the FBY Habitability conference.

Structural Modeling of the Treponema pallidum Outer Membrane Protein Repertoire: a Road Map for Deconvolution of Syphilis Pathogenesis and Development of a Syphilis Vaccine

Treponema pallidum, an obligate human pathogen, has an outer membrane (OM) whose physical properties, ultrastructure, and composition differ markedly from those of phylogenetically distant Gram-negative bacteria. We developed structural models for the outer membrane protein (OMP) repertoire (OMPeome) of T. pallidum Nichols using solved Gram-negative structures, computational tools, and small-angle X-ray scattering (SAXS) of selected recombinant periplasmic domains. The T. pallidum "OMPeome" harbors two "stand-alone" proteins (BamA and LptD) involved in OM biogenesis and four paralogous families involved in the influx/efflux of small molecules: 8-stranded β-barrels, long-chain-fatty-acid transporters (FadLs), OM factors (OMFs) for efflux pumps, and T. pallidum repeat proteins (Tprs). BamA (TP0326), the central component of a β-barrel assembly machine (BAM)/translocation and assembly module (TAM) hybrid, possesses a highly flexible polypeptide-transport-associated (POTRA) 1-5 arm predicted to interact with TamB (TP0325). TP0515, an LptD ortholog, contains a novel, unstructured C-terminal domain that models inside the β-barrel. T. pallidum has four 8-stranded β-barrels, each containing positively charged extracellular loops that could contribute to pathogenesis. Three of five FadL-like orthologs have a novel α-helical, presumptively periplasmic C-terminal extension. SAXS and structural modeling further supported the bipartite membrane topology and tridomain architecture of full-length members of the Tpr family. T. pallidum's two efflux pumps presumably extrude noxious small molecules via four coexpressed OMFs with variably charged tunnels. For BamA, LptD, and OMFs, we modeled the molecular machines that deliver their substrates into the OM or external milieu. The spirochete's extended families of OM transporters collectively confer a broad capacity for nutrient uptake. The models also furnish a structural road map for vaccine development. IMPORTANCE The unusual outer membrane (OM) of T. pallidum, the syphilis spirochete, is the ultrastructural basis for its well-recognized capacity for invasiveness, immune evasion, and persistence. In recent years, we have made considerable progress in identifying T. pallidum's repertoire of OMPs. Here, we developed three-dimensional (3D) models for the T. pallidum Nichols OMPeome using structural modeling, bioinformatics, and solution scattering. The OM contains three families of OMP transporters, an OMP family involved in the extrusion of noxious molecules, and two "stand-alone" proteins involved in OM biogenesis. This work represents a major advance toward elucidating host-pathogen interactions during syphilis; understanding how T. pallidum, an extreme auxotroph, obtains a wide array of biomolecules from its obligate human host; and developing a vaccine with global efficacy.

Probing the Role of Chirality in Phospholipid Membranes

Nucleotides, amino acids, sugars, and lipids are almost ubiquitously homochiral within individual cells on Earth. While oligonucleotides and proteins exist as one natural chirality throughout the tree of life, two stereoisomers of phospholipids have separately emerged in archaea and bacteria, an evolutionary divergence known as "the lipid divide". Within this review, we focus on the emergence of phospholipid homochirality and compare the stability of synthetic homochiral and heterochiral membranes in vitro. We discuss chemical probes designed to study the stereospecific interactions of lipid membranes in vitro. Overall, we aim to highlight studies that help elucidate the determinants of stereospecific interactions between lipids, peptides, and small molecule ligands. Continued work in understanding the drivers of favorable interactions between chiral molecules and biological membranes will lead to the design of increasingly selective chemical tools for bioorthogonal labeling of lipid membranes and safer membrane-associating pharmaceuticals.

Advancing knowledge of the plant nuclear periphery and its application for crop science

In this review, we explore recent advances in knowledge of the structure and dynamics of the plant nuclear envelope. As a paradigm, we focused our attention on the Linker of Nucleoskeleton and Cytoskeleton (LINC) complex, a structurally conserved bridging complex comprising SUN domain proteins in the inner nuclear membrane and KASH domain proteins in the outer nuclear membrane. Studies have revealed that this bridging complex has multiple functions with structural roles in positioning the nucleus within the cell, conveying signals across the membrane and organizing chromatin in the 3D nuclear space with impact on gene transcription. We also provide an up-to-date survey in nuclear dynamics research achieved so far in the model plant Arabidopsis thaliana that highlights its potential impact on several key plant functions such as growth, seed maturation and germination, reproduction and response to biotic and abiotic stress. Finally, we bring evidences that most of the constituents of the LINC Complex and associated components are, with some specificities, conserved in monocot and dicot crop species and are displaying very similar functions to those described for Arabidopsis. This leads us to suggest that a better knowledge of this system and a better account of its potential applications will in the future enhance the resilience and productivity of crop plants.

Phylogenetic analysis of mutational robustness based on codon usage supports that the standard genetic code does not prefer extreme environments

The mutational robustness of the genetic code is rarely discussed in the context of biological diversity, such as codon usage and related factors, often considered as independent of the actual organism's proteome. Here we put the living beings back to picture and use distortion as a metric of mutational robustness. Distortion estimates the expected severities of non-synonymous mutations measuring it by amino acid physicochemical properties and weighting for codon usage. Using the biological variance of codon frequencies, we interpret the mutational robustness of the standard genetic code with regards to their corresponding environments and genomic compositions (GC-content). Employing phylogenetic analyses, we show that coding fidelity in physicochemical properties can deteriorate with codon usages adapted to extreme environments and these putative effects are not the artefacts of phylogenetic bias. High temperature environments select for codon usages with decreased mutational robustness of hydrophobic, volumetric, and isoelectric properties. Selection at high saline concentrations also leads to reduced fidelity in polar and isoelectric patterns. These show that the genetic code performs best with mesophilic codon usages, strengthening the view that LUCA or its ancestors preferred lower temperature environments. Taxonomic implications, such as rooting the tree of life, are also discussed.

Early Evolution of Transcription Systems and Divergence of Archaea and Bacteria

DNA template-dependent multi-subunit RNA polymerases (RNAPs) found in all three domains of life and some viruses are of the two-double-Ψ-β-barrel (DPBB) type. The 2-DPBB protein format is also found in some RNA template-dependent RNAPs and a major replicative DNA template-dependent DNA polymerase (DNAP) from Archaea (PolD). The 2-DPBB family of RNAPs and DNAPs probably evolved prior to the last universal common cellular ancestor (LUCA). Archaeal Transcription Factor B (TFB) and bacterial σ factors include homologous strings of helix-turn-helix units. The consequences of TFB-σ homology are discussed in terms of the evolution of archaeal and bacterial core promoters. Domain-specific DPBB loop inserts functionally connect general transcription factors to the RNAP active site. Archaea appear to be more similar to LUCA than Bacteria. Evolution of bacterial σ factors from TFB appears to have driven divergence of Bacteria from Archaea, splitting the prokaryotic domains.

Expanded diversity of Asgard archaea and their relationships with eukaryotes

Asgard is a recently discovered superphylum of archaea that appears to include the closest archaeal relatives of eukaryotes1-5. Debate continues as to whether the archaeal ancestor of eukaryotes belongs within the Asgard superphylum or whether this ancestor is a sister group to all other archaea (that is, a two-domain versus a three-domain tree of life)6-8. Here we present a comparative analysis of 162 complete or nearly complete genomes of Asgard archaea, including 75 metagenome-assembled genomes that-to our knowledge-have not previously been reported. Our results substantially expand the phylogenetic diversity of Asgard and lead us to propose six additional phyla that include a deep branch that we have provisionally named Wukongarchaeota. Our phylogenomic analysis does not resolve unequivocally the evolutionary relationship between eukaryotes and Asgard archaea, but instead-depending on the choice of species and conserved genes used to build the phylogeny-supports either the origin of eukaryotes from within Asgard (as a sister group to the expanded Heimdallarchaeota-Wukongarchaeota branch) or a deeper branch for the eukaryote ancestor within archaea. Our comprehensive protein domain analysis using the 162 Asgard genomes results in a major expansion of the set of eukaryotic signature proteins. The Asgard eukaryotic signature proteins show variable phyletic distributions and domain architectures, which is suggestive of dynamic evolution through horizontal gene transfer, gene loss, gene duplication and domain shuffling. The phylogenomics of the Asgard archaea points to the accumulation of the components of the mobile archaeal 'eukaryome' in the archaeal ancestor of eukaryotes (within or outside Asgard) through extensive horizontal gene transfer.

The Asgard Archaeal-Unique Contribution to Protein Families of the Eukaryotic Common Ancestor Was 0.3

The identification of the asgard archaea has fueled speculations regarding the nature of the archaeal host in eukaryogenesis and its level of complexity prior to endosymbiosis. Here, we analyzed the coding capacity of 150 eukaryotes, 1,000 bacteria, and 226 archaea, including the only cultured member of the asgard archaea. Clustering methods that consistently recover endosymbiotic contributions to eukaryotic genomes recover an asgard archaeal-unique contribution of a mere 0.3% to protein families present in the last eukaryotic common ancestor, while simultaneously suggesting that this group's diversity rivals that of all other archaea combined. The number of homologs shared exclusively between asgard archaea and eukaryotes is only 27 on average. This tiny asgard archaeal-unique contribution to the root of eukaryotic protein families questions claims that archaea evolved complexity prior to eukaryogenesis. Genomic and cellular complexity remains a eukaryote-specific feature and is best understood as the archaeal host's solution to housing an endosymbiont.

The past, present and future of the tree of life

The advent of comparative genomics in the late 1990s led to the discovery of extensive lateral gene transfer in prokaryotes. The resulting debate over whether life as a whole is best represented as a tree or a network has since given way to a general consensus in which trees and networks co-exist rather than stand in opposition. Embracing this consensus allows us to move beyond the question of which is true or false. The future of the tree of life debate lies in asking what trees and networks can, and should, do for science.

Reconstructing the Evolutionary History of a Highly Conserved Operon Cluster in Gammaproteobacteria and Bacilli

The evolution of gene order rearrangements within bacterial chromosomes is a fast process. Closely related species can have almost no conservation in long-range gene order. A prominent exception to this rule is a >40 kb long cluster of five core operons (secE-rpoBC-str-S10-spc-alpha) and three variable adjacent operons (cysS, tufB, and ecf) that together contain 57 genes of the transcriptional and translational machinery. Previous studies have indicated that at least part of this operon cluster might have been present in the last common ancestor of bacteria and archaea. Using 204 whole genome sequences, ∼2 Gy of evolution of the operon cluster were reconstructed back to the last common ancestors of the Gammaproteobacteria and of the Bacilli. A total of 163 independent evolutionary events were identified in which the operon cluster was altered. Further examination showed that the process of disconnecting two operons generally follows the same pattern. Initially, a small number of genes is inserted between the operons breaking the concatenation followed by a second event that fully disconnects the operons. While there is a general trend for loss of gene synteny over time, there are examples of increased alteration rates at specific branch points or within specific bacterial orders. This indicates the recurrence of relaxed selection on the gene order within bacterial chromosomes. The analysis of the alternation events indicates that segmental genome duplications and/or transposon-directed recombination play a crucial role in rearrangements of the operon cluster.

Evolution of ribosomal protein network architectures

To perform an accurate protein synthesis, ribosomes accomplish complex tasks involving the long-range communication between its functional centres such as the peptidyl transfer centre, the tRNA bindings sites and the peptide exit tunnel. How information is transmitted between these sites remains one of the major challenges in current ribosome research. Many experimental studies have revealed that some r-proteins play essential roles in remote communication and the possible involvement of r-protein networks in these processes have been recently proposed. Our phylogenetic, structural and mathematical study reveals that of the three kingdom's r-protein networks converged towards non-random graphs where r-proteins collectively coevolved to optimize interconnection between functional centres. The massive acquisition of conserved aromatic residues at the interfaces and along the extensions of the newly connected eukaryotic r-proteins also highlights that a strong selective pressure acts on their sequences probably for the formation of new allosteric pathways in the network.

Comparative patterns of modified nucleotides in individual tRNA species from a mesophilic and two thermophilic archaea

To improve and complete our knowledge of archaeal tRNA modification patterns, we have identified and compared the modification pattern (type and location) in tRNAs of three very different archaeal species, Methanococcus maripaludis (a mesophilic methanogen), Pyrococcus furiosus (a hyperthermophile thermococcale), and Sulfolobus acidocaldarius (an acidophilic thermophilic sulfolobale). Most abundant isoacceptor tRNAs (79 in total) for each of the 20 amino acids were isolated by two-dimensional gel electrophoresis followed by in-gel RNase digestions. The resulting oligonucleotide fragments were separated by nanoLC and their nucleotide content analyzed by mass spectrometry (MS/MS). Analysis of total modified nucleosides obtained from complete digestion of bulk tRNAs was also performed. Distinct base- and/or ribose-methylations, cytidine acetylations, and thiolated pyrimidines were identified, some at new positions in tRNAs. Novel, some tentatively identified, modifications were also found. The least diversified modification landscape is observed in the mesophilic Methanococcus maripaludis and the most complex one in Sulfolobus acidocaldarius Notable observations are the frequent occurrence of ac⁴C nucleotides in thermophilic archaeal tRNAs, the presence of m⁷G at positions 1 and 10 in Pyrococcus furiosus tRNAs, and the use of wyosine derivatives at position 37 of tRNAs, especially those decoding U1- and C1-starting codons. These results complete those already obtained by others with sets of archaeal tRNAs from Methanocaldococcus jannaschii and Haloferax volcanii.

G-Quadruplexes in the Archaea Domain

The importance of unusual DNA structures in the regulation of basic cellular processes is an emerging field of research. Amongst local non-B DNA structures, G-quadruplexes (G4s) have gained in popularity during the last decade, and their presence and functional relevance at the DNA and RNA level has been demonstrated in a number of viral, bacterial, and eukaryotic genomes, including humans. Here, we performed the first systematic search of G4-forming sequences in all archaeal genomes available in the NCBI database. In this article, we investigate the presence and locations of G-quadruplex forming sequences using the G4Hunter algorithm. G-quadruplex-prone sequences were identified in all archaeal species, with highly significant differences in frequency, from 0.037 to 15.31 potential quadruplex sequences per kb. While G4 forming sequences were extremely abundant in Hadesarchaea archeon (strikingly, more than 50% of the Hadesarchaea archaeon isolate WYZ-LMO6 genome is a potential part of a G4-motif), they were very rare in the Parvarchaeota phylum. The presence of G-quadruplex forming sequences does not follow a random distribution with an over-representation in non-coding RNA, suggesting possible roles for ncRNA regulation. These data illustrate the unique and non-random localization of G-quadruplexes in Archaea.

A distinct lineage of Caudovirales that encodes a deeply branching multi-subunit RNA polymerase

Bacteriophages play critical roles in the biosphere, but their vast genomic diversity has obscured their evolutionary origins, and phylogenetic analyses have traditionally been hindered by their lack of universal phylogenetic marker genes. In this study we mine metagenomic data and identify a clade of Caudovirales that encodes the β and β' subunits of multi-subunit RNA polymerase (RNAP), a high-resolution phylogenetic marker which enables detailed evolutionary analyses. Our RNAP phylogeny revealed that the Caudovirales RNAP forms a clade distinct from cellular homologs, suggesting an ancient acquisition of this enzyme. Within these multimeric RNAP-encoding Caudovirales (mReC), we find that the similarity of major capsid proteins and terminase large subunits further suggests they form a distinct clade with common evolutionary origin. Our study characterizes a clade of RNAP-encoding Caudovirales and suggests the ancient origin of this enzyme in this group, underscoring the important role of viruses in the early evolution of life on Earth.

Plant-archaea relationships: a potential means to improve crop production in arid and semi-arid regions

Crop production in arid and semi-arid regions of the world is limited by several abiotic factors, including water stress, temperature extremes, low soil fertility, high soil pH, low soil water-holding capacity, and low soil organic matter. Moreover, arid and semi-arid areas experience low levels of rainfall with high spatial and temporal variability. Also, the indiscriminate use of chemicals, a practice that characterizes current agricultural practice, promotes crop and soil pollution potentially resulting in serious human health and environmental hazards. A reliable and sustainable alternative to current farming practice is, therefore, a necessity. One such option includes the use of plant growth-promoting microbes that can help to ameliorate some of the adverse effects of these multiple stresses. In this regard, archaea, functional components of the plant microbiome that are found both in the rhizosphere and the endosphere may contribute to the promotion of plant growth. Archaea can survive in extreme habitats such as areas with high temperatures and hypersaline water. No cases of archaea pathogenicity towards plants have been reported. Archaea appear to have the potential to promote plant growth, improve nutrient supply and protect plants against various abiotic stresses. A better understanding of recent developments in archaea functional diversity, plant colonizing ability, and modes of action could facilitate their eventual usage as reliable components of sustainable agricultural systems. The research discussed herein, therefore, addresses the potential role of archaea to improve sustainable crop production in arid and semi-arid areas.

Descent of Bacteria and Eukarya From an Archaeal Root of Life

The 3 biological domains delineated based on small subunit ribosomal RNAs (SSU rRNAs) are confronted by uncertainties regarding the relationship between Archaea and Bacteria, and the origin of Eukarya. The similarities between the paralogous valyl-tRNA and isoleucyl-tRNA synthetases in 5398 species estimated by BLASTP, which decreased from Archaea to Bacteria and further to Eukarya, were consistent with vertical gene transmission from an archaeal root of life close to Methanopyrus kandleri through a Primitive Archaea Cluster to an Ancestral Bacteria Cluster, and to Eukarya. The predominant similarities of the ribosomal proteins (rProts) of eukaryotes toward archaeal rProts relative to bacterial rProts established that an archaeal parent rather than a bacterial parent underwent genome merger with bacteria to generate eukaryotes with mitochondria. Eukaryogenesis benefited from the predominantly archaeal accelerated gene adoption (AGA) phenotype pertaining to horizontally transferred genes from other prokaryotes and expedited genome evolution via both gene-content mutations and nucleotidyl mutations. Archaeons endowed with substantial AGA activity were accordingly favored as candidate archaeal parents. Based on the top similarity bitscores displayed by their proteomes toward the eukaryotic proteomes of Giardia and Trichomonas, and high AGA activity, the Aciduliprofundum archaea were identified as leading candidates of the archaeal parent. The Asgard archaeons and a number of bacterial species were among the foremost potential contributors of eukaryotic-like proteins to Eukarya.

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.

PaSiT: a novel approach based on short-oligonucleotide frequencies for efficient bacterial identification and typing

MOTIVATION

One of the most widespread methods used in taxonomy studies to distinguish between strains or taxa is the calculation of average nucleotide identity. It requires a computationally expensive alignment step and is therefore not suitable for large-scale comparisons. Short oligonucleotide-based methods do offer a faster alternative but at the expense of accuracy. Here, we aim to address this shortcoming by providing a software that implements a novel method based on short-oligonucleotide frequencies to compute inter-genomic distances.

RESULTS

Our tetranucleotide and hexanucleotide implementations, which were optimized based on a taxonomically well-defined set of over 200 newly sequenced bacterial genomes, are as accurate as the short oligonucleotide-based method TETRA and average nucleotide identity, for identifying bacterial species and strains, respectively. Moreover, the lightweight nature of this method makes it applicable for large-scale analyses.

AVAILABILITY AND IMPLEMENTATION

The method introduced here was implemented, together with other existing methods, in a dependency-free software written in C, GenDisCal, available as source code from https://github.com/LM-UGent/GenDisCal. The software supports multithreading and has been tested on Windows and Linux (CentOS). In addition, a Java-based graphical user interface that acts as a wrapper for the software is also available.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

The Evolutionary History of Ephs and Ephrins: Toward Multicellular Organisms

Eph receptor (Eph) and ephrin signaling regulate fundamental developmental processes through both forward and reverse signaling triggered upon cell-cell contact. In vertebrates, they are both classified into classes A and B, and some representatives have been identified in many metazoan groups, where their expression and functions have been well studied. We have extended previous phylogenetic analyses and examined the presence of Eph and ephrins in the tree of life to determine their origin and evolution. We have found that 1) premetazoan choanoflagellates may already have rudimental Eph/ephrin signaling as they have an Eph-/ephrin-like pair and homologs of downstream-signaling genes; 2) both forward- and reverse-downstream signaling might already occur in Porifera since sponges have most genes involved in these types of signaling; 3) the nonvertebrate metazoan Eph is a type-B receptor that can bind ephrins regardless of their membrane-anchoring structure, glycosylphosphatidylinositol, or transmembrane; 4) Eph/ephrin cross-class binding is specific to Gnathostomata; and 5) kinase-dead Eph receptors can be traced back to Gnathostomata. We conclude that Eph/ephrin signaling is of older origin than previously believed. We also examined the presence of protein domains associated with functional characteristics and the appearance and conservation of downstream-signaling pathways to understand the original and derived functions of Ephs and ephrins. We find that the evolutionary history of these gene families points to an ancestral function in cell-cell interactions that could contribute to the emergence of multicellularity and, in particular, to the required segregation of cell populations.

PhoH2 proteins couple RNA helicase and RNAse activities

PhoH2 proteins are found in a very diverse range of microorganisms that span bacteria and archaea. These proteins are composed of two domains: an N-terminal PIN-domain fused with a C-terminal PhoH domain. Collectively this fusion functions as an RNA helicase and ribonuclease. In other genomic contexts, PINdomains and PhoHdomains are separate but adjacent suggesting association to achieve similar function. Exclusively among the mycobacteria, PhoH2 proteins are encoded in the genome with an upstream gene, phoAT, which is thought to play the role of an antitoxin (in place of the traditional VapB antitoxin that lies upstream of the 47 other PINdomains in the mycobacterial genome). This review examines PhoH2 proteins as a whole and describes the bioinformatics, biochemical, structural, and biological properties of the two domains that make up PhoH2: PIN and PhoH. We review the transcriptional regulators of phoH2 from two mycobacterial species and speculate on the function of PhoH2 proteins in the context of a Type II toxin-antitoxin system which are thought to play a role in the stress response in bacteria.

Somatic multicellularity as a satisficing solution to the prediction-error minimization problem

Adaptive success in the biosphere requires the dynamic ability to adjust physiological, transcriptional, and behavioral responses to environmental conditions. From chemical networks to organisms to whole communities, biological entities at all levels of organization seek to optimize their predictive power. Here, we argue that this fundamental drive provides a novel perspective on the origin of multicellularity. One way for unicellular organisms to minimize surprise with respect to external inputs is to be surrounded by reproductively-disabled, i.e. somatic copies of themselves - highly predictable agents which in effect reduce uncertainty in their microenvironments. We show that the transition to multicellularity can be modeled as a phase transition driven by environmental threats. We present modeling results showing how multicellular bodies can arise if non-reproductive somatic cells protect their reproductive parents from environmental lethality. We discuss how a somatic body can be interpreted as a Markov blanket around one or more reproductive cells, and how the transition to somatic multicellularity can be represented as a transition from exposure of reproductive cells to a high-uncertainty environment to their protection from environmental uncertainty by this Markov blanket. This is, effectively, a transition by the Markov blanket from transparency to opacity for the variational free energy of the environment. We suggest that the ability to arrest the cell cycle of daughter cells and redirect their resource utilization from division to environmental threat amelioration is the key innovation of obligate multicellular eukaryotes, that the nervous system evolved to exercise this control over long distances, and that cancer is an escape by somatic cells from the control of reproductive cells. Our quantitative model illustrates the evolutionary dynamics of this system, provides a novel hypothesis for the origin of multicellular animal bodies, and suggests a fundamental link between the architectures of complex organisms and information processing in proto-cognitive cellular agents.