Unveiling the unknown viral world in groundwater.

Viruses as the prevailing biological entities are poorly understood in underground realms. Here, we establish the first metagenomic Groundwater Virome Catalogue (GWVC) comprising 280,420 viral species ( ≥ 5 kb) detected from 607 monitored wells in seven geo-environmental zones throughout China. In expanding ~10-fold the global portfolio of known groundwater viruses, we uncover over 99% novel viruses and about 95% novel viral clusters. By linking viruses to hosts from 119 prokaryotic phyla, we double the number of microbial phyla known to be virus-infected in groundwater. As keystone ultrasmall symbionts in aquifers, CPR bacteria and DPANN archaea are susceptible to virulent viruses. Certain complete CPR viruses even likely infect non-CPR bacteria, while partial CPR/DPANN viruses harbor cell-surface modification genes that assist symbiont cell adhesion to free-living microbes. This study reveals the unknown viral world and auxiliary metabolism associated with methane, nitrogen, sulfur, and phosphorus cycling in groundwater, and highlights the importance of subsurface virosphere in viral ecology.

Decomposing a San Francisco estuary microbiome using long-read metagenomics reveals species- and strain-level dominance from picoeukaryotes to viruses.

Although long-read sequencing has enabled obtaining high-quality and complete genomes from metagenomes, many challenges still remain to completely decompose a metagenome into its constituent prokaryotic and viral genomes. This study focuses on decomposing an estuarine metagenome to obtain a more accurate estimate of microbial diversity. To achieve this, we developed a new bead-based DNA extraction method, a novel bin refinement method, and obtained 150 Gbp of Nanopore sequencing. We estimate that there are ~500 bacterial and archaeal species in our sample and obtained 68 high-quality bins (>90% complete, <5% contamination, ≤5 contigs, contig length of >100 kbp, and all ribosomal and tRNA genes). We also obtained many contigs of picoeukaryotes, environmental DNA of larger eukaryotes such as mammals, and complete mitochondrial and chloroplast genomes and detected ~40,000 viral populations. Our analysis indicates that there are only a few strains that comprise most of the species abundances. IMPORTANCE: Ocean and estuarine microbiomes play critical roles in global element cycling and ecosystem function. Despite the importance of these microbial communities, many species still have not been cultured in the lab. Environmental sequencing is the primary way the function and population dynamics of these communities can be studied. Long-read sequencing provides an avenue to overcome limitations of short-read technologies to obtain complete microbial genomes but comes with its own technical challenges, such as needed sequencing depth and obtaining high-quality DNA. We present here new sampling and bioinformatics methods to attempt decomposing an estuarine microbiome into its constituent genomes. Our results suggest there are only a few strains that comprise most of the species abundances from viruses to picoeukaryotes, and to fully decompose a metagenome of this diversity requires 1 Tbp of long-read sequencing. We anticipate that as long-read sequencing technologies continue to improve, less sequencing will be needed.

Asgard archaea defense systems and their roles in the origin of eukaryotic immunity.

Dozens of new antiviral systems have been recently characterized in bacteria. Some of these systems are present in eukaryotes and appear to have originated in prokaryotes, but little is known about these defense mechanisms in archaea. Here, we explore the diversity and distribution of defense systems in archaea and identify 2610 complete systems in Asgardarchaeota, a group of archaea related to eukaryotes. The Asgard defense systems comprise 89 unique systems, including argonaute, NLR, Mokosh, viperin, Lassamu, and CBASS. Asgard viperin and argonaute proteins have structural homology to eukaryotic proteins, and phylogenetic analyses suggest that eukaryotic viperin proteins were derived from Asgard viperins. We show that Asgard viperins display anti-phage activity when heterologously expressed in bacteria. Eukaryotic and bacterial argonaute proteins appear to have originated in Asgardarchaeota, and Asgard argonaute proteins have argonaute-PIWI domains, key components of eukaryotic RNA interference systems. Our results support that Asgardarchaeota played important roles in the origin of antiviral defense systems in eukaryotes.

Unveiling the microbial realm with VEBA 2.0: a modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic and viral multi-omics from either short- or long-read sequencing.

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible software suite that bridges the gap between genomics and biotechnological solutions.

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.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes.

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

PacBio sequencing of human fecal samples uncovers the DNA methylation landscape of 22 673 gut phages.

Gut phages have an important impact on human health. Methylation plays key roles in DNA recognition, gene expression regulation and replication for phages. However, the DNA methylation landscape of gut phages is largely unknown. Here, with PacBio sequencing (2120×, 4785 Gb), we detected gut phage methylation landscape based on 22 673 gut phage genomes, and presented diverse methylation motifs and methylation differences in genomic elements. Moreover, the methylation rate of phages was associated with taxonomy and host, and N6-methyladenine methylation rate was higher in temperate phages than in virulent phages, suggesting an important role for methylation in phage-host interaction. In particular, 3543 (15.63%) phage genomes contained restriction-modification system, which could aid in evading clearance by the host. This study revealed the DNA methylation landscape of gut phage and its potential roles, which will advance the understanding of gut phage survival and human health.

Diverse Mechanisms of CRISPR-Cas9 Inhibition by Type II Anti-CRISPR Proteins.

Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas (CRISPR-associated proteins) systems provide bacteria and archaea with an adaptive immune response against invasion by mobile genetic elements like phages, plasmids, and transposons. These systems have been repurposed as very powerful biotechnological tools for gene editing applications in both bacterial and eukaryotic systems. The discovery of natural off-switches for CRISPR-Cas systems, known as anti-CRISPR proteins, provided a mechanism for controlling CRISPR-Cas activity and opened avenues for the development of more precise editing tools. In this review, we focus on the inhibitory mechanisms of anti-CRISPRs that are active against type II CRISPR-Cas systems and briefly discuss their biotechnological applications.

In Silico Approaches for Prediction of Anti-CRISPR Proteins.

Numerous viruses infecting bacteria and archaea encode CRISPR-Cas system inhibitors, known as anti-CRISPR proteins (Acr). The Acrs typically are highly specific for particular CRISPR variants, resulting in remarkable sequence and structural diversity and complicating accurate prediction and identification of Acrs. In addition to their intrinsic interest for understanding the coevolution of defense and counter-defense systems in prokaryotes, Acrs could be natural, potent on-off switches for CRISPR-based biotechnological tools, so their discovery, characterization and application are of major importance. Here we discuss the computational approaches for Acr prediction. Due to the enormous diversity and likely multiple origins of the Acrs, sequence similarity searches are of limited use. However, multiple features of protein and gene organization have been successfully harnessed to this end including small protein size and distinct amino acid compositions of the Acrs, association of acr genes in virus genomes with genes encoding helix-turn-helix proteins that regulate Acr expression (Acr-associated proteins, Aca), and presence of self-targeting CRISPR spacers in bacterial and archaeal genomes containing Acr-encoding proviruses. Productive approaches for Acr prediction also involve genome comparison of closely related viruses, of which one is resistant and the other one is sensitive to a particular CRISPR variant, and "guilt by association" whereby genes adjacent to a homolog of a known Aca are identified as candidate Acrs. The distinctive features of Acrs are employed for Acr prediction both by developing dedicated search algorithms and through machine learning. New approaches will be needed to identify novel types of Acrs that are likely to exist.

Non-canonical inhibition strategies and structural basis of anti-CRISPR proteins targeting type I CRISPR-Cas systems.

Mobile genetic elements (MGEs) such as bacteriophages and their host prokaryotes are trapped in an eternal battle against each other. To cope with foreign infection, bacteria and archaea have evolved multiple immune strategies, out of which CRISPR-Cas system is up to now the only discovered adaptive system in prokaryotes. Despite the fact that CRISPR-Cas system provides powerful and delicate protection against MGEs, MGEs have also evolved anti-CRISPR proteins (Acrs) to counteract the CRISPR-Cas immune defenses. To date, 46 families of Acrs targeting type I CRISPR-Cas system have been characterized, out of which structure information of 21 families have provided insights on their inhibition strategies. Here, we review the non-canonical inhibition strategies adopted by Acrs targeting type I CRISPR-Cas systems based on their structure information by incorporating the most recent advances in this field, and discuss our current understanding and future perspectives. The delicate interplay between type I CRISPR-Cas systems and their Acrs provides us with important insights into the ongoing fierce arms race between prokaryotic hosts and their predators.

Anti-CRISPR Discovery: Using Magnets to Find Needles in Haystacks.

CRISPR-Cas immune systems in bacteria and archaea protect against viral infection, which has spurred viruses to develop dedicated inhibitors of these systems called anti-CRISPRs (Acrs). Like most host-virus arms races, many diverse examples of these immune and counter-immune proteins are encoded by the genomes of bacteria, archaea, and their viruses. For the case of Acrs, it is almost certain that just a small minority of nature's true diversity has been described. In this review, I discuss the various approaches used to identify these Acrs and speculate on the future for Acr discovery. Because Acrs can determine infection outcomes in nature and regulate CRISPR-Cas activities in applied settings, they have a dual importance to both host-virus conflicts and emerging biotechnologies. Thus, revealing the largely hidden world of Acrs should provide important lessons in microbiology that have the potential to ripple far beyond the field.

Diverse viruses of marine archaea discovered using metagenomics.

During the past decade, metagenomics became a method of choice for the discovery of novel viruses. However, host assignment for uncultured viruses remains challenging, especially for archaeal viruses, which are grossly undersampled compared to viruses of bacteria and eukaryotes. Here, we assessed the utility of CRISPR spacer targeting, tRNA gene matching and homology searches for viral signature proteins, such as major capsid proteins, for the assignment of archaeal hosts and validated these approaches on metaviromes from Yangshan Harbor (YSH). We report 35 new genomes of viruses which could be confidently assigned to hosts representing diverse lineages of marine archaea. We show that the archaeal YSH virome is highly diverse, with some viruses enriching the previously described virus groups, such as magroviruses of Marine Group II Archaea (Poseidoniales), and others representing novel groups of marine archaeal viruses. Metagenomic recruitment of Tara Oceans datasets on the YSH viral genomes demonstrated the presence of YSH Poseidoniales and Nitrososphaeria viruses in the global oceans, but also revealed the endemic YSH-specific viral lineages. Furthermore, our results highlight the relationship between the soil and marine thaumarchaeal viruses. We propose three new families within the class Caudoviricetes for the classification of the five complete viral genomes predicted to replicate in marine Poseidoniales and Nitrososphaeria, two ecologically important and widespread archaeal groups. This study illustrates the utility of viral metagenomics in exploring the archaeal virome and provides new insights into the diversity, distribution and evolution of marine archaeal viruses.

Prokaryotic innate immunity through pattern recognition of conserved viral proteins.

Many organisms have evolved specialized immune pattern-recognition receptors, including nucleotide-binding oligomerization domain-like receptors (NLRs) of the STAND superfamily that are ubiquitous in plants, animals, and fungi. Although the roles of NLRs in eukaryotic immunity are well established, it is unknown whether prokaryotes use similar defense mechanisms. Here, we show that antiviral STAND (Avs) homologs in bacteria and archaea detect hallmark viral proteins, triggering Avs tetramerization and the activation of diverse N-terminal effector domains, including DNA endonucleases, to abrogate infection. Cryo-electron microscopy reveals that Avs sensor domains recognize conserved folds, active-site residues, and enzyme ligands, allowing a single Avs receptor to detect a wide variety of viruses. These findings extend the paradigm of pattern recognition of pathogen-specific proteins across all three domains of life.

PADLOC: a web server for the identification of antiviral defence systems in microbial genomes.

Most bacteria and archaea possess multiple antiviral defence systems that protect against infection by phages, archaeal viruses and mobile genetic elements. Our understanding of the diversity of defence systems has increased greatly in the last few years, and many more systems likely await discovery. To identify defence-related genes, we recently developed the Prokaryotic Antiviral Defence LOCator (PADLOC) bioinformatics tool. To increase the accessibility of PADLOC, we describe here the PADLOC web server (freely available at https://padloc.otago.ac.nz), allowing users to analyse whole genomes, metagenomic contigs, plasmids, phages and archaeal viruses. The web server includes a more than 5-fold increase in defence system types detected (since the first release) and expanded functionality enabling detection of CRISPR arrays and retron ncRNAs. Here, we provide user information such as input options, description of the multiple outputs, limitations and considerations for interpretation of the results, and guidance for subsequent analyses. The PADLOC web server also houses a precomputed database of the defence systems in > 230,000 RefSeq genomes. These data reveal two taxa, Campylobacterota and Spriochaetota, with unusual defence system diversity and abundance. Overall, the PADLOC web server provides a convenient and accessible resource for the detection of antiviral defence systems.

Characterization of Blf4, an Archaeal Lytic Virus Targeting a Member of the Methanomicrobiales.

Today, the number of known viruses infecting methanogenic archaea is limited. Here, we report on a novel lytic virus, designated Blf4, and its host strain Methanoculleus bourgensis E02.3, a methanogenic archaeon belonging to the Methanomicrobiales, both isolated from a commercial biogas plant in Germany. The virus consists of an icosahedral head 60 nm in diameter and a long non-contractile tail of 125 nm in length, which is consistent with the new isolate belonging to the Siphoviridae family. Electron microscopy revealed that Blf4 attaches to the vegetative cells of M. bourgensis E02.3 as well as to cellular appendages. Apart from M. bourgensis E02.3, none of the tested Methanoculleus strains were lysed by Blf4, indicating a narrow host range. The complete 37 kb dsDNA genome of Blf4 contains 63 open reading frames (ORFs), all organized in the same transcriptional direction. For most of the ORFs, potential functions were predicted. In addition, the genome of the host M. bourgensis E02.3 was sequenced and assembled, resulting in a 2.6 Mbp draft genome consisting of nine contigs. All genes required for a hydrogenotrophic lifestyle were predicted. A CRISPR/Cas system (type I-U) was identified with six spacers directed against Blf4, indicating that this defense system might not be very efficient in fending off invading Blf4 virus.

Identification and classification of antiviral defence systems in bacteria and archaea with PADLOC reveals new system types.

To provide protection against viral infection and limit the uptake of mobile genetic elements, bacteria and archaea have evolved many diverse defence systems. The discovery and application of CRISPR-Cas adaptive immune systems has spurred recent interest in the identification and classification of new types of defence systems. Many new defence systems have recently been reported but there is a lack of accessible tools available to identify homologs of these systems in different genomes. Here, we report the Prokaryotic Antiviral Defence LOCator (PADLOC), a flexible and scalable open-source tool for defence system identification. With PADLOC, defence system genes are identified using HMM-based homologue searches, followed by validation of system completeness using gene presence/absence and synteny criteria specified by customisable system classifications. We show that PADLOC identifies defence systems with high accuracy and sensitivity. Our modular approach to organising the HMMs and system classifications allows additional defence systems to be easily integrated into the PADLOC database. To demonstrate application of PADLOC to biological questions, we used PADLOC to identify six new subtypes of known defence systems and a putative novel defence system comprised of a helicase, methylase and ATPase. PADLOC is available as a standalone package (https://github.com/padlocbio/padloc) and as a webserver (https://padloc.otago.ac.nz).

Prokaryotic viruses impact functional microorganisms in nutrient removal and carbon cycle in wastewater treatment plants.

As one of the largest biotechnological applications, activated sludge (AS) systems in wastewater treatment plants (WWTPs) harbor enormous viruses, with 10-1,000-fold higher concentrations than in natural environments. However, the compositional variation and host-connections of AS viruses remain poorly explored. Here, we report a catalogue of ~50,000 prokaryotic viruses from six WWTPs, increasing the number of described viral species of AS by 23-fold, and showing the very high viral diversity which is largely unknown (98.4-99.6% of total viral contigs). Most viral genera are represented in more than one AS system with 53 identified across all. Viral infection widely spans 8 archaeal and 58 bacterial phyla, linking viruses with aerobic/anaerobic heterotrophs, and other functional microorganisms controlling nitrogen/phosphorous removal. Notably, Mycobacterium, notorious for causing AS foaming, is associated with 402 viral genera. Our findings expand the current AS virus catalogue and provide reference for the phage treatment to control undesired microorganisms in WWTPs.

Bacterial Viruses Subcommittee and Archaeal Viruses Subcommittee of the ICTV: update of taxonomy changes in 2021.

In this article, we - the Bacterial Viruses Subcommittee and the Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) - summarise the results of our activities for the period March 2020 - March 2021. We report the division of the former Bacterial and Archaeal Viruses Subcommittee in two separate Subcommittees, welcome new members, a new Subcommittee Chair and Vice Chair, and give an overview of the new taxa that were proposed in 2020, approved by the Executive Committee and ratified by vote in 2021. In particular, a new realm, three orders, 15 families, 31 subfamilies, 734 genera and 1845 species were newly created or redefined (moved/promoted).