Towards functional characterization of archaeal genomic dark matter

Regulatory sequence-based discovery of anti-defense genes in archaeal viruses

In silico identification of viral anti-CRISPR proteins (Acrs) has relied largely on the guilt-by-association method using known Acrs or anti-CRISPR associated proteins (Acas) as the bait. However, the low number and limited spread of the characterized archaeal Acrs and Aca hinders our ability to identify Acrs using guilt-by-association. Here, based on the observation that the few characterized archaeal Acrs and Aca are transcribed immediately post viral infection, we hypothesize that these genes, and many other unidentified anti-defense genes (ADG), are under the control of conserved regulatory sequences including a strong promoter, which can be used to predict anti-defense genes in archaeal viruses. Using this consensus sequence based method, we identify 354 potential ADGs in 57 archaeal viruses and 6 metagenome-assembled genomes. Experimental validation identified a CRISPR subtype I-A inhibitor and the first virally encoded inhibitor of an archaeal toxin-antitoxin based immune system. We also identify regulatory proteins potentially akin to Acas that can facilitate further identification of ADGs combined with the guilt-by-association approach. These results demonstrate the potential of regulatory sequence analysis for extensive identification of ADGs in viruses of archaea and bacteria.

Diversity, origin, and evolution of the ESCRT systems

Endosomal sorting complexes required for transport (ESCRT) play key roles in protein sorting between membrane-bounded compartments of eukaryotic cells. Homologs of many ESCRT components are identifiable in various groups of archaea, especially in Asgardarchaeota, the archaeal phylum that is currently considered to include the closest relatives of eukaryotes, but not in bacteria. We performed a comprehensive search for ESCRT protein homologs in archaea and reconstructed ESCRT evolution using the phylogenetic tree of Vps4 ATPase (ESCRT IV) as a scaffold and using sensitive protein sequence analysis and comparison of structural models to identify previously unknown ESCRT proteins. Several distinct groups of ESCRT systems in archaea outside of Asgard were identified, including proteins structurally similar to ESCRT-I and ESCRT-II, and several other domains involved in protein sorting in eukaryotes, suggesting an early origin of these components. Additionally, distant homologs of CdvA proteins were identified in Thermoproteales which are likely components of the uncharacterized cell division system in these archaea. We propose an evolutionary scenario for the origin of eukaryotic and Asgard ESCRT complexes from ancestral building blocks, namely, the Vps4 ATPase, ESCRT-III components, wH (winged helix-turn-helix fold) and possibly also coiled-coil, and Vps28-like domains. The last archaeal common ancestor likely encompassed a complex ESCRT system that was involved in protein sorting. Subsequent evolution involved either simplification, as in the TACK superphylum, where ESCRT was co-opted for cell division, or complexification as in Asgardarchaeota. In Asgardarchaeota, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was already established.IMPORTANCEAll eukaryotic cells possess complex intracellular membrane organization. Endosomal sorting complexes required for transport (ESCRT) play a central role in membrane remodeling which is essential for cellular functionality in eukaryotes. Recently, it has been shown that Asgard archaea, the archaeal phylum that includes the closest known relatives of eukaryotes, encode homologs of many components of the ESCRT systems. We employed protein sequence and structure comparisons to reconstruct the evolution of ESCRT systems in archaea and identified several previously unknown homologs of ESCRT subunits, some of which can be predicted to participate in cell division. The results of this reconstruction indicate that the last archaeal common ancestor already encoded a complex ESCRT system that was involved in protein sorting. In Asgard archaea, ESCRT systems evolved toward greater complexity, and in particular, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was established.

Widespread photosynthesis reaction centre barrel proteins are necessary for haloarchaeal cell division

Cell division is fundamental to all cellular life. Most archaea depend on either the prokaryotic tubulin homologue FtsZ or the endosomal sorting complex required for transport for division but neither system has been robustly characterized. Here, we show that three of the four photosynthesis reaction centre barrel domain proteins of Haloferax volcanii (renamed cell division proteins B1/2/3 (CdpB1/2/3)) play important roles in cell division. CdpB1 interacts directly with the FtsZ membrane anchor SepF and is essential for cell division, whereas deletion of cdpB2 and cdpB3 causes a major and a minor division defect, respectively. Orthologues of CdpB proteins are also involved in cell division in other haloarchaea, indicating a conserved function of these proteins. Phylogenetic analysis shows that photosynthetic reaction centre barrel proteins are widely distributed among archaea and appear to be central to cell division in most if not all archaea.

Computational analysis of genes with lethal knockout phenotype and prediction of essential genes in archaea

The identification of microbial genes essential for survival as those with lethal knockout phenotype (LKP) is a common strategy for functional interrogation of genomes. However, interpretation of the LKP is complicated because a substantial fraction of the genes with this phenotype remains poorly functionally characterized. Furthermore, many genes can exhibit LKP not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes (conditionally essential genes). We analyzed the sets of LKP genes for two archaea, Methanococcus maripaludis and Sulfolobus islandicus, using a variety of computational approaches aiming to differentiate between essential and conditionally essential genes and to predict at least a general function for as many of the proteins encoded by these genes as possible. This analysis allowed us to predict the functions of several LKP genes including previously uncharacterized subunit of the GINS protein complex with an essential function in genome replication and of the KEOPS complex that is responsible for an essential tRNA modification as well as GRP protease implicated in protein quality control. Additionally, several novel antitoxins (conditionally essential genes) were predicted, and this prediction was experimentally validated by showing that the deletion of these genes together with the adjacent genes apparently encoding the cognate toxins caused no growth defect. We applied principal component analysis based on sequence and comparative genomic features showing that this approach can separate essential genes from conditionally essential ones and used it to predict essential genes in other archaeal genomes.IMPORTANCEOnly a relatively small fraction of the genes in any bacterium or archaeon is essential for survival as demonstrated by the lethal effect of their disruption. The identification of essential genes and their functions is crucial for understanding fundamental cell biology. However, many of the genes with a lethal knockout phenotype remain poorly functionally characterized, and furthermore, many genes can exhibit this phenotype not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes. We applied state-of-the-art computational methods to predict the functions of a number of uncharacterized genes with the lethal knockout phenotype in two archaeal species and developed a computational approach to predict genes involved in essential functions. These findings advance the current understanding of key functionalities of archaeal cells.

Diversity, Origin and Evolution of the ESCRT Systems

Endosomal Sorting Complexes Required for Transport (ESCRT) play key roles in protein sorting between membrane-bounded compartments of eukaryotic cells. Homologs of many ESCRT components are identifiable in various groups of archaea, especially in Asgardarchaeota, the archaeal phylum that is currently considered to include the closest relatives of eukaryotes, but not in bacteria. We performed a comprehensive search for ESCRT protein homologs in archaea and reconstructed ESCRT evolution using the phylogenetic tree of Vps4 ATPase (ESCRT IV) as a scaffold, using sensitive protein sequence analysis and comparison of structural models to identify previously unknown ESCRT proteins. Several distinct groups of ESCRT systems in archaea outside of Asgard were identified, including proteins structurally similar to ESCRT-I and ESCRT-II, and several other domains involved in protein sorting in eukaryotes, suggesting an early origin of these components. Additionally, distant homologs of CdvA proteins were identified in Thermoproteales which are likely components of the uncharacterized cell division system in these archaea. We propose an evolutionary scenario for the origin of eukaryotic and Asgard ESCRT complexes from ancestral building blocks, namely, the Vps4 ATPase, ESCRT-III components, wH (winged helix-turn-helix fold) and possibly also coiled-coil, and Vps28-like domains. The Last Archaeal Common Ancestor likely encompassed a complex ESCRT system that was involved in protein sorting. Subsequent evolution involved either simplification, as in the TACK superphylum, where ESCRT was co-opted for cell division, or complexification as in Asgardarchaeota. In Asgardarchaeota, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was already established.

Functional and evolutionary significance of unknown genes from uncultivated taxa

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

A unified compendium of prokaryotic and viral genomes from over 300 anaerobic digestion microbiomes

BACKGROUND

The anaerobic digestion process degrades organic matter into simpler compounds and occurs in strictly anaerobic and microaerophilic environments. The process is carried out by a diverse community of microorganisms where each species has a unique role and it has relevant biotechnological applications since it is used for biogas production. Some aspects of the microbiome, including its interaction with phages, remains still unclear: a better comprehension of the community composition and role of each species is crucial for a cured understanding of the carbon cycle in anaerobic systems and improving biogas production.

RESULTS

The primary objective of this study was to expand our understanding on the anaerobic digestion microbiome by jointly analyzing its prokaryotic and viral components. By integrating 192 additional datasets into a previous metagenomic database, the binning process generated 11,831 metagenome-assembled genomes from 314 metagenome samples published between 2014 and 2022, belonging to 4,568 non-redundant species based on ANI calculation and quality verification. CRISPR analysis on these genomes identified 76 archaeal genomes with active phage interactions. Moreover, single-nucleotide variants further pointed to archaea as the most critical members of the community. Among the MAGs, two methanogenic archaea, Methanothrix sp. 43zhSC_152 and Methanoculleus sp. 52maCN_3230, had the highest number of SNVs, with the latter having almost double the density of most other MAGs.

CONCLUSIONS

This study offers a more comprehensive understanding of microbial community structures that thrive at different temperatures. The findings revealed that the fraction of archaeal species characterized at the genome level and reported in public databases is higher than that of bacteria, although still quite limited. The identification of shared spacers between phages and microbes implies a history of phage-bacterial interactions, and specifically lysogenic infections. A significant number of SNVs were identified, primarily comprising synonymous and nonsynonymous variants. Together, the findings indicate that methanogenic archaea are subject to intense selective pressure and suggest that genomic variants play a critical role in the anaerobic digestion process. Overall, this study provides a more balanced and diverse representation of the anaerobic digestion microbiota in terms of geographic location, temperature range and feedstock utilization.

Functional prediction of proteins from the human gut archaeome

The human gastrointestinal tract contains diverse microbial communities, including archaea. Among them, Methanobrevibacter smithii represents a highly active and clinically relevant methanogenic archaeon, being involved in gastrointestinal disorders, such as inflammatory bowel disease and obesity. Herein, we present an integrated approach using sequence and structure information to improve the annotation of M. smithii proteins using advanced protein structure prediction and annotation tools, such as AlphaFold2, trRosetta, ProFunc, and DeepFri. Of an initial set of 873 481 archaeal proteins, we found 707 754 proteins exclusively present in the human gut. Having analysed archaeal proteins together with 87 282 994 bacterial proteins, we identified unique archaeal proteins and archaeal-bacterial homologs. We then predicted and characterized functional domains and structures of 73 unique and homologous archaeal protein clusters linked the human gut and M. smithii. We refined annotations based on the predicted structures, extending existing sequence similarity-based annotations. We identified gut-specific archaeal proteins that may be involved in defense mechanisms, virulence, adhesion, and the degradation of toxic substances. Interestingly, we identified potential glycosyltransferases that could be associated with N-linked and O-glycosylation. Additionally, we found preliminary evidence for interdomain horizontal gene transfer between Clostridia species and M. smithii, which includes sporulation Stage V proteins AE and AD. Our study broadens the understanding of archaeal biology, particularly M. smithii, and highlights the importance of considering both sequence and structure for the prediction of protein function.

Challenges and Approaches of Culturing the Unculturable Archaea

Since Carl Woese's discovery of archaea as a third domain of life, numerous archaeal species have been discovered, yet archaeal diversity is poorly characterized. Culturing archaea is complicated, but several queries about archaeal cell biology, evolution, physiology, and diversity need to be solved by culturing and culture-dependent techniques. Increasing interest in demand for innovative culturing methods has led to various technological and methodological advances. The current review explains frequent hurdles hindering uncultured archaea isolation and discusses features for more archaeal cultivation. This review also discusses successful strategies and available media for archaeal culturing, which might be helpful for future culturing practices.

The Archaeome's Role in Colorectal Cancer: Unveiling the DPANN Group and Investigating Archaeal Functional Signatures

BACKGROUND AND AIMS

Gut microbial imbalances are linked to colorectal cancer (CRC), but archaea's role remains underexplored. Here, using previously published metagenomic data from different populations including Austria, Germany, Italy, Japan, China, and India, we performed bioinformatic and statistical analysis to identify archaeal taxonomic and functional signatures related to CRC.

METHODS

We analyzed published fecal metagenomic data from 390 subjects, comparing the archaeomes of CRC and healthy individuals. We conducted a biostatistical analysis to investigate the relationship between Candidatus Mancarchaeum acidiphilum (DPANN superphylum) and other archaeal species associated with CRC. Using the Prokka tool, we annotated the data focusing on archaeal genes, subsequently linking them to CRC and mapping them against UniprotKB and GO databases for specific archaeal gene functions.

RESULTS

Our analysis identified enrichment of methanogenic archaea in healthy subjects, with an exception for Methanobrevibacter smithii, which correlated with CRC. Notably, CRC showed a strong association with archaeal species, particularly Natrinema sp. J7-2, Ferroglobus placidus, and Candidatus Mancarchaeum acidiphilum. Furthermore, the DPANN archaeon exhibited a significant correlation with other CRC-associated archaea (p < 0.001). Functionally, we found a marked association between MvhB-type polyferredoxin and colorectal cancer. We also highlighted the association of archaeal proteins involved in the biosynthesis of leucine and the galactose metabolism process with the healthy phenotype.

CONCLUSIONS

The archaeomes of CRC patients show identifiable alterations, including a decline in methanogens and an increase in Halobacteria species. MvhB-type polyferredoxin, linked with CRC and species like Candidatus Mancarchaeum acidiphilum, Natrinema sp. J7-2, and Ferroglobus placidus emerge as potential archaeal biomarkers. Archaeal proteins may also offer gut protection, underscoring archaea's role in CRC dynamics.

Widespread PRC barrel proteins play critical roles in archaeal cell division

Cell division is fundamental to all cellular life. Most of the archaea employ one of two alternative division machineries, one centered around the prokaryotic tubulin homolog FtsZ and the other around the endosomal sorting complex required for transport (ESCRT). However, neither of these mechanisms has been thoroughly characterized in archaea. Here, we show that three of the four PRC (Photosynthetic Reaction Center) barrel domain proteins of Haloferax volcanii (renamed Cell division proteins B1/2/3 (CdpB1/2/3)), play important roles in division. CdpB1 interacts directly with the FtsZ membrane anchor SepF and is essential for division, whereas deletion of cdpB2 and cdpB3 causes a major and a minor division defect, respectively. Orthologs of CdpB proteins are also involved in cell division in other haloarchaea. Phylogenetic analysis shows that PRC barrel proteins are widely distributed among archaea, including the highly conserved CdvA protein of the crenarchaeal ESCRT-based division system. Thus, diverse PRC barrel proteins appear to be central to cell division in most if not all archaea. Further study of these proteins is expected to elucidate the division mechanisms in archaea and their evolution.

Genomics, Transcriptomics, and Proteomics of SSV1 and Related Fusellovirus: A Minireview

Saccharolobus spindle-shaped virus 1 (SSV1) was one of the first viruses identified in the archaeal kingdom. Originally isolated from a Japanese species of Saccharolobus back in 1984, it has been extensively used as a model system for genomic, transcriptomic, and proteomic studies, as well as to unveil the molecular mechanisms governing the host–virus interaction. The purpose of this mini review is to supply a compendium of four decades of research on the SSV1 virus.

A roadmap for the functional annotation of protein families: a community perspective

Over the last 25 years, biology has entered the genomic era and is becoming a science of ‘big data’. Most interpretations of genomic analyses rely on accurate functional annotations of the proteins encoded by more than 500 000 genomes sequenced to date. By different estimates, only half the predicted sequenced proteins carry an accurate functional annotation, and this percentage varies drastically between different organismal lineages. Such a large gap in knowledge hampers all aspects of biological enterprise and, thereby, is standing in the way of genomic biology reaching its full potential. A brainstorming meeting to address this issue funded by the National Science Foundation was held during 3–4 February 2022. Bringing together data scientists, biocurators, computational biologists and experimentalists within the same venue allowed for a comprehensive assessment of the current state of functional annotations of protein families. Further, major issues that were obstructing the field were identified and discussed, which ultimately allowed for the proposal of solutions on how to move forward.

Functional characterization of prokaryotic dark matter: the road so far and what lies ahead

Eight-hundred thousand to one trillion prokaryotic species may inhabit our planet. Yet, fewer than two-hundred thousand prokaryotic species have been described. This uncharted fraction of microbial diversity, and its undisclosed coding potential, is known as the "microbial dark matter" (MDM). Next-generation sequencing has allowed to collect a massive amount of genome sequence data, leading to unprecedented advances in the field of genomics. Still, harnessing new functional information from the genomes of uncultured prokaryotes is often limited by standard classification methods. These methods often rely on sequence similarity searches against reference genomes from cultured species. This hinders the discovery of unique genetic elements that are missing from the cultivated realm. It also contributes to the accumulation of prokaryotic gene products of unknown function among public sequence data repositories, highlighting the need for new approaches for sequencing data analysis and classification. Increasing evidence indicates that these proteins of unknown function might be a treasure trove of biotechnological potential. Here, we outline the challenges, opportunities, and the potential hidden within the functional dark matter (FDM) of prokaryotes. We also discuss the pitfalls surrounding molecular and computational approaches currently used to probe these uncharted waters, and discuss future opportunities for research and applications.

The tRNA discriminator base defines the mutual orthogonality of two distinct pyrrolysyl-tRNA synthetase/tRNAPyl pairs in the same organism

Site-specific incorporation of distinct non-canonical amino acids into proteins via genetic code expansion requires mutually orthogonal aminoacyl-tRNA synthetase/tRNA pairs. Pyrrolysyl-tRNA synthetase (PylRS)/tRNAPyl pairs are ideal for genetic code expansion and have been extensively engineered for developing mutually orthogonal pairs. Here, we identify two novel wild-type PylRS/tRNAPyl pairs simultaneously present in the deep-rooted extremely halophilic euryarchaeal methanogen Candidatus Methanohalarchaeum thermophilum HMET1, and show that both pairs are functional in the model halophilic archaeon Haloferax volcanii. These pairs consist of two different PylRS enzymes and two distinct tRNAs with dissimilar discriminator bases. Surprisingly, these two PylRS/tRNAPyl pairs display mutual orthogonality enabled by two unique features, the A73 discriminator base of tRNAPyl2 and a shorter motif 2 loop in PylRS2. In vivo translation experiments show that tRNAPyl2 charging by PylRS2 is defined by the enzyme's shortened motif 2 loop. Finally, we demonstrate that the two HMET1 PylRS/tRNAPyl pairs can simultaneously decode UAG and UAA codons for incorporation of two distinct noncanonical amino acids into protein. This example of a single base change in a tRNA leading to additional coding capacity suggests that the growth of the genetic code is not yet limited by the number of identity elements fitting into the tRNA structure.

Capillary Electrophoresis-Based Functional Genomics Screening to Discover Novel Archaeal DNA Modifying Enzymes

It has been predicted that 30 to 80% of archaeal genomes remain annotated as hypothetical proteins with no assigned gene function. Further, many archaeal organisms are difficult to grow or are unculturable. To overcome these technical and experimental hurdles, we developed a high-throughput functional genomics screen that utilizes capillary electrophoresis (CE) to identify nucleic acid modifying enzymes based on activity rather than sequence homology. Here, we describe a functional genomics screening workflow to find DNA modifying enzyme activities encoded by the hyperthermophile Thermococcus kodakarensis (T. kodakarensis). Large DNA insert fosmid libraries representing an ∼5-fold average coverage of the T. kodakarensis genome were prepared in Escherichia coli. RNA-seq showed a high fraction (84%) of T. kodakarensis genes were transcribed in E. coli despite differences in promoter structure and translational machinery. Our high-throughput screening workflow used fluorescently labeled DNA substrates directly in heat-treated lysates of fosmid clones with capillary electrophoresis detection of reaction products. Using this method, we identified both a new DNA endonuclease activity for a previously described RNA endonuclease (Nob1) and a novel AP lyase DNA repair enzyme family (termed 'TK0353') that is found only in a small subset of Thermococcales. The screening methodology described provides a fast and efficient way to explore the T. kodakarensis genome for a variety of nucleic acid modifying activities and may have implications for similar exploration of enzymes and pathways that underlie core cellular processes in other Archaea. IMPORTANCE This study provides a rapid, simple, high-throughput method to discover novel archaeal nucleic acid modifying enzymes by utilizing a fosmid genomic library, next-generation sequencing, and capillary electrophoresis. The method described here provides the details necessary to create 384-well fosmid library plates from Thermococcus kodakarensis genomic DNA, sequence 384-well fosmids plates using Illumina next-generation sequencing, and perform high-throughput functional read-out assays using capillary electrophoresis to identify a variety of nucleic acid modifying activities, including DNA cleavage and ligation. We used this approach to identify a new DNA endonuclease activity for a previously described RNA endonuclease (Nob1) and identify a novel AP lyase enzyme (TK0353) that lacks sequence homology to known nucleic acid modifying enzymes.

Small Proteins in Archaea, a Mainly Unexplored World

In recent years, increasing numbers of small proteins have moved into the focus of science. Small proteins have been identified and characterized in all three domains of life, but the majority remains functionally uncharacterized, lack secondary structure, and exhibit limited evolutionary conservation. While quite a few have already been described for bacteria and eukaryotic organisms, the amount of known and functionally analyzed archaeal small proteins is still very limited. In this review, we compile the current state of research, show strategies for systematic approaches for global identification of small archaeal proteins, and address selected functionally characterized examples. Besides, we document exemplarily for one archaeon the tool development and optimization to identify small proteins using genome-wide approaches.

OUP accepted manuscript

Ammonia-oxidising archaea (AOA) are environmentally important microorganisms involved in the biogeochemical cycling of nitrogen. Routine cultivation of AOA is exclusively performed in liquid cultures and reports on their growth on solid medium are scarce. The ability to grow AOA on solid medium would be beneficial for not only the purification of enrichment cultures but also for developing genetic tools. The aim of this study was to develop a reliable method for growing individual colonies from AOA cultures on solid medium. Three phylogenetically distinct AOA strains were tested: ‘Candidatus Nitrosocosmicus franklandus C13’, Nitrososphaera viennensis EN76 and ‘Candidatus Nitrosotalea sinensis Nd2’. Of the gelling agents tested, agar and Bacto-agar severely inhibited growth of all three strains. In contrast, both ‘Ca. N. franklandus C13’ and N. viennensis EN76 tolerated Phytagel™ while the acidophilic ‘Ca. N. sinensis Nd2’ was completely inhibited. Based on these observations, we developed a Liquid-Solid (LS) method that involves immobilising cells in Phytagel™ and overlaying with liquid medium. This approach resulted in the development of visible distinct colonies from ‘Ca. N. franklandus C13’ and N. viennensis EN76 cultures and lays the groundwork for the genetic manipulation of this group of microorganisms.

Assessment of the plasmidome of an extremophilic microbial community from the Diamante Lake, Argentina

Diamante Lake located at 4589 m.a.s.l. in the Andean Puna constitutes an extreme environment. It is exposed to multiple extreme conditions such as an unusually high concentration of arsenic (over 300 mg L^-1) and low oxygen pressure. Microorganisms thriving in the lake display specific genotypes that facilitate survival, which include at least a multitude of plasmid-encoded resistance traits. Hence, the genetic information provided by the plasmids essentially contributes to understand adaptation to different stressors. Though plasmids from cultivable organisms have already been analyzed to the sequence level, the impact of the entire plasmid-borne genetic information on such microbial ecosystem is not known. This study aims at assessing the plasmidome from Diamante Lake, which facilitates the identification of potential hosts and prediction of gene functions as well as the ecological impact of mobile genetic elements. The deep-sequencing analysis revealed a large fraction of previously unknown DNA sequences of which the majority encoded putative proteins of unknown function. Remarkably, functions related to the oxidative stress response, DNA repair, as well as arsenic- and antibiotic resistances were annotated. Additionally, all necessary capacities related to plasmid replication, mobilization and maintenance were detected. Sequences characteristic for megaplasmids and other already known plasmid-associated genes were identified as well. The study highlights the potential of the deep-sequencing approach specifically targeting plasmid populations as it allows to evaluate the ecological impact of plasmids from (cultivable and non-cultivable) microorganisms, thereby contributing to the understanding of the distribution of resistance factors within an extremophilic microbial community.

Unprecedented Diversity of Unique CRISPR-Cas-Related Systems and Cas1 Homologs in Asgard Archaea

The principal function of archaeal and bacterial CRISPR-Cas systems is antivirus adaptive immunity. However, recent genome analyses identified a variety of derived CRISPR-Cas variants at least some of which appear to perform different functions. Here, we describe a unique repertoire of CRISPR-Cas-related systems that we discovered by searching archaeal metagenome-assemble genomes of the Asgard superphylum. Several of these variants contain extremely diverged homologs of Cas1, the integrase involved in CRISPR adaptation as well as casposon transposition. Strikingly, the diversity of Cas1 in Asgard archaea alone is greater than that detected so far among the rest of archaea and bacteria. The Asgard CRISPR-Cas derivatives also encode distinct forms of Cas4, Cas5, and Cas7 proteins, and/or additional nucleases. Some of these systems are predicted to perform defense functions, but possibly not programmable ones, whereas others are likely to represent previously unknown mobile genetic elements.

Living to the High Extreme: Unraveling the Composition, Structure, and Functional Insights of Bacterial Communities Thriving in the Arsenic-Rich Salar de Huasco Altiplanic Ecosystem

Microbial communities inhabiting extreme environments such as Salar de Huasco (SH) in northern Chile are adapted to thrive while exposed to several abiotic pressures and the presence of toxic elements such as arsenic (As). Hence, we aimed to uncover the role of As in shaping bacterial composition, structure, and functional potential in five different sites in this altiplanic wetland using a shotgun metagenomic approach. The sites exhibit wide gradients of As (9 to 321 mg/kg), and our results showed highly diverse communities and a clear dominance exerted by the Proteobacteria and Bacteroidetes phyla. Functional potential analyses show broadly convergent patterns, contrasting with their great taxonomic variability. As-related metabolism, as well as other functional categories such as those related to the CH4 and S cycles, differs among the five communities. Particularly, we found that the distribution and abundance of As-related genes increase as the As concentration rises. Approximately 75% of the detected genes for As metabolism belong to expulsion mechanisms; arsJ and arsP pumps are related to sites with higher As concentrations and are present almost exclusively in Proteobacteria. Furthermore, taxonomic diversity and functional potential are reflected in the 12 reconstructed high-quality metagenome assembled genomes (MAGs) belonging to the Bacteroidetes (5), Proteobacteria (5), Cyanobacteria (1), and Gemmatimonadetes (1) phyla. We conclude that SH microbial communities are diverse and possess a broad genetic repertoire to thrive under extreme conditions, including increasing concentrations of highly toxic As. Finally, this environment represents a reservoir of unknown and undescribed microorganisms, with great metabolic versatility, which needs further study. IMPORTANCE As microbial communities inhabiting extreme environments are fundamental for maintaining ecosystems, many studies concerning composition, functionality, and interactions have been carried out. However, much is still unknown. Here, we sampled microbial communities in the Salar de Huasco, an extreme environment subjected to several abiotic stresses (high UV radiation, salinity and arsenic; low pressure and temperatures). We found that although microbes are taxonomically diverse, functional potential seems to have an important degree of convergence, suggesting high levels of adaptation. Particularly, arsenic metabolism showed differences associated with increasing concentrations of the metalloid throughout the area, and it effectively exerts a significant pressure over these organisms. Thus, the significance of this research is that we describe highly specialized communities thriving in little-explored environments subjected to several pressures, considered analogous of early Earth and other planets, that have the potential for unraveling technologies to face the repercussions of climate change in many areas of interest.

Fervidicoccus fontis Strain 3639Fd, the First Crenarchaeon Capable of Growth on Lipids

Abstract—

Up to now, ability of prokaryotes to grow on lipids has been shown only for bacteria. Thermococcus sibiricus, member of the phylum Euryarchaeota isolated from a high-temperature oil well and capable of growth on olive oil, is the only exception. The present work reports isolation of a pure culture of a strictly anaerobic archaeon, strain 3639Fd (=VKM B-3509, =KCTC 25228) from a Kamchatka thermal spring, capable of growth on various lipids (tributyrin, triolein, and sesame, cottonseed, and sunflower oil) at 70°C and pH 5.5–6.0. Growth on tributyrin resulted in formation of butyrate, CO2, and hydrogen. According to the results of the 16S rRNA gene sequence analysis and in silico DNA–DNA hybridization, the isolate was classified as a strain of Fervidicoccus fontis, an archaeon of the phylum Crenarchaeota. The closest characterized homologs of the α/β-hydrolases, encoded in the genomes of F. fontis 3639Fd and of the type strain of this species, Kam940T, were various carboxyl esterases (EC 3.1.1), the enzymes responsible for lipid hydrolysis. Thus, F. fontis is the first crenarchaeon able to obtain energy by hydrolysis of lipid substrates.

Microbial Arsenal of Antiviral Defenses - Part I

Bacteriophages or phages are viruses that infect bacterial cells (for the scope of this review we will also consider viruses that infect Archaea). Constant threat of phage infection is a major force that shapes evolution of the microbial genomes. To withstand infection, bacteria had evolved numerous strategies to avoid recognition by phages or to directly interfere with phage propagation inside the cell. Classical molecular biology and genetic engineering have been deeply intertwined with the study of phages and host defenses. Nowadays, owing to the rise of phage therapy, broad application of CRISPR-Cas technologies, and development of bioinformatics approaches that facilitate discovery of new systems, phage biology experiences a revival. This review describes variety of strategies employed by microbes to counter phage infection, with a focus on novel systems discovered in recent years. First chapter covers defense associated with cell surface, role of small molecules, and innate immunity systems relying on DNA modification.

GCsnap: Interactive Snapshots for the Comparison of Protein-Coding Genomic Contexts

The biological function and evolutionary history of protein-coding genes are not only written in their nucleotide sequences. The comparison of genomic contexts throughout different lineages may highlight genomic mechanisms in the generation of new protein families, while the conservation of gene clusters may unravel, for instance, metabolic pathways. Various tools and databases exist that allow for the analysis and comparison of genomic contexts, but each has its own limitations. Online databases allow for quick comparisons, but only for those genomes for which data were pre-calculated. More advanced tools may allow for the comparison of any genome, but are often limited to a given phylogenetic kingdom or provide only a snapshot of the genomic contexts without further information about the genes involved. Here, we introduce GCsnap, a flexible Python-based tool that allows for the interactive comparison of the genomic contexts of protein-coding genes from any genome at any taxonomic level, integrating them with functional and structural information. By connecting the output to different protein databases, users can navigate through the different genomic contexts from a simple interactive platform, facilitating the further analysis of the contexts found. GCsnap is not limited to a single input format, can perform batch jobs and accepts protein classification maps. Results are stored in detailed, human and machine-readable files, and customizable, publication-ready figures. GCsnap is freely available from https://github.com/JoanaMPereira/GCsnap and can be set up easily on any computer.

Microbial dark matter filling the niche in hypersaline microbial mats

BACKGROUND

Shark Bay, Australia, harbours one of the most extensive and diverse systems of living microbial mats that are proposed to be analogs of some of the earliest ecosystems on Earth. These ecosystems have been shown to possess a substantial abundance of uncultivable microorganisms. These enigmatic microbes, jointly coined as 'microbial dark matter' (MDM), are hypothesised to play key roles in modern microbial mats.

RESULTS

We reconstructed 115 metagenome-assembled genomes (MAGs) affiliated to MDM, spanning 42 phyla. This study reports for the first time novel microorganisms (Zixibacterial order GN15) putatively taking part in dissimilatory sulfate reduction in surface hypersaline settings, as well as novel eukaryote signature proteins in the Asgard archaea. Despite possessing reduced-size genomes, the MDM MAGs are capable of fermenting and degrading organic carbon, suggesting a role in recycling organic carbon. Several forms of RuBisCo were identified, allowing putative CO₂ incorporation into nucleotide salvaging pathways, which may act as an alternative carbon and phosphorus source. High capacity of hydrogen production was found among Shark Bay MDM. Putative schizorhodopsins were also identified in Parcubacteria, Asgard archaea, DPANN archaea, and Bathyarchaeota, allowing these members to potentially capture light energy. Diversity-generating retroelements were prominent in DPANN archaea that likely facilitate the adaptation to a dynamic, host-dependent lifestyle.

CONCLUSIONS

This is the first study to reconstruct and describe in detail metagenome-assembled genomes (MAGs) affiliated with microbial dark matter in hypersaline microbial mats. Our data suggests that these microbial groups are major players in these systems. In light of our findings, we propose H₂, ribose and CO/CO₂ as the main energy currencies of the MDM community in these mat systems. Video Abstract.

Evolution: Two Domains of Life or Three?

That Bacteria, Archaea and Eukarya (eukaryotes) represent three separate domains of Life, no one having evolved from within any other, has been taken as fact for three decades. Recent work shows this to be untrue. Eukarya arose from well within Archaea and are specifically related to newly discovered archaeal species with eukaryote-like features.

The Archaeal Proteome Project advances knowledge about archaeal cell biology through comprehensive proteomics

While many aspects of archaeal cell biology remain relatively unexplored, systems biology approaches like mass spectrometry (MS) based proteomics offer an opportunity for rapid advances. Unfortunately, the enormous amount of MS data generated often remains incompletely analyzed due to a lack of sophisticated bioinformatic tools and field-specific biological expertise for data interpretation. Here we present the initiation of the Archaeal Proteome Project (ArcPP), a community-based effort to comprehensively analyze archaeal proteomes. Starting with the model archaeon Haloferax volcanii, we reanalyze MS datasets from various strains and culture conditions. Optimized peptide spectrum matching, with strict control of false discovery rates, facilitates identifying > 72% of the reference proteome, with a median protein sequence coverage of 51%. These analyses, together with expert knowledge in diverse aspects of cell biology, provide meaningful insights into processes such as N-terminal protein maturation, N-glycosylation, and metabolism. Altogether, ArcPP serves as an invaluable blueprint for comprehensive prokaryotic proteomics.

Genome Analyses and Genome-Centered Metatranscriptomics of Methanothermobacter wolfeii Strain SIV6, Isolated from a Thermophilic Production-Scale Biogas Fermenter

In the thermophilic biogas-producing microbial community, the genus Methanothermobacter was previously described to be frequently abundant. The aim of this study was to establish and analyze the genome sequence of the archaeal strain Methanothermobacter wolfeii SIV6 originating from a thermophilic industrial-scale biogas fermenter and compare it to related reference genomes. The circular chromosome has a size of 1,686,891 bases, featuring a GC content of 48.89%. Comparative analyses considering three completely sequenced Methanothermobacter strains revealed a core genome of 1494 coding sequences and 16 strain specific genes for M. wolfeii SIV6, which include glycosyltransferases and CRISPR/cas associated genes. Moreover, M. wolfeii SIV6 harbors all genes for the hydrogenotrophic methanogenesis pathway and genome-centered metatranscriptomics indicates the high metabolic activity of this strain, with 25.18% of all transcripts per million (TPM) belong to the hydrogenotrophic methanogenesis pathway and 18.02% of these TPM exclusively belonging to the mcr operon. This operon encodes the different subunits of the enzyme methyl-coenzyme M reductase (EC: 2.8.4.1), which catalyzes the final and rate-limiting step during methanogenesis. Finally, fragment recruitment of metagenomic reads from the thermophilic biogas fermenter on the SIV6 genome showed that the strain is abundant (1.2%) within the indigenous microbial community. Detailed analysis of the archaeal isolate M. wolfeii SIV6 indicates its role and function within the microbial community of the thermophilic biogas fermenter, towards a better understanding of the biogas production process and a microbial-based management of this complex process.

Current Status and Potential Applications of Underexplored Prokaryotes

Thousands of prokaryotic genera have been published, but methodological bias in the study of prokaryotes is noted. Prokaryotes that are relatively easy to isolate have been well-studied from multiple aspects. Massive quantities of experimental findings and knowledge generated from the well-known prokaryotic strains are inundating scientific publications. However, researchers may neglect or pay little attention to the uncommon prokaryotes and hard-to-cultivate microorganisms. In this review, we provide a systematic update on the discovery of underexplored culturable and unculturable prokaryotes and discuss the insights accumulated from various research efforts. Examining these neglected prokaryotes may elucidate their novelties and functions and pave the way for their industrial applications. In addition, we hope that this review will prompt the scientific community to reconsider these untapped pragmatic resources.

Diversity of "Ca. Micrarchaeota" in Two Distinct Types of Acidic Environments and Their Associations with Thermoplasmatales

“Candidatus Micrarchaeota” are widely distributed in acidic environments; however, their cultivability and our understanding of their interactions with potential hosts are very limited. Their habitats were so far attributed with acidic sites, soils, peats, freshwater systems, and hypersaline mats. Using cultivation and culture-independent approaches (16S rRNA gene clonal libraries, high-throughput amplicon sequencing of V3-V4 region of 16S rRNA genes), we surveyed the occurrence of these archaea in geothermal areas on Kamchatka Peninsula and Kunashir Island and assessed their taxonomic diversity in relation with another type of low-pH environment, acid mine drainage stream (Wales, UK). We detected “Ca. Micrarchaeota” in thermophilic heterotrophic enrichment cultures of Kunashir and Kamchatka that appeared as two different phylotypes, namely “Ca. Mancarchaeum acidiphilum”-, and ARMAN-2-related, alongside their potential hosts, Cuniculiplasma spp. and other Thermoplasmatales archaea without defined taxonomic position. These clusters of “Ca. Micrarchaeota” together with three other groups were also present in mesophilic acid mine drainage community. Present work expands our knowledge on the diversity of “Ca. Micrarchaeota” in thermophilic and mesophilic acidic environments, suggests cultivability patterns of acidophilic archaea and establishes potential links between low-abundance species of thermophilic “Ca. Micrarchaeota” and certain Thermoplasmatales, such as Cuniculiplasma spp. in situ.