Evolution of an archaeal virus nucleocapsid protein from the CRISPR-associated Cas4 nuclease

Diverse viral cas genes antagonize CRISPR immunity

Prokaryotic CRISPR-Cas immunity is subverted by anti-CRISPRs (Acrs), which inhibit Cas protein activities when expressed during the phage lytic cycle or from resident prophages or plasmids1. Acrs often bind to specific cognate Cas proteins, and hence inhibition is typically limited to a single CRISPR-Cas subtype2. Furthermore, although acr genes are frequently organized together in phage-associated gene clusters3, how such inhibitors initially evolve has remained unclear. Here we investigated the Acr content and inhibition specificity of diverse Listeria isolates, which naturally harbour four CRISPR-Cas systems (types I-B, II-A, II-C and VI-A). We observed widespread antagonism of CRISPR, which we traced to 11 previously unknown and 4 known acr gene families encoded by endogenous mobile elements. Among these were two Acrs that possess sequence homology to type I-B Cas proteins, one of which assembles into a defective interference complex. Surprisingly, an additional type I-B Cas homologue did not affect type I immunity, but instead inhibited the RNA-targeting type VI CRISPR system by means of CRISPR RNA (crRNA) degradation. By probing viral sequence databases, we detected abundant orphan cas genes located within putative anti-defence gene clusters. Among them, we verified the activity of a particularly broad-spectrum cas3 homologue that inhibits type I-B, II-A and VI-A CRISPR immunity. Our observations provide direct evidence of Acr evolution by cas gene co-option, and new genes with potential for broad-spectrum control of genome editing technologies.

Phage Genome Diversity in a Biogas-Producing Microbiome Analyzed by Illumina and Nanopore GridION Sequencing

The microbial biogas network is complex and intertwined, and therefore relatively stable in its overall functionality. However, if key functional groups of microorganisms are affected by biotic or abiotic factors, the entire efficacy may be impaired. Bacteriophages are hypothesized to alter the steering process of the microbial network. In this study, an enriched fraction of virus-like particles was extracted from a mesophilic biogas reactor and sequenced on the Illumina MiSeq and Nanopore GridION sequencing platforms. Metagenome data analysis resulted in identifying 375 metagenome-assembled viral genomes (MAVGs). Two-thirds of the classified sequences were only assigned to the superkingdom Viruses and the remaining third to the family Siphoviridae, followed by Myoviridae, Podoviridae, Tectiviridae, and Inoviridae. The metavirome showed a close relationship to the phage genomes that infect members of the classes Clostridia and Bacilli. Using publicly available biogas metagenomic data, a fragment recruitment approach showed the widespread distribution of the MAVGs studied in other biogas microbiomes. In particular, phage sequences from mesophilic microbiomes were highly similar to the phage sequences of this study. Accordingly, the virus particle enrichment approach and metavirome sequencing provided additional genome sequence information for novel virome members, thus expanding the current knowledge of viral genetic diversity in biogas reactors.

Evolution of a major virion protein of the giant pandoraviruses from an inactivated bacterial glycoside hydrolase

The diverse viruses in the phylum Nucleocytoviricota (also known as NLCDVs, Nucleo-cytoplasmic Large DNA Viruses) typically possess large icosahedral virions. However, in several families of Nucleocytoviricota, the icosahedral capsid was replaced by irregular particle shapes, most notably, the amphora-like virions of pandoraviruses and pithoviruses, the largest known virus particles in the entire virosphere. Pandoraviruses appear to be the most highly derived viruses in this phylum because their evolution involved not only the change in the virion shape, but also, the actual loss of the gene encoding double-jelly roll major capsid protein (DJR MCP), the main building block of icosahedral capsids in this virus assemblage. Instead, pandoravirus virions are built of unrelated abundant proteins. Here we show that the second most abundant virion protein of pandoraviruses, major virion protein 2 (MVP2), evolved from an inactivated derivative of a bacterial glycoside hydrolase of the GH16 family. The ancestral form of MVP2 was apparently acquired early in the evolution of the Nucleocytoviricota, to become a minor virion protein. After a duplication in the common ancestor of pandoraviruses and molliviruses, one of the paralogs displaces DJR MCP in pandoraviruses, conceivably, opening the way for a major increase in the size of the virion and the genome. Exaptation of a carbohydrate-binding protein for the function of the MVP is a general trend in virus evolution and might underlie the transformation of the virion shape in other groups of the Nucleocytoviricota as well.

CRISPR-Cas Systems and the Paradox of Self-Targeting Spacers

CRISPR-Cas immune systems in bacteria and archaea record prior infections as spacers within each system’s CRISPR arrays. Spacers are normally derived from invasive genetic material and direct the immune system to complementary targets as part of future infections. However, not all spacers appear to be derived from foreign genetic material and instead can originate from the host genome. Their presence poses a paradox, as self-targeting spacers would be expected to induce an autoimmune response and cell death. In this review, we discuss the known frequency of self-targeting spacers in natural CRISPR-Cas systems, how these spacers can be incorporated into CRISPR arrays, and how the host can evade lethal attack. We also discuss how self-targeting spacers can become the basis for alternative functions performed by CRISPR-Cas systems that extend beyond adaptive immunity. Overall, the acquisition of genome-targeting spacers poses a substantial risk but can aid in the host’s evolution and potentially lead to or support new functionalities.

Structural conservation in a membrane-enveloped filamentous virus infecting a hyperthermophilic acidophile

Different forms of viruses that infect archaea inhabiting extreme environments continue to be discovered at a surprising rate, suggesting that the current sampling of these viruses is sparse. We describe here Sulfolobus filamentous virus 1 (SFV1), a membrane-enveloped virus infecting Sulfolobus shibatae. The virus encodes two major coat proteins which display no apparent sequence similarity with each other or with any other proteins in databases. We have used cryo-electron microscopy at 3.7 Å resolution to show that these two proteins form a nearly symmetrical heterodimer, which wraps around A-form DNA, similar to what has been shown for SIRV2 and AFV1, two other archaeal filamentous viruses. The thin (∼ 20 Å) membrane of SFV1 is mainly archaeol, a lipid species that accounts for only 1% of the host lipids. Our results show how relatively conserved structural features can be maintained across evolution by both proteins and lipids that have diverged considerably.

Only a few archaeal filamentous viruses have been structurally characterized. Here the authors describe the membrane-enveloped Sulfolobus filamentous virus 1 that infects Sulfolobus shibatae and present its 3.7 Å resolution cryo-EM structure, which reveals that major coat proteins are structurally conserved among archaeal filamentous viruses.

Comparative genetic and genomic analysis of the novel fusellovirus Sulfolobus spindle-shaped virus 10

Viruses that infect thermophilic Archaea are unique in both their structure and genetic makeup. The lemon-shaped fuselloviruses—which infect members of the order Sulfolobales, growing optimally at 80 °C and pH 3—are some of the most ubiquitous and best studied viruses of the thermoacidophilic Archaea. Nonetheless, much remains to be learned about these viruses. In order to investigate fusellovirus evolution, we have isolated and characterized a novel fusellovirus, Sulfolobus spindle-shaped virus 10 (formerly SSV-L1). Comparative genomic analyses highlight significant similarity with both SSV8 and SSV9, as well as conservation of promoter elements within the Fuselloviridae. SSV10 encodes five ORFs with no homology within or outside of the Fuselloviridae, as well as a putatively functional Cas4-like ORF, which may play a role in evading CRISPR-mediated host defenses. Moreover, we demonstrate the ability of SSV10 to withstand mutation in a fashion consistent with mutagenesis in SSV1.

Using bioinformatic and phylogenetic approaches to classify transposable elements and understand their complex evolutionary histories

In recent years, much attention has been paid to comparative genomic studies of transposable elements (TEs) and the ensuing problems of their identification, classification, and annotation. Different approaches and diverse automated pipelines are being used to catalogue and categorize mobile genetic elements in the ever-increasing number of prokaryotic and eukaryotic genomes, with little or no connectivity between different domains of life. Here, an overview of the current picture of TE classification and evolutionary relationships is presented, updating the diversity of TE types uncovered in sequenced genomes. A tripartite TE classification scheme is proposed to account for their replicative, integrative, and structural components, and the need to expand in vitro and in vivo studies of their structural and biological properties is emphasized. Bioinformatic studies have now become front and center of novel TE discovery, and experimental pursuits of these discoveries hold great promise for both basic and applied science.

Phylogenomics of Cas4 family nucleases

Background

The Cas4 family endonuclease is a component of the adaptation module in many variants of CRISPR-Cas adaptive immunity systems. Unlike most of the other Cas proteins, Cas4 is often encoded outside CRISPR-cas loci (solo-Cas4) and is also found in mobile genetic elements (MGE-Cas4).

Results

As part of our ongoing investigation of CRISPR-Cas evolution, we explored the phylogenomics of the Cas4 family. About 90% of the archaeal genomes encode Cas4 compared to only about 20% of the bacterial genomes. Many archaea encode both the CRISPR-associated form (CAS-Cas4) and solo-Cas4, whereas in bacteria, this combination is extremely rare. The solo-cas4 genes are over-represented in environmental bacteria and archaea with small genomes that typically lack CRISPR-Cas, suggesting that Cas4 could perform uncharacterized defense or repair functions in these microbes. Phylogenomic analysis indicates that both the CRISPR-associated cas4 genes are often transferred horizontally but almost exclusively, as part of the adaptation module. The evolutionary integrity of the adaptation module sharply contrasts the rampant shuffling of CRISPR-cas modules whereby a given variant of the adaptation module can combine with virtually any effector module. The solo-cas4 genes evolve primarily via vertical inheritance and are subject only to occasional horizontal transfer. The selection pressure on cas4 genes does not substantially differ between CAS-Cas4 and solo-cas4, and is close to the genomic median. Thus, cas4 genes, similarly to cas1 and cas2, evolve similarly to ‘regular’ microbial genes involved in various cellular functions, showing no evidence of direct involvement in virus-host arms races. A notable feature of the Cas4 family evolution is the frequent recruitment of cas4 genes by various mobile genetic elements (MGE), particularly, archaeal viruses. The functions of Cas4 in these elements are unknown and potentially might involve anti-defense roles.

Conclusions

Unlike most of the other Cas proteins, Cas4 family members are as often encoded by stand-alone genes as they are incorporated in CRISPR-Cas systems. In addition, cas4 genes were repeatedly recruited by MGE, perhaps, for anti-defense functions. Experimental characterization of the solo and MGE-encoded Cas4 nucleases is expected to reveal currently uncharacterized defense and anti-defense systems and their interactions with CRISPR-Cas systems.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1081-1) contains supplementary material, which is available to authorized users.

Viruses of archaea: Structural, functional, environmental and evolutionary genomics

Viruses of archaea represent one of the most enigmatic parts of the virosphere. Most of the characterized archaeal viruses infect extremophilic hosts and display remarkable diversity of virion morphotypes, many of which have never been observed among viruses of bacteria or eukaryotes. The uniqueness of the virion morphologies is matched by the distinctiveness of the genomes of these viruses, with ∼75% of genes encoding unique proteins, refractory to functional annotation based on sequence analyses. In this review, we summarize the state-of-the-art knowledge on various aspects of archaeal virus genomics. First, we outline how structural and functional genomics efforts provided valuable insights into the functions of viral proteins and revealed intricate details of the archaeal virus-host interactions. We then highlight recent metagenomics studies, which provided a glimpse at the diversity of uncultivated viruses associated with the ubiquitous archaea in the oceans, including Thaumarchaeota, Marine Group II Euryarchaeota, and others. These findings, combined with the recent discovery that archaeal viruses mediate a rapid turnover of thaumarchaea in the deep sea ecosystems, illuminate the prominent role of these viruses in the biosphere. Finally, we discuss the origins and evolution of archaeal viruses and emphasize the evolutionary relationships between viruses and non-viral mobile genetic elements. Further exploration of the archaeal virus diversity as well as functional studies on diverse virus-host systems are bound to uncover novel, unexpected facets of the archaeal virome.

The enigmatic archaeal virosphere

One of the most prominent features of archaea is the extraordinary diversity of their DNA viruses. Many archaeal viruses differ substantially in morphology from bacterial and eukaryotic viruses and represent unique virus families. The distinct nature of archaeal viruses also extends to the gene composition and architectures of their genomes and the properties of the proteins that they encode. Environmental research has revealed prominent roles of archaeal viruses in influencing microbial communities in ocean ecosystems, and recent metagenomic studies have uncovered new groups of archaeal viruses that infect extremophiles and mesophiles in diverse habitats. In this Review, we summarize recent advances in our understanding of the genomic and morphological diversity of archaeal viruses and the molecular biology of their life cycles and virus-host interactions, including interactions with archaeal CRISPR-Cas systems. We also examine the potential origins and evolution of archaeal viruses and discuss their place in the global virosphere.

Mobile Genetic Elements and Evolution of CRISPR-Cas Systems: All the Way There and Back

The Clustered Regularly Interspaced Palindromic Repeats (CRISPR)-CRISPR-associated proteins (Cas) systems of bacterial and archaeal adaptive immunity show multifaceted evolutionary relationships with at least five classes of mobile genetic elements (MGE). First, the adaptation module of CRISPR-Cas that is responsible for the formation of the immune memory apparently evolved from a Casposon, a self-synthesizing transposon that employs the Cas1 protein as the integrase and might have brought additional cas genes to the emerging immunity loci. Second, a large subset of type III CRISPR-Cas systems recruited a reverse transcriptase from a Group II intron, providing for spacer acquisition from RNA. Third, effector nucleases of Class 2 CRISPR-Cas systems that are responsible for the recognition and cleavage of the target DNA were derived from transposon-encoded TnpB nucleases, most likely, on several independent occasions. Fourth, accessory nucleases in some variants of types I and III toxin and type VI effectors RNases appear to be ultimately derived from toxin nucleases of microbial toxin-antitoxin modules. Fifth, the opposite direction of evolution is manifested in the recruitment of CRISPR-Cas systems by a distinct family of Tn7-like transposons that probably exploit the capacity of CRISPR-Cas to recognize unique DNA sites to facilitate transposition as well as by bacteriophages that employ them to cope with host defense. Additionally, individual Cas proteins, such as the Cas4 nuclease, were recruited by bacteriophages and transposons. The two-sided evolutionary connection between CRISPR-Cas and MGE fits the "guns for hire" paradigm whereby homologous enzymatic machineries, in particular nucleases, are shuttled between MGE and defense systems and are used alternately as means of offense or defense.

Extreme Mutation Tolerance: Nearly Half of the Archaeal Fusellovirus Sulfolobus Spindle-Shaped Virus 1 Genes Are Not Required for Virus Function, Including the Minor Capsid Protein Gene vp3

Viruses infecting the Archaea harbor a tremendous amount of genetic diversity. This is especially true for the spindle-shaped viruses of the family Fuselloviridae, where >90% of the viral genes do not have detectable homologs in public databases. This significantly limits our ability to elucidate the role of viral proteins in the infection cycle. To address this, we have developed genetic techniques to study the well-characterized fusellovirus Sulfolobus spindle-shaped virus 1 (SSV1), which infects Sulfolobus solfataricus in volcanic hot springs at 80°C and pH 3. Here, we present a new comparative genome analysis and a thorough genetic analysis of SSV1 using both specific and random mutagenesis and thereby generate mutations in all open reading frames. We demonstrate that almost half of the SSV1 genes are not essential for infectivity, and the requirement for a particular gene correlates well with its degree of conservation within the Fuselloviridae The major capsid gene vp1 is essential for SSV1 infectivity. However, the universally conserved minor capsid gene vp3 could be deleted without a loss in infectivity and results in virions with abnormal morphology.IMPORTANCE Most of the putative genes in the spindle-shaped archaeal hyperthermophile fuselloviruses have no sequences that are clearly similar to characterized genes. In order to determine which of these SSV genes are important for function, we disrupted all of the putative genes in the prototypical fusellovirus, SSV1. Surprisingly, about half of the genes could be disrupted without destroying virus function. Even deletions of one of the known structural protein genes that is present in all known fuselloviruses, vp3, allows the production of infectious viruses. However, viruses lacking vp3 have abnormal shapes, indicating that the vp3 gene is important for virus structure. Identification of essential genes will allow focused research on minimal SSV genomes and further understanding of the structure of these unique, ubiquitous, and extremely stable archaeal viruses.

Multiple origins of viral capsid proteins from cellular ancestors

Viruses are the most abundant biological entities on earth and show remarkable diversity of genome sequences, replication and expression strategies, and virion structures. Evolutionary genomics of viruses revealed many unexpected connections but the general scenario(s) for the evolution of the virosphere remains a matter of intense debate among proponents of the cellular regression, escaped genes, and primordial virus world hypotheses. A comprehensive sequence and structure analysis of major virion proteins indicates that they evolved on about 20 independent occasions, and in some of these cases likely ancestors are identifiable among the proteins of cellular organisms. Virus genomes typically consist of distinct structural and replication modules that recombine frequently and can have different evolutionary trajectories. The present analysis suggests that, although the replication modules of at least some classes of viruses might descend from primordial selfish genetic elements, bona fide viruses evolved on multiple, independent occasions throughout the course of evolution by the recruitment of diverse host proteins that became major virion components.

Bipartite Network Analysis of the Archaeal Virosphere: Evolutionary Connections between Viruses and Capsidless Mobile Elements

Archaea and particularly hyperthermophilic crenarchaea are hosts to many unusual viruses with diverse virion shapes and distinct gene compositions. As is typical of viruses in general, there are no universal genes in the archaeal virosphere. Therefore, to obtain a comprehensive picture of the evolutionary relationships between viruses, network analysis methods are more productive than traditional phylogenetic approaches. Here we present a comprehensive comparative analysis of genomes and proteomes from all currently known taxonomically classified and unclassified, cultivated and uncultivated archaeal viruses. We constructed a bipartite network of archaeal viruses that includes two classes of nodes, the genomes and gene families that connect them. Dissection of this network using formal community detection methods reveals strong modularity, with 10 distinct modules and 3 putative supermodules. However, compared to similar previously analyzed networks of eukaryotic and bacterial viruses, the archaeal virus network is sparsely connected. With the exception of the tailed viruses related to bacteriophages of the order Caudovirales and the families Turriviridae and Sphaerolipoviridae that are linked to a distinct supermodule of eukaryotic and bacterial viruses, there are few connector genes shared by different archaeal virus modules. In contrast, most of these modules include, in addition to viruses, capsidless mobile elements, emphasizing tight evolutionary connections between the two types of entities in archaea. The relative contributions of distinct evolutionary origins, in particular from nonviral elements, and insufficient sampling to the sparsity of the archaeal virus network remain to be determined by further exploration of the archaeal virosphere.

IMPORTANCE

Viruses infecting archaea are among the most mysterious denizens of the virosphere. Many of these viruses display no genetic or even morphological relationship to viruses of bacteria and eukaryotes, raising questions regarding their origins and position in the global virosphere. Analysis of 5,740 protein sequences from 116 genomes allowed dissection of the archaeal virus network and showed that most groups of archaeal viruses are evolutionarily connected to capsidless mobile genetic elements, including various plasmids and transposons. This finding could reflect actual independent origins of the distinct groups of archaeal viruses from different nonviral elements, providing important insights into the emergence and evolution of the archaeal virome.

A virus of hyperthermophilic archaea with a unique architecture among DNA viruses

Viruses package their genetic material in diverse ways. Most known strategies include encapsulation of nucleic acids into spherical or filamentous virions with icosahedral or helical symmetry, respectively. Filamentous viruses with dsDNA genomes are currently associated exclusively with Archaea. Here, we describe a filamentous hyperthermophilic archaeal virus, Pyrobaculum filamentous virus 1 (PFV1), with a type of virion organization not previously observed in DNA viruses. The PFV1 virion, 400 ± 20 × 32 ± 3 nm, contains an envelope and an inner core consisting of two structural units: a rod-shaped helical nucleocapsid formed of two 14-kDa major virion proteins and a nucleocapsid-encompassing protein sheath composed of a single major virion protein of 18 kDa. The virion organization of PFV1 is superficially similar to that of negative-sense RNA viruses of the family Filoviridae, including Ebola virus and Marburg virus. The linear dsDNA of PFV1 carries 17,714 bp, including 60-bp-long terminal inverted repeats, and contains 39 predicted ORFs, most of which do not show similarities to sequences in public databases. PFV1 is a lytic virus that completely disrupts the host cell membrane at the end of the infection cycle.