Inference and reconstruction of the heimdallarchaeial ancestry of eukaryotes.

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

Expanding Archaeal Diversity and Phylogeny: Past, Present, and Future.

The discovery of the Archaea is a major scientific hallmark of the twentieth century. Since then, important features of their cell biology, physiology, ecology, and diversity have been revealed. Over the course of some 40 years, the diversity of known archaea has expanded from 2 to about 30 phyla comprising over 20,000 species. Most of this archaeal diversity has been revealed by environmental 16S rRNA gene amplicon sequencing surveys using a broad range of universal and targeted primers. Of the few primers that target a large fraction of known archaeal diversity, all display a bias against recently discovered lineages, which limits studies aiming to survey overall archaeal diversity. Induced by genomic exploration of archaeal diversity, and improved phylogenomics approaches, archaeal taxonomic classification has been frequently revised. Due to computational limitations and continued discovery of new lineages, a stable archaeal phylogeny is not yet within reach. Obtaining phylogenetic and taxonomic consensus of archaea should be a high priority for the archaeal research community.

Adaptation of the late ISC pathway in the anaerobic mitochondrial organelles of Giardia intestinalis.

Mitochondrial metabolism is entirely dependent on the biosynthesis of the [4Fe-4S] clusters, which are part of the subunits of the respiratory chain. The mitochondrial late ISC pathway mediates the formation of these clusters from simpler [2Fe-2S] molecules and transfers them to client proteins. Here, we characterized the late ISC pathway in one of the simplest mitochondria, mitosomes, of the anaerobic protist Giardia intestinalis that lost the respiratory chain and other hallmarks of mitochondria. In addition to IscA2, Nfu1 and Grx5 we identified a novel BolA1 homologue in G. intestinalis mitosomes. It specifically interacts with Grx5 and according to the high-affinity pulldown also with other core mitosomal components. Using CRISPR/Cas9 we were able to establish full bolA1 knock out, the first cell line lacking a mitosomal protein. Despite the ISC pathway being the only metabolic role of the mitosome no significant changes in the mitosome biology could be observed as neither the number of the mitosomes or their capability to form [2Fe-2S] clusters in vitro was affected. We failed to identify natural client proteins that would require the [2Fe-2S] or [4Fe-4S] cluster within the mitosomes, with the exception of [2Fe-2S] ferredoxin, which is itself part of the ISC pathway. The overall uptake of iron into the cellular proteins remained unchanged as also observed for the grx5 knock out cell line. The pull-downs of all late ISC components were used to build the interactome of the pathway showing specific position of IscA2 due to its interaction with the outer mitosomal membrane proteins. Finally, the comparative analysis across Metamonada species suggested that the adaptation of the late ISC pathway identified in G. intestinalis occurred early in the evolution of this supergroup of eukaryotes.

Deep mitochondrial origin outside the sampled alphaproteobacteria.

Mitochondria are ATP-generating organelles, the endosymbiotic origin of which was a key event in the evolution of eukaryotic cells 1 . Despite strong phylogenetic evidence that mitochondria had an alphaproteobacterial ancestry 2 , efforts to pinpoint their closest relatives among sampled alphaproteobacteria have generated conflicting results, complicating detailed inferences about the identity and nature of the mitochondrial ancestor. While most studies support the idea that mitochondria evolved from an ancestor related to Rickettsiales3-9, an order that includes several host-associated pathogenic and endosymbiotic lineages10,11, others have suggested that mitochondria evolved from a free-living group12-14. Here we re-evaluate the phylogenetic placement of mitochondria. We used genome-resolved binning of oceanic metagenome datasets and increased the genomic sampling of Alphaproteobacteria with twelve divergent clades, and one clade representing a sister group to all Alphaproteobacteria. Subsequent phylogenomic analyses that specifically address long branch attraction and compositional bias artefacts suggest that mitochondria did not evolve from Rickettsiales or any other currently recognized alphaproteobacterial lineage. Rather, our analyses indicate that mitochondria evolved from a proteobacterial lineage that branched off before the divergence of all sampled alphaproteobacteria. In light of this new result, previous hypotheses on the nature of the mitochondrial ancestor6,15,16 should be re-evaluated.

Asgard archaea illuminate the origin of eukaryotic cellular complexity.

The origin and cellular complexity of eukaryotes represent a major enigma in biology. Current data support scenarios in which an archaeal host cell and an alphaproteobacterial (mitochondrial) endosymbiont merged together, resulting in the first eukaryotic cell. The host cell is related to Lokiarchaeota, an archaeal phylum with many eukaryotic features. The emergence of the structural complexity that characterizes eukaryotic cells remains unclear. Here we describe the 'Asgard' superphylum, a group of uncultivated archaea that, as well as Lokiarchaeota, includes Thor-, Odin- and Heimdallarchaeota. Asgard archaea affiliate with eukaryotes in phylogenomic analyses, and their genomes are enriched for proteins formerly considered specific to eukaryotes. Notably, thorarchaeal genomes encode several homologues of eukaryotic membrane-trafficking machinery components, including Sec23/24 and TRAPP domains. Furthermore, we identify thorarchaeal proteins with similar features to eukaryotic coat proteins involved in vesicle biogenesis. Our results expand the known repertoire of 'eukaryote-specific' proteins in Archaea, indicating that the archaeal host cell already contained many key components that govern eukaryotic cellular complexity.

Convergent Evolution of Hydrogenosomes from Mitochondria by Gene Transfer and Loss.

Hydrogenosomes are H2-producing mitochondrial homologs found in some anaerobic microbial eukaryotes that provide a rare intracellular niche for H2-utilizing endosymbiotic archaea. Among ciliates, anaerobic and aerobic lineages are interspersed, demonstrating that the switch to an anaerobic lifestyle with hydrogenosomes has occurred repeatedly and independently. To investigate the molecular details of this transition, we generated genomic and transcriptomic data sets from anaerobic ciliates representing three distinct lineages. Our data demonstrate that hydrogenosomes have evolved from ancestral mitochondria in each case and reveal different degrees of independent mitochondrial genome and proteome reductive evolution, including the first example of complete mitochondrial genome loss in ciliates. Intriguingly, the FeFe-hydrogenase used for generating H2 has a unique domain structure among eukaryotes and appears to have been present, potentially through a single lateral gene transfer from an unknown donor, in the common aerobic ancestor of all three lineages. The early acquisition and retention of FeFe-hydrogenase helps to explain the facility whereby mitochondrial function can be so radically modified within this diverse and ecologically important group of microbial eukaryotes.

PhyloMagnet: fast and accurate screening of short-read meta-omics data using gene-centric phylogenetics.

MOTIVATION: Metagenomic and metatranscriptomic sequencing have become increasingly popular tools for producing massive amounts of short-read data, often used for the reconstruction of draft genomes or the detection of (active) genes in microbial communities. Unfortunately, sequence assemblies of such datasets generally remain a computationally challenging task. Frequently, researchers are only interested in a specific group of organisms or genes; yet, the assembly of multiple datasets only to identify candidate sequences for a specific question is sometimes prohibitively slow, forcing researchers to select a subset of available datasets to address their question. Here, we present PhyloMagnet, a workflow to screen meta-omics datasets for taxa and genes of interest using gene-centric assembly and phylogenetic placement of sequences. RESULTS: Using PhyloMagnet, we could identify up to 87% of the genera in an in vitro mock community with variable abundances, while the false positive predictions per single gene tree ranged from 0 to 23%. When applied to a group of metagenomes for which a set of metagenome assembled genomes (MAGs) have been published, we could detect the majority of the taxonomic labels that the MAGs had been annotated with. In a metatranscriptomic setting, the phylogenetic placement of assembled contigs corresponds to that of transcripts obtained from transcriptome assembly. AVAILABILITY AND IMPLEMENTATION: PhyloMagnet is built using Nextflow, available at github.com/maxemil/PhyloMagnet and is developed and tested on Linux. It is released under the open source GNU GPL licence and documentation is available at phylomagnet.readthedocs.io. Version 0.5 of PhyloMagnet was used for all benchmarking experiments. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Complex archaea that bridge the gap between prokaryotes and eukaryotes.

The origin of the eukaryotic cell remains one of the most contentious puzzles in modern biology. Recent studies have provided support for the emergence of the eukaryotic host cell from within the archaeal domain of life, but the identity and nature of the putative archaeal ancestor remain a subject of debate. Here we describe the discovery of 'Lokiarchaeota', a novel candidate archaeal phylum, which forms a monophyletic group with eukaryotes in phylogenomic analyses, and whose genomes encode an expanded repertoire of eukaryotic signature proteins that are suggestive of sophisticated membrane remodelling capabilities. Our results provide strong support for hypotheses in which the eukaryotic host evolved from a bona fide archaeon, and demonstrate that many components that underpin eukaryote-specific features were already present in that ancestor. This provided the host with a rich genomic 'starter-kit' to support the increase in the cellular and genomic complexity that is characteristic of eukaryotes.

An efficient single-cell transcriptomics workflow for microbial eukaryotes benchmarked on Giardia intestinalis cells.

BACKGROUND: Most diversity in the eukaryotic tree of life is represented by microbial eukaryotes, which is a polyphyletic group also referred to as protists. Among the protists, currently sequenced genomes and transcriptomes give a biased view of the actual diversity. This biased view is partly caused by the scientific community, which has prioritized certain microbes of biomedical and agricultural importance. Additionally, some protists remain difficult to maintain in cultures, which further influences what has been studied. It is now possible to bypass the time-consuming process of cultivation and directly analyze the gene content of single protist cells. Single-cell genomics was used in the first experiments where individual protists cells were genomically explored. Unfortunately, single-cell genomics for protists is often associated with low genome recovery and the assembly process can be complicated because of repetitive intergenic regions. Sequencing repetitive sequences can be avoided if single-cell transcriptomics is used, which only targets the part of the genome that is transcribed. RESULTS: In this study we test different modifications of Smart-seq2, a single-cell RNA sequencing protocol originally developed for mammalian cells, to establish a robust and more cost-efficient workflow for protists. The diplomonad Giardia intestinalis was used in all experiments and the available genome for this species allowed us to benchmark our results. We could observe increased transcript recovery when freeze-thaw cycles were added as an extra step to the Smart-seq2 protocol. Further we reduced the reaction volume and purified the amplified cDNA with alternative beads to test different cost-reducing changes of Smart-seq2. Neither improved the procedure, and reducing the volumes by half led to significantly fewer genes detected. We also added a 5' biotin modification to our primers and reduced the concentration of oligo-dT, to potentially reduce generation of artifacts. Except adding freeze-thaw cycles and reducing the volume, no other modifications lead to a significant change in gene detection. Therefore, we suggest adding freeze-thaw cycles to Smart-seq2 when working with protists and further consider our other modification described to improve cost and time-efficiency. CONCLUSIONS: The presented single-cell RNA sequencing workflow represents an efficient method to explore the diversity and cell biology of individual protist cells.

Microbial community assembly and evolution in subseafloor sediment.

Bacterial and archaeal communities inhabiting the subsurface seabed live under strong energy limitation and have growth rates that are orders of magnitude slower than laboratory-grown cultures. It is not understood how subsurface microbial communities are assembled and whether populations undergo adaptive evolution or accumulate mutations as a result of impaired DNA repair under such energy-limited conditions. Here we use amplicon sequencing to explore changes of microbial communities during burial and isolation from the surface to the >5,000-y-old subsurface of marine sediment and identify a small core set of mostly uncultured bacteria and archaea that is present throughout the sediment column. These persisting populations constitute a small fraction of the entire community at the surface but become predominant in the subsurface. We followed patterns of genome diversity with depth in four dominant lineages of the persisting populations by mapping metagenomic sequence reads onto single-cell genomes. Nucleotide sequence diversity was uniformly low and did not change with age and depth of the sediment. Likewise, there was no detectable change in mutation rates and efficacy of selection. Our results indicate that subsurface microbial communities predominantly assemble by selective survival of taxa able to persist under extreme energy limitation.

Integrative modeling of gene and genome evolution roots the archaeal tree of life.

A root for the archaeal tree is essential for reconstructing the metabolism and ecology of early cells and for testing hypotheses that propose that the eukaryotic nuclear lineage originated from within the Archaea; however, published studies based on outgroup rooting disagree regarding the position of the archaeal root. Here we constructed a consensus unrooted archaeal topology using protein concatenation and a multigene supertree method based on 3,242 single gene trees, and then rooted this tree using a recently developed model of genome evolution. This model uses evidence from gene duplications, horizontal transfers, and gene losses contained in 31,236 archaeal gene families to identify the most likely root for the tree. Our analyses support the monophyly of DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, Nanohaloarchaea), a recently discovered cosmopolitan and genetically diverse lineage, and, in contrast to previous work, place the tree root between DPANN and all other Archaea. The sister group to DPANN comprises the Euryarchaeota and the TACK Archaea, including Lokiarchaeum, which our analyses suggest are monophyletic sister lineages. Metabolic reconstructions on the rooted tree suggest that early Archaea were anaerobes that may have had the ability to reduce CO2 to acetate via the Wood-Ljungdahl pathway. In contrast to proposals suggesting that genome reduction has been the predominant mode of archaeal evolution, our analyses infer a relatively small-genomed archaeal ancestor that subsequently increased in complexity via gene duplication and horizontal gene transfer.

Oxygen induces the expression of invasion and stress response genes in the anaerobic salmon parasite Spironucleus salmonicida.

BACKGROUND: Spironucleus salmonicida is an anaerobic parasite that can cause systemic infections in Atlantic salmon. Unlike other diplomonad parasites, such as the human pathogen Giardia intestinalis, Spironucleus species can infiltrate the blood stream of their hosts eventually colonizing organs, skin and gills. How this presumed anaerobe can persist and invade oxygenated tissues, despite having a strictly anaerobic metabolism, remains elusive. RESULTS: To investigate how S. salmonicida response to oxygen stress, we performed RNAseq transcriptomic analyses of cells grown in the presence of oxygen or antioxidant-free medium. We found that over 20% of the transcriptome is differentially regulated in oxygen (1705 genes) and antioxidant-depleted (2280 genes) conditions. These differentially regulated transcripts encode proteins related to anaerobic metabolism, cysteine and Fe-S cluster biosynthesis, as well as a large number of proteins of unknown function. S. salmonicida does not encode genes involved in the classical elements of oxygen metabolism (e.g., catalases, superoxide dismutase, glutathione biosynthesis, oxidative phosphorylation). Instead, we found that genes encoding bacterial-like oxidoreductases were upregulated in response to oxygen stress. Phylogenetic analysis revealed some of these oxygen-responsive genes (e.g., nadh oxidase, rubrerythrin, superoxide reductase) are rare in eukaryotes and likely derived from lateral gene transfer (LGT) events into diplomonads from prokaryotes. Unexpectedly, we observed that many host evasion- and invasion-related genes were also upregulated under oxidative stress suggesting that oxygen might be an important signal for pathogenesis. CONCLUSION: While oxygen is toxic for related organisms, such as G. intestinalis, we find that oxygen is likely a gene induction signal for host invasion- and evasion-related pathways in S. salmonicida. These data provide the first molecular evidence for how S. salmonicida could tolerate oxic host environments and demonstrate how LGT can have a profound impact on the biology of anaerobic parasites.

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.

Vast diversity of prokaryotic virus genomes encoding double jelly-roll major capsid proteins uncovered by genomic and metagenomic sequence analysis.

BACKGROUND: Analysis of metagenomic sequences has become the principal approach for the study of the diversity of viruses. Many recent, extensive metagenomic studies on several classes of viruses have dramatically expanded the visible part of the virosphere, showing that previously undetected viruses, or those that have been considered rare, actually are important components of the global virome. RESULTS: We investigated the provenance of viruses related to tail-less bacteriophages of the family Tectiviridae by searching genomic and metagenomics sequence databases for distant homologs of the tectivirus-like Double Jelly-Roll major capsid proteins (DJR MCP). These searches resulted in the identification of numerous genomes of virus-like elements that are similar in size to tectiviruses (10-15 kilobases) and have diverse gene compositions. By comparison of the gene repertoires, the DJR MCP-encoding genomes were classified into 6 distinct groups that can be predicted to differ in reproduction strategies and host ranges. Only the DJR MCP gene that is present by design is shared by all these genomes, and most also encode a predicted DNA-packaging ATPase; the rest of the genes are present only in subgroups of this unexpectedly diverse collection of DJR MCP-encoding genomes. Only a minority encode a DNA polymerase which is a hallmark of the family Tectiviridae and the putative family "Autolykiviridae". Notably, one of the identified putative DJR MCP viruses encodes a homolog of Cas1 endonuclease, the integrase involved in CRISPR-Cas adaptation and integration of transposon-like elements called casposons. This is the first detected occurrence of Cas1 in a virus. Many of the identified elements are individual contigs flanked by inverted or direct repeats and appear to represent complete, extrachromosomal viral genomes, whereas others are flanked by bacterial genes and thus can be considered as proviruses. These contigs come from metagenomes of widely different environments, some dominated by archaea and others by bacteria, suggesting that collectively, the DJR MCP-encoding elements have a broad host range among prokaryotes. CONCLUSIONS: The findings reported here greatly expand the known host range of (putative) viruses of bacteria and archaea that encode a DJR MCP. They also demonstrate the extreme diversity of genome architectures in these viruses that encode no universal proteins other than the capsid protein that was used as the marker for their identification. From a supposedly minor group of bacterial and archaeal viruses, these viruses are emerging as a substantial component of the prokaryotic virome.

Archaeal Actin-Family Filament Systems.

Actin represents one of the most abundant and conserved eukaryotic proteins over time, and has an important role in many different cellular processes such as cell shape determination, motility, force generation, cytokinesis, amongst many others. Eukaryotic actin has been studied for decades and was for a long time considered a eukaryote-specific trait. However, in the early 2000s a bacterial actin homolog, MreB, was identified, characterized and found to have a cytoskeletal function and group within the superfamily of actin proteins. More recently, an actin cytoskeleton was also identified in archaea. The genome of the hyperthermophilic crenarchaeon Pyrobaculum calidifontis contains a five-gene cluster named Arcade encoding for an actin homolog, Crenactin, polymerizing into helical filaments spanning the whole length of the cell. Phylogenetic and structural studies place Crenactin closer to the eukaryotic actin than to the bacterial homologues. A significant difference, however, is that Crenactin can form single helical filaments in addition to filaments containing two intertwined proto filaments. The genome of the recently discovered Lokiarchaeota encodes several different actin homologues, termed Lokiactins, which are even more closely related to the eukaryotic actin than Crenactin. A primitive, dynamic actin-based cytoskeleton in archaea could have enabled the engulfment of the alphaproteobacterial progenitor of the mitochondria, a key-event in the evolution of eukaryotes.

Tracing the Archaeal Origins of Eukaryotic Membrane-Trafficking System Building Blocks.

In contrast to prokaryotes, eukaryotic cells are characterized by a complex set of internal membrane-bound compartments. A subset of these, and the protein machineries that move material between them, define the membrane-trafficking system (MTS), the emergence of which represents a landmark in eukaryotic evolution. Unlike mitochondria and plastids, MTS organelles have autogenous origins. Much of the MTS machinery is composed of building blocks, including small GTPase, coiled-coil, beta-propeller + alpha-solenoid, and longin domains. Despite the identification of prokaryotic proteins containing these domains, only few represent direct orthologues, leaving the origins and early evolution of the MTS poorly understood. Here, we present an in-depth analysis of MTS building block homologues in the composite genome of Lokiarchaeum, the recently discovered archaeal sister clade of eukaryotes, yielding several key insights. We identify two previously unreported Eukaryotic Signature Proteins; orthologues of the Gtr/Rag family GTPases, involved in target of rapamycin complex signaling, and of the RLC7 dynein component. We could not identify golgin or SNARE (coiled-coil) or beta-propeller + alpha-solenoid orthologues, nor typical MTS domain fusions, suggesting that these either were lost from Lokiarchaeum or emerged later in eukaryotic evolution. Furthermore, our phylogenetic analyses of lokiarchaeal GTPases support a split into Ras-like and Arf-like superfamilies, with different prokaryotic antecedents, before the advent of eukaryotes. While no GTPase activating proteins or exchange factors were identified, we show that Lokiarchaeum encodes numerous roadblock domain proteins and putative longin domain proteins, confirming the latter's origin from Archaea. Altogether, our study provides new insights into the emergence and early evolution of the eukaryotic membrane-trafficking system.

A Rickettsiales symbiont of amoebae with ancient features.

The Rickettsiae comprise intracellular bacterial symbionts and pathogens infecting diverse eukaryotes. Here, we provide a detailed characterization of 'Candidatus Jidaibacter acanthamoeba', a rickettsial symbiont of Acanthamoeba. The bacterium establishes the infection in its amoeba host within 2 h where it replicates within vacuoles. Higher bacterial loads and accelerated spread of infection at elevated temperatures were observed. The infection had a negative impact on host growth rate, although no increased levels of host cell lysis were seen. Phylogenomic analysis identified this bacterium as member of the Midichloriaceae. Its 2.4 Mb genome represents the largest among Rickettsiales and is characterized by a moderate degree of pseudogenization and a high coding density. We found an unusually large number of genes encoding proteins with eukaryotic-like domains such as ankyrins, leucine-rich repeats and tetratricopeptide repeats, which likely function in host interaction. There are a total of three divergent, independently acquired type IV secretion systems, and 35 flagellar genes representing the most complete set found in an obligate intracellular Alphaproteobacterium. The deeply branching phylogenetic position of 'Candidatus Jidaibacter acanthamoeba' together with its ancient features place it closely to the rickettsial ancestor and helps to better understand the transition from a free-living to an intracellular lifestyle.