More than the sum of its parts: uncovering emerging effects of microbial interactions in complex communities

Microbial communities are not only shaped by the diversity of microorganisms and their individual metabolic potential, but also by the vast amount of intra- and interspecies interactions that can occur pairwise interactions among microorganisms, we suggest that more attention should be drawn towards the effects on the entire microbiome that emerge from individual interactions between community members. The production of certain metabolites that can be tied to a specific microbe-microbe interaction might subsequently influence the physicochemical parameters of the habitat, stimulate a change in the trophic network of the community or create new micro-habitats through the formation of biofilms, similar to the production of antimicrobial substances which might negatively affect only one microorganism but cause a ripple effect on the abundance of other community members. Here, we argue that combining established as well as innovative laboratory and computational methods is needed to predict novel interactions and assess their secondary effects. Such efforts will enable future microbiome studies to expand our knowledge on the dynamics of complex microbial communities.

Type III-B CRISPR-Cas cascade of proteolytic cleavages

The generation of cyclic oligoadenylates and subsequent allosteric activation of proteins that carry sensory domains is a distinctive feature of type III CRISPR-Cas systems. In this work, we characterize a set of associated genes of a type III-B system from Haliangium ochraceum that contains two caspase-like proteases, SAVED-CHAT and PCaspase (prokaryotic caspase), co-opted from a cyclic oligonucleotide-based antiphage signaling system (CBASS). Cyclic tri-adenosine monophosphate (AMP)-induced oligomerization of SAVED-CHAT activates proteolytic activity of the CHAT domains, which specifically cleave and activate PCaspase. Subsequently, activated PCaspase cleaves a multitude of proteins, which results in a strong interference phenotype in vivo in Escherichia coli. Taken together, our findings reveal how a CRISPR-Cas-based detection of a target RNA triggers a cascade of caspase-associated proteolytic activities.

A genus in the bacterial phylum Aquificota appears to be endemic to Aotearoa-New Zealand

Allopatric speciation has been difficult to examine among microorganisms, with prior reports of endemism restricted to sub-genus level taxa. Previous microbial community analysis via 16S rRNA gene sequencing of 925 geothermal springs from the Taupō Volcanic Zone (TVZ), Aotearoa-New Zealand, revealed widespread distribution and abundance of a single bacterial genus across 686 of these ecosystems (pH 1.2-9.6 and 17.4-99.8 °C). Here, we present evidence to suggest that this genus, Venenivibrio (phylum Aquificota), is endemic to Aotearoa-New Zealand. A specific environmental niche that increases habitat isolation was identified, with maximal read abundance of Venenivibrio occurring at pH 4-6, 50-70 °C, and low oxidation-reduction potentials. This was further highlighted by genomic and culture-based analyses of the only characterised species for the genus, Venenivibrio stagnispumantis CP.B2T, which confirmed a chemolithoautotrophic metabolism dependent on hydrogen oxidation. While similarity between Venenivibrio populations illustrated that dispersal is not limited across the TVZ, extensive amplicon, metagenomic, and phylogenomic analyses of global microbial communities from DNA sequence databases indicates Venenivibrio is geographically restricted to the Aotearoa-New Zealand archipelago. We conclude that geographic isolation, complemented by physicochemical constraints, has resulted in the establishment of an endemic bacterial genus.

Candidatus Nemesobacterales is a sponge-specific clade of the candidate phylum Desulfobacterota adapted to a symbiotic lifestyle

Members of the candidate phylum Dadabacteria, recently reassigned to the phylum Candidatus Desulfobacterota, are cosmopolitan in the marine environment found both free-living and associated with hosts that are mainly marine sponges. Yet, these microorganisms are poorly characterized, with no cultured representatives and an ambiguous phylogenetic position in the tree of life. Here, we performed genome-centric metagenomics to elucidate their phylogenomic placement and predict the metabolism of the sponge-associated members of this lineage. Rank-based phylogenomics revealed several new species and a novel family (Candidatus Spongomicrobiaceae) within a sponge-specific order, named here Candidatus Nemesobacterales. Metabolic reconstruction suggests that Ca. Nemesobacterales are aerobic heterotrophs, capable of synthesizing most amino acids, vitamins and cofactors and degrading complex carbohydrates. We also report functional divergence between sponge- and seawater-associated metagenome-assembled genomes. Niche-specific adaptations to the sponge holobiont were evident from significantly enriched genes involved in defense mechanisms against foreign DNA and environmental stressors, host-symbiont interactions and secondary metabolite production. Fluorescence in situ hybridization gave a first glimpse of the morphology and lifestyle of a member of Ca. Desulfobacterota. Candidatus Nemesobacterales spp. were found both inside sponge cells centred around sponge nuclei and in the mesohyl of the sponge Geodia barretti. This study sheds light on the enigmatic group Ca. Nemesobacterales and their functional characteristics that reflect a symbiotic lifestyle.

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.

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.

The AEGEAN-169 clade of bacterioplankton is synonymous with SAR11 subclade V (HIMB59) and metabolically distinct

Bacterioplankton of the SAR11 clade are the most abundant marine microorganisms and consist of numerous subclades spanning order-level divergence (Pelagibacterales). The assignment of the earliest diverging subclade V (a.k.a. HIMB59) to the Pelagibacterales is highly controversial, with multiple recent phylogenetic studies placing them completely separate from SAR11. Other than through phylogenomics, subclade V has not received detailed examination due to limited genomes from this group. Here, we assessed the ecogenomic characteristics of subclade V to better understand the role of this group in comparison to the Pelagibacterales. We used a new isolate genome, recently released single-amplified genomes and metagenome-assembled genomes, and previously established SAR11 genomes to perform a comprehensive comparative genomics analysis. We paired this analysis with the recruitment of metagenomes spanning the open ocean, coastal, and brackish systems. Phylogenomics, average amino acid identity, and 16S rRNA gene phylogeny indicate that SAR11 subclade V is synonymous with the ubiquitous AEGEAN-169 clade and support the contention that this group represents a taxonomic family. AEGEAN-169 shared many bulk genome qualities with SAR11, such as streamlining and low GC content, but genomes were generally larger. AEGEAN-169 had overlapping distributions with SAR11 but was metabolically distinct from SAR11 in its potential to transport and utilize a broader range of sugars as well as in the transport of trace metals and thiamin. Thus, regardless of the ultimate phylogenetic placement of AEGEAN-169, these organisms have distinct metabolic capacities that likely allow them to differentiate their niche from canonical SAR11 taxa.IMPORTANCEOne goal of marine microbiologists is to uncover the roles various microorganisms are playing in biogeochemical cycles. Success in this endeavor relies on differentiating groups of microbes and circumscribing their relationships. An early-diverging group (subclade V) of the most abundant bacterioplankton, SAR11, has recently been proposed as a separate lineage that does not share a most recent common ancestor. But beyond phylogenetics, little has been done to evaluate how these organisms compare with SAR11. Our work leverages dozens of new genomes to demonstrate the similarities and differences between subclade V and SAR11. In our analysis, we also establish that subclade V is synonymous with a group of bacteria established from 16S rRNA gene sequences, AEGEAN-169. Subclade V/AEGEAN-169 has clear metabolic distinctions from SAR11 and their shared traits point to remarkable convergent evolution if they do not share a most recent common ancestor.

A Mitosome With Distinct Metabolism in the Uncultured Protist Parasite Paramikrocytos canceri (Rhizaria, Ascetosporea)

Ascetosporea are endoparasites of marine invertebrates that include economically important pathogens of aquaculture species. Owing to their often-minuscule cell sizes, strict intracellular lifestyle, lack of cultured representatives and minimal availability of molecular data, these unicellular parasites remain poorly studied. Here, we sequenced and assembled the genome and transcriptome of Paramikrocytos canceri, an endoparasite isolated from the European edible crab Cancer pagurus. Using bioinformatic predictions, we show that P. canceri likely possesses a mitochondrion-related organelle (MRO) with highly reduced metabolism, resembling the mitosomes of other parasites but with key differences. Like other mitosomes, this MRO is predicted to have reduced metabolic capacity and lack an organellar genome and function in iron-sulfur cluster (ISC) pathway-mediated Fe-S cluster biosynthesis. However, the MRO in P. canceri is uniquely predicted to produce ATP via a partial glycolytic pathway and synthesize phospholipids de novo through the CDP-DAG pathway. Heterologous gene expression confirmed that proteins from the ISC and CDP-DAG pathways retain mitochondrial targeting sequences that are recognized by yeast mitochondria. This represents a unique combination of metabolic pathways in an MRO, including the first reported case of a mitosome-like organelle able to synthesize phospholipids de novo. Some of these phospholipids, such as phosphatidylserine, are vital in other protist endoparasites that invade their host through apoptotic mimicry.

Gene gain facilitated endosymbiotic evolution of Chlamydiae

Chlamydiae is a bacterial phylum composed of obligate animal and protist endosymbionts. However, other members of the Planctomycetes-Verrucomicrobia-Chlamydiae superphylum are primarily free living. How Chlamydiae transitioned to an endosymbiotic lifestyle is still largely unresolved. Here we reconstructed Planctomycetes-Verrucomicrobia-Chlamydiae species relationships and modelled superphylum genome evolution. Gene content reconstruction from 11,996 gene families suggests a motile and facultatively anaerobic last common Chlamydiae ancestor that had already gained characteristic endosymbiont genes. Counter to expectations for genome streamlining in strict endosymbionts, we detected substantial gene gain within Chlamydiae. We found that divergence in energy metabolism and aerobiosis observed in extant lineages emerged later during chlamydial evolution. In particular, metabolic and aerobic genes characteristic of the more metabolically versatile protist-infecting chlamydiae were gained, such as respiratory chain complexes. Our results show that metabolic complexity can increase during endosymbiont evolution, adding an additional perspective for understanding symbiont evolutionary trajectories across the tree of life.

The evolutionary origin of host association in the Rickettsiales

The evolution of obligate host-association of bacterial symbionts and pathogens remains poorly understood. The Rickettsiales are an alphaproteobacterial order of obligate endosymbionts and parasites that infect a wide variety of eukaryotic hosts, including humans, livestock, insects and protists. Induced by their host-associated lifestyle, Rickettsiales genomes have undergone reductive evolution, leading to small, AT-rich genomes with limited metabolic capacities. Here we uncover eleven deep-branching alphaproteobacterial metagenome assembled genomes from aquatic environments, including data from the Tara Oceans initiative and other publicly available datasets, distributed over three previously undescribed Rickettsiales-related clades. Phylogenomic analyses reveal that two of these clades, Mitibacteraceae and Athabascaceae, branch sister to all previously sampled Rickettsiales. The third clade, Gamibacteraceae, branch sister to the recently identified ectosymbiotic ‘Candidatus Deianiraea vastatrix’. Comparative analyses indicate that the gene complement of Mitibacteraceae and Athabascaceae is reminiscent of that of free-living and biofilm-associated bacteria. Ancestral genome content reconstruction across the Rickettsiales species tree further suggests that the evolution of host association in Rickettsiales was a gradual process that may have involved the repurposing of a type IV secretion system.

Phylogenomic analyses reveal novel environmental clades of Rickettsiales providing insights into their evolution from free-living to host-associated lifestyle.

Asgard archaea shed light on the evolutionary origins of the eukaryotic ubiquitin-ESCRT machinery

The ESCRT machinery, comprising of multiple proteins and subcomplexes, is crucial for membrane remodelling in eukaryotic cells, in processes that include ubiquitin-mediated multivesicular body formation, membrane repair, cytokinetic abscission, and virus exit from host cells. This ESCRT system appears to have simpler, ancient origins, since many archaeal species possess homologues of ESCRT-III and Vps4, the components that execute the final membrane scission reaction, where they have been shown to play roles in cytokinesis, extracellular vesicle formation and viral egress. Remarkably, metagenome assemblies of Asgard archaea, the closest known living relatives of eukaryotes, were recently shown to encode homologues of the entire cascade involved in ubiquitin-mediated membrane remodelling, including ubiquitin itself, components of the ESCRT-I and ESCRT-II subcomplexes, and ESCRT-III and Vps4. Here, we explore the phylogeny, structure, and biochemistry of Asgard homologues of the ESCRT machinery and the associated ubiquitylation system. We provide evidence for the ESCRT-I and ESCRT-II subcomplexes being involved in ubiquitin-directed recruitment of ESCRT-III, as it is in eukaryotes. Taken together, our analyses suggest a pre-eukaryotic origin for the ubiquitin-coupled ESCRT system and a likely path of ESCRT evolution via a series of gene duplication and diversification events.

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 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. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Innovations to culturing the uncultured microbial majority

Despite the surge of microbial genome data, experimental testing is important to confirm inferences about the cell biology, ecological roles and evolution of microorganisms. As the majority of archaeal and bacterial diversity remains uncultured and poorly characterized, culturing is a priority. The growing interest in and need for efficient cultivation strategies has led to many rapid methodological and technological advances. In this Review, we discuss common barriers that can hamper the isolation and culturing of novel microorganisms and review emerging, innovative methods for targeted or high-throughput cultivation. We also highlight recent examples of successful cultivation of novel archaea and bacteria, and suggest key microorganisms for future cultivation attempts.

Bubble biofilm: Bacterial colonization of air-air interface

Microbial mats or biofilms are known to colonize a wide range of substrates in aquatic environments. These dense benthic communities efficiently recycle nutrients and often exhibit high tolerance to environmental stressors, characteristics that enable them to inhabit harsh ecological niches. In some special cases, floating biofilms form at the air-water interface residing on top of a hydrophobic microlayer. Here, we describe biofilms that reside at the air-air interface by forming gas bubbles (bubble biofilms) in the former Ytterby mine, Sweden. The bubbles are built by micrometer thick membrane-like biofilm that holds enough water to sustain microbial activity. Molecular identification shows that the biofilm communities are dominated by the neuston bacterium Nevskia. Gas bubbles contain mostly air with a slightly elevated concentration of carbon dioxide. Biofilm formation and development was monitored in situ using a time-lapse camera over one year, taking one image every second hour. The bubbles were stable over long periods of time (weeks, even months) and gas build-up occurred in pulses as if the bedrock suddenly exhaled. The result was however not a passive inflation of a dying biofilm becoming more fragile with time (as a result of overstretching of the organic material). To the contrary, microbial growth lead to a more robust, hydrophobic bubble biofilm that kept the bubbles inflated for extended periods (several weeks, and in some cases even months).

Author Correction: Roadmap for naming uncultivated Archaea and Bacteria

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Hikarchaeia demonstrate an intermediate stage in the methanogen-to-halophile transition

Halobacteria (henceforth: Haloarchaea) are predominantly aerobic halophiles that are thought to have evolved from anaerobic methanogens. This remarkable transformation most likely involved an extensive influx of bacterial genes. Whether it entailed a single massive transfer event or a gradual stream of transfers remains a matter of debate. To address this, genomes that descend from methanogen-to-halophile intermediates are necessary. Here, we present five such near-complete genomes of Marine Group IV archaea (Hikarchaeia), the closest known relatives of Haloarchaea. Their inclusion in gene tree-aware ancestral reconstructions reveals an intermediate stage that had already lost a large number of genes, including nearly all of those involved in methanogenesis and the Wood-Ljungdahl pathway. In contrast, the last Haloarchaea common ancestor gained a large number of genes and expanded its aerobic respiration and salt/UV resistance gene repertoire. Our results suggest that complex and gradual patterns of gain and loss shaped the methanogen-to-halophile transition.

Microbiomes in a manganese oxide producing ecosystem in the Ytterby mine, Sweden: impact on metal mobility

Microbe-mediated precipitation of Mn-oxides enriched in rare earth elements (REE) and other trace elements was discovered in tunnels leading to the main shaft of the Ytterby mine, Sweden. Defining the spatial distribution of microorganisms and elements in this ecosystem provide a better understanding of specific niches and parameters driving the emergence of these communities and associated mineral precipitates. Along with elemental analyses, high-throughput sequencing of the following four subsystems were conducted: (i) water seeping from a rock fracture into the tunnel, (ii) Mn-oxides and associated biofilm; referred to as the Ytterby Black Substance (YBS) biofilm (iii) biofilm forming bubbles on the Mn-oxides; referred to as the bubble biofilm and (iv) fracture water that has passed through the biofilms. Each subsystem hosts a specific collection of microorganisms. Differentially abundant bacteria in the YBS biofilm were identified within the Rhizobiales (e.g. Pedomicrobium), PLTA13 Gammaproteobacteria, Pirellulaceae, Hyphomonadaceae, Blastocatellia and Nitrospira. These taxa, likely driving the Mn-oxide production, were not detected in the fracture water. This biofilm binds Mn, REE and other trace elements in an efficient, dynamic process, as indicated by substantial depletion of these metals from the fracture water as it passes through the Mn deposit zone. Microbe-mediated oxidation of Mn(II) and formation of Mn(III/IV)-oxides can thus have considerable local environmental impact by removing metals from aquatic environments.

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.

Chlamydial contribution to anaerobic metabolism during eukaryotic evolution

The origin of eukaryotes is a major open question in evolutionary biology. Multiple hypotheses posit that eukaryotes likely evolved from a syntrophic relationship between an archaeon and an alphaproteobacterium based on H₂ exchange. However, there are no strong indications that modern eukaryotic H₂ metabolism originated from archaea or alphaproteobacteria. Here, we present evidence for the origin of H₂ metabolism genes in eukaryotes from an ancestor of the Anoxychlamydiales-a group of anaerobic chlamydiae, newly described here, from marine sediments. Among Chlamydiae, these bacteria uniquely encode genes for H₂ metabolism and other anaerobiosis-associated pathways. Phylogenetic analyses of several components of H₂ metabolism reveal that Anoxychlamydiales homologs are the closest relatives to eukaryotic sequences. We propose that an ancestor of the Anoxychlamydiales contributed these key genes during the evolution of eukaryotes, supporting a mosaic evolutionary origin of eukaryotic metabolism.

A novel alphaproteobacterium with a small genome identified from the digestive gland of multiple species of abalone

We identified an alphaproteobacterium in the digestive gland of the abalone species Haliotis discus hannai. This phylotype dominated our 16S rRNA clone libraries from the digestive gland of H. discus hannai. Diversity surveys revealed that this phylotype was associated with H. discus hannai and also in another host species, H. gigantea. Whole genome phylogenies placed this bacterium as a new member affiliated with the family Rhodospirillaceae in Alphaproteobacteria. Gene annotation revealed a nearly complete glycolysis pathway but no TCA cycle, but the presence of anaerobic ribonucleoside-triphosphate reductase and oxygen-insensitive NAD(P)H-dependent nitroreductase, which show the genomic potential for anaerobic metabolism. A large cluster of genes encoding ankyrin repeat proteins (ANK) of eukaryotic-like repeat domains and a large gene set for the flagellar system were also detected. Alginate-binding periplasmic proteins and key genes responsible for alginate assimilation were found in the genome, which could potentially contribute to the breakdown of the host's alginate-rich macroalgal diet. These results raise the possibility that this novel alphaproteobacterium is a widespread member of the abalone microbiome that may use polysaccharides derived from its host's macroalgal diet.

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.

Roadmap for naming uncultivated Archaea and Bacteria

The assembly of single-amplified genomes (SAGs) and metagenome-assembled genomes (MAGs) has led to a surge in genome-based discoveries of members affiliated with Archaea and Bacteria, bringing with it a need to develop guidelines for nomenclature of uncultivated microorganisms. The International Code of Nomenclature of Prokaryotes (ICNP) only recognizes cultures as 'type material', thereby preventing the naming of uncultivated organisms. In this Consensus Statement, we propose two potential paths to solve this nomenclatural conundrum. One option is the adoption of previously proposed modifications to the ICNP to recognize DNA sequences as acceptable type material; the other option creates a nomenclatural code for uncultivated Archaea and Bacteria that could eventually be merged with the ICNP in the future. Regardless of the path taken, we believe that action is needed now within the scientific community to develop consistent rules for nomenclature of uncultivated taxa in order to provide clarity and stability, and to effectively communicate microbial diversity.

Marine Sediments Illuminate Chlamydiae Diversity and Evolution

The bacterial phylum Chlamydiae is so far composed of obligate symbionts of eukaryotic hosts. Well known for Chlamydiaceae, pathogens of humans and other animals, Chlamydiae also include so-called environmental lineages that primarily infect microbial eukaryotes. Environmental surveys indicate that Chlamydiae are found in a wider range of environments than anticipated previously. However, the vast majority of this chlamydial diversity has been underexplored, biasing our current understanding of their biology, ecological importance, and evolution. Here, we report that previously undetected and active chlamydial lineages dominate microbial communities in deep anoxic marine sediments taken from the Arctic Mid-Ocean Ridge. Reaching relative abundances of up to 43% of the bacterial community, and a maximum diversity of 163 different species-level taxonomic units, these Chlamydiae represent important community members. Using genome-resolved metagenomics, we reconstructed 24 draft chlamydial genomes, expanding by over a third the known genomic diversity in this phylum. Phylogenomic analyses revealed several novel clades across the phylum, including a previously unknown sister lineage of the Chlamydiaceae, providing new insights into the origin of pathogenicity in this family. We were unable to identify putative eukaryotic hosts for these marine sediment chlamydiae, despite identifying genomic features that may be indicative of host-association. The high abundance and genomic diversity of Chlamydiae in these anoxic marine sediments indicate that some members could play an important, and thus far overlooked, ecological role in such environments and may indicate alternate lifestyle strategies.

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.

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.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0634-8) contains supplementary material, which is available to authorized users.

Genomes from uncultivated prokaryotes: a comparison of metagenome-assembled and single-amplified genomes

BACKGROUND

Prokaryotes dominate the biosphere and regulate biogeochemical processes essential to all life. Yet, our knowledge about their biology is for the most part limited to the minority that has been successfully cultured. Molecular techniques now allow for obtaining genome sequences of uncultivated prokaryotic taxa, facilitating in-depth analyses that may ultimately improve our understanding of these key organisms.

RESULTS

We compared results from two culture-independent strategies for recovering bacterial genomes: single-amplified genomes and metagenome-assembled genomes. Single-amplified genomes were obtained from samples collected at an offshore station in the Baltic Sea Proper and compared to previously obtained metagenome-assembled genomes from a time series at the same station. Among 16 single-amplified genomes analyzed, seven were found to match metagenome-assembled genomes, affiliated with a diverse set of taxa. Notably, genome pairs between the two approaches were nearly identical (average 99.51% sequence identity; range 98.77-99.84%) across overlapping regions (30-80% of each genome). Within matching pairs, the single-amplified genomes were consistently smaller and less complete, whereas the genetic functional profiles were maintained. For the metagenome-assembled genomes, only on average 3.6% of the bases were estimated to be missing from the genomes due to wrongly binned contigs.

CONCLUSIONS

The strong agreement between the single-amplified and metagenome-assembled genomes emphasizes that both methods generate accurate genome information from uncultivated bacteria. Importantly, this implies that the research questions and the available resources are allowed to determine the selection of genomics approach for microbiome studies.

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

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

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.

RNA Sequencing of Stentor Cell Fragments Reveals Transcriptional Changes during Cellular Regeneration

While ciliates of the genus Stentor are known for their ability to regenerate when their cells are damaged or even fragmented, the physical and molecular mechanisms underlying this process are poorly understood. To identify genes involved in the regenerative capability of Stentor cells, RNA sequencing of individual Stentor polymorphus cell fragments was performed. After splitting a cell over the anterior-posterior axis, the posterior fragment has to regenerate the oral apparatus, while the anterior part needs to regenerate the hold fast. Altogether, differential expression analysis of both posterior and anterior S. polymorphus cell fragments for four different post-split time points revealed over 10,000 upregulated genes throughout the regeneration process. Among these, genes involved in cell signaling, microtubule-based movement, and cell cycle regulation seemed to be particularly important during cellular regeneration. We identified roughly nine times as many upregulated genes in regenerating S. polymorphus posterior fragments as compared to anterior fragments, indicating that regeneration of the anterior oral apparatus is a complex process that involves many genes. Our analyses identified several expanded groups of genes, such as dual-specific tyrosine-(Y)-phosphorylation-regulated kinases and MORN domain-containing proteins that seemingly act as key regulators of cellular regeneration. In agreement with earlier morphological and cell biological studies [1, 2], our differential expression analyses indicate that cellular regeneration and vegetative division share many similarities.

Functional reconstruction of a eukaryotic-like E1/E2/(RING) E3 ubiquitylation cascade from an uncultured archaeon

The covalent modification of protein substrates by ubiquitin regulates a diverse range of critical biological functions. Although it has been established that ubiquitin-like modifiers evolved from prokaryotic sulphur transfer proteins it is less clear how complex eukaryotic ubiquitylation system arose and diversified from these prokaryotic antecedents. The discovery of ubiquitin, E1-like, E2-like and small-RING finger (srfp) protein components in the Aigarchaeota and the Asgard archaea superphyla has provided a substantive step toward addressing this evolutionary question. Encoded in operons, these components are likely representative of the progenitor apparatus that founded the modern eukaryotic ubiquitin modification systems. Here we report that these proteins from the archaeon Candidatus ‘Caldiarchaeum subterraneum’ operate together as a bona fide ubiquitin modification system, mediating a sequential ubiquitylation cascade reminiscent of the eukaryotic process. Our observations support the hypothesis that complex eukaryotic ubiquitylation signalling pathways have developed from compact systems originally inherited from an archaeal ancestor.

In eukaryotic cells, the ubiquitylation system regulates several cellular processes central to protein homoeostasis. Here the authors demonstrate the existence of an eukaryotic-like ubiquitylation cascade requiring E1, E2 and E3-like enzymes in the archaeon C. subterraneum, shedding light on the evolution of the ubiquitin-proteasome system.

Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea

Standards for sequencing the microbial 'uncultivated majority', namely bacterial and archaeal single-cell genome sequences, and genome sequences from metagenomic datasets, are proposed.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.

The novel bacterial phylum Calditrichaeota is diverse, widespread and abundant in marine sediments and has the capacity to degrade detrital proteins

Calditrichaeota is a recently recognized bacterial phylum with three cultured representatives, isolated from hydrothermal vents. Here we expand the phylogeny and ecology of this novel phylum with metagenome-derived and single-cell genomes from six uncultivated bacteria previously not recognized as members of Calditrichaeota. Using 16S rRNA gene sequences from these genomes, we then identified 322 16S rRNA gene sequences from cultivation-independent studies that can now be classified as Calditrichaeota for the first time. This dataset was used to re-analyse a collection of 16S rRNA gene amplicon datasets from marine sediments showing that the Calditrichaeota are globally distributed in the seabed at high abundance, making up to 6.7% of the total bacterial community. This wide distribution and high abundance of Calditrichaeota in cold marine sediment has gone unrecognized until now. All Calditrichaeota genomes show indications of a chemoorganoheterotrophic metabolism with the potential to degrade detrital proteins through the use of extracellular peptidases. Most of the genomes contain genes encoding proteins that confer O₂ tolerance, consistent with the relatively high abundance of Calditrichaeota in surficial bioturbated part of the seabed and, together with the genes encoding extracellular peptidases, suggestive of a general ecophysiological niche for this newly recognized phylum in marine sediment.

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.

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.

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.

Picodroplet partitioned whole genome amplification of low biomass samples preserves genomic diversity for metagenomic analysis

BACKGROUND

Whole genome amplification (WGA) is a challenging, key step in metagenomic studies of samples containing minute amounts of DNA, such as samples from low biomass environments. It is well known that multiple displacement amplification (MDA), the most commonly used WGA method for microbial samples, skews the genomic representation in the sample. We have combined MDA with droplet microfluidics to perform the reaction in a homogeneous emulsion. Each droplet in this emulsion can be considered an individual reaction chamber, allowing partitioning of the MDA reaction into millions of parallel reactions with only one or very few template molecules per droplet.

RESULTS

As a proof-of-concept, we amplified genomic DNA from a synthetic metagenome by MDA either in one bulk reaction or in emulsion and found that after sequencing, the species distribution was better preserved and the coverage depth was more evenly distributed across the genomes when the MDA reaction had been performed in emulsion.

CONCLUSIONS

Partitioning MDA reactions into millions of reactions by droplet microfluidics is a straightforward way to improve the uniformity of MDA reactions for amplifying complex samples with limited amounts of DNA.

Exploring microbial dark matter to resolve the deep archaeal ancestry of eukaryotes

The origin of eukaryotes represents an enigmatic puzzle, which is still lacking a number of essential pieces. Whereas it is currently accepted that the process of eukaryogenesis involved an interplay between a host cell and an alphaproteobacterial endosymbiont, we currently lack detailed information regarding the identity and nature of these players. A number of studies have provided increasing support for the emergence of the eukaryotic host cell from within the archaeal domain of life, displaying a specific affiliation with the archaeal TACK superphylum. Recent studies have shown that genomic exploration of yet-uncultivated archaea, the so-called archaeal ‘dark matter’, is able to provide unprecedented insights into the process of eukaryogenesis. Here, we provide an overview of state-of-the-art cultivation-independent approaches, and demonstrate how these methods were used to obtain draft genome sequences of several novel members of the TACK superphylum, including Lokiarchaeum, two representatives of the Miscellaneous Crenarchaeotal Group (Bathyarchaeota), and a Korarchaeum-related lineage. The maturation of cultivation-independent genomics approaches, as well as future developments in next-generation sequencing technologies, will revolutionize our current view of microbial evolution and diversity, and provide profound new insights into the early evolution of life, including the enigmatic origin of the eukaryotic cell.

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.

Genomic inference of the metabolism of cosmopolitan subsurface Archaea, Hadesarchaea

The subsurface biosphere is largely unexplored and contains a broad diversity of uncultured microbes(1). Despite being one of the few prokaryotic lineages that is cosmopolitan in both the terrestrial and marine subsurface(2-4), the physiological and ecological roles of SAGMEG (South-African Gold Mine Miscellaneous Euryarchaeal Group) Archaea are unknown. Here, we report the metabolic capabilities of this enigmatic group as inferred from genomic reconstructions. Four high-quality (63-90% complete) genomes were obtained from White Oak River estuary and Yellowstone National Park hot spring sediment metagenomes. Phylogenomic analyses place SAGMEG Archaea as a deeply rooting sister clade of the Thermococci, leading us to propose the name Hadesarchaea for this new Archaeal class. With an estimated genome size of around 1.5 Mbp, the genomes of Hadesarchaea are distinctly streamlined, yet metabolically versatile. They share several physiological mechanisms with strict anaerobic Euryarchaeota. Several metabolic characteristics make them successful in the subsurface, including genes involved in CO and H2 oxidation (or H2 production), with potential coupling to nitrite reduction to ammonia (DNRA). This first glimpse into the metabolic capabilities of these cosmopolitan Archaea suggests they are mediating key geochemical processes and are specialized for survival in the subsurface biosphere.