On the Origin of Reverse Transcriptase-Using CRISPR-Cas Systems and Their Hyperdiverse, Enigmatic Spacer Repertoires

Unity among the diverse RNA-guided CRISPR-Cas interference mechanisms

CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated) systems are adaptive immune systems that protect bacteria and archaea from invading mobile genetic elements (MGEs). The Cas protein-CRISPR RNA (crRNA) complex uses complementarity of the crRNA "guide" region to specifically recognize the invader genome. CRISPR effectors that perform targeted destruction of the foreign genome have emerged independently as multi-subunit protein complexes (Class 1 systems) and as single multi-domain proteins (Class 2). These different CRISPR-Cas systems can cleave RNA, DNA, and protein in an RNA-guided manner to eliminate the invader, and in some cases, they initiate programmed cell death/dormancy. The versatile mechanisms of the different CRISPR-Cas systems to target and destroy nucleic acids have been adapted to develop various programmable-RNA-guided tools and have revolutionized the development of fast, accurate, and accessible genomic applications. In this review, we present the structure and interference mechanisms of different CRISPR-Cas systems and an analysis of their unified features. The three types of Class 1 systems (I, III, and IV) have a conserved right-handed helical filamentous structure that provides a backbone for sequence-specific targeting while using unique proteins with distinct mechanisms to destroy the invader. Similarly, all three Class 2 types (II, V, and VI) have a bilobed architecture that binds the RNA-DNA/RNA hybrid and uses different nuclease domains to cleave invading MGEs. Additionally, we highlight the mechanistic similarities of CRISPR-Cas enzymes with other RNA-cleaving enzymes and briefly present the evolutionary routes of the different CRISPR-Cas systems.

Mechanisms used for cDNA synthesis and site-specific integration of RNA into DNA genomes by a reverse transcriptase-Cas1 fusion protein

Reverse transcriptase-Cas1 (RT-Cas1) fusion proteins found in some CRISPR systems enable spacer acquisition from both RNA and DNA, but the mechanism of RNA spacer acquisition has remained unclear. Here, we found that Marinomonas mediterranea RT-Cas1/Cas2 adds short 3'-DNA (dN) tails to RNA protospacers, enabling their direct integration into CRISPR arrays as 3'-dN-RNAs or 3'-dN-RNA/cDNA duplexes at rates comparable to similarly configured DNAs. Reverse transcription of RNA protospacers is initiated at 3' proximal sites by multiple mechanisms, including recently described de novo initiation, protein priming with any dNTP, and use of short exogenous or synthesized DNA oligomer primers, enabling synthesis of near full-length cDNAs of diverse RNAs without fixed sequence requirements. The integration of 3'-dN-RNAs or single-stranded DNAs (ssDNAs) is favored over duplexes at higher protospacer concentrations, potentially relevant to spacer acquisition from abundant pathogen RNAs or ssDNA fragments generated by phage defense nucleases. Our findings reveal mechanisms for site-specifically integrating RNA into DNA genomes with potential biotechnological applications.

CRISPR-Cas systems of lactic acid bacteria and applications in food science

CRISPR-Cas (Clustered regularly interspaced short palindromic repeats-CRISPR associated proteins) systems are widely distributed in lactic acid bacteria (LAB), contributing to their RNA-mediated adaptive defense immunity. The CRISPR-Cas-based genetic tools have exhibited powerful capability. It has been highly utilized in different organisms, accelerating the development of life science. The review summarized the components, adaptive immunity mechanisms, and classification of CRISPR-Cas systems; analyzed the distribution and characteristics of CRISPR-Cas system in LAB. The review focuses on the development of CRISPR-Cas-based genetic tools in LAB for providing latest development and future trend. The diverse and broad applications of CRISPR-Cas systems in food/probiotic industry are introduced. LAB harbor a plenty of CRISPR-Cas systems, which contribute to generate safer and more robust strains with increased resistance against bacteriophage and prevent the dissemination of plasmids carrying antibiotic-resistance markers. Furthermore, the CRISPR-Cas system from LAB could be used to exploit novel, flexible, programmable genome editing tools of native host and other organisms, resolving the limitation of genetic operation of some LAB species, increasing the important biological functions of probiotics, improving the adaptation of probiotics in complex environments, and inhibiting the growth of foodborne pathogens. The development of the genetic tools based on CRISPR-Cas system in LAB, especially the endogenous CRISPR-Cas system, will open new avenues for precise regulation, rational design, and flexible application of LAB.

Characterization of Limnospira platensis PCC 9108 R-M and CRISPR-Cas systems

The filamentous cyanobacterium Limnospira platensis, formerly known as Arthrospira platensis or spirulina, is one of the most commercially important species of microalgae. Due to its high nutritional value, pharmacological and industrial applications it is extensively cultivated on a large commercial scale. Despite its widespread use, its precise manipulation is still under development due to the lack of effective genetic protocols. Genetic transformation of Limnospira has been attempted but the methods reported have not been generally reproducible in other laboratories. Knowledge of the transformation defense mechanisms is essential for understanding its physiology and for broadening their applications. With the aim to understand more about the genetic defenses of L. platensis, in this work we have identified the restriction-modification and CRISPR-Cas systems and we have cloned and characterized thirteen methylases. In parallel, we have also characterized the methylome and orphan methyltransferases using genome-wide analysis of DNA methylation patterns and RNA-seq. The identification and characterization of these enzymes will be a valuable resource to know how this strain avoids being genetically manipulated and for further genomics studies.

Harnessing CRISPR-Cas adaptation for RNA recording and beyond

Prokaryotes encode clustered regularly interspaced short palindromic repeat (CRISPR) arrays and CRISPR-associated (Cas) genes as an adaptive immune machinery. CRISPR-Cas systems effectively protect hosts from the invasion of foreign enemies, such as bacteriophages and plasmids. During a process called 'adaptation', non-self-nucleic acid fragments are acquired as spacers between repeats in the host CRISPR array, to establish immunological memory. The highly conserved Cas1-Cas2 complexes function as molecular recorders to integrate spacers in a time course manner, which can subsequently be expressed as crRNAs complexed with Cas effector proteins for the RNAguided interference pathways. In some of the RNA-targeting type III systems, Cas1 proteins are fused with reverse transcriptase (RT), indicating that RT-Cas1-Cas2 complexes can acquire RNA transcripts for spacer acquisition. In this review, we summarize current studies that focus on the molecular structure and function of the RT-fused Cas1-Cas2 integrase, and its potential applications as a directional RNA-recording tool in cells. Furthermore, we highlight outstanding questions for RT-Cas1-Cas2 studies and future directions for RNA-recording CRISPR technologies. [BMB Reports 2024; 57(1): 40-49].

Mechanisms used for cDNA synthesis and site-specific integration of RNA into DNA genomes by a reverse transcriptase-Cas1 fusion protein

Reverse transcriptase-Cas1 (RT-Cas1) fusion proteins found in some CRISPR systems enable spacer acquisition from both RNA and DNA, but the mechanism of RNA spacer acquisition has remained unclear. Here, we found Marinomonas mediterranea RT-Cas1/Cas2 adds short 3'-DNA (dN) tails to RNA protospacers enabling their direct integration into CRISPR arrays as 3'-dN-RNA/cDNA duplexes or 3'-dN-RNAs at rates comparable to similarly configured DNAs. Reverse transcription of RNA protospacers occurs by multiple mechanisms, including recently described de novo initiation, protein priming with any dNTP, and use of short exogenous or synthesized DNA oligomer primers, enabling synthesis of cDNAs from diverse RNAs without fixed sequence requirements. The integration of 3'-dN-RNAs or single-stranded (ss) DNAs is favored over duplexes at higher protospacer concentrations, potentially relevant to spacer acquisition from abundant pathogen RNAs or ssDNA fragments generated by phage-defense nucleases. Our findings reveal novel mechanisms for site-specifically integrating RNA into DNA genomes with potential biotechnological applications.

Reverse transcriptases prime DNA synthesis

The discovery of reverse transcriptases (RTs) challenged the central dogma by establishing that genetic information can also flow from RNA to DNA. Although they act as DNA polymerases, RTs are distantly related to replicases that also possess de novo primase activity. Here we identify that CRISPR associated RTs (CARTs) directly prime DNA synthesis on both RNA and DNA. We demonstrate that RT-dependent priming is utilized by some CRISPR-Cas complexes to synthesise new spacers and integrate these into CRISPR arrays. Expanding our analyses, we show that primer synthesis activity is conserved in representatives of other major RT classes, including group II intron RT, telomerase and retroviruses. Together, these findings establish a conserved innate ability of RTs to catalyse de novo DNA primer synthesis, independently of accessory domains or alternative priming mechanisms, which likely plays important roles in a wide variety of biological pathways.

Global phylogenomic novelty of the Cas1 gene from hot spring microbial communities

The Cas1 protein is essential for the functioning of CRISPR-Cas adaptive systems. However, despite the high prevalence of CRISPR-Cas systems in thermophilic microorganisms, few studies have investigated the occurrence and diversity of Cas1 across hot spring microbial communities. Phylogenomic analysis of 2,150 Cas1 sequences recovered from 48 metagenomes representing hot springs (42–80°C, pH 6–9) from three continents, revealed similar ecological diversity of Cas1 and 16S rRNA associated with geographic location. Furthermore, phylogenetic analysis of the Cas1 sequences exposed a broad taxonomic distribution in thermophilic bacteria, with new clades of Cas1 homologs branching at the root of the tree or at the root of known clades harboring reference Cas1 types. Additionally, a new family of casposases was identified from hot springs, which further completes the evolutionary landscape of the Cas1 superfamily. This ecological study contributes new Cas1 sequences from known and novel locations worldwide, mainly focusing on under-sampled hot spring microbial mat taxa. Results herein show that circumneutral hot springs are environments harboring high diversity and novelty related to adaptive immunity systems.

Structural biology of CRISPR-Cas immunity and genome editing enzymes

CRISPR-Cas systems provide resistance against foreign mobile genetic elements and have a wide range of genome editing and biotechnological applications. In this Review, we examine recent advances in understanding the molecular structures and mechanisms of enzymes comprising bacterial RNA-guided CRISPR-Cas immune systems and deployed for wide-ranging genome editing applications. We explore the adaptive and interference aspects of CRISPR-Cas function as well as open questions about the molecular mechanisms responsible for genome targeting. These structural insights reflect close evolutionary links between CRISPR-Cas systems and mobile genetic elements, including the origins and evolution of CRISPR-Cas systems from DNA transposons, retrotransposons and toxin-antitoxin modules. We discuss how the evolution and structural diversity of CRISPR-Cas systems explain their functional complexity and utility as genome editing tools.

Expansion of the global RNA virome reveals diverse clades of bacteriophages

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

Adaptation by Type III CRISPR-Cas Systems: Breakthrough Findings and Open Questions

CRISPR-Cas systems acquire heritable defense memory against invading nucleic acids through adaptation. Type III CRISPR-Cas systems have unique and intriguing features of defense and are important in method development for Genetics research. We started to understand the common and unique properties of type III CRISPR-Cas adaptation in recent years. This review summarizes our knowledge regarding CRISPR-Cas adaptation with the emphasis on type III systems and discusses open questions for type III adaptation studies.

Recruitment of Mobile Genetic Elements for Diverse Cellular Functions in Prokaryotes

Prokaryotic genomes are replete with mobile genetic elements (MGE) that span a continuum of replication autonomy. On numerous occasions during microbial evolution, diverse MGE lose their autonomy altogether but, rather than being quickly purged from the host genome, assume a new function that benefits the host, rendering the immobilized MGE subject to purifying selection, and resulting in its vertical inheritance. This mini-review highlights the diversity of the repurposed (exapted) MGE as well as the plethora of cellular functions that they perform. The principal contribution of the exaptation of MGE and their components is to the prokaryotic functional systems involved in biological conflicts, and in particular, defense against viruses and other MGE. This evolutionary entanglement between MGE and defense systems appears to stem both from mechanistic similarities and from similar evolutionary predicaments whereby both MGEs and defense systems tend to incur fitness costs to the hosts and thereby evolve mechanisms for survival including horizontal mobility, causing host addiction, and exaptation for functions beneficial to the host. The examples discussed demonstrate that the identity of an MGE, overall mobility and relationship with the host cell (mutualistic, symbiotic, commensal, or parasitic) are all factors that affect exaptation.

Different modes of spacer acquisition by the Staphylococcus epidermidis type III-A CRISPR-Cas system

CRISPR-Cas systems provide prokaryotic organisms with an adaptive defense mechanism that acquires immunological memories of infections. This is accomplished by integration of short fragments from the genome of invaders such as phages and plasmids, called ‘spacers’, into the CRISPR locus of the host. Depending on their genetic composition, CRISPR-Cas systems can be classified into six types, I-VI, however spacer acquisition has been extensively studied only in type I and II systems. Here, we used an inducible spacer acquisition assay to study this process in the type III-A CRISPR-Cas system of Staphylococcus epidermidis, in the absence of phage selection. Similarly to type I and II spacer acquisition, this type III system uses Cas1 and Cas2 to preferentially integrate spacers from the chromosomal terminus and free dsDNA ends produced after DNA breaks, in a manner that is enhanced by the AddAB DNA repair complex. Surprisingly, a different mode of spacer acquisition from rRNA and tRNA loci, which spans only the transcribed sequences of these genes and is not enhanced by AddAB, was also detected. Therefore, our findings reveal both common mechanistic principles that may be conserved in all CRISPR-Cas systems, as well as unique and intriguing features of type III spacer acquisition.

Evolutionary plasticity and functional versatility of CRISPR systems

The principal biological function of bacterial and archaeal CRISPR systems is RNA-guided adaptive immunity against viruses and other mobile genetic elements (MGEs). These systems show remarkable evolutionary plasticity and functional versatility at multiple levels, including both the defense mechanisms that lead to direct, specific elimination of the target DNA or RNA and those that cause programmed cell death (PCD) or induction of dormancy. This flexibility is also evident in the recruitment of CRISPR systems for nondefense functions. Defective CRISPR systems or individual CRISPR components have been recruited by transposons for RNA-guided transposition, by plasmids for interplasmid competition, and by viruses for antidefense and interviral conflicts. Additionally, multiple highly derived CRISPR variants of yet unknown functions have been discovered. A major route of innovation in CRISPR evolution is the repurposing of diverged repeat variants encoded outside CRISPR arrays for various structural and regulatory functions. The evolutionary plasticity and functional versatility of CRISPR systems are striking manifestations of the ubiquitous interplay between defense and “normal” cellular functions.

The CRISPR systems show remarkable functional versatility beyond their principal function as an adaptive immune mechanism. This Essay discusses how derived CRISPR systems have been recruited by transposons on multiple occasions and mediate RNA-guided transposition; derived CRISPR RNAs are frequently recruited for regulatory functions.

Identification and classification of reverse transcriptases in bacterial genomes and metagenomes

Reverse transcriptases (RTs) are found in different systems including group II introns, Diversity Generating Retroelements (DGRs), retrons, CRISPR-Cas systems, and Abortive Infection (Abi) systems in prokaryotes. Different classes of RTs can play different roles, such as template switching and mobility in group II introns, spacer acquisition in CRISPR-Cas systems, mutagenic retrohoming in DGRs, programmed cell suicide in Abi systems, and recently discovered phage defense in retrons. While some classes of RTs have been studied extensively, others remain to be characterized. There is a lack of computational tools for identifying and characterizing various classes of RTs. In this study, we built a tool (called myRT) for identification and classification of prokaryotic RTs. In addition, our tool provides information about the genomic neighborhood of each RT, providing potential functional clues. We applied our tool to predict RTs in all complete and draft bacterial genomes, and created a collection that can be used for exploration of putative RTs and their associated protein domains. Application of myRT to metagenomes showed that gut metagenomes encode proportionally more RTs related to DGRs, outnumbering retron-related RTs, as compared to the collection of reference genomes. MyRT is both available as a standalone software (https://github.com/mgtools/myRT) and also through a website (https://omics.informatics.indiana.edu/myRT/).

Prokaryotic reverse transcriptases: from retroelements to specialized defense systems

Reverse transcriptases (RTs) catalyze the polymerization of DNA from an RNA template. These enzymes were first discovered in RNA tumor viruses in 1970, but it was not until 1989 that they were found in prokaryotes as a key component of retrons. Apart from RTs encoded by the 'selfish' mobile retroelements known as group II introns, prokaryotic RTs are extraordinarily diverse, but their function has remained elusive. However, recent studies have revealed that different lineages of prokaryotic RTs, including retrons, those associated with CRISPR-Cas systems, Abi-like RTs and other yet uncharacterized RTs, are key components of different lines of defense against phages and other mobile genetic elements. Prokaryotic RTs participate in various antiviral strategies, including abortive infection (Abi), in which the infected cell is induced to commit suicide to protect the host population, adaptive immunity, in which a memory of previous infection is used to build an efficient defense, and other as yet unidentified mechanisms. These prokaryotic enzymes are attracting considerable attention, both for use in cutting-edge technologies, such as genome editing, and as an emerging research topic. In this review, we discuss what is known about prokaryotic RTs, and the exciting evidence for their domestication from retroelements to create specialized defense systems.

Characterization of 67 Confirmed Clustered Regularly Interspaced Short Palindromic Repeats Loci in 52 Strains of Staphylococci

Staphylococcus aureus (S. aureus), which is one of the most important species of Staphylococci, poses a great threat to public health. Clustered regularly interspaced short palindromic repeats (CRISPR) and their CRISPR-associated proteins (Cas) are an adaptive immune platform to combat foreign mobile genetic elements (MGEs) such as plasmids and phages. The aim of this study is to describe the distribution and structure of CRISPR-Cas system in S. aureus, and to explore the relationship between CRISPR and horizontal gene transfer (HGT). Here, we analyzed 67 confirmed CRISPR loci and 15 companion Cas proteins in 52 strains of Staphylococci with bioinformatics methods. Comparing with the orphan CRISPR loci in Staphylococci, the strains harboring complete CRISPR-Cas systems contained multiple CRISPR loci, direct repeat sequences (DR) forming stable RNA secondary structures with lower minimum free energy (MFE), and variable spacers with detectable protospacers. In S. aureus, unlike the orphan CRISPRs away from Staphylococcal cassette chromosome mec (SCCmec), the complete CRISPR-Cas systems were in J1 region of SCCmec. In addition, we found a conserved motif 5'-TTCTCGT-3' that may protect their downstream sequences from DNA interference. In general, orphan CRISPR locus in S. aureus differed greatly from the structural characteristics of the CRISPR-Cas system. Collectively, our results provided new insight into the diversity and characterization of the CRISPR-Cas system in S. aureus.

Type III CRISPR-Cas Systems: Deciphering the Most Complex Prokaryotic Immune System

The emergence and persistence of selfish genetic elements is an intrinsic feature of all living systems. Cellular organisms have evolved a plethora of elaborate defense systems that limit the spread of such genetic parasites. CRISPR-Cas are RNA-guided defense systems used by prokaryotes to recognize and destroy foreign nucleic acids. These systems acquire and store fragments of foreign nucleic acids and utilize the stored sequences as guides to recognize and destroy genetic invaders. CRISPR-Cas systems have been extensively studied, as some of them are used in various genome editing technologies. Although Type III CRISPR-Cas systems are among the most common CRISPR-Cas systems, they are also some of the least investigated ones, mostly due to the complexity of their action compared to other CRISPR-Cas system types. Type III effector complexes specifically recognize and cleave RNA molecules. The recognition of the target RNA activates the effector large subunit - the so-called CRISPR polymerase - which cleaves DNA and produces small cyclic oligonucleotides that act as signaling molecules to activate auxiliary effectors, notably non-specific RNases. In this review, we provide a historical overview of the sometimes meandering pathway of the Type III CRISPR research. We also review the current data on the structures and activities of Type III CRISPR-Cas systems components, their biological roles, and evolutionary history. Finally, using structural modeling with AlphaFold2, we show that the archaeal HRAMP signature protein, which heretofore has had no assigned function, is a degenerate relative of Type III CRISPR-Cas signature protein Cas10, suggesting that HRAMP systems have descended from Type III CRISPR-Cas systems or their ancestors.

CRISPR-Cas System: The Powerful Modulator of Accessory Genomes in Prokaryotes

Being frequently exposed to foreign nucleic acids, bacteria and archaea have developed an ingenious adaptive defense system, called CRISPR-Cas. The system is composed of the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) array, together with CRISPR (cas)-associated genes. This system consists of a complex machinery that integrates fragments of foreign nucleic acids from viruses and mobile genetic elements (MGEs), into CRISPR arrays. The inserted segments (spacers) are transcribed and then used by cas proteins as guide RNAs for recognition and inactivation of the targets. Different types and families of CRISPR-Cas systems consist of distinct adaptation and effector modules with evolutionary trajectories, partially independent. The origin of the effector modules and the mechanism of spacer integration/deletion is far less clear. A review of the most recent data regarding the structure, ecology, and evolution of CRISPR-Cas systems and their role in the modulation of accessory genomes in prokaryotes is proposed in this article. The CRISPR-Cas system's impact on the physiology and ecology of prokaryotes, modulation of horizontal gene transfer events, is also discussed here. This system gained popularity after it was proposed as a tool for plant and animal embryo editing, in cancer therapy, as antimicrobial against pathogenic bacteria, and even for combating the novel coronavirus - SARS-CoV-2; thus, the newest and promising applications are reviewed as well.

Structural coordination between active sites of a CRISPR reverse transcriptase-integrase complex

CRISPR-Cas systems provide adaptive immunity in bacteria and archaea, beginning with integration of foreign sequences into the host CRISPR genomic locus and followed by transcription and maturation of CRISPR RNAs (crRNAs). In some CRISPR systems, a reverse transcriptase (RT) fusion to the Cas1 integrase and Cas6 maturase creates a single protein that enables concerted sequence integration and crRNA production. To elucidate how the RT-integrase organizes distinct enzymatic activities, we present the cryo-EM structure of a Cas6-RT-Cas1-Cas2 CRISPR integrase complex. The structure reveals a heterohexamer in which the RT directly contacts the integrase and maturase domains, suggesting functional coordination between all three active sites. Together with biochemical experiments, our data support a model of sequential enzymatic activities that enable CRISPR sequence acquisition from RNA and DNA substrates. These findings highlight an expanded capacity of some CRISPR systems to acquire diverse sequences that direct CRISPR-mediated interference.

Spacer acquisition by Type III CRISPR-Cas system during bacteriophage infection of Thermus thermophilus

Type III CRISPR–Cas systems provide immunity to foreign DNA by targeting its transcripts. Target recognition activates RNases and DNases that may either destroy foreign DNA directly or elicit collateral damage inducing death of infected cells. While some Type III systems encode a reverse transcriptase to acquire spacers from foreign transcripts, most contain conventional spacer acquisition machinery found in DNA-targeting systems. We studied Type III spacer acquisition in phage-infected Thermus thermophilus, a bacterium that lacks either a standalone reverse transcriptase or its fusion to spacer integrase Cas1. Cells with spacers targeting a subset of phage transcripts survived the infection, indicating that Type III immunity does not operate through altruistic suicide. In the absence of selection spacers were acquired from both strands of phage DNA, indicating that no mechanism ensuring acquisition of RNA-targeting spacers exists. Spacers that protect the host from the phage demonstrate a very strong strand bias due to positive selection during infection. Phages that escaped Type III interference accumulated deletions of integral number of codons in an essential gene and much longer deletions in a non-essential gene. This and the fact that Type III immunity can be provided by plasmid-borne mini-arrays open ways for genomic manipulation of Thermus phages.

IMG/VR v.2.0: an integrated data management and analysis system for cultivated and environmental viral genomes

The Integrated Microbial Genome/Virus (IMG/VR) system v.2.0 (https://img.jgi.doe.gov/vr/) is the largest publicly available data management and analysis platform dedicated to viral genomics. Since the last report published in the 2016, NAR Database Issue, the data has tripled in size and currently contains genomes of 8389 cultivated reference viruses, 12 498 previously published curated prophages derived from cultivated microbial isolates, and 735 112 viral genomic fragments computationally predicted from assembled shotgun metagenomes. Nearly 60% of the viral genomes and genome fragments are clustered into 110 384 viral Operational Taxonomic Units (vOTUs) with two or more members. To improve data quality and predictions of host specificity, IMG/VR v.2.0 now separates prokaryotic and eukaryotic viruses, utilizes known prophage sequences to improve taxonomic assignments, and provides viral genome quality scores based on the estimated genome completeness. New features also include enhanced BLAST search capabilities for external queries. Finally, geographic map visualization to locate user-selected viral genomes or genome fragments has been implemented and download options have been extended. All of these features make IMG/VR v.2.0 a key resource for the study of viruses.

Real-time observation of CRISPR spacer acquisition by Cas1-Cas2 integrase

Cas1 integrase associates with Cas2 to insert short DNA fragments into a CRISPR array, establishing nucleic acid memory in prokaryotes. Here we applied single-molecule FRET methods to the Enterococcus faecalis (Efa) Cas1–Cas2 system to establish a kinetic framework describing target-searching, integration, and post-synapsis events. EfaCas1–Cas2 on its own is not able to find the CRISPR repeat in the CRISPR array; it only does so after prespacer loading. The leader sequence adjacent to the repeat further stabilizes EfaCas1–Cas2 contacts, enabling leader-side integration and subsequent spacer-side integration. The resulting post-synaptic complex has a surprisingly short mean lifetime. Remarkably, transcription efficiently resolves the postsynaptic complex and we predict that this is a conserved mechanism that ensures efficient and directional spacer integration in many CRISPR systems. Overall, our study provides a complete model of spacer acquisition, which can be harnessed for DNA-based information storage and cell lineage tracing technologies.

Origins and evolution of CRISPR-Cas systems

CRISPR-Cas, the bacterial and archaeal adaptive immunity systems, encompass a complex machinery that integrates fragments of foreign nucleic acids, mostly from mobile genetic elements (MGE), into CRISPR arrays embedded in microbial genomes. Transcripts of the inserted segments (spacers) are employed by CRISPR-Cas systems as guide (g)RNAs for recognition and inactivation of the cognate targets. The CRISPR-Cas systems consist of distinct adaptation and effector modules whose evolutionary trajectories appear to be at least partially independent. Comparative genome analysis reveals the origin of the adaptation module from casposons, a distinct type of transposons, which employ a homologue of Cas1 protein, the integrase responsible for the spacer incorporation into CRISPR arrays, as the transposase. The origin of the effector module(s) is far less clear. The CRISPR-Cas systems are partitioned into two classes, class 1 with multisubunit effectors, and class 2 in which the effector consists of a single, large protein. The class 2 effectors originate from nucleases encoded by different MGE, whereas the origin of the class 1 effector complexes remains murky. However, the recent discovery of a signalling pathway built into the type III systems of class 1 might offer a clue, suggesting that type III effector modules could have evolved from a signal transduction system involved in stress-induced programmed cell death. The subsequent evolution of the class 1 effector complexes through serial gene duplication and displacement, primarily of genes for proteins containing RNA recognition motif domains, can be hypothetically reconstructed. In addition to the multiple contributions of MGE to the evolution of CRISPR-Cas, the reverse flow of information is notable, namely, recruitment of minimalist variants of CRISPR-Cas systems by MGE for functions that remain to be elucidated. Here, we attempt a synthesis of the diverse threads that shed light on CRISPR-Cas origins and evolution. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Expansion of known ssRNA phage genomes: From tens to over a thousand

The first sequenced genome was that of the 3569-nucleotide single-stranded RNA (ssRNA) bacteriophage MS2. Despite the recent accumulation of vast amounts of DNA and RNA sequence data, only 12 representative ssRNA phage genome sequences are available from the NCBI Genome database (June 2019). The difficulty in detecting RNA phages in metagenomic datasets raises questions as to their abundance, taxonomic structure, and ecological importance. In this study, we iteratively applied profile hidden Markov models to detect conserved ssRNA phage proteins in 82 publicly available metatranscriptomic datasets generated from activated sludge and aquatic environments. We identified 15,611 nonredundant ssRNA phage sequences, including 1015 near-complete genomes. This expansion in the number of known sequences enabled us to complete a phylogenetic assessment of both sequences identified in this study and known ssRNA phage genomes. Our expansion of these viruses from two environments suggests that they have been overlooked within microbiome studies.

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.

Molecular mechanisms of CRISPR-Cas spacer acquisition

Many bacteria and archaea have the unique ability to heritably alter their genomes by incorporating small fragments of foreign DNA, called spacers, into CRISPR loci. Once transcribed and processed into individual CRISPR RNAs, spacer sequences guide Cas effector nucleases to destroy complementary, invading nucleic acids. Collectively, these two processes are known as the CRISPR-Cas immune response. In this Progress article, we review recent studies that have advanced our understanding of the molecular mechanisms underlying spacer acquisition and that have revealed a fundamental link between the two phases of CRISPR immunity that ensures optimal immunity from newly acquired spacers. Finally, we highlight important open questions and discuss the potential basic and applied impact of spacer acquisition research.

Targeted transposition with Tn7 elements: safe sites, mobile plasmids, CRISPR/Cas and beyond

Transposon Tn7 is notable for the control it exercises over where transposition events are directed. One Tn7 integration pathways recognizes a highly conserved attachment (att) site in the chromosome, while a second pathway specifically recognizes mobile plasmids that facilitate transfer of the element to new hosts. In this review, I discuss newly discovered families of Tn7-like elements with different targeting pathways. Perhaps the most exciting examples are multiple instances where Tn7-like elements have repurposed CRISPR/Cas systems. In these cases, the CRISPR/Cas systems have lost their canonical defensive function to destroy incoming mobile elements; instead, the systems have been naturally adapted to use guide RNAs to specifically direct transposition into these mobile elements. The new families of Tn7-like elements also include a variety of novel att sites in bacterial chromosomes where genome islands can form. Interesting families have also been revealed where proteins described in the prototypic Tn7 element are fused or otherwise repurposed for the new dual activities. This expanded understanding of Tn7-like elements broadens our view of how genetic systems are repurposed and provides potentially exciting new tools for genome modification and genomics. Future opportunities and challenges to understanding the impact of the new families of Tn7-like elements are discussed.

Multiple origins of reverse transcriptases linked to CRISPR-Cas systems

Prokaryotic genomes harbour a plethora of uncharacterized reverse transcriptases (RTs). RTs phylogenetically related to those encoded by group-II introns have been found associated with type III CRISPR-Cas systems, adjacent or fused at the C-terminus to Cas1. It is thought that these RTs may have a relevant function in the CRISPR immune response mediating spacer acquisition from RNA molecules. The origin and relationships of these RTs and the ways in which the various protein domains evolved remain matters of debate. We carried out a large survey of annotated RTs in databases (198,760 sequences) and constructed a large dataset of unique representative sequences (9,141). The combined phylogenetic reconstruction and identification of the RTs and their various protein domains in the vicinity of CRISPR adaptation and effector modules revealed three different origins for these RTs, consistent with their emergence on multiple occasions: a larger group that have evolved from group-II intron RTs, and two minor lineages that may have arisen more recently from Retron/retron-like sequences and Abi-P2 RTs, the latter associated with type I-C systems. We also identified a particular group of RTs associated with CRISPR-cas loci in clade 12, fused C-terminally to an archaeo-eukaryotic primase (AEP), a protein domain (AE-Prim_S_like) forming a particular family within the AEP proper clade. Together, these data provide new insight into the evolution of CRISPR-Cas/RT systems.

CRISPR-Cas systems are present predominantly on mobile genetic elements in Vibrio species

Background

Bacteria are prey for many viruses that hijack the bacterial cell in order to propagate, which can result in bacterial cell lysis and death. Bacteria have developed diverse strategies to counteract virus predation, one of which is the clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR associated (Cas) proteins immune defense system. Species within the bacterial family Vibrionaceae are marine organisms that encounter large numbers of phages. Our goal was to determine the significance of CRISPR-Cas systems as a mechanism of defense in this group by investigating their prevalence, phylogenetic distribution, and genome context.

Results

Herein, we describe all the CRISPR-Cas system types and their distribution within the family Vibrionaceae. In Vibrio cholerae genomes, we identified multiple variant type I-F systems, which were also present in 41 additional species. In a large number of Vibrio species, we identified a mini type I-F system comprised of tniQcas5cas7cas6f, which was always associated with Tn7-like transposons. The Tn7-like elements, in addition to the CRISPR-Cas system, also contained additional cargo genes such as restriction modification systems and type three secretion systems. A putative hybrid CRISPR-Cas system was identified containing type III-B genes followed by a type I-F cas6f and a type I-F CRISPR that was associated with a prophage in V. cholerae and V. metoecus strains. Our analysis identified CRISPR-Cas types I-C, I-E, I-F, II-B, III-A, III-B, III-D, and the rare type IV systems as well as cas loci architectural variants among 70 species. All systems described contained a CRISPR array that ranged in size from 3 to 179 spacers. The systems identified were present predominantly within mobile genetic elements (MGEs) such as genomic islands, plasmids, and transposon-like elements. Phylogenetic analysis of Cas proteins indicated that the CRISPR-Cas systems were acquired by horizontal gene transfer.

Conclusions

Our data show that CRISPR-Cas systems are phylogenetically widespread but sporadic in occurrence, actively evolving, and present on MGEs within Vibrionaceae.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5439-1) contains supplementary material, which is available to authorized users.

A Reverse Transcriptase-Cas1 Fusion Protein Contains a Cas6 Domain Required for Both CRISPR RNA Biogenesis and RNA Spacer Acquisition

Prokaryotic CRISPR-Cas systems provide adaptive immunity by integrating portions of foreign nucleic acids (spacers) into genomic CRISPR arrays. Cas6 proteins then process CRISPR array transcripts into spacer-derived RNAs (CRISPR RNAs; crRNAs) that target Cas nucleases to matching invaders. We find that a Marinomonas mediterranea fusion protein combines three enzymatic domains (Cas6, reverse transcriptase [RT], and Cas1), which function in both crRNA biogenesis and spacer acquisition from RNA and DNA. We report a crystal structure of this divergent Cas6, identify amino acids required for Cas6 activity, show that the Cas6 domain is required for RT activity and RNA spacer acquisition, and demonstrate that CRISPR-repeat binding to Cas6 regulates RT activity. Co-evolution of putative interacting surfaces suggests a specific structural interaction between the Cas6 and RT domains, and phylogenetic analysis reveals repeated, stable association of free-standing Cas6s with CRISPR RTs in multiple microbial lineages, indicating that a functional interaction between these proteins preceded evolution of the fusion.

Classification and Nomenclature of CRISPR-Cas Systems: Where from Here?

As befits an immune mechanism, CRISPR-Cas systems are highly variable with respect to Cas protein sequences, gene composition, and organization of the genomic loci. Optimal classification of CRISPR-Cas systems and rational nomenclature for CRISPR-associated genes are essential for further progress of CRISPR research. These are highly challenging tasks because of the complexity of CRISPR-Cas and their fast evolution, including frequent module shuffling, as well as the lack of universal markers for a consistent evolutionary classification. The complexity and variability of CRISPR-Cas systems necessitate a multipronged approach to classification and nomenclature. We present a brief summary of the current state of the art and discuss further directions in this area.

CRISPR-Cas systems in multicellular cyanobacteria

Novel CRISPR-Cas systems possess substantial potential for genome editing and manipulation of gene expression. The types and numbers of CRISPR-Cas systems vary substantially between different organisms. Some filamentous cyanobacteria harbor > 40 different putative CRISPR repeat-spacer cassettes, while the number of cas gene instances is much lower. Here we addressed the types and diversity of CRISPR-Cas systems and of CRISPR-like repeat-spacer arrays in 171 publicly available genomes of multicellular cyanobacteria. The number of 1328 repeat-spacer arrays exceeded the total of 391 encoded Cas1 proteins suggesting a tendency for fragmentation or the involvement of alternative adaptation factors. The model cyanobacterium Anabaena sp. PCC 7120 contains only three cas1 genes but hosts three Class 1, possibly one Class 2 and five orphan repeat-spacer arrays, all of which exhibit crRNA-typical expression patterns suggesting active transcription, maturation and incorporation into CRISPR complexes. The CRISPR-Cas system within the element interrupting the Anabaena sp. PCC 7120 fdxN gene, as well as analogous arrangements in other strains, occupy the genetic elements that become excised during the differentiation-related programmed site-specific recombination. This fact indicates the propensity of these elements for the integration of CRISPR-cas systems and points to a previously not recognized connection. The gene all3613 resembling a possible Class 2 effector protein is linked to a short repeat-spacer array and a single tRNA gene, similar to its homologs in other cyanobacteria. The diversity and presence of numerous CRISPR-Cas systems in DNA elements that are programmed for homologous recombination make filamentous cyanobacteria a prolific resource for their study. Abbreviations: Cas: CRISPR associated sequences; CRISPR: Clustered Regularly Interspaced Short Palindromic Repeats; C2c: Class 2 candidate; SDR: small dispersed repeat; TSS: transcriptional start site; UTR: untranslated region.

On the Origin and Evolutionary Relationships of the Reverse Transcriptases Associated With Type III CRISPR-Cas Systems

Reverse transcriptases (RTs) closely related to those encoded by group II introns but lacking the intron RNA structure have been found associated with type III clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems, a prokaryotic immune system against invading viruses and foreign genetic elements. Two models have been proposed to explain the origin and evolutionary relationships of these RTs: (i) the “single point of origin” model, according to which these RTs originated from a single acquisition event in bacterial, with the various protein domains (RT, RT-Cas1, and Cas6-RT-Cas1 fusions) corresponding to single points in evolution; and (ii) the “various origins” model, according to which, independent acquisition events in different evolutionary episodes led to these fusions. We tested these alternative hypotheses, by analyzing and integrating published datasets of RT sequences associated with CRISPR-Cas systems and inferring phylogenetic trees by maximum likelihood (ML) methods. The RTs studied could be grouped into 13 clades, mostly in bacteria, in which they probably evolved. The various clades appear to form three independent lineages in bacteria and a recent lineage in archaea. Our data show that the Cas6 domain was acquired twice, independently, through RT-Cas1 fusion, in the bacterial lineages. Taken together, there more evidence to support the “various origins” hypothesis.

Systematic prediction of genes functionally linked to CRISPR-Cas systems by gene neighborhood analysis

The CRISPR-Cas systems of bacterial and archaeal adaptive immunity consist of direct repeat arrays separated by unique spacers and multiple CRISPR-associated (cas) genes encoding proteins that mediate all stages of the CRISPR response. In addition to the relatively small set of core cas genes that are typically present in all CRISPR-Cas systems of a given (sub)type and are essential for the defense function, numerous genes occur in CRISPR-cas loci only sporadically. Some of these have been shown to perform various ancillary roles in CRISPR response, but the functional relevance of most remains unknown. We developed a computational strategy for systematically detecting genes that are likely to be functionally linked to CRISPR-Cas. The approach is based on a "CRISPRicity" metric that measures the strength of CRISPR association for all protein-coding genes from sequenced bacterial and archaeal genomes. Uncharacterized genes with CRISPRicity values comparable to those of cas genes are considered candidate CRISPR-linked genes. We describe additional criteria to predict functionally relevance for genes in the candidate set and identify 79 genes as strong candidates for functional association with CRISPR-Cas systems. A substantial majority of these CRISPR-linked genes reside in type III CRISPR-cas loci, which implies exceptional functional versatility of type III systems. Numerous candidate CRISPR-linked genes encode integral membrane proteins suggestive of tight membrane association of CRISPR-Cas systems, whereas many others encode proteins implicated in various signal transduction pathways. These predictions provide ample material for improving annotation of CRISPR-cas loci and experimental characterization of previously unsuspected aspects of CRISPR-Cas system functionality.

A Small RNA Isolation and Sequencing Protocol and Its Application to Assay CRISPR RNA Biogenesis in Bacteria

Next generation high-throughput sequencing has enabled sensitive and unambiguous analysis of RNA populations in cells. Here, we describe a method for isolation and strand-specific sequencing of small RNA pools from bacteria that can be multiplexed to accommodate multiple biological samples in a single experiment. Small RNAs are isolated by polyacrylamide gel electrophoresis and treated with T4 polynucleotide kinase. This allows for 3' adapter ligation to CRISPR RNAs, which don't have pre-existing 3'-OH ends. Pre-adenylated adapters are then ligated using T4 RNA ligase 1 in the absence of ATP and with a high concentration of polyethylene glycol (PEG). The 3' capture step enables precise determination of the 3' ends of diverse RNA molecules. Additionally, a random hexamer in the ligated adapter helps control for potential downstream amplification bias. Following reverse-transcription, the cDNA product is circularized and libraries are prepared by PCR. We show that the amplified library need not be visible by gel electrophoresis for efficient sequencing of the desired product. Using this method, we routinely prepare RNA sequencing libraries from minute amounts of purified small RNA. This protocol is tailored to assay for CRISPR RNA biogenesis in bacteria through sequencing of mature CRISPR RNAs, but can be used to sequence diverse classes of small RNAs. We also provide a fully worked example of our data processing pipeline, with instructions for running the provided scripts.

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.

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.