Optimal number of spacers in CRISPR arrays

Assessing spacer acquisition rates in E. coli type I-E CRISPR arrays

CRISPR/Cas is an adaptive defense mechanism protecting prokaryotes from viruses and other potentially harmful genetic elements. Through an adaptation process, short "spacer" sequences, captured from these elements and incorporated into a CRISPR array, provide target specificity for the immune response. CRISPR arrays and array expansion are also central to many emerging biotechnologies. The rates at which spacers integrate into native arrays within bacterial populations have not been quantified. Here, we measure naïve spacer acquisition rates in Escherichia coli Type I-E CRISPR, identify factors that affect these rates, and model this process fundamental to CRISPR/Cas defense. Prolonged Cas1-Cas2 expression produced fewer new spacers per cell on average than predicted by the model. Subsequent experiments revealed that this was due to a mean fitness reduction linked to array-expanded populations. In addition, the expression of heterologous non-homologous end-joining DNA-repair genes was found to augment spacer acquisition rates, translating to enhanced phage infection defense. Together, these results demonstrate the impact of intracellular factors that modulate spacer acquisition and identify an intrinsic fitness effect associated with array-expanded populations.

Ecological drivers of CRISPR immune systems

CRISPR-Cas is the only known adaptive immune system of prokaryotes. It is a powerful defense system against mobile genetic elements such as bacteriophages. While CRISPR-Cas systems can be found throughout the prokaryotic tree of life, they are distributed unevenly across taxa and environments. Since adaptive immunity is more useful in environments where pathogens persist or reoccur, the density and/or diversity of the host/pathogen community may drive the uneven distribution of CRISPR systems. We directly tested hypotheses connecting CRISPR incidence with prokaryotic density/diversity by analyzing 16S rRNA and metagenomic data from publicly available environmental sequencing projects. In terms of density, we found that CRISPR systems are significantly favored in lower abundance (less dense) taxa and disfavored in higher abundance taxa, at least in marine environments. When we extended this work to compare taxonomic diversity between samples, we found CRISPR system incidence strongly correlated with diversity in human oral environments. Together, these observations confirm that, at least in certain types of environments, the prokaryotic ecological context indeed plays a key role in selecting for CRISPR immunity.

IMPORTANCE

Microbes must constantly defend themselves against viral pathogens, and a large proportion of prokaryotes do so using the highly effective CRISPR-Cas adaptive immune system. However, many prokaryotes do not. We investigated the ecological factors behind this uneven distribution of CRISPR-Cas immune systems in natural microbial populations. We found strong patterns linking CRISPR-Cas systems to prokaryotic density within ocean environments and to prokaryotic diversity within human oral environments. Our study validates previous within-lab experimental results that suggested these factors might be important and confirms that local environment and ecological context interact to select for CRISPR immunity.

SpacerPlacer: ancestral reconstruction of CRISPR arrays reveals the evolutionary dynamics of spacer deletions

Bacteria employ CRISPR-Cas systems for defense by integrating invader-derived sequences, termed spacers, into the CRISPR array, which constitutes an immunity memory. While spacer deletions occur randomly across the array, newly acquired spacers are predominantly integrated at the leader end. Consequently, spacer arrays can be used to derive the chronology of spacer insertions. Reconstruction of ancestral spacer acquisitions and deletions could help unravel the coevolution of phages and bacteria, the evolutionary dynamics in microbiomes, or track pathogens. However, standard reconstruction methods produce misleading results by overlooking insertion order and joint deletions of spacers. Here, we present SpacerPlacer, a maximum likelihood-based ancestral reconstruction approach for CRISPR array evolution. We used SpacerPlacer to reconstruct and investigate ancestral deletion events of 4565 CRISPR arrays, revealing that spacer deletions occur 374 times more frequently than mutations and are regularly deleted jointly, with an average of 2.7 spacers. Surprisingly, we observed a decrease in the spacer deletion frequency towards both ends of the reconstructed arrays. While the resulting trailer-end conservation is commonly observed, a reduced deletion frequency is now also detectable towards the variable leader end. Finally, our results point to the hypothesis that frequent loss of recently acquired spacers may provide a selective advantage.

CRISPR Tools for Engineering Prokaryotic Systems: Recent Advances and New Applications

In the past decades, the broad selection of CRISPR-Cas systems has revolutionized biotechnology by enabling multimodal genetic manipulation in diverse organisms. Rooted in a molecular engineering perspective, we recapitulate the different CRISPR components and how they can be designed for specific genetic engineering applications. We first introduce the repertoire of Cas proteins and tethered effectors used to program new biological functions through gene editing and gene regulation. We review current guide RNA (gRNA) design strategies and computational tools and how CRISPR-based genetic circuits can be constructed through regulated gRNA expression. Then, we present recent advances in CRISPR-based biosensing, bioproduction, and biotherapeutics across in vitro and in vivo prokaryotic systems. Finally, we discuss forthcoming applications in prokaryotic CRISPR technology that will transform synthetic biology principles in the near future.

Ecological drivers of CRISPR immune systems

CRISPR-Cas is the only known adaptive immune system of prokaryotes. It is a powerful defense system against mobile genetic elements such as bacteriophages. While CRISPR-Cas systems can be found throughout the prokaryotic tree of life, they are distributed unevenly across taxa and environments. Since adaptive immunity is more useful in environments where pathogens persist or reoccur, the density and/or diversity of the host/pathogen community may drive the uneven distribution of CRISPR system. We directly tested hypotheses connecting CRISPR incidence with prokaryotic density/diversity by analyzing 16S rRNA and metagenomic data from publicly available environmental sequencing projects. In terms of density, we found that CRISPR systems are significantly favored in lower abundance (less dense) taxa and disfavored in higher abundance taxa, at least in marine environments. When we extended this work to compare taxonomic diversity between samples, we found CRISPR system incidence strongly correlated with diversity in human oral environments. Together, these observations confirm that, at least in certain types of environments, the prokaryotic ecological context indeed plays a key role in selecting for CRISPR immunity.

Importance

2Microbes must constantly defend themselves against viral pathogens, and a large proportion of prokaryotes do so using the highly effective CRISPR-Cas adaptive immune system. However, many prokaryotes do not. We investigated the ecological factors behind this uneven distribution of CRISPR-Cas immune systems in natural microbial populations. We found strong patterns linking CRISPR-Cas systems to prokaryotic density within ocean environments and to prokaryotic diversity within human oral environments. Our study validates previous within-lab experimental results that suggested these factors might be important and confirms that local environment and ecological context interact to select for CRISPR immunity.

Insights into the inhibition of protospacer integration via direct interaction between Cas2 and AcrVA5

Spacer acquisition step in CRISPR-Cas system involves the recognition and subsequent integration of protospacer by the Cas1-Cas2 complex in CRISPR-Cas systems. Here we report an anti-CRISPR protein, AcrVA5, and reveal the mechanisms by which it strongly inhibits protospacer integration. Our biochemical data shows that the integration by Cas1-Cas2 was abrogated in the presence of AcrVA5. AcrVA5 exhibits low binding affinity towards Cas2 and acetylates Cas2 at Lys55 on the binding interface of the Cas2 and AcrVA5 N-terminal peptide complex to inhibit the Cas2-mediated endonuclease activity. Moreover, a detailed structural comparison between our crystal structure and homolog structure shows that binding of AcrVA5 to Cas2 causes steric hindrance to the neighboring protospacer resulting in the partial disassembly of the Cas1-Cas2 and protospacer complex, as demonstrated by electrophoretic mobility shift assay. Our study focuses on this mechanism of spacer acquisition inhibition and provides insights into the biology of CRISPR-Cas systems.

Cbp1 and Cren7 form chromatin-like structures that ensure efficient transcription of long CRISPR arrays

CRISPR arrays form the physical memory of CRISPR adaptive immune systems by incorporating foreign DNA as spacers that are often AT-rich and derived from viruses. As promoter elements such as the TATA-box are AT-rich, CRISPR arrays are prone to harbouring cryptic promoters. Sulfolobales harbour extremely long CRISPR arrays spanning several kilobases, a feature that is accompanied by the CRISPR-specific transcription factor Cbp1. Aberrant Cbp1 expression modulates CRISPR array transcription, but the molecular mechanisms underlying this regulation are unknown. Here, we characterise the genome-wide Cbp1 binding at nucleotide resolution and characterise the binding motifs on distinct CRISPR arrays, as well as on unexpected non-canonical binding sites associated with transposons. Cbp1 recruits Cren7 forming together 'chimeric' chromatin-like structures at CRISPR arrays. We dissect Cbp1 function in vitro and in vivo and show that the third helix-turn-helix domain is responsible for Cren7 recruitment, and that Cbp1-Cren7 chromatinization plays a dual role in the transcription of CRISPR arrays. It suppresses spurious transcription from cryptic promoters within CRISPR arrays but enhances CRISPR RNA transcription directed from their cognate promoters in their leader region. Our results show that Cbp1-Cren7 chromatinization drives the productive expression of long CRISPR arrays.

Widespread CRISPR-derived RNA regulatory elements in CRISPR-Cas systems

CRISPR-cas loci typically contain CRISPR arrays with unique spacers separating direct repeats. Spacers along with portions of adjacent repeats are transcribed and processed into CRISPR(cr) RNAs that target complementary sequences (protospacers) in mobile genetic elements, resulting in cleavage of the target DNA or RNA. Additional, standalone repeats in some CRISPR-cas loci produce distinct cr-like RNAs implicated in regulatory or other functions. We developed a computational pipeline to systematically predict crRNA-like elements by scanning for standalone repeat sequences that are conserved in closely related CRISPR-cas loci. Numerous crRNA-like elements were detected in diverse CRISPR-Cas systems, mostly, of type I, but also subtype V-A. Standalone repeats often form mini-arrays containing two repeat-like sequence separated by a spacer that is partially complementary to promoter regions of cas genes, in particular cas8, or cargo genes located within CRISPR-Cas loci, such as toxins-antitoxins. We show experimentally that a mini-array from a type I-F1 CRISPR-Cas system functions as a regulatory guide. We also identified mini-arrays in bacteriophages that could abrogate CRISPR immunity by inhibiting effector expression. Thus, recruitment of CRISPR effectors for regulatory functions via spacers with partial complementarity to the target is a common feature of diverse CRISPR-Cas systems.

Long-term CRISPR locus dynamics and stable host-virus co-existence in subsurface fractured shales

Viruses are the most ubiquitous biological entities on Earth. Even so, elucidating the impact of viruses on microbial communities and associated ecosystem processes often requires identification of unambiguous host-virus linkages-an undeniable challenge in many ecosystems. Subsurface fractured shales present a unique opportunity to first make these strong linkages via spacers in CRISPR-Cas arrays and subsequently reveal complex long-term host-virus dynamics. Here, we sampled two replicated sets of fractured shale wells for nearly 800 days, resulting in 78 metagenomes from temporal sampling of six wells in the Denver-Julesburg Basin (Colorado, USA). At the community level, there was strong evidence for CRISPR-Cas defense systems being used through time and likely in response to viral interactions. Within our host genomes, represented by 202 unique MAGs, we also saw that CRISPR-Cas systems were widely encoded. Together, spacers from host CRISPR loci facilitated 2,110 CRISPR-based viral linkages across 90 host MAGs spanning 25 phyla. We observed less redundancy in host-viral linkages and fewer spacers associated with hosts from the older, more established wells, possibly reflecting enrichment of more beneficial spacers through time. Leveraging temporal patterns of host-virus linkages across differing well ages, we report how host-virus co-existence dynamics develop and converge through time, possibly reflecting selection for viruses that can evade host CRISPR-Cas systems. Together, our findings shed light on the complexities of host-virus interactions as well as long-term dynamics of CRISPR-Cas defense among diverse microbial populations.

Widespread CRISPR repeat-like RNA regulatory elements in CRISPR-Cas systems

CRISPR- cas loci typically contain CRISPR arrays with unique spacers separating direct repeats. Spacers along with portions of adjacent repeats are transcribed and processed into CRISPR(cr) RNAs that target complementary sequences (protospacers) in mobile genetic elements, resulting in cleavage of the target DNA or RNA. Additional, standalone repeats in some CRISPR- cas loci produce distinct cr-like RNAs implicated in regulatory or other functions. We developed a computational pipeline to systematically predict crRNA-like elements by scanning for standalone repeat sequences that are conserved in closely related CRISPR- cas loci. Numerous crRNA-like elements were detected in diverse CRISPR-Cas systems, mostly, of type I, but also subtype V-A. Standalone repeats often form mini-arrays containing two repeat-like sequence separated by a spacer that is partially complementary to promoter regions of cas genes, in particular cas8 , or cargo genes located within CRISPR-Cas loci, such as toxins-antitoxins. We show experimentally that a mini-array from a type I-F1 CRISPR-Cas system functions as a regulatory guide. We also identified mini-arrays in bacteriophages that could abrogate CRISPR immunity by inhibiting effector expression. Thus, recruitment of CRISPR effectors for regulatory functions via spacers with partial complementarity to the target is a common feature of diverse CRISPR-Cas systems.

Genome-centric view of the microbiome in a new deep-sea glass sponge species Bathydorus sp

Sponges are widely distributed in the global ocean and harbor diverse symbiotic microbes with mutualistic relationships. However, sponge symbionts in the deep sea remain poorly studied at the genome level. Here, we report a new glass sponge species of the genus Bathydorus and provide a genome-centric view of its microbiome. We obtained 14 high-quality prokaryotic metagenome-assembled genomes (MAGs) affiliated with the phyla Nitrososphaerota, Pseudomonadota, Nitrospirota, Bdellovibrionota, SAR324, Bacteroidota, and Patescibacteria. In total, 13 of these MAGs probably represent new species, suggesting the high novelty of the deep-sea glass sponge microbiome. An ammonia-oxidizing Nitrososphaerota MAG B01, which accounted for up to 70% of the metagenome reads, dominated the sponge microbiomes. The B01 genome had a highly complex CRISPR array, which likely represents an advantageous evolution toward a symbiotic lifestyle and forceful ability to defend against phages. A sulfur-oxidizing Gammaproteobacteria species was the second most dominant symbiont, and a nitrite-oxidizing Nitrospirota species could also be detected, but with lower relative abundance. Bdellovibrio species represented by two MAGs, B11 and B12, were first reported as potential predatory symbionts in deep-sea glass sponges and have undergone dramatic genome reduction. Comprehensive functional analysis indicated that most of the sponge symbionts encoded CRISPR-Cas systems and eukaryotic-like proteins for symbiotic interactions with the host. Metabolic reconstruction further illustrated their essential roles in carbon, nitrogen, and sulfur cycles. In addition, diverse putative phages were identified from the sponge metagenomes. Our study expands the knowledge of microbial diversity, evolutionary adaption, and metabolic complementarity in deep-sea glass sponges.

The coordination of anti-phage immunity mechanisms in bacterial cells

Bacterial cells are equipped with a variety of immune strategies to fight bacteriophage infections. Such strategies include unspecific mechanisms directed against any phage infecting the cell, ranging from the identification and cleavage of the viral DNA by restriction nucleases (restriction-modification systems) to the suicidal death of infected host cells (abortive infection, Abi). In addition, CRISPR-Cas systems generate an immune memory that targets specific phages in case of reinfection. However, the timing and coordination of different antiviral systems in bacterial cells are poorly understood. Here, we use simple mathematical models of immune responses in individual bacterial cells to propose that the intracellular dynamics of phage infections are key to addressing these questions. Our models suggest that the rates of viral DNA replication and cleavage inside host cells define functional categories of phages that differ in their susceptibility to bacterial anti-phage mechanisms, which could give raise to alternative phage strategies to escape bacterial immunity. From this viewpoint, the combined action of diverse bacterial defenses would be necessary to reduce the chances of phage immune evasion. The decision of individual infected cells to undergo suicidal cell death or to incorporate new phage sequences into their immune memory would be determined by dynamic interactions between the host's immune mechanisms and the phage DNA. Our work highlights the importance of within-cell dynamics to understand bacterial immunity, and formulates hypotheses that may inspire future research in this area.

Fitness costs of CRISPR-Cas systems in bacteria

CRISPR-Cas systems provide bacteria with both specificity and adaptability in defence against invading genetic elements. From a theoretical perspective, CRISPR-Cas systems confer many benefits. However, they are observed at an unexpectedly low prevalence across the bacterial domain. While these defence systems can be gained horizontally, fitness costs may lead to selection against their carriage. Understanding the source of CRISPR-related fitness costs will help us to understand the evolutionary dynamics of CRISPR-Cas systems and their role in shaping bacterial genome evolution. Here, we review our current understanding of the potential fitness costs associated with CRISPR-Cas systems. In addition to potentially restricting the acquisition of genetic material that could confer fitness benefits, we explore five alternative biological factors that from a theoretical perspective may influence the fitness costs associated with CRISPR-Cas system carriage: (1) the repertoire of defence mechanisms a bacterium has available to it, (2) the potential for a metabolic burden, (3) larger-scale population and environmental factors, (4) the phenomenon of self-targeting spacers, and (5) alternative non-defence roles for CRISPR-Cas.

A scaling law in CRISPR repertoire sizes arises from the avoidance of autoimmunity

Some prokaryotes possess CRISPR-Cas systems that use DNA segments called spacers, which are acquired from invading phages, to guide immune defense. Here, we propose that cross-reactive CRISPR targeting can, however, lead to "heterologous autoimmunity," whereby foreign spacers guide self-targeting in a spacer-length-dependent fashion. Balancing antiviral defense against autoimmunity predicts a scaling relation between spacer length and CRISPR repertoire size. We find evidence for this scaling through a comparative analysis of sequenced prokaryotic genomes and show that this association also holds at the level of CRISPR types. By contrast, the scaling is absent in strains with nonfunctional CRISPR loci. Finally, we demonstrate that stochastic spacer loss can explain variations around the scaling relation, even between strains of the same species. Our results suggest that heterologous autoimmunity is a selective factor shaping the evolution of CRISPR-Cas systems, analogous to the trade-offs between immune specificity, breadth, and autoimmunity that constrain the diversity of adaptive immune systems in vertebrates.

Single cell variability of CRISPR-Cas interference and adaptation

While CRISPR-Cas defence mechanisms have been studied on a population level, their temporal dynamics and variability in individual cells have remained unknown. Using a microfluidic device, time-lapse microscopy and mathematical modelling, we studied invader clearance in Escherichia coli across multiple generations. We observed that CRISPR interference is fast with a narrow distribution of clearance times. In contrast, for invaders with escaping PAM mutations we found large cell-to-cell variability, which originates from primed CRISPR adaptation. Faster growth and cell division and higher levels of Cascade increase the chance of clearance by interference, while slower growth is associated with increased chances of clearance by priming. Our findings suggest that Cascade binding to the mutated invader DNA, rather than spacer integration, is the main source of priming heterogeneity. The highly stochastic nature of primed CRISPR adaptation implies that only subpopulations of bacteria are able to respond quickly to invading threats. We conjecture that CRISPR-Cas dynamics and heterogeneity at the cellular level are crucial to understanding the strategy of bacteria in their competition with other species and phages.

Spacer prioritization in CRISPR-Cas9 immunity is enabled by the leader RNA

CRISPR-Cas systems store fragments of foreign DNA called spacers as immunological recordings used to combat future infections. Of the many spacers stored in a CRISPR array, the newest spacers are known to be prioritized for immune defense. However, the underlying mechanism remains unclear. Here we show that the leader region upstream of CRISPR arrays in CRISPR-Cas9 systems enhances CRISPR RNA (crRNA) processing from the newest spacer, prioritizing defense against the matching invader. Using the CRISPR-Cas9 system from Streptococcus pyogenes as a model, we found that the transcribed leader interacts with the conserved repeats bordering the newest spacer. The resulting interaction promotes tracrRNA hybridization with the second repeat, accelerating crRNA processing. Accordingly, disrupting this structure reduces the abundance of the associated crRNA and immune defense against targeted plasmids and bacteriophages. Beyond the S. pyogenes system, bioinformatics analyses revealed that leader-repeat structures appear across CRISPR-Cas9 systems. CRISPR-Cas systems thus possess an RNA-based mechanism to prioritize defense against the most recently encountered invaders.

CRISPR and KRAS: a match yet to be made

CRISPR (clustered regularly interspaced short palindromic repeats) systems are one of the most fascinating tools of the current era in molecular biotechnology. With the ease that they provide in genome editing, CRISPR systems generate broad opportunities for targeting mutations. Specifically in recent years, disease-causing mutations targeted by the CRISPR systems have been of main research interest; particularly for those diseases where there is no current cure, including cancer. KRAS mutations remain untargetable in cancer. Mutations in this oncogene are main drivers in common cancers, including lung, colorectal and pancreatic cancers, which are severe causes of public health burden and mortality worldwide, with no cure at hand. CRISPR systems provide an opportunity for targeting cancer causing mutations. In this review, we highlight the work published on CRISPR applications targeting KRAS mutations directly, as well as CRISPR applications targeting mutations in KRAS-related molecules. In specific, we focus on lung, colorectal and pancreatic cancers. To date, the limited literature on CRISPR applications targeting KRAS, reflect promising results. Namely, direct targeting of mutant KRAS variants using various CRISPR systems resulted in significant decrease in cell viability and proliferation in vitro, as well as tumor growth inhibition in vivo. In addition, the effect of mutant KRAS knockdown, via CRISPR, has been observed to exert regulatory effects on the downstream molecules including PI3K, ERK, Akt, Stat3, and c-myc. Molecules in the KRAS pathway have been subjected to CRISPR applications more often than KRAS itself. The aim of using CRISPR systems in these studies was mainly to analyze the therapeutic potential of possible downstream and upstream effectors of KRAS, as well as to discover further potential molecules. Although there have been molecules identified to have such potential in treatment of KRAS-driven cancers, a substantial amount of effort is still needed to establish treatment strategies based on these discoveries. We conclude that, at this point in time, despite being such a powerful directed genome editing tool, CRISPR remains to be underutilized for targeting KRAS mutations in cancer. Efforts channelled in this direction, might pave the way in solving the long-standing challenge of targeting the KRAS mutations in cancers.

Divalent-Metal-Ion Selectivity of the CRISPR-Cas System-Associated Cas1 Protein: Insights from Classical Molecular Dynamics Simulations and Electronic Structure Calculations

CRISPR-associated protein 1 (Cas1) is a universally conserved essential metalloenzyme of the clustered regularly interspaced short palindromic repeat (CRISPR) immune system of prokaryotes (bacteria, archaea) that can cut and integrate a part of viral DNA to its host genome with the help of other proteins. The integrated DNA acts as a memory of viral infection, which can be transcribed to RNA and stop future infection by recognition (based on the RNA/DNA complementarity principle) followed by protein-mediated degradation of the viral DNA. It has been proposed that the presence of a single manganese (Mn²⁺) ion in a conserved divalent-metal-ion binding pocket (key residues: E190, H254, D265, D268) of Cas1 is crucial for its function. Cas1-mediated DNA degradation was proposed to be hindered by metal substitution, metal chelation, or mutation of the binding pocket residues. Cas1 is active toward dsDNA degradation with both Mn²⁺ and Mg²⁺. X-ray structures of Cas1 revealed an intricate atomic interaction network of the divalent-metal-ion binding pocket and opened up the possibility of modeling related metal ions (viz., Mg²⁺, Ca²⁺) in the binding pocket of wild-type (WT) and mutated Cas1 proteins for computational analysis, which includes (1) quantitative estimation of the energetics of the divalent-metal-ion preference and (2) exploring the structural and dynamical aspects of the protein in response to divalent-metal-ion substitution or amino acid mutation. Using the X-ray structure of the Cas1 protein from Pseudomonas aeruginosa as a template (PDB 3GOD), we performed (∼2.23 μs) classical molecular dynamics (MD) simulations to compare structural and dynamical differences between Mg²⁺- and Ca²⁺-bound binding pockets of wild-type (WT) and mutant (E190A, H254A, D265A, D268A) Cas1. Furthermore, reduced binding pocket models were generated from X-ray and molecular dynamics (MD) trajectories, and the resulting structures were subjected to quantum chemical calculations. Results suggest that Cas1 prefers Mg²⁺ binding relative to Ca²⁺ and the preference is the strongest for WT and the weakest for the D268A mutant. Quantum chemical calculations indicate that Mn²⁺ is the most preferred relative to both Mg²⁺ and Ca²⁺ in the wild-type and mutant Cas1. Substitution of Mg²⁺ by Ca²⁺ does not alter the interaction network between Cas1 and the divalent metal ion but increases the wetness of the binding pocket by introducing a single water molecule in the first coordination shell of the latter. The strength of metal-ion preference (Mg²⁺ versus Ca²⁺) seems to be dependent on the solvent accessibility of the divalent-metal-ion binding pocket, strongest for wild-type Cas1 (in which the metal-ion binding pocket is dry, which includes two water molecules) and the weakest for the D268A mutant (in which the metal-ion binding pocket is wet, which includes four water molecules).

SCRAMBLER: A Tool for De Novo CRISPR Array Reconstruction and Its Application for Analysis of the Structure of Prokaryotic Populations

CRISPR arrays are prokaryotic genomic loci consisting of repeat sequences alternating with unique spacers acquired from foreign nucleic acids. As one of the fastest-evolving parts of the genome, CRISPR arrays can be used to differentiate closely related prokaryotic lineages and track individual strains in prokaryotic communities. However, the assembly of full-length CRISPR arrays sequences remains a problem. Here, we developed SCRAMBLER, a tool that includes several pipelines for assembling CRISPR arrays from high-throughput short-read sequencing data. We assessed its performance with model data sets (Escherichia coli strains containing different CRISPR arrays and imitating prokaryotic communities of different complexities) and intestinal microbiomes of extant and extinct pachyderms. Evaluation of SCRAMBLER's performance using model data sets demonstrated its ability to assemble CRISPR arrays correctly from reads containing pairs of spacers, yielding a precision rate of >80% and a recall rate of 60-85% when checked against ground-truth data. Likewise, SCRAMBLER successfully assembled CRISPR arrays from the environmental samples, as attested by their matching with database entries. SCRAMBLER, an open-source software (github.com/biolab-tools/SCRAMBLER), can facilitate analysis of the composition and dynamics of CRISPR arrays in complex communities.

Coevolution between bacterial CRISPR-Cas systems and their bacteriophages

CRISPR-Cas systems provide bacteria and archaea with adaptive, heritable immunity against their viruses (bacteriophages and phages) and other parasitic genetic elements. CRISPR-Cas systems are highly diverse, and we are only beginning to understand their relative importance in phage defense. In this review, we will discuss when and why CRISPR-Cas immunity against phages evolves, and how this, in turn, selects for the evolution of immune evasion by phages. Finally, we will discuss our current understanding of if, and when, we observe coevolution between CRISPR-Cas systems and phages, and how this may be influenced by the mechanism of CRISPR-Cas immunity.

Pruning and Tending Immune Memories: Spacer Dynamics in the CRISPR Array

CRISPR-Cas (Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR-associated genes) is a type of prokaryotic immune system that is unique in its ability to provide sequence-specific adaptive protection, which can be updated in response to new threats. CRISPR-Cas does this by storing fragments of DNA from invading genetic elements in an array interspersed with short repeats. The CRISPR array can be continuously updated through integration of new DNA fragments (termed spacers) at one end, but over time existing spacers become obsolete. To optimize immunity, spacer uptake, residency, and loss must be regulated. This mini-review summarizes what is known about how spacers are organized, maintained, and lost from CRISPR arrays.

Searching for fat tails in CRISPR-Cas systems: Data analysis and mathematical modeling

Understanding CRISPR-Cas systems-the adaptive defence mechanism that about half of bacterial species and most of archaea use to neutralise viral attacks-is important for explaining the biodiversity observed in the microbial world as well as for editing animal and plant genomes effectively. The CRISPR-Cas system learns from previous viral infections and integrates small pieces from phage genomes called spacers into the microbial genome. The resulting library of spacers collected in CRISPR arrays is then compared with the DNA of potential invaders. One of the most intriguing and least well understood questions about CRISPR-Cas systems is the distribution of spacers across the microbial population. Here, using empirical data, we show that the global distribution of spacer numbers in CRISPR arrays across multiple biomes worldwide typically exhibits scale-invariant power law behaviour, and the standard deviation is greater than the sample mean. We develop a mathematical model of spacer loss and acquisition dynamics which fits observed data from almost four thousand metagenomes well. In analogy to the classical 'rich-get-richer' mechanism of power law emergence, the rate of spacer acquisition is proportional to the CRISPR array size, which allows a small proportion of CRISPRs within the population to possess a significant number of spacers. Our study provides an alternative explanation for the rarity of all-resistant super microbes in nature and why proliferation of phages can be highly successful despite the effectiveness of CRISPR-Cas systems.

Guide RNA Categorization Enables Target Site Choice in Tn7-CRISPR-Cas Transposons

CRISPR-Cas defense systems have been coopted multiple times in nature for guide RNA-directed transposition by Tn7-like elements. Prototypic Tn7 uses dedicated proteins for two targeting pathways: one targeting a neutral and conserved attachment site in the chromosome and a second directing transposition into mobile plasmids facilitating cell-to-cell transfer. We show that Tn7-CRISPR-Cas elements evolved a system of guide RNA categorization to accomplish the same two-pathway lifestyle. Multiple mechanisms allow functionally distinct guide RNAs for transposition: a conventional system capable of acquiring guide RNAs to new plasmid and phage targets and a second providing long-term memory for access to chromosomal sites upon entry into a new host. Guide RNAs are privatized to be recognized only by the transposon-adapted system via sequence specialization, mismatch tolerance, and selective regulation to avoid toxic self-targeting by endogenous CRISPR-Cas defense systems. This information reveals promising avenues to engineer guide RNAs for enhanced CRISPR-Cas functionality for genome modification.

Transcription termination and antitermination of bacterial CRISPR arrays

A hallmark of CRISPR-Cas immunity systems is the CRISPR array, a genomic locus consisting of short, repeated sequences ('repeats') interspersed with short, variable sequences ('spacers'). CRISPR arrays are transcribed and processed into individual CRISPR RNAs that each include a single spacer, and direct Cas proteins to complementary sequences in invading nucleic acid. Most bacterial CRISPR array transcripts are unusually long for untranslated RNA, suggesting the existence of mechanisms to prevent premature transcription termination by Rho, a conserved bacterial transcription termination factor that rapidly terminates untranslated RNA. We show that Rho can prematurely terminate transcription of bacterial CRISPR arrays, and we identify a widespread antitermination mechanism that antagonizes Rho to facilitate complete transcription of CRISPR arrays. Thus, our data highlight the importance of transcription termination and antitermination in the evolution of bacterial CRISPR-Cas systems.

Comparative Genomics Uncovers the Genetic Diversity and Characters of Veillonella atypica and Provides Insights Into Its Potential Applications

Veillonella atypica is a bacterium that is present in the gut and the oral cavity of mammals and plays diverse roles in different niches. A recent study demonstrated that Veillonella is highly associated with marathon running and approved that V. atypica gavage improves treadmill run time in mice, revealing that V. atypica has a high biotechnological potential in improving athlete performance. However, a comprehensive analysis of the genetic diversity, function traits, and genome editing method of V. atypica remains elusive. In the present study, we conducted a systemically comparative analysis of the genetic datasets of nine V. atypica strains. The pan-genome of V. atypica consisted of 2,065 homologous clusters and exhibited an open pan-genome structure. A phylogenetic analysis of V. atypica with two different categories revealed that V. atypica OK5 was the most distant from the other eight V. atypica strains. A total of 43 orthologous genes were identified as CAZyme genes and grouped into 23 CAZyme families. The CAZyme components derived from accessory clusters contributed to the differences in the ability of the nine V. atypica strains to utilize carbohydrates. An integrated analysis of the metabolic pathways of V. atypica suggested that V. atypica strains harbored vancomycin resistance and were involved in several biosynthesis pathways of secondary metabolites. The V. atypica strains harbored four main Cas proteins, namely, CAS-Type IIIA, CAS-Type IIA, CAS-Type IIC, and CAS-Type IIID. This pilot study provides an in-depth understanding of and a fundamental knowledge about the biology of V. atypica that allow the possibility to increase the biotechnological potential of this bacterium.

The size of the immune repertoire of bacteria

Some bacteria possess an adaptive immune system that maintains a memory of past viral infections in the CRISPR loci of their genomes. This memory is used to mount targeted responses against later threats but is remarkably shallow: it remembers only a few dozen to a few hundred viruses. We present a statistical theory of CRISPR-based immunity that quantitatively predicts the depth of bacterial immune memory in terms of a tradeoff with fundamental constraints of the cellular biochemical machinery.

Some bacteria and archaea possess an immune system, based on the CRISPR-Cas mechanism, that confers adaptive immunity against viruses. In such species, individual prokaryotes maintain cassettes of viral DNA elements called spacers as a memory of past infections. Typically, the cassettes contain several dozen expressed spacers. Given that bacteria can have very large genomes and since having more spacers should confer a better memory, it is puzzling that so little genetic space would be devoted by prokaryotes to their adaptive immune systems. Here, assuming that CRISPR functions as a long-term memory-based defense against a diverse landscape of viral species, we identify a fundamental tradeoff between the amount of immune memory and effectiveness of response to a given threat. This tradeoff implies an optimal size for the prokaryotic immune repertoire in the observational range.

Atypical organizations and epistatic interactions of CRISPRs and cas clusters in genomes and their mobile genetic elements

Prokaryotes use CRISPR-Cas systems for adaptive immunity, but the reasons for the frequent existence of multiple CRISPRs and cas clusters remain poorly understood. Here, we analysed the joint distribution of CRISPR and cas genes in a large set of fully sequenced bacterial genomes and their mobile genetic elements. Our analysis suggests few negative and many positive epistatic interactions between Cas subtypes. The latter often result in complex genetic organizations, where a locus has a single adaptation module and diverse interference mechanisms that might provide more effective immunity. We typed CRISPRs that could not be unambiguously associated with a cas cluster and found that such complex loci tend to have unique type I repeats in multiple CRISPRs. Many chromosomal CRISPRs lack a neighboring Cas system and they often have repeats compatible with the Cas systems encoded in trans. Phages and 25 000 prophages were almost devoid of CRISPR-Cas systems, whereas 3% of plasmids had CRISPR-Cas systems or isolated CRISPRs. The latter were often compatible with the chromosomal cas clusters, suggesting that plasmids can co-opt the latter. These results highlight the importance of interactions between CRISPRs and cas present in multiple copies and in distinct genomic locations in the function and evolution of bacterial immunity.

Direct Visualization of Native CRISPR Target Search in Live Bacteria Reveals Cascade DNA Surveillance Mechanism

CRISPR-Cas systems encode RNA-guided surveillance complexes to find and cleave invading DNA elements. While it is thought that invaders are neutralized minutes after cell entry, the mechanism and kinetics of target search and its impact on CRISPR protection levels have remained unknown. Here, we visualize individual Cascade complexes in a native type I CRISPR-Cas system. We uncover an exponential relation between Cascade copy number and CRISPR interference levels, pointing to a time-driven arms race between invader replication and target search, in which 20 Cascade complexes provide 50% protection. Driven by PAM-interacting subunit Cas8e, Cascade spends half its search time rapidly probing DNA (∼30 ms) in the nucleoid. We further demonstrate that target DNA transcription and CRISPR arrays affect the integrity of Cascade and affect CRISPR interference. Our work establishes the mechanism of cellular DNA surveillance by Cascade that allows the timely detection of invading DNA in a crowded, DNA-packed environment.

Long reads reveal the diversification and dynamics of CRISPR reservoir in microbiomes

BACKGROUND

Sequencing of microbiomes has accelerated the characterization of the diversity of CRISPR-Cas immune systems. However, the utilization of next generation short read sequences for the characterization of CRISPR-Cas dynamics remains limited due to the repetitive nature of CRISPR arrays. CRISPR arrays are comprised of short spacer segments (derived from invaders' genomes) interspaced between flanking repeat sequences. The repetitive structure of CRISPR arrays poses a computational challenge for the accurate assembly of CRISPR arrays from short reads. In this paper we evaluate the use of long read sequences for the analysis of CRISPR-Cas system dynamics in microbiomes.

RESULTS

We analyzed a dataset of Illumina's TruSeq Synthetic Long-Reads (SLR) derived from a gut microbiome. We showed that long reads captured CRISPR spacers at a high degree of redundancy, which highlights the spacer conservation of spacer sharing CRISPR variants, enabling the study of CRISPR array dynamics in ways difficult to achieve though short read sequences. We introduce compressed spacer graphs, a visual abstraction of spacer sharing CRISPR arrays, to provide a simplified view of complex organizational structures present within CRISPR array dynamics. Utilizing compressed spacer graphs, several key defining characteristics of CRISPR-Cas system dynamics were observed including spacer acquisition and loss events, conservation of the trailer end spacers, and CRISPR arrays' directionality (transcription orientation). Other result highlights include the observation of intense array contraction and expansion events, and reconstruction of a full-length genome for a potential invader (Faecalibacterium phage) based on identified spacers.

CONCLUSION

We demonstrate in an in silico system that long reads provide the necessary context for characterizing the organization of CRISPR arrays in a microbiome, and reveal dynamic and evolutionary features of CRISPR-Cas systems in a microbial population.

Probing the Behaviour of Cas1-Cas2 upon Protospacer Binding in CRISPR-Cas Systems using Molecular Dynamics Simulations

Adaptation in CRISPR-Cas systems enables the generation of an immunological memory to defend against invading viruses. This process is driven by foreign DNA spacer (termed protospacer) selection and integration mediated by Cas1-Cas2 protein. Recently, different states of Cas1-Cas2, in its free form and in complex with protospacer DNAs, were solved by X-ray crystallography. In this paper, molecular dynamics (MD) simulations are employed to study crystal structures of one free and two protospacer-bound Cas1-Cas2 complexes. The simulated results indicate that the protospacer binding markedly increases the system stability, in particular when the protospacer containing the PAM-complementary sequence. The hydrogen bond and binding free energy calculations explain that PAM recognition introduces more specific interactions to increase the cleavage activity of Cas1. By using principal component analysis (PCA) and intramolecular angle calculation, this study observes two dominant slow motions associated with the binding of Ca1-Cas2 to the protospacer and potential target DNAs respectively. The comparison of DNA structural deformation further implies a cooperative conformational change of Cas1-Cas2 and protospacer for the target DNA capture. We propose that this cooperativity is the intrinsic requirement of the CRISPR integration complex formation. This study provides some new insights into the understanding of CRISPR-Cas adaptation.

Imprecise Spacer Acquisition Generates CRISPR-Cas Immune Diversity through Primed Adaptation

Many prokaryotes possess CRISPR-Cas adaptive immune systems to defend against viruses and invading mobile genetic elements. CRISPR-Cas immunity relies on genetic memories, termed spacers, for sequence-specific recognition of infections. The diversity of spacers within host populations is important for immune resilience, but we have limited understanding of how CRISPR diversity is generated. Type I CRISPR-Cas systems use existing spacers to enhance the acquisition of new spacers through primed CRISPR adaptation (priming). Here, we present a pathway to priming that is stimulated by imprecisely acquired (slipped) spacers. Slipped spacers are less effective for immunity but increase priming compared with canonical spacers. The benefits of slipping depend on the relative rates of phage mutation and adaptation during defense. We propose that slipped spacers provide a route to increase population-level spacer diversity that pre-empts phage escape mutant proliferation and that the trade-off between adaptation and immunity is important in diverse CRISPR-Cas systems.

Selective Maintenance of Multiple CRISPR Arrays Across Prokaryotes

Prokaryotes are under nearly constant attack by viral pathogens. To protect against this threat of infection, bacteria and archaea have evolved a wide array of defense mechanisms, singly and in combination. While immune diversity in a single organism likely reduces the chance of pathogen evolutionary escape, it remains puzzling why many prokaryotes also have multiple, seemingly redundant, copies of the same type of immune system. Here, we focus on the highly flexible CRISPR adaptive immune system, which is present in multiple copies in a surprising 28% of the prokaryotic genomes in RefSeq. We use a comparative genomics approach looking across all prokaryotes to demonstrate that on average, organisms are under selection to maintain more than one CRISPR array. Given this surprising conclusion, we consider several hypotheses concerning the source of selection and include a theoretical analysis of the possibility that a trade-off between memory span and learning speed could select for both "long-term memory" and "short-term memory" CRISPR arrays.