Heterogeneous diversity of spacers within CRISPR (clustered regularly interspaced short palindromic repeats)

Characterization of CRISPR-Cas Systems in Shewanella algae and Shewanella haliotis: Insights into the Adaptation and Survival of Marine Pathogens

CRISPR-Cas systems are adaptive immune mechanisms present in most prokaryotes that play an important role in the adaptation of bacteria and archaea to new environments. Shewanella algae is a marine zoonotic pathogen with worldwide distribution, which accounts for the majority of clinical cases of Shewanella infections. However, the characterization of Shewanella algae CRISPR-Cas systems has not been well investigated yet. Through whole genome sequence analysis, we characterized the CRISPR-Cas systems in S. algae. Our results indicate that CRISPR-Cas systems are prevalent in S. algae, with the majority of strains containing the Type I-F system. This study provides new insights into the diversity and function of CRISPR-Cas systems in S. algae and highlights their potential role in the adaptation and survival of these marine pathogens.

Transient eco-evolutionary dynamics early in a phage epidemic have strong and lasting impact on the long-term evolution of bacterial defences

Organisms have evolved a range of constitutive (always active) and inducible (elicited by parasites) defence mechanisms, but we have limited understanding of what drives the evolution of these orthogonal defence strategies. Bacteria and their phages offer a tractable system to study this: Bacteria can acquire constitutive resistance by mutation of the phage receptor (surface mutation, sm) or induced resistance through their CRISPR-Cas adaptive immune system. Using a combination of theory and experiments, we demonstrate that the mechanism that establishes first has a strong advantage because it weakens selection for the alternative resistance mechanism. As a consequence, ecological factors that alter the relative frequencies at which the different resistances are acquired have a strong and lasting impact: High growth conditions promote the evolution of sm resistance by increasing the influx of receptor mutation events during the early stages of the epidemic, whereas a high infection risk during this stage of the epidemic promotes the evolution of CRISPR immunity, since it fuels the (infection-dependent) acquisition of CRISPR immunity. This work highlights the strong and lasting impact of the transient evolutionary dynamics during the early stages of an epidemic on the long-term evolution of constitutive and induced defences, which may be leveraged to manipulate phage resistance evolution in clinical and applied settings.

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.

CRISPR-Cas9: A History of Its Discovery and Ethical Considerations of Its Use in Genome Editing

The development of a method for genome editing based on CRISPR-Cas9 technology was awarded The Nobel Prize in Chemistry in 2020, less than a decade after the discovery of all principal molecular components of the system. For the first time in history a Nobel prize was awarded to two women, Emmanuelle Charpentier and Jennifer Doudna, who made key discoveries in the field of DNA manipulation with the CRISPR-Cas9 system, so-called "genetic scissors". It is difficult to overestimate the importance of the technique as it enables one not only to manipulate genomes of model organisms in scientific experiments, and modify characteristics of important crops and animals, but also has the potential of introducing revolutionary changes in medicine, especially in treatment of genetic diseases. The original biological function of CRISPR-Cas9 system is the protection of prokaryotes from mobile genetic elements, in particular viruses. Currently, CRISPR-Cas9 and related technologies have been successfully used to cure life-threatening diseases, make coronavirus detection tests, and even to modify human embryo cells with the consequent birth of babies carrying the introduced modifications. This intervention with human germplasm cells resulted in wide disapproval in the scientific community due to ethical concerns, and calls for a moratorium on inheritable genomic manipulations. This review focuses on the history of the discovery of the CRISPR-Cas9 system with some aspects of its current applications, including ethical concerns about its use in humans.

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.

CRISPRidentify: identification of CRISPR arrays using machine learning approach

CRISPR–Cas are adaptive immune systems that degrade foreign genetic elements in archaea and bacteria. In carrying out their immune functions, CRISPR–Cas systems heavily rely on RNA components. These CRISPR (cr) RNAs are repeat-spacer units that are produced by processing of pre-crRNA, the transcript of CRISPR arrays, and guide Cas protein(s) to the cognate invading nucleic acids, enabling their destruction. Several bioinformatics tools have been developed to detect CRISPR arrays based solely on DNA sequences, but all these tools employ the same strategy of looking for repetitive patterns, which might correspond to CRISPR array repeats. The identified patterns are evaluated using a fixed, built-in scoring function, and arrays exceeding a cut-off value are reported. Here, we instead introduce a data-driven approach that uses machine learning to detect and differentiate true CRISPR arrays from false ones based on several features. Our CRISPR detection tool, CRISPRidentify, performs three steps: detection, feature extraction and classification based on manually curated sets of positive and negative examples of CRISPR arrays. The identified CRISPR arrays are then reported to the user accompanied by detailed annotation. We demonstrate that our approach identifies not only previously detected CRISPR arrays, but also CRISPR array candidates not detected by other tools. Compared to other methods, our tool has a drastically reduced false positive rate. In contrast to the existing tools, our approach not only provides the user with the basic statistics on the identified CRISPR arrays but also produces a certainty score as a practical measure of the likelihood that a given genomic region is a CRISPR array.

Cost and benefits of clustered regularly interspaced short palindromic repeats spacer acquisition

Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas-mediated immunity in bacteria allows bacterial populations to protect themselves against pathogens. However, it also exposes them to the dangers of auto-immunity by developing protection that targets its own genome. Using a simple model of the coupled dynamics of phage and bacterial populations, we explore how acquisition rates affect the probability of the bacterial colony going extinct. We find that the optimal strategy depends on the initial population sizes of both viruses and bacteria. Additionally, certain combinations of acquisition and dynamical rates and initial population sizes guarantee protection, owing to a dynamical balance between the evolving population sizes, without relying on acquisition of viral spacers. Outside this regime, the high cost of auto-immunity limits the acquisition rate. We discuss these optimal strategies that minimize the probability of the colony going extinct in terms of recent experiments.

This article is part of a discussion meeting issue ‘The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems’.

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.

Non-conventional therapeutic technique to replace CRISPR bacteria from biofilm by inducible lysogen

Bacteriophage can be an effective means of regulating bacterial populations when conditions allow phage invasion of bacterial colonies. Phage can either infect and lyse a host cell, or insert their DNA into the host cell genome; the latter process is called lysogeny. The clustered regularly interspaced short palindromic repeat (CRISPR) system, linked with CRISPR-associated (Cas) genes, is a regulatory system present in a variety of bacteria which confers immunity against bacteriophage. Studies of the group behaviour of bacteria with CRISPR/Cas systems have provided evidence that CRISPR in lysogenized bacteria can cause an inability to form biofilm. This allows CRISPR-immune bacteria in biofilms to effectively resist phage therapy. Our recent work has described a potential therapeutic technique to eradicate CRISPR-immune bacteria from a biofilm by a continuous influx of lysogens carrying an identical phage sequence. However, this model predicted that the CRISPR-immune population could persist for long times before eradication. Our current focus is on the use of diverse lysogens against CRISPR-capable bacterial populations. The goal of this work is to find a suitable strategy which can eradicate bacteria with a CRISPR system through the influx of finite amounts of distinct lysogens over fixed intervals.

How adaptive immunity constrains the composition and fate of large bacterial populations

Features of the CRISPR-Cas system, in which bacteria integrate small segments of phage genome (spacers) into their DNA to neutralize future attacks, suggest that its effect is not limited to individual bacteria but may control the fate and structure of whole populations. Emphasizing the population-level impact of the CRISPR-Cas system, recent experiments show that some bacteria regulate CRISPR-associated genes via the quorum sensing (QS) pathway. Here we present a model that shows that from the highly stochastic dynamics of individual spacers under QS control emerges a rank-abundance distribution of spacers that is time invariant, a surprising prediction that we test with dynamic spacer-tracking data from literature. This distribution depends on the state of the competing phage-bacteria population, which due to QS-based regulation may coexist in multiple stable states that vary significantly in their phage-to-bacterium ratio, a widely used ecological measure to characterize microbial systems.

The physicist's guide to one of biotechnology's hottest new topics: CRISPR-Cas

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated proteins (Cas) constitute a multi-functional, constantly evolving immune system in bacteria and archaea cells. A heritable, molecular memory is generated of phage, plasmids, or other mobile genetic elements that attempt to attack the cell. This memory is used to recognize and interfere with subsequent invasions from the same genetic elements. This versatile prokaryotic tool has also been used to advance applications in biotechnology. Here we review a large body of CRISPR-Cas research to explore themes of evolution and selection, population dynamics, horizontal gene transfer, specific and cross-reactive interactions, cost and regulation, non-immunological CRISPR functions that boost host cell robustness, as well as applicable mechanisms for efficient and specific genetic engineering. We offer future directions that can be addressed by the physics community. Physical understanding of the CRISPR-Cas system will advance uses in biotechnology, such as developing cell lines and animal models, cell labeling and information storage, combatting antibiotic resistance, and human therapeutics.

Bacterial CRISPR Regions: General Features and their Potential for Epidemiological Molecular Typing Studies

Introduction

CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) loci as novel and applicable regions in prokaryotic genomes have gained great attraction in the post genomics era.

Methods

These unique regions are diverse in number and sequence composition in different pathogenic bacteria and thereby can be a suitable candidate for molecular epidemiology and genotyping studies. Results:Furthermore, the arrayed structure of CRISPR loci (several unique repeats spaced with the variable sequence) and associated cas genes act as an active prokaryotic immune system against viral replication and conjugative elements. This property can be used as a tool for RNA editing in bioengineering studies.

Conclusion

The aim of this review was to survey some details about the history, nature, and potential applications of CRISPR arrays in both genetic engineering and bacterial genotyping studies.

Optimal number of spacers in CRISPR arrays

Prokaryotic organisms survive under constant pressure of viruses. CRISPR-Cas system provides its prokaryotic host with an adaptive immune defense against viruses that have been previously encountered. It consists of two components: Cas-proteins that cleave the foreign DNA and CRISPR array that suits as a virus recognition key. CRISPR array consists of a series of spacers, short pieces of DNA that originate from and match the corresponding parts of viral DNA called protospacers. Here we estimate the number of spacers in a CRISPR array of a prokaryotic cell which maximizes its protection against a viral attack. The optimality follows from a competition between two trends: too few distinct spacers make host vulnerable to an attack by a virus with mutated corresponding protospacers, while an excessive variety of spacers dilutes the number of the CRISPR complexes armed with the most recent and thus most useful spacers. We first evaluate the optimal number of spacers in a simple scenario of an infection by a single viral species and later consider a more general case of multiple viral species. We find that depending on such parameters as the concentration of CRISPR-Cas interference complexes and its preference to arm with more recently acquired spacers, the rate of viral mutation, and the number of viral species, the predicted optimal number of spacers lies within a range that agrees with experimentally-observed values.

CRISPR-Cas systems provide adaptive immunity defense in bacteria and archaea against viruses. They function by accumulating in prokaryotic genome an array of spacers, or fragments of virus DNA from previous attacks. By matching spacers to corresponding parts of viral DNA called protospacers, a CRISPR-Cas system identifies and destroys intruder DNA. Here we theoretically estimate the number of spacers that maximizes prokaryotic host cell survival. This optimum emerges from a competition between two trends: More spacers allow a prokaryotic cell to hedge against mutations in viral protospacers. However, the older spacers loose efficiency as corresponding protospacers mutate. For a limited pool of CRISPR-Cas molecular machines, keeping too many spacers leaves fewer of such machines armed with more efficient young (most recently acquired) spacers. We have shown that a higher efficiency of CRISPR-Cas system allows a prokaryotic cell to utilize more spacers, increasing the optimal number of spacers. On contrary, a higher viral mutation rate makes older spacers useless and favors shorter arrays. A higher diversity in viral species reduces the efficiency of CRISPR-Cas but does not necessary lead to longer arrays. Our study provides a new viewpoint at a variety of the observed array spacer number and could be used as a base for evolutionary models of host-phage coexistence.

Non-classical phase diagram for virus bacterial coevolution mediated by clustered regularly interspaced short palindromic repeats

CRISPR is a newly discovered prokaryotic immune system. Bacteria and archaea with this system incorporate genetic material from invading viruses into their genomes, providing protection against future infection by similar viruses. The condition for coexistence of prokaryots and viruses is an interesting problem in evolutionary biology. In this work, we show an intriguing phase diagram of the virus extinction probability, which is more complex than that of the classical predator-prey model. As the CRISPR incorporates genetic material, viruses are under pressure to evolve to escape recognition by CRISPR. When bacteria have a small rate of deleting spacers, a new parameter region in which bacteria and viruses can coexist arises, and it leads to a more complex coexistence patten for bacteria and viruses. For example, when the virus mutation rate is low, the virus extinction probability changes non-montonically with the bacterial exposure rate. The virus and bacteria coevolution not only alters the virus extinction probability, but also changes the bacterial population structure. Additionally, we show that recombination is a successful strategy for viruses to escape from CRISPR recognition when viruses have multiple proto-spacers, providing support for a recombination-mediated escape mechanism suggested experimentally. Finally, we suggest that the re-entrant phase diagram, in which phages can progress through three phases of extinction and two phases of abundance at low spacer deletion rates as a function of exposure rate to bacteria, is an experimentally testable phenomenon.

Mathematical modelling of CRISPR-Cas system effects on biofilm formation

Clustered regularly interspaced short palindromic repeats (CRISPR), linked with CRISPR associated (Cas) genes, can confer adaptive immunity to bacteria, against bacteriophage infections. Thus from a therapeutic standpoint, CRISPR immunity increases biofilm resistance to phage therapy. Recently, however, CRISPR-Cas genes have been implicated in reducing biofilm formation in lysogenized cells. Thus CRISPR immunity can have complex effects on phage-host-lysogen interactions, particularly in a biofilm. In this contribution, we develop and analyse a series of dynamical systems to elucidate and disentangle these interactions. Two competition models are used to study the effects of lysogens (first model) and CRISPR-immune bacteria (second model) in the biofilm. In the third model, the effect of delivering lysogens to a CRISPR-immune biofilm is investigated. Using standard analyses of equilibria, stability and bifurcations, our models predict that lysogens may be able to displace CRISPR-immune bacteria in a biofilm, and thus suggest strategies to eliminate phage-resistant biofilms.

Dynamics of adaptive immunity against phage in bacterial populations

The CRISPR (clustered regularly interspaced short palindromic repeats) mechanism allows bacteria to adaptively defend against phages by acquiring short genomic sequences (spacers) that target specific sequences in the viral genome. We propose a population dynamical model where immunity can be both acquired and lost. The model predicts regimes where bacterial and phage populations can co-exist, others where the populations exhibit damped oscillations, and still others where one population is driven to extinction. Our model considers two key parameters: (1) ease of acquisition and (2) spacer effectiveness in conferring immunity. Analytical calculations and numerical simulations show that if spacers differ mainly in ease of acquisition, or if the probability of acquiring them is sufficiently high, bacteria develop a diverse population of spacers. On the other hand, if spacers differ mainly in their effectiveness, their final distribution will be highly peaked, akin to a “winner-take-all” scenario, leading to a specialized spacer distribution. Bacteria can interpolate between these limiting behaviors by actively tuning their overall acquisition probability.

The CRISPR system in bacteria and archaea provides adaptive immunity by incorporating foreign DNA (spacers) into the genome, and later targeting DNA sequences that match these spacers. The way in which bacteria choose spacer sequences from a clonal phage population is not understood. Our model considers competing effects of ease of acquisition and effectiveness against infections in shaping the spacer distribution. The model suggests that a diverse spacer population results when the acquisition rate is high, or when spacers are similarly effective. At moderate acquisition rates, the spacer distribution becomes highly sensitive to spacer effectiveness. There is a rich landscape of behaviors including bacteria-phage coexistence and oscillations in the populations.

Reversals and collisions optimize protein exchange in bacterial swarms

Swarming groups of bacteria coordinate their behavior by self-organizing as a population to move over surfaces in search of nutrients and optimal niches for colonization. Many open questions remain about the cues used by swarming bacteria to achieve this self-organization. While chemical cue signaling known as quorum sensing is well-described, swarming bacteria often act and coordinate on time scales that could not be achieved via these extracellular quorum sensing cues. Here, cell-cell contact-dependent protein exchange is explored as a mechanism of intercellular signaling for the bacterium Myxococcus xanthus. A detailed biologically calibrated computational model is used to study how M. xanthus optimizes the connection rate between cells and maximizes the spread of an extracellular protein within the population. The maximum rate of protein spreading is observed for cells that reverse direction optimally for swarming. Cells that reverse too slowly or too fast fail to spread extracellular protein efficiently. In particular, a specific range of cell reversal frequencies was observed to maximize the cell-cell connection rate and minimize the time of protein spreading. Furthermore, our findings suggest that predesigned motion reversal can be employed to enhance the collective behavior of biological synthetic active systems.

Computational models of populations of bacteria and lytic phage

The use of phages to control and reduce numbers of unwanted bacteria can be traced back to the early 1900s, when phages were explored as a tool to treat infections before the wide scale use of antibiotics. Recently, phage therapy has received renewed interest as a method to treat multiresistant bacteria. Phages are also widely used in the food industry to prevent the growth of certain bacteria in foods, and are currently being explored as a tool for use in bioremediation and wastewater treatment. Despite the large body of biological research on phages, relatively little attention has been given to computational modeling of the population dynamics of phage and bacterial interactions. The earliest model was described by Campbell in the 1960s. Subsequent modifications to this model include partial or complete resistance, multiple phage binding sites, and spatial heterogeneity. This review provides a general introduction to modeling of the population dynamics of bacteria and phage. The review introduces the basic model and relevant concepts and evaluates more complex variations of the basic model published to date, including a model of disease epidemics caused by infectious bacteria. Finally, the shortcomings and potential ways to improve the models are discussed.

CRISPR-Based Typing and Next-Generation Tracking Technologies

Bacteria occur ubiquitously in nature and are broadly relevant throughout the food supply chain, with diverse and variable tolerance levels depending on their origin, biological role, and impact on the quality and safety of the product as well as on the health of the consumer. With increasing knowledge of and accessibility to the microbial composition of our environments, food supply, and host-associated microbiota, our understanding of and appreciation for the ratio of beneficial to undesirable bacteria are rapidly evolving. Therefore, there is a need for tools and technologies that allow definite, accurate, and high-resolution identification and typing of various groups of bacteria that include beneficial microbes such as starter cultures and probiotics, innocuous commensals, and undesirable pathogens and spoilage organisms. During the transition from the current molecular biology-based PFGE (pulsed-field gel electrophoresis) gold standard to the increasingly accessible omics-level whole-genome sequencing (WGS) N-gen standard, high-resolution technologies such as CRISPR-based genotyping constitute practical and powerful alternatives that provide valuable insights into genome microevolution and evolutionary trajectories. Indeed, several studies have shown potential for CRISPR-based typing of industrial starter cultures, health-promoting probiotic strains, animal commensal species, and problematic pathogens. Emerging CRISPR-based typing methods open new avenues for high-resolution typing of a broad range of bacteria and constitute a practical means for rapid tracking of a diversity of food-associated microbes.

Comparative genomics of Streptococcus pyogenes M1 isolates differing in virulence and propensity to cause systemic infection in mice

Streptococcus pyogenes serotype M1 is a frequent cause of severe infections in humans. Some M1 isolates are pathogenic in mice and used in studies on infection pathogenesis. We observed marked differences in murine infections caused by M1 strain SF370, 5448, 5448AP or AP1 which prompted us to sequence the whole genome of isolates 5448 and AP1 for comparative analysis. Strain 5448 is known to acquire inactivating mutations in the CovRS two-component system during mouse infection, producing hypervirulent progeny such as 5448AP. Isolates AP1 and 5448AP, more than 5448, caused disseminating infections that became systemic and lethal. SF370 was not pathogenic. Phages caused gross genetic differences and increased the gene content of AP1 by 8% as compared to 5448 and SF370. Each of six examined M1 genomes contained two CRISPR-Cas systems. Phage insertion destroyed a type II CRISPR-Cas system in AP1 and other strains of serotypes M1, M3, M6 and M24, but not in M1 strains 5448, SF370, MGAS5005, A20 or M1 476. A resulting impaired defence against invading genetic elements could have led to the wealth of phages in AP1. AP1 lacks genetic features of the MGAS5005-like clonal complex including the streptodornase that drives selection for hypervirulent clones with inactivated CovRS system. Still, inactivating mutations in covS were a common genetic feature of AP1 and the MGAS5005-like isolate 5448AP. Abolished expression of the cysteine proteinase SpeB, due to CovRS inactivation could be a common cause for hypervirulence of the two isolates. Moreover, an additional protein H-coding gene and a mutation in the regulator gene rofA distinguished AP1 form other M1 isolates. In conclusion, hypervirulence of S. pyogenes M1 in mice is not limited to the MGAS5005-like genotype.

How a well-adapted immune system is organized

The repertoire of lymphocyte receptors in the adaptive immune system protects organisms from diverse pathogens. A well-adapted repertoire should be tuned to the pathogenic environment to reduce the cost of infections. We develop a general framework for predicting the optimal repertoire that minimizes the cost of infections contracted from a given distribution of pathogens. The theory predicts that the immune system will have more receptors for rare antigens than expected from the frequency of encounters; individuals exposed to the same infections will have sparse repertoires that are largely different, but nevertheless exploit cross-reactivity to provide the same coverage of antigens; and the optimal repertoires can be reached via the dynamics of competitive binding of antigens by receptors and selective amplification of stimulated receptors. Our results follow from a tension between the statistics of pathogen detection, which favor a broader receptor distribution, and the effects of cross-reactivity, which tend to concentrate the optimal repertoire onto a few highly abundant clones. Our predictions can be tested in high-throughput surveys of receptor and pathogen diversity.

Evolution of the CRISPR-Cas adaptive immunity systems in prokaryotes: models and observations on virus-host coevolution

CRISPR-Cas is an adaptive immunity system in prokaryotes that functions via a unique mechanism which involves incorporation of foreign DNA fragments into CRISPR arrays and subsequent utilization of transcripts of these inserts (known as spacers) as guide RNAs to cleave the cognate selfish element genome. Multiple attempts have been undertaken to explore the coevolution of viruses and microbial hosts carrying CRISPR-Cas using mathematical models that employ either systems of differential equations or an agent-based approach, or combinations thereof. Analysis of these models reveals highly complex co-evolutionary dynamics that ensues from the combination of the heritability of the CRISPR-mediated adaptive immunity with the existence of different degrees of immunity depending on the number of cognate spacers and the cost of carrying a CRISPR-Cas locus. Depending on the details of the models, a variety of testable, sometimes conflicting predictions have been made on the dependence of the degree of immunity and the benefit of maintaining CRISPR-Cas on the abundance and diversity of hosts and viruses. Some of these predictions have already been directly validated experimentally. In particular, both the reality of the virus-host arms race, with viruses escaping resistance and hosts reacquiring it through the capture of new spacers, and the fitness cost of CRISPR-Cas due to the curtailment of beneficial HGT have been reproduced in the laboratory. However, to test the predictions of the models more specifically, detailed studies of coevolving populations of microbes and viruses both in nature and in the laboratory are essential. Such analyses are expected to yield disagreements with the predictions of the current, oversimplified models and to trigger a new round of theoretical developments.

Evolutionary causes and consequences of diversified CRISPR immune profiles in natural populations

Host-pathogen co-evolution is a significant force which shapes the ecology and evolution of all types of organisms, and such interactions are driven by resistance and immunity mechanisms of the host. Diversity of resistance and immunity can affect the co-evolutionary trajectory of both host and pathogen. The microbial CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) system is one host immunity mechanism which offers a tractable model for examining the dynamics of diversity in an immune system. In the present article, we review CRISPR variation observed in a variety of natural populations, examine the forces which can push CRISPRs towards high or low diversity, and investigate the consequences of various levels of diversity on microbial populations.

CRISPR-induced distributed immunity in microbial populations

In bacteria and archaea, viruses are the primary infectious agents, acting as virulent, often deadly pathogens. A form of adaptive immune defense known as CRISPR-Cas enables microbial cells to acquire immunity to viral pathogens by recognizing specific sequences encoded in viral genomes. The unique biology of this system results in evolutionary dynamics of host and viral diversity that cannot be fully explained by the traditional models used to describe microbe-virus coevolutionary dynamics. Here, we show how the CRISPR-mediated adaptive immune response of hosts to invading viruses facilitates the emergence of an evolutionary mode we call distributed immunity - the coexistence of multiple, equally-fit immune alleles among individuals in a microbial population. We use an eco-evolutionary modeling framework to quantify distributed immunity and demonstrate how it emerges and fluctuates in multi-strain communities of hosts and viruses as a consequence of CRISPR-induced coevolution under conditions of low viral mutation and high relative numbers of viral protospacers. We demonstrate that distributed immunity promotes sustained diversity and stability in host communities and decreased viral population density that can lead to viral extinction. We analyze sequence diversity of experimentally coevolving populations of Streptococcus thermophilus and their viruses where CRISPR-Cas is active, and find the rapid emergence of distributed immunity in the host population, demonstrating the importance of this emergent phenomenon in evolving microbial communities.

Accurate computational prediction of the transcribed strand of CRISPR non-coding RNAs

MOTIVATION

CRISPR RNAs (crRNAs) are a type of small non-coding RNA that form a key part of an acquired immune system in prokaryotes. Specific prediction methods find crRNA-encoding loci in nearly half of sequenced bacterial, and three quarters of archaeal, species. These Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) arrays consist of repeat elements alternating with specific spacers. Generally one strand is transcribed, producing long pre-crRNAs, which are processed to short crRNAs that base pair with invading nucleic acids to facilitate their destruction. No current software for the discovery of CRISPR loci predicts the direction of crRNA transcription.

RESULTS

We have developed an algorithm that accurately predicts the strand of the resulting crRNAs. The method uses as input CRISPR repeat predictions. CRISPRDirection uses parameters that are calculated from the CRISPR repeat predictions and flanking sequences, which are combined by weighted voting. The prediction may use prior coding sequence annotation but this is not required. CRISPRDirection correctly predicted the orientation of 94% of a reference set of arrays.

AVAILABILITY AND IMPLEMENTATION

The Perl source code is freely available from http://bioanalysis.otago.ac.nz/CRISPRDirection.

Evolutionary dynamics of the prokaryotic adaptive immunity system CRISPR-Cas in an explicit ecological context

A stochastic, agent-based mathematical model of the coevolution of the archaeal and bacterial adaptive immunity system, CRISPR-Cas, and lytic viruses shows that CRISPR-Cas immunity can stabilize the virus-host coexistence rather than leading to the extinction of the virus. In the model, CRISPR-Cas immunity does not specifically promote viral diversity, presumably because the selection pressure on each single proto-spacer is too weak. However, the overall virus diversity in the presence of CRISPR-Cas grows due to the increase of the host and, accordingly, the virus population size. Above a threshold value of total viral diversity, which is proportional to the viral mutation rate and population size, the CRISPR-Cas system becomes ineffective and is lost due to the associated fitness cost. Our previous modeling study has suggested that the ubiquity of CRISPR-Cas in hyperthermophiles, which contrasts its comparative low prevalence in mesophiles, is due to lower rates of mutation fixation in thermal habitats. The present findings offer a complementary, simpler perspective on this contrast through the larger population sizes of mesophiles compared to hyperthermophiles, because of which CRISPR-Cas can become ineffective in mesophiles. The efficacy of CRISPR-Cas sharply increases with the number of proto-spacers per viral genome, potentially explaining the low information content of the proto-spacer-associated motif (PAM) that is required for spacer acquisition by CRISPR-Cas because a higher specificity would restrict the number of spacers available to CRISPR-Cas, thus hampering immunity. The very existence of the PAM might reflect the tradeoff between the requirement of diverse spacers for efficient immunity and avoidance of autoimmunity.

Reassortment of CRISPR repeat-spacer loci in Sulfolobus islandicus

Virus-host interactions are a key factor shaping population dynamics of microbial species. The CRISPR-Cas adaptive immune system confers sequence-specific immunity to viral infection and has the potential to dramatically shape coevolutionary interactions between viruses and their microbial hosts. To assess evolutionary dynamics of CRISPR loci, we have sampled a population of closely related Sulfolobus islandicus strains from Kamchatka, Russia at two time points, 10 years apart. Sequence analysis of the conserved trailer sequences reveals that alleles are reassorted among three CRISPR spacer loci into combinatorial genotypes. Reassortment provides the evolutionary independence of CRISPR loci from one another as demonstrated by the differential change in allele frequencies between two time points. Genome sequences of 12 strains from this population also reveal very recent horizontal gene transfer of novel, divergent cas gene cassettes. The evolutionary independence of CRISPR loci from each other and of the cas genes that control their function are consistent with the evolutionary expectation that reassortment increases the efficiency of adaptation at these loci that are likely under strong selection by lytic viruses.

Physical model of the immune response of bacteria against bacteriophage through the adaptive CRISPR-Cas immune system

Bacteria and archaea have evolved an adaptive, heritable immune system that recognizes and protects against viruses or plasmids. This system, known as the CRISPR-Cas system, allows the host to recognize and incorporate short foreign DNA or RNA sequences, called 'spacers' into its CRISPR system. Spacers in the CRISPR system provide a record of the history of bacteria and phage coevolution. We use a physical model to study the dynamics of this coevolution as it evolves stochastically over time. We focus on the impact of mutation and recombination on bacteria and phage evolution and evasion. We discuss the effect of different spacer deletion mechanisms on the coevolutionary dynamics. We make predictions about bacteria and phage population growth, spacer diversity within the CRISPR locus, and spacer protection against the phage population.

CRISPR-mediated adaptive immune systems in bacteria and archaea

Effective clearance of an infection requires that the immune system rapidly detects and neutralizes invading parasites while strictly avoiding self-antigens that would result in autoimmunity. The cellular machinery and complex signaling pathways that coordinate an effective immune response have generally been considered properties of the eukaryotic immune system. However, a surprisingly sophisticated adaptive immune system that relies on small RNAs for sequence-specific targeting of foreign nucleic acids was recently discovered in bacteria and archaea. Molecular vaccination in prokaryotes is achieved by integrating short fragments of foreign nucleic acids into a repetitive locus in the host chromosome known as a CRISPR (clustered regularly interspaced short palindromic repeat). Here we review the mechanisms of CRISPR-mediated immunity and discuss the ecological and evolutionary implications of these adaptive defense systems.

Viral diversity threshold for adaptive immunity in prokaryotes

Bacteria and archaea face continual onslaughts of rapidly diversifying viruses and plasmids. Many prokaryotes maintain adaptive immune systems known as clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated genes (Cas). CRISPR-Cas systems are genomic sensors that serially acquire viral and plasmid DNA fragments (spacers) that are utilized to target and cleave matching viral and plasmid DNA in subsequent genomic invasions, offering critical immunological memory. Only 50% of sequenced bacteria possess CRISPR-Cas immunity, in contrast to over 90% of sequenced archaea. To probe why half of bacteria lack CRISPR-Cas immunity, we combined comparative genomics and mathematical modeling. Analysis of hundreds of diverse prokaryotic genomes shows that CRISPR-Cas systems are substantially more prevalent in thermophiles than in mesophiles. With sequenced bacteria disproportionately mesophilic and sequenced archaea mostly thermophilic, the presence of CRISPR-Cas appears to depend more on environmental temperature than on bacterial-archaeal taxonomy. Mutation rates are typically severalfold higher in mesophilic prokaryotes than in thermophilic prokaryotes. To quantitatively test whether accelerated viral mutation leads microbes to lose CRISPR-Cas systems, we developed a stochastic model of virus-CRISPR coevolution. The model competes CRISPR-Cas-positive (CRISPR-Cas+) prokaryotes against CRISPR-Cas-negative (CRISPR-Cas−) prokaryotes, continually weighing the antiviral benefits conferred by CRISPR-Cas immunity against its fitness costs. Tracking this cost-benefit analysis across parameter space reveals viral mutation rate thresholds beyond which CRISPR-Cas cannot provide sufficient immunity and is purged from host populations. These results offer a simple, testable viral diversity hypothesis to explain why mesophilic bacteria disproportionately lack CRISPR-Cas immunity. More generally, fundamental limits on the adaptability of biological sensors (Lamarckian evolution) are predicted.

A remarkable recent discovery in microbiology is that bacteria and archaea possess systems conferring immunological memory and adaptive immunity. Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated genes (CRISPR-Cas) are genomic sensors that allow prokaryotes to acquire DNA fragments from invading viruses and plasmids. Providing immunological memory, these stored fragments destroy matching DNA in future viral and plasmid invasions. CRISPR-Cas systems also provide adaptive immunity, keeping up with mutating viruses and plasmids by continually acquiring new DNA fragments. Surprisingly, less than 50% of mesophilic bacteria, in contrast to almost 90% of thermophilic bacteria and Archaea, maintain CRISPR-Cas immunity. Using mathematical modeling, we probe this dichotomy, showing how increased viral mutation rates can explain the reduced prevalence of CRISPR-Cas systems in mesophiles. Rapidly mutating viruses outrun CRISPR-Cas immune systems, likely decreasing their prevalence in bacterial populations. Thus, viral adaptability may select against, rather than for, immune adaptability in prokaryotes.

Comparative genomic structures of Mycobacterium CRISPR-Cas

Clustered regularly interspaced short palindromic repeats (CRISPR) are inheritable genetic elements of many archaea and bacteria, conferring acquired immunity against invading nucleic acids. CRISPR might be indicative of the bacterial niche adaptation and evolutionary. Mycobacterium is an important genus occupying diverse niches with profound medical and environmental significance. To present a comparative genomic landscape of the Mycobacterium CRISPR, the feature of mycobacterium CRISPR structures with sequenced complete genomes were bioinformatically analyzed. The results show that CRISPR structures can be found among 14 mycobacteria, and all loci are chromosomally located. Long CRISPRs present in three species, namely M. tuberculosis, M. bovis, and M. avium. Integrated CRISPR-Cas system can only be found in M. tuberculosis and M. bovis, with highly conserved repeat sequences, very short leaders, and promoterless. M. tuberculosis and M. bovis repeat sequences cannot form stable RNA secondary structure, consistent with a Cas6-binding sequence. M. avium repeat sequences can form classical stem-loop structure. A three-step model of M. tuberculosis CRISPR-Cas system action was put forward based on the composition and function of cas genes cluster. M. tuberculosis and M. bovis CRISPRs might interfere with the invading nucleic acids, but have somehow lost the capacity to incorporate new spacers and co-evolve with corresponding mycobacteriophages.

Spatial structure and Lamarckian adaptation explain extreme genetic diversity at CRISPR locus

Even within similar bacterial strains, it has been found that the clustered, regularly interspaced short palindromic repeat (CRISPR) shows a large variability of spacers. Modeling bacterial strains with different levels of immunity to infection by a single virulent phage, we find that coexistence in a well-mixed environment is possible only when these levels are distinctly different. When bacterial strains are similar, one subpopulation collapses. In the case of bacteria with various levels of CRISPR immunity to a range of phages, small differences in spacer composition will accordingly be suppressed under well-mixed conditions. Using a numerical model of populations spreading in space, we predict that it is the Lamarckian nature of CRISPR evolution that combines with spatial correlations to sustain the experimentally observed distribution of spacer diversity.

Evolution of animal Piwi-interacting RNAs and prokaryotic CRISPRs

Piwi-interacting RNAs (piRNAs) and CRISPR RNAs (crRNAs) are two recently discovered classes of small noncoding RNA that are found in animals and prokaryotes, respectively. Both of these novel RNA species function as components of adaptive immune systems that protect their hosts from foreign nucleic acids-piRNAs repress transposable elements in animal germlines, whereas crRNAs protect their bacterial hosts from phage and plasmids. The piRNA and CRISPR systems are nonhomologous but rather have independently evolved into logically similar defense mechanisms based on the specificity of targeting via nucleic acid base complementarity. Here we review what is known about the piRNA and CRISPR systems with a focus on comparing their evolutionary properties. In particular, we highlight the importance of several factors on the pattern of piRNA and CRISPR evolution, including the population genetic environment, the role of alternate defense systems and the mechanisms of acquisition of new piRNAs and CRISPRs.

Persisting viral sequences shape microbial CRISPR-based immunity

Well-studied innate immune systems exist throughout bacteria and archaea, but a more recently discovered genomic locus may offer prokaryotes surprising immunological adaptability. Mediated by a cassette-like genomic locus termed Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), the microbial adaptive immune system differs from its eukaryotic immune analogues by incorporating new immunities unidirectionally. CRISPR thus stores genomically recoverable timelines of virus-host coevolution in natural organisms refractory to laboratory cultivation. Here we combined a population genetic mathematical model of CRISPR-virus coevolution with six years of metagenomic sequencing to link the recoverable genomic dynamics of CRISPR loci to the unknown population dynamics of virus and host in natural communities. Metagenomic reconstructions in an acid-mine drainage system document CRISPR loci conserving ancestral immune elements to the base-pair across thousands of microbial generations. This 'trailer-end conservation' occurs despite rapid viral mutation and despite rapid prokaryotic genomic deletion. The trailer-ends of many reconstructed CRISPR loci are also largely identical across a population. 'Trailer-end clonality' occurs despite predictions of host immunological diversity due to negative frequency dependent selection (kill the winner dynamics). Statistical clustering and model simulations explain this lack of diversity by capturing rapid selective sweeps by highly immune CRISPR lineages. Potentially explaining 'trailer-end conservation,' we record the first example of a viral bloom overwhelming a CRISPR system. The polyclonal viruses bloom even though they share sequences previously targeted by host CRISPR loci. Simulations show how increasing random genomic deletions in CRISPR loci purges immunological controls on long-lived viral sequences, allowing polyclonal viruses to bloom and depressing host fitness. Our results thus link documented patterns of genomic conservation in CRISPR loci to an evolutionary advantage against persistent viruses. By maintaining old immunities, selection may be tuning CRISPR-mediated immunity against viruses reemerging from lysogeny or migration.

Multiscale model of CRISPR-induced coevolutionary dynamics: diversification at the interface of Lamarck and Darwin

The CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) system is a recently discovered type of adaptive immune defense in bacteria and archaea that functions via directed incorporation of viral and plasmid DNA into host genomes. Here, we introduce a multiscale model of dynamic coevolution between hosts and viruses in an ecological context that incorporates CRISPR immunity principles. We analyze the model to test whether and how CRISPR immunity induces host and viral diversification and the maintenance of many coexisting strains. We show that hosts and viruses coevolve to form highly diverse communities. We observe the punctuated replacement of existent strains, such that populations have very low similarity compared over the long term. However, in the short term, we observe evolutionary dynamics consistent with both incomplete selective sweeps of novel strains (as single strains and coalitions) and the recurrence of previously rare strains. Coalitions of multiple dominant host strains are predicted to arise because host strains can have nearly identical immune phenotypes mediated by CRISPR defense albeit with different genotypes. We close by discussing how our explicit eco-evolutionary model of CRISPR immunity can help guide efforts to understand the drivers of diversity seen in microbial communities where CRISPR systems are active.

CRISPR-Cas systems in bacteria and archaea: versatile small RNAs for adaptive defense and regulation

Bacteria and archaea have evolved defense and regulatory mechanisms to cope with various environmental stressors, including virus attack. This arsenal has been expanded by the recent discovery of the versatile CRISPR-Cas system, which has two novel features. First, the host can specifically incorporate short sequences from invading genetic elements (virus or plasmid) into a region of its genome that is distinguished by clustered regularly interspaced short palindromic repeats (CRISPRs). Second, when these sequences are transcribed and precisely processed into small RNAs, they guide a multifunctional protein complex (Cas proteins) to recognize and cleave incoming foreign genetic material. This adaptive immunity system, which uses a library of small noncoding RNAs as a potent weapon against fast-evolving viruses, is also used as a regulatory system by the host. Exciting breakthroughs in understanding the mechanisms of the CRISPR-Cas system and its potential for biotechnological applications and understanding evolutionary dynamics are discussed.

CRISPR: new horizons in phage resistance and strain identification

Bacteria have been widely used as starter cultures in the food industry, notably for the fermentation of milk into dairy products such as cheese and yogurt. Lactic acid bacteria used in food manufacturing, such as lactobacilli, lactococci, streptococci, Leuconostoc, pediococci, and bifidobacteria, are selectively formulated based on functional characteristics that provide idiosyncratic flavor and texture attributes, as well as their ability to withstand processing and manufacturing conditions. Unfortunately, given frequent viral exposure in industrial environments, starter culture selection and development rely on defense systems that provide resistance against bacteriophage predation, including restriction-modification, abortive infection, and recently discovered CRISPRs (clustered regularly interspaced short palindromic repeats). CRISPRs, together with CRISPR-associated genes (cas), form the CRISPR/Cas immune system, which provides adaptive immunity against phages and invasive genetic elements. The immunization process is based on the incorporation of short DNA sequences from virulent phages into the CRISPR locus. Subsequently, CRISPR transcripts are processed into small interfering RNAs that guide a multifunctional protein complex to recognize and cleave matching foreign DNA. Hypervariable CRISPR loci provide insights into the phage and host population dynamics, and new avenues for enhanced phage resistance and genetic typing and tagging of industrial strains.

Targeted bacterial immunity buffers phage diversity

Bacteria have evolved diverse defense mechanisms that allow them to fight viral attacks. One such mechanism, the clustered, regularly interspaced, short palindromic repeat (CRISPR) system, is an adaptive immune system consisting of genetic loci that can take up genetic material from invasive elements (viruses and plasmids) and later use them to reject the returning invaders. It remains an open question how, despite the ongoing evolution of attack and defense mechanisms, bacteria and viral phages manage to coexist. Using a simple mathematical model and a two-dimensional numerical simulation, we found that CRISPR adaptive immunity allows for robust phage-bacterium coexistence even when the number of virus species far exceeds the capacity of CRISPR-encoded genetic memory. Coexistence is predicted to be a consequence of the presence of many interdependent species that stress but do not overrun the bacterial defense system.

CRISPR inhibition of prophage acquisition in Streptococcus pyogenes

Streptococcus pyogenes, one of the major human pathogens, is a unique species since it has acquired diverse strain-specific virulence properties mainly through the acquisition of streptococcal prophages. In addition, S. pyogenes possesses clustered regularly interspaced short palindromic repeats (CRISPR)/Cas systems that can restrict horizontal gene transfer (HGT) including phage insertion. Therefore, it was of interest to examine the relationship between CRISPR and acquisition of prophages in S. pyogenes. Although two distinct CRISPR loci were found in S. pyogenes, some strains lacked CRISPR and these strains possess significantly more prophages than CRISPR harboring strains. We also found that the number of spacers of S. pyogenes CRISPR was less than for other streptococci. The demonstrated spacer contents, however, suggested that the CRISPR appear to limit phage insertions. In addition, we found a significant inverse correlation between the number of spacers and prophages in S. pyogenes. It was therefore suggested that S. pyogenes CRISPR have permitted phage insertion by lacking its own spacers. Interestingly, in two closely related S. pyogenes strains (SSI-1 and MGAS315), CRISPR activity appeared to be impaired following the insertion of phage genomes into the repeat sequences. Detailed analysis of this prophage insertion site suggested that MGAS315 is the ancestral strain of SSI-1. As a result of analysis of 35 additional streptococcal genomes, it was suggested that the influences of the CRISPR on the phage insertion vary among species even within the same genus. Our results suggested that limitations in CRISPR content could explain the characteristic acquisition of prophages and might contribute to strain-specific pathogenesis in S. pyogenes.

Clustered regularly interspaced short palindromic repeats (CRISPRs): the hallmark of an ingenious antiviral defense mechanism in prokaryotes

Many prokaryotes contain the recently discovered defense system against mobile genetic elements. This defense system contains a unique type of repetitive DNA stretches, termed Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs). CRISPRs consist of identical repeated DNA sequences (repeats), interspaced by highly variable sequences referred to as spacers. The spacers originate from either phages or plasmids and comprise the prokaryotes' ‘immunological memory’. CRISPR-associated (cas) genes encode conserved proteins that together with CRISPRs make-up the CRISPR/Cas system, responsible for defending the prokaryotic cell against invaders. CRISPR-mediated resistance has been proposed to involve three stages: (i) CRISPR-Adaptation, the invader DNA is encountered by the CRISPR/Cas machinery and an invader-derived short DNA fragment is incorporated in the CRISPR array. (ii) CRISPR-Expression, the CRISPR array is transcribed and the transcript is processed by Cas proteins. (iii) CRISPR-Interference, the invaders' nucleic acid is recognized by complementarity to the crRNA and neutralized. An application of the CRISPR/Cas system is the immunization of industry-relevant prokaryotes (or eukaryotes) against mobile-genetic invasion. In addition, the high variability of the CRISPR spacer content can be exploited for phylogenetic and evolutionary studies. Despite impressive progress during the last couple of years, the elucidation of several fundamental details will be a major challenge in future research.