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

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.

Characterization, complete genome sequencing, and CRISPR/Cas9 system-based decontamination of a novel Escherichia coli phage TR1 from fermentation substrates

Phage contamination has become a major concern for industrial bacteria, such as Escherichia coli BL21(DE3), used in fermentation processes. Herein, we report a CRISPR/Cas9 defense system-based strategy to precisely prey and degrade phage DNA to decontaminate target phages. First, we isolated a novel phage from fermentation substrates with BL21(DE3) as the host, named TR1. It showed a typical podovirus morphology with a head diameter of 51.46 ± 2.04 nm and a tail length of 9.31 ± 2.77 nm. The burst size of phage TR1 was 151 PFU/cell, suggesting its strong fecundity in the fermentation system. Additionally, whole-genome sequencing revealed that phage TR1 has a DNA genome of 44,099 bp in length with a 43.8% GC content, encoding a total of 68 open reading frames. Comparative genomics and phylogenetic analysis designated this phage to be a new species of the genus Christensenvirus. To counteract phage TR1, we employed the CRISPR/Cas9 system-based strategy and constructed two phage-resistant E. coli strains, BL21-C and BL21-T, based on conserved genes. Both EOP assays and growth curves indicated strong phage resistance of the recombinant strains, without affecting cell growth. Therefore, this study aimed to provide a resilient strategy to respond to ever-changing phages and ongoing phage-host arm race in industrial fermentation environments by the personalized design of spacers in the recombinant CRISPR/Cas system-containing plasmid. More importantly, our research sparks the use of phage defense mechanism to prevent phage contamination in extensive biotechnological applications.

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

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

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.

Viral Emerging Pathogen Evolution

The emergence of emerging and reemerging infectious diseases throughout history has been the product of the interaction of infectious agents, immunity, selective pressure, and environmental factors. The origin of the viral ancestors remains controversial, the debate remains whether the viruses existed before their host cells, evolved as molecular structures that adapted to the intracellular development or of the cells that infected or housed them as hosts. The nature of the viruses could propose a hypothetical evolutionary model based on the way they parasitize cells and the way they fuse with their host cells. So far there are no “universal” viral genes that have not allowed the hypotheses of the “common viral ancestor” to be conclusive. Insects have a special role in understanding viral evolution and how they became pathogens, suggesting a coevolutionary scenario of colonization that would initially allow the survival of viruses and later an ecological niche required for the cycle of virus replication. The coevolution of viruses and defense systems of hosts represents a fundamental role in the evolution of both viruses and host cells, a kind of arms race between viruses and hosts, and the deescalation of that war could explain the coexistence of viruses and hosts. Viruses are playing an important role in the understanding of cell evolution and have earned the right to be included in the tree of life. The aim of this chapter is to review the viral evolution and the insects as host.

Collateral damage and CRISPR genome editing

The simplicity and the versatility of clustered regularly interspaced short palindromic repeats/CRISPR-associated protein (CRISPR-Cas) systems have enabled the genetic modification of virtually every organism and offer immense therapeutic potential for the treatment of human disease. Although these systems may function efficiently within eukaryotic cells, there remain concerns about the accuracy of Cas endonuclease effectors and their use for precise gene editing. Recently, two independent reports investigating the editing accuracy of the CRISPR-Cas9 system were published by separate groups at the Wellcome Sanger Institute; our study-Iyer and colleagues [1]-defined the landscape of off-target mutations, whereas the other by Kosicki and colleagues [2] detailed the existence of on-target, potentially deleterious deletions. Although both studies found evidence of large on-target CRISPR-induced deletions, they reached seemingly very different conclusions.

CRISPRs for Strain Tracking and Their Application to Microbiota Transplantation Data Analysis

CRISPR-Cas systems are adaptive immune systems naturally found in bacteria and archaea. Prokaryotes use these immune systems to defend against invaders, which include phages, plasmids, and other mobile genetic elements. Relying on the integration of spacers derived from invader sequences (protospacers) into CRISPR loci (forming spacers flanked by repeats), CRISPR-Cas systems are able to store the memory of past immunological encounters. While CRISPR-Cas systems have evolved in response to invading mobile genetic elements, invaders have also developed mechanisms to avoid detection. As a result of an arms race between CRISPR-Cas systems and their targets, CRISPR arrays typically undergo rapid turnover of spacers through the acquisition and loss events. Additionally, microbiomes of different individuals rarely share spacers. Here, we present a computational pipeline, CRISPRtrack, for strain tracking based on CRISPR spacer content, and we applied it to fecal transplantation microbiome data to study the retention of donor strains in recipients. Our results demonstrate the potential use of CRISPRs as a simple yet effective tool for donor-strain tracking in fecal transplantation and as a general purpose tool for quantifying microbiome similarity.

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.

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.

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.

Regulated CRISPR Modules Exploit a Dual Defense Strategy of Restriction and Abortive Infection in a Model of Prokaryote-Phage Coevolution

CRISPRs offer adaptive immunity in prokaryotes by acquiring genomic fragments from infecting phage and subsequently exploiting them for phage restriction via an RNAi-like mechanism. Here, we develop and analyze a dynamical model of CRISPR-mediated prokaryote-phage coevolution that incorporates classical CRISPR kinetics along with the recently discovered infection-induced activation and autoimmunity side effects. Our analyses reveal two striking characteristics of the CRISPR defense strategy: that both restriction and abortive infections operate during coevolution with phages, driving phages to much lower densities than possible with restriction alone, and that CRISPR maintenance is determined by a key dimensionless combination of parameters, which upper bounds the activation level of CRISPRs in uninfected populations. We contrast these qualitative observations with experimental data on CRISPR kinetics, which offer insight into the spacer deletion mechanism and the observed low CRISPR prevalence in clinical isolates. More generally, we exploit numerical simulations to delineate four regimes of CRISPR dynamics in terms of its host, kinetic, and regulatory parameters.

To counteract viral infections, bacteria and archaea have evolved a variety of defense systems. These can broadly be classified into either restriction or suicide mechanisms. The former enforces nicks in the invading DNA making it unusable for production of further infectious particles; the latter, by contrast, induces cell death whereby an infected cell activates specific host suicidal pathways that are otherwise strongly repressed, thus inhibiting further infection. Examples of the former class include restriction-modification (R-M) and the recently discovered CRISPR systems, while the latter class includes a variety of toxin/anti-toxin systems. CRISPRs, in contrast to R-Ms, adapt to target viral genomes by updating the database of target sites they recognize. The adverse side effect of such a mechanism, however, is that CRISPRs can target the host genome itself resulting in undesirable cell death (autoimmunity). The recent discovery of infection-induced activation of CRISPR systems suggests that these negative side effects may be limited to periods of infection. This led us to hypothesize that such regulatory control—similar to abortive infection mechanisms—can be advantageous by limiting viral spread through suicide of infected cells. To test this hypothesis, we mathematically model CRISPR induced prokaryote-phage coevolutionary dynamics in the presence of infection-regulated CRISPR activity. Our results indicate that, except in limited growth rates, regulated CRISPRs exploit both autoimmunity and target restriction and can therefore be considered a hybrid class that leverages both restriction and suicide mechanisms to limit phage infection.

The Contribution of Genetic Recombination to CRISPR Array Evolution

CRISPR (clustered regularly interspaced short palindromic repeats) is a microbial immune system against foreign DNA. Recognition sequences (spacers) encoded within the CRISPR array mediate the immune reaction in a sequence-specific manner. The known mechanisms for the evolution of CRISPR arrays include spacer acquisition from foreign DNA elements at the time of invasion and array erosion through spacer deletion. Here, we consider the contribution of genetic recombination between homologous CRISPR arrays to the evolution of spacer repertoire. Acquisition of spacers from exogenic arrays via recombination may confer the recipient with immunity against unencountered antagonists. For this purpose, we develop a novel method for the detection of recombination in CRISPR arrays by modeling the spacer order in arrays from multiple strains from the same species. Because the evolutionary signal of spacer recombination may be similar to that of pervasive spacer deletions or independent spacer acquisition, our method entails a robustness analysis of the recombination inference by a statistical comparison to resampled and perturbed data sets. We analyze CRISPR data sets from four bacterial species: two Gammaproteobacteria species harboring CRISPR type I and two Streptococcus species harboring CRISPR type II loci. We find that CRISPR array evolution in Escherichia coli and Streptococcus agalactiae can be explained solely by vertical inheritance and differential spacer deletion. In Pseudomonas aeruginosa, we find an excess of single spacers potentially incorporated into the CRISPR locus during independent acquisition events. In Streptococcus thermophilus, evidence for spacer acquisition by recombination is present in 5 out of 70 strains. Genetic recombination has been proposed to accelerate adaptation by combining beneficial mutations that arose in independent lineages. However, for most species under study, we find that CRISPR evolution is shaped mainly by spacer acquisition and loss rather than recombination. Since the evolution of spacer content is characterized by a rapid turnover, it is likely that recombination is not beneficial for improving phage resistance in the strains under study, or that it cannot be detected in the resolution of intraspecies comparisons.

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.

Novel strategies to prevent or exploit phages in fermentations, insights from phage-host interactions

Phages infecting lactic acid bacteria (LAB) provide some of the most advanced model systems for (tailed) phage-host interactions. In particular the identification of receptor molecules of representative lactococcal phages combined with the elucidation of the structure of the receptor-binding protein has permitted crucial insights into the early stages of infection. Dairy and biotechnological fermentations are persistently marred by the destructive activities of phages. Here, we discuss how recent advances in our knowledge on LAB phage-host interactions have provided a basis for the next generation anti-phage strategies. Furthermore, the significant increase in genomic data has furthered our understanding of the genetics of these phages, thereby permitting the exploitation of phage-derived components for food safety and biotechnological applications.

Pseudo-chaotic oscillations in CRISPR-virus coevolution predicted by bifurcation analysis

BACKGROUND

The CRISPR-Cas systems of adaptive antivirus immunity are present in most archaea and many bacteria, and provide resistance to specific viruses or plasmids by inserting fragments of foreign DNA into the host genome and then utilizing transcripts of these spacers to inactivate the cognate foreign genome. The recent development of powerful genome engineering tools on the basis of CRISPR-Cas has sharply increased the interest in the diversity and evolution of these systems. Comparative genomic data indicate that during evolution of prokaryotes CRISPR-Cas loci are lost and acquired via horizontal gene transfer at high rates. Mathematical modeling and initial experimental studies of CRISPR-carrying microbes and viruses reveal complex coevolutionary dynamics.

RESULTS

We performed a bifurcation analysis of models of coevolution of viruses and microbial host that possess CRISPR-Cas hereditary adaptive immunity systems. The analyzed Malthusian and logistic models display complex, and in particular, quasi-chaotic oscillation regimes that have not been previously observed experimentally or in agent-based models of the CRISPR-mediated immunity. The key factors for the appearance of the quasi-chaotic oscillations are the non-linear dependence of the host immunity on the virus load and the partitioning of the hosts into the immune and susceptible populations, so that the system consists of three components.

CONCLUSIONS

Bifurcation analysis of CRISPR-host coevolution model predicts complex regimes including quasi-chaotic oscillations. The quasi-chaotic regimes of virus-host coevolution are likely to be biologically relevant given the evolutionary instability of the CRISPR-Cas loci revealed by comparative genomics. The results of this analysis might have implications beyond the CRISPR-Cas systems, i.e. could describe the behavior of any adaptive immunity system with a heritable component, be it genetic or epigenetic. These predictions are experimentally testable.

REVIEWERS' REPORTS

This manuscript was reviewed by Sandor Pongor, Sergei Maslov and Marek Kimmel. For the complete reports, go to the Reviewers' Reports section.

Involvement of eukaryotic small RNA pathways in host defense and viral pathogenesis

Post-transcriptional gene regulation by small RNAs is now established as an important branch of the gene regulatory system. Many different classes of small RNAs have been discovered; among these are short interfering RNAs (siRNAs) and microRNA (miRNAs). Though differences in the processing and function of small RNAs exist between plants and animals, both groups utilize small RNA-mediated gene regulation in response to pathogens. Host encoded miRNAs and siRNAs are generated from viral RNA function in host defense and pathogenic resistance in plants. In animals, miRNAs are key regulators in both immune system development and in immune function. Pathogens, in particular viruses, have evolved mechanisms to usurp the host’s small RNA-mediated regulatory system. Overall, small RNAs are a major component of host defense and immunity in eukaryotes. The goal of this review is to summarize our current knowledge of the involvement of eukaryotic small RNA pathways in host defense and viral pathogenesis.

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.