Asgard archaea illuminate the origin of eukaryotic cellular complexity

Isolation of an archaeon at the prokaryote-eukaryote interface

The origin of eukaryotes remains unclear1-4. Current data suggest that eukaryotes may have emerged from an archaeal lineage known as 'Asgard' archaea5,6. Despite the eukaryote-like genomic features that are found in these archaea, the evolutionary transition from archaea to eukaryotes remains unclear, owing to the lack of cultured representatives and corresponding physiological insights. Here we report the decade-long isolation of an Asgard archaeon related to Lokiarchaeota from deep marine sediment. The archaeon-'Candidatus Prometheoarchaeum syntrophicum' strain MK-D1-is an anaerobic, extremely slow-growing, small coccus (around 550 nm in diameter) that degrades amino acids through syntrophy. Although eukaryote-like intracellular complexes have been proposed for Asgard archaea6, the isolate has no visible organelle-like structure. Instead, Ca. P. syntrophicum is morphologically complex and has unique protrusions that are long and often branching. On the basis of the available data obtained from cultivation and genomics, and reasoned interpretations of the existing literature, we propose a hypothetical model for eukaryogenesis, termed the entangle-engulf-endogenize (also known as E3) model.

The molecular foundations of zygosis

Zygosis is the generation of new biological individuals by the sexual fusion of gamete cells. Our current understanding of eukaryotic phylogeny indicates that sex is ancestral to all extant eukaryotes. Although sexual development is extremely diverse, common molecular elements have been retained. HAP2-GCS1, a protein that promotes the fusion of gamete cell membranes that is related in structure to certain viral fusogens, is conserved in many eukaryotic lineages, even though gametes vary considerably in form and behaviour between species. Similarly, although zygotes have dramatically different forms and fates in different organisms, diverse eukaryotes share a common developmental programme in which homeodomain-containing transcription factors play a central role. These common mechanistic elements suggest possible common evolutionary histories that, if correct, would have profound implications for our understanding of eukaryogenesis.

An Introduction to Phagocytosis

Phagocytosis is usually defined as the cellular process by which cells internalise particulate matter larger than about 0.5 μm in diameter. It is an endocytic process, distinct from pinocytosis and macropinocytosis. These latter processes may internalise small particles suspended the extracellular fluid, but this is a by-product of internalising the fluid, and is not phagocytosis per se. In contrast, phagocytosis is targeted at solid particulates, usually microbes, which are internalised and "digested" either to provide food, or as part of the immune system of higher animals. The mechanism of phagocytosis may have, at its core, many primitive elements, but it is a highly complex and coordinated series of cell biological and molecular events which together result in the uptake of a particle. In this introduction, the basis of phagocytosis and some ideas of its origin are discussed.

Phylogenomics provides robust support for a two-domains tree of life

Hypotheses about the origin of eukaryotic cells are classically framed within the context of a universal 'tree of life' based on conserved core genes. Vigorous ongoing debate about eukaryote origins is based on assertions that the topology of the tree of life depends on the taxa included and the choice and quality of genomic data analysed. Here we have reanalysed the evidence underpinning those claims and apply more data to the question by using supertree and coalescent methods to interrogate >3,000 gene families in archaea and eukaryotes. We find that eukaryotes consistently originate from within the archaea in a two-domains tree when due consideration is given to the fit between model and data. Our analyses support a close relationship between eukaryotes and Asgard archaea and identify the Heimdallarchaeota as the current best candidate for the closest archaeal relatives of the eukaryotic nuclear lineage.

Metabolic activity analyses demonstrate that Lokiarchaeon exhibits homoacetogenesis in sulfidic marine sediments

The genomes of the Asgard superphylum of Archaea hold clues pertaining to the nature of the host cell that acquired the mitochondrion at the origin of eukaryotes^1-4. Representatives of the Asgard candidate phylum Candidatus Lokiarchaeota (Lokiarchaeon) have the capacity for acetogenesis and fermentation^5-7, but how their metabolic activity responds to environmental conditions is poorly understood. Here, we show that in anoxic Namibian shelf sediments, Lokiarchaeon gene expression levels are higher than those of bacterial phyla and increase with depth below the seafloor. Lokiarchaeon gene expression was significantly different across a hypoxic-sulfidic redox gradient, whereby genes involved in growth, fermentation and H₂-dependent carbon fixation had the highest expression under the most reducing (sulfidic) conditions. Quantitative stable isotope probing revealed that anaerobic utilization of CO₂ and diatomaceous extracellular polymeric substances by Lokiarchaeon was higher than the bacterial average, consistent with higher expression of Lokiarchaeon genes, including those involved in transport and fermentation of sugars and amino acids. The quantitative stable isotope probing and gene expression data demonstrate homoacetogenic activity of Candidatus Lokiarchaeota, whereby fermentative H₂ production from organic substrates is coupled with the Wood-Ljungdahl carbon fixation pathway⁸. The high energetic efficiency provided by homoacetogenesis⁸ helps to explain the elevated metabolic activity of Lokiarchaeon in this anoxic, energy-limited setting.

Phagocytosis-like cell engulfment by a planctomycete bacterium

Phagocytosis is a key eukaryotic feature, conserved from unicellular protists to animals, that enabled eukaryotes to feed on other organisms. It could also be a driving force behind endosymbiosis, a process by which α-proteobacteria and cyanobacteria evolved into mitochondria and plastids, respectively. Here we describe a planctomycete bacterium, 'Candidatus Uab amorphum', which is able to engulf other bacteria and small eukaryotic cells through a phagocytosis-like mechanism. Observations via light and electron microscopy suggest that this bacterium digests prey cells in specific compartments. With the possible exception of a gene encoding an actin-like protein, analysis of the 'Ca. Uab amorphum' genomic sequence does not reveal any genes homologous to eukaryotic phagocytosis genes, suggesting that cell engulfment in this microorganism is probably not homologous to eukaryotic phagocytosis. The discovery of this "phagotrophic" bacterium expands our understanding of the cellular complexity of prokaryotes, and may be relevant to the origin of eukaryotic cells.

PhySpeTree: an automated pipeline for reconstructing phylogenetic species trees

Background

Phylogenetic species trees are widely used in inferring evolutionary relationships. Existing software and algorithms mainly focus on phylogenetic inference. However, less attention has been paid to intermediate steps, such as processing extremely large sequences and preparing configure files to connect multiple software. When the species number is large, the intermediate steps become a bottleneck that may seriously affect the efficiency of tree building.

Results

Here, we present an easy-to-use pipeline named PhySpeTree to facilitate the reconstruction of species trees across bacterial, archaeal, and eukaryotic organisms. Users need only to input the abbreviations of species names; PhySpeTree prepares complex configure files for different software, then automatically downloads genomic data, cleans sequences, and builds trees. PhySpeTree allows users to perform critical steps such as sequence alignment and tree construction by adjusting advanced options. PhySpeTree provides two parallel pipelines based on concatenated highly conserved proteins and small subunit ribosomal RNA sequences, respectively. Accessory modules, such as those for inserting new species, generating visualization configurations, and combining trees, are distributed along with PhySpeTree.

Conclusions

Together with accessory modules, PhySpeTree significantly simplifies tree reconstruction. PhySpeTree is implemented in Python running on modern operating systems (Linux, macOS, and Windows). The source code is freely available with detailed documentation (https://github.com/yangfangs/physpetools).

Phylogenomics of 10,575 genomes reveals evolutionary proximity between domains Bacteria and Archaea

Rapid growth of genome data provides opportunities for updating microbial evolutionary relationships, but this is challenged by the discordant evolution of individual genes. Here we build a reference phylogeny of 10,575 evenly-sampled bacterial and archaeal genomes, based on a comprehensive set of 381 markers, using multiple strategies. Our trees indicate remarkably closer evolutionary proximity between Archaea and Bacteria than previous estimates that were limited to fewer "core" genes, such as the ribosomal proteins. The robustness of the results was tested with respect to several variables, including taxon and site sampling, amino acid substitution heterogeneity and saturation, non-vertical evolution, and the impact of exclusion of candidate phyla radiation (CPR) taxa. Our results provide an updated view of domain-level relationships.

On the Origin of Iron/Sulfur Cluster Biosynthesis in Eukaryotes

Iron and sulfur are indispensable elements of every living cell, but on their own these elements are toxic and require dedicated machineries for the formation of iron/sulfur (Fe/S) clusters. In eukaryotes, proteins requiring Fe/S clusters (Fe/S proteins) are found in or associated with various organelles including the mitochondrion, endoplasmic reticulum, cytosol, and the nucleus. These proteins are involved in several pathways indispensable for the viability of each living cell including DNA maintenance, protein translation and metabolic pathways. Thus, the formation of Fe/S clusters and their delivery to these proteins has a fundamental role in the functions and the evolution of the eukaryotic cell. Currently, most eukaryotes harbor two (located in cytosol and mitochondrion) or three (located in plastid) machineries for the assembly of Fe/S clusters, but certain anaerobic microbial eukaryotes contain sulfur mobilization (SUF) machineries that were previously thought to be present only in archaeal linages. These machineries could not only stipulate which pathway was present in the last eukaryotic common ancestor (LECA), but they could also provide clues regarding presence of an Fe/S cluster machinery in the proto-eukaryote and evolution of Fe/S cluster assembly machineries in all eukaryotes.

The Known, the New, and a Possible Surprise: A Re-Evaluation of the Nucleomorph-Encoded Proteome of Cryptophytes

Nucleomorphs are small nuclei that evolved from the nucleus of former eukaryotic endosymbionts of cryptophytes and chlorarachniophytes. These enigmatic organelles reside in their complex plastids and harbor the smallest and most compacted eukaryotic genomes investigated so far. Although the coding capacity of the nucleomorph genomes is small, a significant percentage of the encoded proteins (predicted nucleomorph-encoded proteins, pNMPs) is still not functionally annotated. We have analyzed pNMPs with unknown functions via Phyre2, a bioinformatic tool for prediction and modeling of protein structure, resulting in a functional annotation of 215 pNMPs out of 826 uncharacterized open reading frames of cryptophytes. The newly annotated proteins are predicted to participate in nucleomorph-specific functions such as chromosome organization and expression, as well as in modification and degradation of nucleomorph-encoded proteins. Additionally, we have functionally assigned nucleomorph-encoded, putatively plastid-targeted proteins among the reinvestigated pNMPs. Hints for a putative function in the periplastid compartment, the cytoplasm surrounding the nucleomorphs, emerge from the identification of pNMPs that might be homologs of endomembrane system-related proteins. These proteins are discussed in respect to their putative functions.

The E. coli MinCDE system in the regulation of protein patterns and gradients

Molecular self-organziation, also regarded as pattern formation, is crucial for the correct distribution of cellular content. The processes leading to spatiotemporal patterns often involve a multitude of molecules interacting in complex networks, so that only very few cellular pattern-forming systems can be regarded as well understood. Due to its compositional simplicity, the Escherichia coli MinCDE system has, thus, become a paradigm for protein pattern formation. This biological reaction diffusion system spatiotemporally positions the division machinery in E. coli and is closely related to ParA-type ATPases involved in most aspects of spatiotemporal organization in bacteria. The ATPase MinD and the ATPase-activating protein MinE self-organize on the membrane as a reaction matrix. In vivo, these two proteins typically oscillate from pole-to-pole, while in vitro they can form a variety of distinct patterns. MinC is a passenger protein supposedly operating as a downstream cue of the system, coupling it to the division machinery. The MinCDE system has helped to extract not only the principles underlying intracellular patterns, but also how they are shaped by cellular boundaries. Moreover, it serves as a model to investigate how patterns can confer information through specific and non-specific interactions with other molecules. Here, we review how the three Min proteins self-organize to form patterns, their response to geometric boundaries, and how these patterns can in turn induce patterns of other molecules, focusing primarily on experimental approaches and developments.

Mitochondria and endomembrane origins

In this Guest Editorial, Heidi McBride introduces our special issue on membranes with a discussion of the contribution of mitochondria to the emergence of the endomembrane system.

Contributions of single-cell genomics to our understanding of planktonic marine archaea

Single-cell genomics has transformed many fields of biology, marine microbiology included. Here, we consider the impact of single-cell genomics on a specific group of marine microbes-the planktonic marine archaea. Despite single-cell enabled discoveries of novel metabolic function in the marine thaumarchaea, population-level investigations are hindered by an overall lower than expected recovery of thaumarchaea in single-cell studies. Metagenome-assembled genomes have so far been a more useful method for accessing genome-resolved insights into the Marine Group II euryarchaea. Future progress in the application of single-cell genomics to archaeal biology in the ocean would benefit from more targeted sorting approaches, and a more systematic investigation of potential biases against archaea in single-cell workflows including cell lysis, genome amplification and genome screening. This article is part of a discussion meeting issue 'Single cell ecology'.

Targeted isolation and cultivation of uncultivated bacteria by reverse genomics

Most microorganisms from all taxonomic levels are uncultured. Single-cell genomes and metagenomes continue to increase the known diversity of Bacteria and Archaea, but while ‘omics can be used to infer physiological or ecological roles for species in a community, most of those hypothetical roles remain unvalidated. Here we report an approach to capture specific microorganisms from complex communities into pure cultures using genome-informed antibody engineering. We apply our reverse genomics approach to isolate and sequence single cells and to cultivate three different species-level lineages of human oral Saccharibacteria/TM7. Using our pure cultures we show that all three saccharibacteria species are epibionts of diverse Actinobacteria. We also isolate and cultivate human oral SR1 bacteria, which are members of a lineage of previously uncultured bacteria. Reverse-genomics-enabled cultivation of microorganisms can be applied to any species from any environment and has the potential to unlock the isolation, cultivation and characterization of species from as-yet-uncultured branches of the microbial tree of life.

Global Transcriptional Programs in Archaea Share Features with the Eukaryotic Environmental Stress Response

The environmental stress response (ESR), a global transcriptional program originally identified in yeast, is characterized by a rapid and transient transcriptional response composed of large, oppositely regulated gene clusters. Genes induced during the ESR encode core components of stress tolerance, macromolecular repair, and maintenance of homeostasis. In this review, we investigate the possibility for conservation of the ESR across the eukaryotic and archaeal domains of life. We first re-analyze existing transcriptomics data sets to illustrate that a similar transcriptional response is identifiable in Halobacterium salinarum, an archaeal model organism. To substantiate the archaeal ESR, we calculated gene-by-gene correlations, gene function enrichment, and comparison of temporal dynamics. We note reported examples of variation in the ESR across fungi, then synthesize high-level trends present in expression data of other archaeal species. In particular, we emphasize the need for additional high-throughput time series expression data to further characterize stress-responsive transcriptional programs in the Archaea. Together, this review explores an open question regarding features of global transcriptional stress response programs shared across domains of life.

A polyamine-independent role for S-adenosylmethionine decarboxylase

The only known function of S-adenosylmethionine decarboxylase (AdoMetDC) is to supply, with its partner aminopropyltransferase enzymes such as spermidine synthase (SpdSyn), the aminopropyl donor for polyamine biosynthesis. Polyamine spermidine is probably essential for the growth of all eukaryotes, most archaea and many bacteria. Two classes of AdoMetDC exist, the prokaryotic class 1a and 1b forms, and the eukaryotic class 2 enzyme, which is derived from an ancient fusion of two prokaryotic class 1b genes. Herein, we show that ‘eukaryotic' class 2 AdoMetDCs are found in bacteria and are enzymatically functional. However, the bacterial AdoMetDC class 2 genes are phylogenetically limited and were likely acquired from a eukaryotic source via transdomain horizontal gene transfer, consistent with the class 2 form of AdoMetDC being a eukaryotic invention. We found that some class 2 and thousands of class 1b AdoMetDC homologues are present in bacterial genomes that also encode a gene fusion of an N-terminal membrane protein of the Major Facilitator Superfamily (MFS) class of transporters and a C-terminal SpdSyn-like domain. Although these AdoMetDCs are enzymatically functional, spermidine is absent, and an entire fusion protein or its SpdSyn-like domain only, does not biochemically complement a SpdSyn deletion strain of E. coli. This suggests that the fusion protein aminopropylates a substrate other than putrescine, and has a role outside of polyamine biosynthesis. Another integral membrane protein found clustered with these genes is DUF350, which is also found in other gene clusters containing a homologue of the glutathionylspermidine synthetase family and occasionally other polyamine biosynthetic enzymes.

The Role of Archaeal Chromatin in Transcription

Genomic organization impacts accessibility and movement of information processing systems along DNA. DNA-bound proteins dynamically dictate gene expression and provide regulatory potential to tune transcription rates to match ever-changing environmental conditions. Archaeal genomes are typically small, circular, gene dense, and organized either by histone proteins that are homologous to their eukaryotic counterparts, or small basic proteins that function analogously to bacterial nucleoid proteins. We review here how archaeal genomes are organized and how such organization impacts archaeal gene expression, focusing on conserved DNA-binding proteins within the clade and the factors that are known to impact transcription initiation and elongation within protein-bound genomes.

Metagenomes from Coastal Marine Sediments Give Insights into the Ecological Role and Cellular Features of Loki- and Thorarchaeota

The genomes of Asgard Archaea, a novel archaeal proposed superphylum, share an enriched repertoire of eukaryotic signature genes and thus promise to provide insights into early eukaryote evolution. However, the distribution, metabolisms, cellular structures, and ecology of the members within this superphylum are not well understood. Here we provide a meta-analysis of the environmental distribution of the Asgard archaea, based on available 16S rRNA gene sequences. Metagenome sequencing of samples from a salt-crusted lagoon on the Baja California Peninsula of Mexico allowed the assembly of a new Thorarchaeota and three Lokiarchaeota genomes. Comparative analyses of all known Lokiarchaeota and Thorarchaeota genomes revealed overlapping genome content, including central carbon metabolism. Members of both groups contained putative reductive dehalogenase genes, suggesting that these organisms might be able to metabolize halogenated organic compounds. Unlike the first report on Lokiarchaeota, we identified genes encoding glycerol-1-phosphate dehydrogenase in all Loki- and Thorarchaeota genomes, suggesting that these organisms are able to synthesize bona fide archaeal lipids with their characteristic glycerol stereochemistry.IMPORTANCE Microorganisms of the superphylum Asgard Archaea are considered to be the closest living prokaryotic relatives of eukaryotes (including plants and animals) and thus promise to give insights into the early evolution of more complex life forms. However, very little is known about their biology as none of the organisms has yet been cultivated in the laboratory. Here we report on the ecological distribution of Asgard Archaea and on four newly sequenced genomes of the Lokiarchaeota and Thorarchaeota lineages that give insight into possible metabolic features that might eventually help to identify these enigmatic groups of archaea in the environment and to culture them.

Prokaryotic Diversity in Mangrove Sediments across Southeastern China Fundamentally Differs from That in Other Biomes

Mangroves, as a blue carbon reservoir, provide an environment for a variety of microorganisms. Mangroves lie in special locations connecting coastal and estuarine areas and experience fluctuating conditions, which are expected to intensify with climate change, creating a need to better understand the relative roles of stochastic and deterministic processes in shaping microbial community assembly. Here, a study of microbial communities inhabiting mangrove sediments across southeastern China, spanning mangroves in six nature reserves, was conducted. We performed high-throughput DNA sequencing of these samples and compared them with data of 1,370 sediment samples collected from the Earth Microbiome Project (EMP) to compare the microbial diversity of mangroves with that of other biomes. Our results showed that prokaryotic alpha diversity in mangroves was significantly higher than that in other biomes and that microbial beta diversity generally clustered according to biome types. The core operational taxonomic units (OTUs) in mangroves were mostly assigned to Gammaproteobacteria, Deltaproteobacteria, Chloroflexi, and Euryarchaeota The majority of beta nearest-taxon index values were higher than 2, indicating that community assembly in mangroves was better explained through a deterministic process than through a stochastic process. Mean annual precipitation (MAP) and total organic carbon (TOC) were main deterministic factors explaining variation in the microbial community. This study fills a gap in addressing the unique microbial diversity of mangrove ecosystems and their microbial community assembly mechanisms.IMPORTANCE Understanding the underlying mechanisms of microbial community assembly patterns is a vital issue in microbial ecology. Mangroves, as an important and special ecosystem, provide a unique environment for examining the relative importance of stochastic and deterministic processes. We made the first global-scale comparison and found that microbial diversity was significantly different in mangrove sediments compared to that of other biomes. Furthermore, our results suggest that a deterministic process is more important in shaping microbial community assembly in mangroves.

Toward completion of the Earth's proteome: an update a decade later

Protein databases are steadily growing driven by the spread of new more efficient sequencing techniques. This growth is dominated by an increase in redundancy (homologous proteins with various degrees of sequence similarity) and by the incapability to process and curate sequence entries as fast as they are created. To understand these trends and aid bioinformatic resources that might be compromised by the increasing size of the protein sequence databases, we have created a less-redundant protein data set. In parallel, we analyzed the evolution of protein sequence databases in terms of size and redundancy. While the SwissProt database has decelerated its growth mostly because of a focus on increasing the level of annotation of its sequences, its counterpart TrEMBL, much less limited by curation steps, is still in a phase of accelerated growth. However, we predict that before 2020, almost all entries deposited in UniProtKB will be homologous to known proteins. We propose that new sequencing projects can be made more useful if they are driven to sequencing voids, parts of the tree of life far from already sequenced species or model organisms. We show these voids are present in the Archaea and Eukarya domains of life. The approach to the certainty of the redundancy of new protein sequence entries leads to the consideration that most of the protein diversity on Earth has already been described, which we estimate to be of around 3.75 million proteins, revising down the prediction we did a decade ago.

The Vulnerability of Microbial Ecosystems in A Changing Climate: Potential Impact in Shark Bay

The potential impact of climate change on eukaryotes, including humans, has been relatively well described. In contrast, the contribution and susceptibility of microorganisms to a changing climate have, until recently, received relatively less attention. In this review, the importance of microorganisms in the climate change discourse is highlighted. Microorganisms are responsible for approximately half of all primary production on earth, support all forms of macroscopic life whether directly or indirectly, and often persist in “extreme” environments where most other life are excluded. In short, microorganisms are the life support system of the biosphere and therefore must be included in decision making regarding climate change. Any effects climate change will have on microorganisms will inevitably impact higher eukaryotes and the activity of microbial communities in turn can contribute to or alleviate the severity of the changing climate. It is of vital importance that unique, fragile, microbial ecosystems are a focus of research efforts so that their resilience to extreme weather events and climate change are thoroughly understood and that conservation efforts can be implemented as a response. One such ecosystem under threat are the evolutionarily significant microbial mats and stromatolites, such as those present in Shark Bay, Western Australia. Climate change models have suggested the duration and severity of extreme weather events in this region will increase, along with rising temperatures, sea levels, and ocean acidification. These changes could upset the delicate balance that fosters the development of microbial mats and stromatolites in Shark Bay. Thus, the challenges facing Shark Bay microbial communities will be presented here as a specific case study.

Division plane placement in pleomorphic archaea is dynamically coupled to cell shape

One mechanism for achieving accurate placement of the cell division machinery is via Turing patterns, where nonlinear molecular interactions spontaneously produce spatiotemporal concentration gradients. The resulting patterns are dictated by cell shape. For example, the Min system of Escherichia coli shows spatiotemporal oscillation between cell poles, leaving a mid-cell zone for division. The universality of pattern-forming mechanisms in divisome placement is currently unclear. We examined the location of the division plane in two pleomorphic archaea, Haloferax volcanii and Haloarcula japonica, and showed that it correlates with the predictions of Turing patterning. Time-lapse analysis of H. volcanii shows that divisome locations after successive rounds of division are dynamically determined by daughter cell shape. For H. volcanii, we show that the location of DNA does not influence division plane location, ruling out nucleoid occlusion. Triangular cells provide a stringent test for Turing patterning, where there is a bifurcation in division plane orientation. For the two archaea examined, most triangular cells divide as predicted by a Turing mechanism; however, in some cases multiple division planes are observed resulting in cells dividing into three viable progeny. Our results suggest that the division site placement is consistent with a Turing patterning system in these archaea.

Group I introns are widespread in archaea

Group I catalytic introns have been found in bacterial, viral, organellar, and some eukaryotic genomes, but not in archaea. All known archaeal introns are bulge-helix-bulge (BHB) introns, with the exception of a few group II introns. It has been proposed that BHB introns arose from extinct group I intron ancestors, much like eukaryotic spliceosomal introns are thought to have descended from group II introns. However, group I introns have little sequence conservation, making them difficult to detect with standard sequence similarity searches. Taking advantage of recent improvements in a computational homology search method that accounts for both conserved sequence and RNA secondary structure, we have identified 39 group I introns in a wide range of archaeal phyla, including examples of group I introns and BHB introns in the same host gene.

Toward a Synthetic Yeast Endosymbiont with a Minimal Genome

Based on the endosymbiotic theory, one of the key events that occurred during mitochondrial evolution was an extensive loss of nonessential genes from the protomitochondrial endosymbiont genome and transfer of some of the essential endosymbiont genes to the host nucleus. We have developed an approach to recapitulate various aspects of endosymbiont genome minimization using a synthetic system consisting of Escherichia coli endosymbionts within host yeast cells. As a first step, we identified a number of E. coli auxotrophs of central metabolites that can form viable endosymbionts within yeast cells. These studies provide a platform to identify nonessential biosynthetic pathways that can be deleted in the E. coli endosymbionts to investigate the evolutionary adaptations in the host and endosymbiont during the evolution of mitochondria.

Origin and evolution of eukaryotic transcription factors

Transcription factors (TFs) have a central role in genome regulation directing gene transcription through binding specific DNA sequences. Eukaryotic genomes encode a large diversity of TF classes, each defined by unique DNA-interaction domains. Recent advances in genome sequencing and phylogenetic placement of diverse eukaryotic and archaeal species are re-defining the evolutionary history of eukaryotic TFs. The emerging view from a comparative genomics perspective is that the Last Eukaryotic Common Ancestor (LECA) had an extensive repertoire of TFs, most of which represent eukaryotic evolutionary novelties. This burst of TF innovation coincides with the emergence of genomic nuclear segregation and complex chromatin organization.

The MEROPS database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database

The MEROPS database (http://www.ebi.ac.uk/merops/) is an integrated source of information about peptidases, their substrates and inhibitors. The hierarchical classification is: protein-species, family, clan, with an identifier at each level. The MEROPS website moved to the EMBL-EBI in 2017, requiring refactoring of the code-base and services provided. The interface to sequence searching has changed and the MEROPS protein sequence libraries can be searched at the EMBL-EBI with HMMER, FastA and BLASTP. Cross-references have been established between MEROPS and the PANTHER database at both the family and protein-species level, which will help to improve curation and coverage between the resources. Because of the increasing size of the MEROPS sequence collection, in future only sequences of characterized proteins, and from completely sequenced genomes of organisms of evolutionary, medical or commercial significance will be added. As an example, peptidase homologues in four proteomes from the Asgard superphylum of Archaea have been identified and compared to other archaean, bacterial and eukaryote proteomes. This has given insights into the origins and evolution of peptidase families, including an expansion in the number of proteasome components in Asgard archaeotes and as organisms increase in complexity. Novel structures for proteasome complexes in archaea are postulated.

Measuring the impact of gene prediction on gene loss estimates in Eukaryotes by quantifying falsely inferred absences

In recent years it became clear that in eukaryotic genome evolution gene loss is prevalent over gene gain. However, the absence of genes in an annotated genome is not always equivalent to the loss of genes. Due to sequencing issues, or incorrect gene prediction, genes can be falsely inferred as absent. This implies that loss estimates are overestimated and, more generally, that falsely inferred absences impact genomic comparative studies. However, reliable estimates of how prevalent this issue is are lacking. Here we quantified the impact of gene prediction on gene loss estimates in eukaryotes by analysing 209 phylogenetically diverse eukaryotic organisms and comparing their predicted proteomes to that of their respective six-frame translated genomes. We observe that 4.61% of domains per species were falsely inferred to be absent for Pfam domains predicted to have been present in the last eukaryotic common ancestor. Between phylogenetically different categories this estimate varies substantially: for clade-specific loss (ancestral loss) we found 1.30% and for species-specific loss 16.88% to be falsely inferred as absent. For BUSCO 1-to-1 orthologous families, 18.30% were falsely inferred to be absent. Finally, we showed that falsely inferred absences indeed impact loss estimates, with the number of losses decreasing by 11.78%. Our work strengthens the increasing number of studies showing that gene loss is an important factor in eukaryotic genome evolution. However, while we demonstrate that on average inferring gene absences from predicted proteomes is reliable, caution is warranted when inferring species-specific absences.

Old cogs, new tricks: the evolution of gene expression in a chromatin context

Sophisticated gene-regulatory mechanisms probably evolved in prokaryotes billions of years before the emergence of modern eukaryotes, which inherited the same basic enzymatic machineries. However, the epigenomic landscapes of eukaryotes are dominated by nucleosomes, which have acquired roles in genome packaging, mitotic condensation and silencing parasitic genomic elements. Although the molecular mechanisms by which nucleosomes are displaced and modified have been described, just how transcription factors, histone variants and modifications and chromatin regulators act on nucleosomes to regulate transcription is the subject of considerable ongoing study. We explore the extent to which these transcriptional regulatory components function in the context of the evolutionarily ancient role of chromatin as a barrier to processes acting on DNA and how chromatin proteins have diversified to carry out evolutionarily recent functions that accompanied the emergence of differentiation and development in multicellular eukaryotes.

Small Genomes and Big Data: Adaptation of Plastid Genomics to the High-Throughput Era

Plastid genome sequences are becoming more readily available with the increase in high-throughput sequencing, and whole-organelle genetic data is available for algae and plants from across the diversity of photosynthetic eukaryotes. This has provided incredible opportunities for studying species which may not be amenable to in vivo study or genetic manipulation or may not yet have been cultured. Research into plastid genomes has pushed the limits of what can be deduced from genomic information, and in particular genomic information obtained from public databases. In this Review, we discuss how research into plastid genomes has benefitted enormously from the explosion of publicly available genome sequence. We describe two case studies in how using publicly available gene data has supported previously held hypotheses about plastid traits from lineage-restricted experiments across algal and plant diversity. We propose how this approach could be used across disciplines for inferring functional and biological characteristics from genomic approaches, including integration of new computational and bioinformatic approaches such as machine learning. We argue that the techniques developed to gain the maximum possible insight from plastid genomes can be applied across the eukaryotic tree of life.

New Microbial Lineages Capable of Carbon Fixation and Nutrient Cycling in Deep-Sea Sediments of the Northern South China Sea

Metagenomics of marine sediments has uncovered a broad diversity of new uncultured taxa and provided insights into their metabolic capabilities. Here, we detected microbial lineages from a sediment core near the Jiulong methane reef of the northern South China Sea (at 1,100-m depth). Assembly and binning of the metagenomes resulted in 11 genomes (>85% complete) that represented nine distinct phyla, including candidate phyla TA06 and LCP-89, Lokiarchaeota, Heimdallarchaeota, and a newly described globally distributed phylum (B38). The genome of LCP-89 has pathways for nitrate, selenate, and sulfate reduction, suggesting that they may be involved in mediating these important processes. B38 are able to participate in the cycling of hydrogen and selenocompounds. Many of these uncultured microbes may also be capable of autotrophic CO₂ fixation, as exemplified by identification of the Wood-Ljungdahl (W-L) pathway. Genes encoding carbohydrate degradation, W-L pathway, Rnf-dependent energy conservation, and Ni/Fe hydrogenases were detected in the transcriptomes of these novel members. Characterization of these new lineages provides insight to the undescribed branches in the tree of life.IMPORTANCE Sedimentary microorganisms in the South China Sea (SCS) remain largely unknown due to the complexity of sediment communities impacted by continent rifting and extension. Distinct geochemical environments may breed special microbial communities including microbes that are still enigmatic. Functional inference of their metabolisms and transcriptional activity provides insight in the ecological roles and substrate-based interactivity of these uncultured Archaea and Bacteria These microorganisms play different roles in utilizing inorganic carbon and scavenging diverse organic compounds involved in the deep-sea carbon cycle. The genomes recovered here contributed undescribed species to the tree of life and laid the foundation for future study on these novel phyla persisting in marginal sediments of the SCS.

Insights into the evolutionary conserved regulation of Rio ATPase activity

Eukaryotic ribosome biogenesis is a complex dynamic process which requires the action of numerous ribosome assembly factors. Among them, the eukaryotic Rio protein family members (Rio1, Rio2 and Rio3) belong to an ancient conserved atypical protein kinase/ ATPase family required for the maturation of the small ribosomal subunit (SSU). Recent structure–function analyses suggested an ATPase-dependent role of the Rio proteins to regulate their dynamic association with the nascent pre-SSU. However, the evolutionary origin of this feature and the detailed molecular mechanism that allows controlled activation of the catalytic activity remained to be determined. In this work we provide functional evidence showing a conserved role of the archaeal Rio proteins for the synthesis of the SSU in archaea. Moreover, we unravel a conserved RNA-dependent regulation of the Rio ATPases, which in the case of Rio2 involves, at least, helix 30 of the SSU rRNA and the P-loop lysine within the shared RIO domain. Together, our study suggests a ribosomal RNA-mediated regulatory mechanism enabling the appropriate stimulation of Rio2 catalytic activity and subsequent release of Rio2 from the nascent pre-40S particle. Based on our findings we propose a unified release mechanism for the Rio proteins.

Confident phylogenetic identification of uncultured prokaryotes through long read amplicon sequencing of the 16S-ITS-23S rRNA operon

Amplicon sequencing of the 16S rRNA gene is the predominant method to quantify microbial compositions and to discover novel lineages. However, traditional short amplicons often do not contain enough information to confidently resolve their phylogeny. Here we present a cost-effective protocol that amplifies a large part of the rRNA operon and sequences the amplicons with PacBio technology. We tested our method on a mock community and developed a read-curation pipeline that reduces the overall read error rate to 0.18%. Applying our method on four environmental samples, we captured near full-length rRNA operon amplicons from a large diversity of prokaryotes. The method operated at moderately high-throughput (22286-37,850 raw ccs reads) and generated a large amount of putative novel archaeal 23S rRNA gene sequences compared to the archaeal SILVA database. These long amplicons allowed for higher resolution during taxonomic classification by means of long (∼1000 bp) 16S rRNA gene fragments and for substantially more confident phylogenies by means of combined near full-length 16S and 23S rRNA gene sequences, compared to shorter traditional amplicons (250 bp of the 16S rRNA gene). We recommend our method to those who wish to cost-effectively and confidently estimate the phylogenetic diversity of prokaryotes in environmental samples at high throughput.

Deciphering the three-domain architecture in schlafens and the structures and roles of human schlafen12 and serpinB12 in transcriptional regulation

Schlafen proteins are important in cell differentiation and defense against viruses, and yet this family of vertebrate proteins is just beginning to be understood at the molecular level. Here, the three-dimensional architecture and molecular interfaces of human schlafen12 (hSLFN12), which promotes intestinal stem cell differentiation, are analyzed by sequence conservation and structural modeling in light of the functions of its homologs and binding partners. Our analysis shows that the schlafen or divergent AAA ATPase domain described in the N-terminal region of schlafens in databases and the literature is a misannotation. This N-terminal region is conclusively an AlbA_2 DNA/RNA binding domain, forming the conserved core of schlafens and their sequence homologs from bacteria through mammals. Group III schlafens additionally contain a AAA NTPase domain in their C-terminal helicase region. In hSLFN12, we have uncovered a domain matching rho GTPases, which directly follows the AlbA_2 domain in all group II-III schlafens. Potential roles for the GTPase-like domain include antiviral activity and cytoskeletal interactions that contribute to nucleocytoplasmic shuttling and cell polarization during differentiation. Based on features conserved with rSlfn13, the AlbA_2 region in hSLFN12 is likely to bind RNA, possibly as a ribonuclease. We hypothesize that RNA binding by hSLFN12 contributes to an RNA-induced transcriptional silencing/E3 ligase complex, given the functions of hSLFN12's partners, SUV39H1, JMJD6, and PDLIM7. hSLFN12's partner hSerpinB12 may contribute to heterochromatin formation, based on its homology to MENT, or directly regulate transcription via its binding to RNA polymerase II. The analysis presented here provides clear architectural and transcriptional regulation hypotheses to guide experimental design for hSLFN12 and the thousands of schlafens that share its motifs.

Investigating the Origins of Membrane Phospholipid Biosynthesis Genes Using Outgroup-Free Rooting

One of the key differences between Bacteria and Archaea is their canonical membrane phospholipids, which are synthesized by distinct biosynthetic pathways with nonhomologous enzymes. This “lipid divide” has important implications for the early evolution of cells and the type of membrane phospholipids present in the last universal common ancestor. One of the main challenges in studies of membrane evolution is that the key biosynthetic genes are ancient and their evolutionary histories are poorly resolved. This poses major challenges for traditional rooting methods because the only available outgroups are distantly related. Here, we address this issue by using the best available substitution models for single-gene trees, by expanding our analyses to the diversity of uncultivated prokaryotes recently revealed by environmental genomics, and by using two complementary approaches to rooting that do not depend on outgroups. Consistent with some previous analyses, our rooted gene trees support extensive interdomain horizontal transfer of membrane phospholipid biosynthetic genes, primarily from Archaea to Bacteria. They also suggest that the capacity to make archaeal-type membrane phospholipids was already present in last universal common ancestor.

Polyploidy in halophilic archaea: regulation, evolutionary advantages, and gene conversion

All analyzed haloarachea are polyploid. In addition, haloarchaea contain more than one type of chromosome, and thus the gene dosage can be regulated independently on different replicons. Haloarchaea and several additional archaea have more than one replication origin on their major chromosome, in stark contrast with bacteria, which have a single replication origin. Two of these replication origins of Haloferax volcanii have been studied in detail and turned out to have very different properties. The chromosome copy number appears to be regulated in response to growth phases and environmental factors. Archaea typically contain about two Origin Recognition Complex (ORC) proteins, which are homologous to eukaryotic ORC proteins. However, haloarchaea are the only archaeal group that contains a multitude of ORC proteins. All 16 ORC protein paralogs from H. volcanii are involved in chromosome copy number regulation. Polyploidy has many evolutionary advantages for haloarchaea, e.g. a high resistance to desiccation, survival over geological times, and the relaxation of cell cycle-specific replication control. A further advantage is the ability to grow in the absence of external phosphate while using the many genome copies as internal phosphate storage polymers. Very efficient gene conversion operates in haloarchaea and results in the unification of genome copies. Taken together, haloarchaea are excellent models to study many aspects of genome biology in prokaryotes, exhibiting properties that have not been found in bacteria.

Carboxydotrophy potential of uncultivated Hydrothermarchaeota from the subseafloor crustal biosphere

The exploration of Earth's terrestrial subsurface biosphere has led to the discovery of several new archaeal lineages of evolutionary significance. Similarly, the deep subseafloor crustal biosphere also harbors many unique, uncultured archaeal taxa, including those belonging to Candidatus Hydrothermarchaeota, formerly known as Marine Benthic Group-E. Recently, Hydrothermarchaeota was identified as an abundant lineage of Juan de Fuca Ridge flank crustal fluids, suggesting its adaptation to this extreme environment. Through the investigation of single-cell and metagenome-assembled genomes, we provide insight into the lineage's evolutionary history and metabolic potential. Phylogenomic analysis reveals the Hydrothermarchaeota to be an early-branching archaeal phylum, branching between the superphylum DPANN, Euryarchaeota, and Asgard lineages. Hydrothermarchaeota genomes suggest a potential for dissimilative and assimilative carbon monoxide oxidation (carboxydotrophy), as well as sulfate and nitrate reduction. There is also a prevalence of chemotaxis and motility genes, indicating adaptive strategies for this nutrient-limited fluid-rock environment. These findings provide the first genomic interpretations of the Hydrothermarchaeota phylum and highlight the anoxic, hot, deep marine crustal biosphere as an important habitat for understanding the evolution of early life.

Phosphorylation of the Archaeal Holliday Junction Resolvase Hjc Inhibits Its Catalytic Activity and Facilitates DNA Repair in Sulfolobus islandicus REY15A

Protein phosphorylation is one of the main protein post-translational modifications and regulates DNA repair in eukaryotes. Archaeal genomes encode eukaryotic-like DNA repair proteins and protein kinases (ePKs), and several proteins involved in homologous recombination repair (HRR) including Hjc, a conserved Holliday junction (HJ) resolvase in Archaea, undergo phosphorylation, indicating that phosphorylation plays important roles in HRR. Herein, we performed phosphorylation analysis of Hjc by various ePKs from Sulfolobus islandicus. It was shown that SiRe_0171, SiRe_2030, and SiRe_2056, were able to phosphorylate Hjc in vitro. These ePKs phosphorylated Hjc at different Ser/Thr residues: SiRe_0171 on S34, SiRe_2030 on both S9 and T138, and SiRe_2056 on T138. The HJ cleavage activity of the phosphorylation-mimic mutants was analyzed and the results showed that the cleavage activity of S34E was completely lost and that of S9E had greatly reduced. S. islandicus strain expressing S34E in replacement of the wild type Hjc was resistant to higher doses of DNA damaging agents. Furthermore, SiRe_0171 deletion mutant exhibited higher sensitivity to DNA damaging agents, suggesting that Hjc phosphorylation by SiRe_0171 enhanced the DNA repair capability. Our results revealed that HJ resolvase is regulated by protein phosphorylation, reminiscent of the regulation of eukaryotic HJ resolvases GEN1 and Yen1.

Mosaic origin of the eukaryotic kinetochore

The emergence of eukaryotes from ancient prokaryotic lineages embodied a remarkable increase in cellular complexity. While prokaryotes operate simple systems to connect DNA to the segregation machinery during cell division, eukaryotes use a highly complex protein assembly known as the kinetochore. Although conceptually similar, prokaryotic segregation systems and the eukaryotic kinetochore are not homologous. Here we investigate the origins of the kinetochore before the last eukaryotic common ancestor (LECA) using phylogenetic trees, sensitive profile-versus-profile homology detection, and structural comparisons of its protein components. We show that LECA's kinetochore proteins share deep evolutionary histories with proteins involved in a few prokaryotic systems and a multitude of eukaryotic processes, including ubiquitination, transcription, and flagellar and vesicular transport systems. We find that gene duplications played a major role in shaping the kinetochore; more than half of LECA's kinetochore proteins have other kinetochore proteins as closest homologs. Some of these have no detectable homology to any other eukaryotic protein, suggesting that they arose as kinetochore-specific folds before LECA. We propose that the primordial kinetochore evolved from proteins involved in various (pre)eukaryotic systems as well as evolutionarily novel folds, after which a subset duplicated to give rise to the complex kinetochore of LECA.