Minimal and hybrid hydrogenases are active from archaea.

Microbial hydrogen (H2) cycling underpins the diversity and functionality of diverse anoxic ecosystems. Among the three evolutionarily distinct hydrogenase superfamilies responsible, [FeFe] hydrogenases were thought to be restricted to bacteria and eukaryotes. Here, we show that anaerobic archaea encode diverse, active, and ancient lineages of [FeFe] hydrogenases through combining analysis of existing and new genomes with extensive biochemical experiments. [FeFe] hydrogenases are encoded by genomes of nine archaeal phyla and expressed by H2-producing Asgard archaeon cultures. We report an ultraminimal hydrogenase in DPANN archaea that binds the catalytic H-cluster and produces H2. Moreover, we identify and characterize remarkable hybrid complexes formed through the fusion of [FeFe] and [NiFe] hydrogenases in ten other archaeal orders. Phylogenetic analysis and structural modeling suggest a deep evolutionary history of hybrid hydrogenases. These findings reveal new metabolic adaptations of archaea, streamlined H2 catalysts for biotechnological development, and a surprisingly intertwined evolutionary history between the two major H2-metabolizing enzymes.

Diverse virus-encoded CRISPR-Cas systems include streamlined genome editors.

CRISPR-Cas systems are host-encoded pathways that protect microbes from viral infection using an adaptive RNA-guided mechanism. Using genome-resolved metagenomics, we find that CRISPR systems are also encoded in diverse bacteriophages, where they occur as divergent and hypercompact anti-viral systems. Bacteriophage-encoded CRISPR systems belong to all six known CRISPR-Cas types, though some lack crucial components, suggesting alternate functional roles or host complementation. We describe multiple new Cas9-like proteins and 44 families related to type V CRISPR-Cas systems, including the Casλ RNA-guided nuclease family. Among the most divergent of the new enzymes identified, Casλ recognizes double-stranded DNA using a uniquely structured CRISPR RNA (crRNA). The Casλ-RNA-DNA structure determined by cryoelectron microscopy reveals a compact bilobed architecture capable of inducing genome editing in mammalian, Arabidopsis, and hexaploid wheat cells. These findings reveal a new source of CRISPR-Cas enzymes in phages and highlight their value as genome editors in plant and human cells.

Major New Microbial Groups Expand Diversity and Alter our Understanding of the Tree of Life.

The recent recovery of genomes for organisms from phyla with no isolated representative (candidate phyla) via cultivation-independent genomics enabled delineation of major new microbial lineages, namely the bacterial candidate phyla radiation (CPR), DPANN archaea, and Asgard archaea. CPR and DPANN organisms are inferred to be mostly symbionts, and some are episymbionts of other microbial community members. Asgard genomes encode typically eukaryotic systems, and their inclusion in phylogenetic analyses results in placement of eukaryotes as a branch within Archaea. Here, we illustrate how new genomes have changed the structure of the tree of life and altered our understanding of biology, evolution, and metabolic roles in biogeochemical processes.

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.

Habitat Transition in the Evolution of Bacteria and Archaea.

Related groups of microbes are widely distributed across Earth's habitats, implying numerous dispersal and adaptation events over evolutionary time. However, relatively little is known about the characteristics and mechanisms of these habitat transitions, particularly for populations that reside in animal microbiomes. Here, we review the literature concerning habitat transitions among a variety of bacterial and archaeal lineages, considering the frequency of migration events, potential environmental barriers, and mechanisms of adaptation to new physicochemical conditions, including the modification of protein inventories and other genomic characteristics. Cells dependent on microbial hosts, particularly bacteria from the Candidate Phyla Radiation, have undergone repeated habitat transitions from environmental sources into animal microbiomes. We compare their trajectories to those of both free-living cells-including the Melainabacteria, Elusimicrobia, and methanogenic archaea-and cellular endosymbionts and bacteriophages, which have made similar transitions. We conclude by highlighting major related topics that may be worthy of future study.

Borgs are giant genetic elements with potential to expand metabolic capacity.

Anaerobic methane oxidation exerts a key control on greenhouse gas emissions1, yet factors that modulate the activity of microorganisms performing this function remain poorly understood. Here we discovered extraordinarily large, diverse DNA sequences that primarily encode hypothetical proteins through studying groundwater, sediments and wetland soil where methane production and oxidation occur. Four curated, complete genomes are linear, up to approximately 1 Mb in length and share genome organization, including replichore structure, long inverted terminal repeats and genome-wide unique perfect tandem direct repeats that are intergenic or generate amino acid repeats. We infer that these are highly divergent archaeal extrachromosomal elements with a distinct evolutionary origin. Gene sequence similarity, phylogeny and local divergence of sequence composition indicate that many of their genes were assimilated from methane-oxidizing Methanoperedens archaea. We refer to these elements as 'Borgs'. We identified at least 19 different Borg types coexisting with Methanoperedens spp. in four distinct ecosystems. Borgs provide methane-oxidizing Methanoperedens archaea access to genes encoding proteins involved in redox reactions and energy conservation (for example, clusters of multihaem cytochromes and methyl coenzyme M reductase). These data suggest that Borgs might have previously unrecognized roles in the metabolism of this group of archaea, which are known to modulate greenhouse gas emissions, but further studies are now needed to establish their functional relevance.

Petabase-scale sequence alignment catalyses viral discovery.

Public databases contain a planetary collection of nucleic acid sequences, but their systematic exploration has been inhibited by a lack of efficient methods for searching this corpus, which (at the time of writing) exceeds 20 petabases and is growing exponentially1. Here we developed a cloud computing infrastructure, Serratus, to enable ultra-high-throughput sequence alignment at the petabase scale. We searched 5.7 million biologically diverse samples (10.2 petabases) for the hallmark gene RNA-dependent RNA polymerase and identified well over 105 novel RNA viruses, thereby expanding the number of known species by roughly an order of magnitude. We characterized novel viruses related to coronaviruses, hepatitis delta virus and huge phages, respectively, and analysed their environmental reservoirs. To catalyse the ongoing revolution of viral discovery, we established a free and comprehensive database of these data and tools. Expanding the known sequence diversity of viruses can reveal the evolutionary origins of emerging pathogens and improve pathogen surveillance for the anticipation and mitigation of future pandemics.

A scoutRNA Is Required for Some Type V CRISPR-Cas Systems.

CRISPR-Cas12c/d proteins share limited homology with Cas12a and Cas9 bacterial CRISPR RNA (crRNA)-guided nucleases used widely for genome editing and DNA detection. However, Cas12c (C2c3)- and Cas12d (CasY)-catalyzed DNA cleavage and genome editing activities have not been directly observed. We show here that a short-complementarity untranslated RNA (scoutRNA), together with crRNA, is required for Cas12d-catalyzed DNA cutting. The scoutRNA differs in secondary structure from previously described tracrRNAs used by CRISPR-Cas9 and some Cas12 enzymes, and in Cas12d-containing systems, scoutRNA includes a conserved five-nucleotide sequence that is essential for activity. In addition to supporting crRNA-directed DNA recognition, biochemical and cell-based experiments establish scoutRNA as an essential cofactor for Cas12c-catalyzed pre-crRNA maturation. These results define scoutRNA as a third type of transcript encoded by a subset of CRISPR-Cas genomic loci and explain how Cas12c/d systems avoid requirements for host factors including ribonuclease III for bacterial RNA-mediated adaptive immunity.

The NIH Somatic Cell Genome Editing program.

The move from reading to writing the human genome offers new opportunities to improve human health. The United States National Institutes of Health (NIH) Somatic Cell Genome Editing (SCGE) Consortium aims to accelerate the development of safer and more-effective methods to edit the genomes of disease-relevant somatic cells in patients, even in tissues that are difficult to reach. Here we discuss the consortium's plans to develop and benchmark approaches to induce and measure genome modifications, and to define downstream functional consequences of genome editing within human cells. Central to this effort is a rigorous and innovative approach that requires validation of the technology through third-party testing in small and large animals. New genome editors, delivery technologies and methods for tracking edited cells in vivo, as well as newly developed animal models and human biological systems, will be assembled-along with validated datasets-into an SCGE Toolkit, which will be disseminated widely to the biomedical research community. We visualize this toolkit-and the knowledge generated by its applications-as a means to accelerate the clinical development of new therapies for a wide range of conditions.

A Functional Mini-Integrase in a Two-Protein-type V-C CRISPR System.

CRISPR-Cas immunity requires integration of short, foreign DNA fragments into the host genome at the CRISPR locus, a site consisting of alternating repeat sequences and foreign-derived spacers. In most CRISPR systems, the proteins Cas1 and Cas2 form the integration complex and are both essential for DNA acquisition. Most type V-C and V-D systems lack the cas2 gene and have unusually short CRISPR repeats and spacers. Here, we show that a mini-integrase comprising the type V-C Cas1 protein alone catalyzes DNA integration with a preference for short (17- to 19-base-pair) DNA fragments. The mini-integrase has weak specificity for the CRISPR array. We present evidence that the Cas1 proteins form a tetramer for integration. Our findings support a model of a minimal integrase with an internal ruler mechanism that favors shorter repeats and spacers. This minimal integrase may represent the function of the ancestral Cas1 prior to Cas2 adoption.

Tandem repeats in giant archaeal Borg elements undergo rapid evolution and create new intrinsically disordered regions in proteins.

Borgs are huge, linear extrachromosomal elements associated with anaerobic methane-oxidizing archaea. Striking features of Borg genomes are pervasive tandem direct repeat (TR) regions. Here, we present six new Borg genomes and investigate the characteristics of TRs in all ten complete Borg genomes. We find that TR regions are rapidly evolving, recently formed, arise independently, and are virtually absent in host Methanoperedens genomes. Flanking partial repeats and A-enriched character constrain the TR formation mechanism. TRs can be in intergenic regions, where they might serve as regulatory RNAs, or in open reading frames (ORFs). TRs in ORFs are under very strong selective pressure, leading to perfect amino acid TRs (aaTRs) that are commonly intrinsically disordered regions. Proteins with aaTRs are often extracellular or membrane proteins, and functionally similar or homologous proteins often have aaTRs composed of the same amino acids. We propose that Borg aaTR-proteins functionally diversify Methanoperedens and all TRs are crucial for specific Borg-host associations and possibly cospeciation.

Clades of huge phages from across Earth's ecosystems.

Bacteriophages typically have small genomes1 and depend on their bacterial hosts for replication2. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems.

Rare ribosomal RNA sequences from archaea stabilize the bacterial ribosome.

The ribosome serves as the universally conserved translator of the genetic code into proteins and supports life across diverse temperatures ranging from below freezing to above 120°C. Ribosomes are capable of functioning across this wide range of temperatures even though the catalytic site for peptide bond formation, the peptidyl transferase center, is nearly universally conserved. Here we find that Thermoproteota, a phylum of thermophilic Archaea, substitute cytidine for uridine at large subunit rRNA positions 2554 and 2555 (Escherichia coli numbering) in the A loop, immediately adjacent to the binding site for the 3'-end of A-site tRNA. We show by cryo-EM that E. coli ribosomes with uridine to cytidine mutations at these positions retain the proper fold and post-transcriptional modification of the A loop. Additionally, these mutations do not affect cellular growth, protect the large ribosomal subunit from thermal denaturation, and increase the mutational robustness of nucleotides in the peptidyl transferase center. This work identifies sequence variation across archaeal ribosomes in the peptidyl transferase center that likely confers stabilization of the ribosome at high temperatures and develops a stable mutant bacterial ribosome that can act as a scaffold for future ribosome engineering efforts.

Eukaryotic RNA-guided endonucleases evolved from a unique clade of bacterial enzymes.

RNA-guided endonucleases form the crux of diverse biological processes and technologies, including adaptive immunity, transposition, and genome editing. Some of these enzymes are components of insertion sequences (IS) in the IS200/IS605 and IS607 transposon families. Both IS families encode a TnpA transposase and a TnpB nuclease, an RNA-guided enzyme ancestral to CRISPR-Cas12s. In eukaryotes, TnpB homologs occur as two distinct types, Fanzor1s and Fanzor2s. We analyzed the evolutionary relationships between prokaryotic TnpBs and eukaryotic Fanzors, which revealed that both Fanzor1s and Fanzor2s stem from a single lineage of IS607 TnpBs with unusual active site arrangement. The widespread nature of Fanzors implies that the properties of this particular lineage of IS607 TnpBs were particularly suited to adaptation in eukaryotes. Biochemical analysis of an IS607 TnpB and Fanzor1s revealed common strategies employed by TnpBs and Fanzors to co-evolve with their cognate transposases. Collectively, our results provide a new model of sequential evolution from IS607 TnpBs to Fanzor2s, and Fanzor2s to Fanzor1s that details how genes of prokaryotic origin evolve to give rise to new protein families in eukaryotes.

Weathered granites and soils harbour microbes with lanthanide-dependent methylotrophic enzymes.

BACKGROUND: Prior to soil formation, phosphate liberated by rock weathering is often sequestered into highly insoluble lanthanide phosphate minerals. Dissolution of these minerals releases phosphate and lanthanides to the biosphere. Currently, the microorganisms involved in phosphate mineral dissolution and the role of lanthanides in microbial metabolism are poorly understood. RESULTS: Although there have been many studies of soil microbiology, very little research has investigated microbiomes of weathered rock. Here, we sampled weathered granite and associated soil to identify the zones of lanthanide phosphate mineral solubilisation and genomically define the organisms implicated in lanthanide utilisation. We reconstructed 136 genomes from 11 bacterial phyla and found that gene clusters implicated in lanthanide-based metabolism of methanol (primarily xoxF3 and xoxF5) are surprisingly common in microbial communities in moderately weathered granite. Notably, xoxF3 systems were found in Verrucomicrobia for the first time, and in Acidobacteria, Gemmatimonadetes and Alphaproteobacteria. The xoxF-containing gene clusters are shared by diverse Acidobacteria and Gemmatimonadetes, and include conserved hypothetical proteins and transporters not associated with the few well studied xoxF systems. Given that siderophore-like molecules that strongly bind lanthanides may be required to solubilise lanthanide phosphates, it is notable that candidate metallophore biosynthesis systems were most prevalent in bacteria in moderately weathered rock, especially in Acidobacteria with lanthanide-based systems. CONCLUSIONS: Phosphate mineral dissolution, putative metallophore production and lanthanide utilisation by enzymes involved in methanol oxidation linked to carbonic acid production co-occur in the zone of moderate granite weathering. In combination, these microbial processes likely accelerate the conversion of granitic rock to soil.

Transporter genes in biosynthetic gene clusters predict metabolite characteristics and siderophore activity.

Biosynthetic gene clusters (BGCs) are operonic sets of microbial genes that synthesize specialized metabolites with diverse functions, including siderophores and antibiotics, which often require export to the extracellular environment. For this reason, genes for transport across cellular membranes are essential for the production of specialized metabolites and are often genomically colocalized with BGCs. Here, we conducted a comprehensive computational analysis of transporters associated with characterized BGCs. In addition to known exporters, in BGCs we found many importer-specific transmembrane domains that co-occur with substrate binding proteins possibly for uptake of siderophores or metabolic precursors. Machine learning models using transporter gene frequencies were predictive of known siderophore activity, molecular weights, and a measure of lipophilicity (log P) for corresponding BGC-synthesized metabolites. Transporter genes associated with BGCs were often equally or more predictive of metabolite features than biosynthetic genes. Given the importance of siderophores as pathogenicity factors, we used transporters specific for siderophore BGCs to identify both known and uncharacterized siderophore-like BGCs in genomes from metagenomes from the infant and adult gut microbiome. We find that 23% of microbial genomes from premature infant guts have siderophore-like BGCs, but only 3% of those assembled from adult gut microbiomes do. Although siderophore-like BGCs from the infant gut are predominantly associated with Enterobacteriaceae and Staphylococcus, siderophore-like BGCs can be identified from taxa in the adult gut microbiome that have rarely been recognized for siderophore production. Taken together, these results show that consideration of BGC-associated transporter genes can inform predictions of specialized metabolite structure and function.

Harnessing Fermentation May Enhance the Performance of Biological Sulfate-Reducing Bioreactors.

Biological sulfate reduction (BSR) represents a promising strategy for bioremediation of sulfate-rich waste streams, yet the impact of metabolic interactions on performance is largely unexplored. Here, genome-resolved metagenomics was used to characterize 17 microbial communities in reactors treating synthetic sulfate-contaminated solutions. Reactors were supplemented with lactate or acetate and a small amount of fermentable substrate. Of the 163 genomes representing all the abundant bacteria, 130 encode 321 NiFe and FeFe hydrogenases and all genomes of the 22 sulfate-reducing microorganisms (SRM) encode genes for H2 uptake. We observed lactate oxidation solely in the first packed bed reactor zone, with propionate and acetate oxidation in the middle and predominantly acetate oxidation in the effluent zone. The energetics of these reactions are very different, yet sulfate reduction kinetics were unaffected by the type of electron donor available. We hypothesize that the comparable rates, despite the typically slow growth of SRM on acetate, are a result of the consumption of H2 generated by fermentation. This is supported by the sustained performance of a predominantly acetate-supplemented stirred tank reactor dominated by diverse fermentative bacteria encoding FeFe hydrogenase genes and SRM capable of acetate and hydrogen consumption and CO2 assimilation. Thus, addition of fermentable substrates to stimulate syntrophic relationships may improve the performance of BSR reactors supplemented with inexpensive acetate.

Novel soil bacteria possess diverse genes for secondary metabolite biosynthesis.

In soil ecosystems, microorganisms produce diverse secondary metabolites such as antibiotics, antifungals and siderophores that mediate communication, competition and interactions with other organisms and the environment1,2. Most known antibiotics are derived from a few culturable microbial taxa 3 , and the biosynthetic potential of the vast majority of bacteria in soil has rarely been investigated 4 . Here we reconstruct hundreds of near-complete genomes from grassland soil metagenomes and identify microorganisms from previously understudied phyla that encode diverse polyketide and nonribosomal peptide biosynthetic gene clusters that are divergent from well-studied clusters. These biosynthetic loci are encoded by newly identified members of the Acidobacteria, Verrucomicobia and Gemmatimonadetes, and the candidate phylum Rokubacteria. Bacteria from these groups are highly abundant in soils5-7, but have not previously been genomically linked to secondary metabolite production with confidence. In particular, large numbers of biosynthetic genes were characterized in newly identified members of the Acidobacteria, which is the most abundant bacterial phylum across soil biomes 5 . We identify two acidobacterial genomes from divergent lineages, each of which encodes an unusually large repertoire of biosynthetic genes with up to fifteen large polyketide and nonribosomal peptide biosynthetic loci per genome. To track gene expression of genes encoding polyketide synthases and nonribosomal peptide synthetases in the soil ecosystem that we studied, we sampled 120 time points in a microcosm manipulation experiment and, using metatranscriptomics, found that gene clusters were differentially co-expressed in response to environmental perturbations. Transcriptional co-expression networks for specific organisms associated biosynthetic genes with two-component systems, transcriptional activation, putative antimicrobial resistance and iron regulation, linking metabolite biosynthesis to processes of environmental sensing and ecological competition. We conclude that the biosynthetic potential of abundant and phylogenetically diverse soil microorganisms has previously been underestimated. These organisms may represent a source of natural products that can address needs for new antibiotics and other pharmaceutical compounds.

Accurate and complete genomes from metagenomes.

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

"Candidatus Nealsonbacteria" Are Likely Biomass Recycling Ectosymbionts of Methanogenic Archaea in a Stable Benzene-Degrading Enrichment Culture.

The Candidate Phyla Radiation (CPR), also referred to as superphylum Patescibacteria, is a very large group of bacteria with no pure culture representatives discovered by 16S rRNA sequencing or genome-resolved metagenomic analyses of environmental samples. Within the CPR, candidate phylum Parcubacteria, previously referred to as OD1, is prevalent in anoxic sediments and groundwater. Previously, we had identified a specific member of the Parcubacteria (referred to as DGGOD1a) as an important member of a methanogenic benzene-degrading consortium. Phylogenetic analyses herein place DGGOD1a within the clade "Candidatus Nealsonbacteria." Because of its persistence over many years, we hypothesized that "Ca. Nealsonbacteria" DGGOD1a must play an important role in sustaining anaerobic benzene metabolism in the consortium. To try to identify its growth substrate, we amended the culture with a variety of defined compounds (pyruvate, acetate, hydrogen, DNA, and phospholipid), as well as crude culture lysate and three subfractions thereof. We observed the greatest (10-fold) increase in the absolute abundance of "Ca. Nealsonbacteria" DGGOD1a only when the consortium was amended with crude cell lysate. These results implicate "Ca. Nealsonbacteria" in biomass recycling. Fluorescence in situ hybridization and cryogenic transmission electron microscope images revealed that "Ca. Nealsonbacteria" DGGOD1a cells were attached to larger archaeal Methanothrix cells. This apparent epibiont lifestyle was supported by metabolic predictions from a manually curated complete genome. This is one of the first examples of bacterial-archaeal episymbiosis and may be a feature of other "Ca. Nealsonbacteria" found in anoxic environments. IMPORTANCE An anaerobic microbial enrichment culture was used to study members of candidate phyla that are difficult to grow in the lab. We were able to visualize tiny "Candidatus Nealsonbacteria" cells attached to a large Methanothrix cell, revealing a novel episymbiosis.