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

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.

Multistep diversification in spatiotemporal bacterial-phage coevolution

The evolutionary arms race between phages and bacteria, where bacteria evolve resistance to phages and phages retaliate with resistance-countering mutations, is a major driving force of molecular innovation and genetic diversification. Yet attempting to reproduce such ongoing retaliation dynamics in the lab has been challenging; laboratory coevolution experiments of phage and bacteria are typically performed in well-mixed environments and often lead to rapid stagnation with little genetic variability. Here, co-culturing motile E. coli with the lytic bacteriophage T7 on swimming plates, we observe complex spatiotemporal dynamics with multiple genetically diversifying adaptive cycles. Systematically quantifying over 10,000 resistance-infectivity phenotypes between evolved bacteria and phage isolates, we observe diversification into multiple coexisting ecotypes showing a complex interaction network with both host-range expansion and host-switch tradeoffs. Whole-genome sequencing of these evolved phage and bacterial isolates revealed a rich set of adaptive mutations in multiple genetic pathways including in genes not previously linked with phage-bacteria interactions. Synthetically reconstructing these new mutations, we discover phage-general and phage-specific resistance phenotypes as well as a strong synergy with the more classically known phage-resistance mutations. These results highlight the importance of spatial structure and migration for driving phage-bacteria coevolution, providing a concrete system for revealing new molecular mechanisms across diverse phage-bacterial systems.

From kill the winner to eliminate the winner in open phage-bacteria systems

Phages and bacteria manage to coexist and sustain ecosystems with a high diversity of strains, despite limited resources and heavy predation. This diversity can be explained by the “kill the winner” model where virulent phages predominantly prey on fast-growing bacteria and thereby suppress the competitive exclusion of slower-growing bacteria. Here we computationally investigate the robustness of these systems against invasions, where new phages or bacteria may interact with more than one of the resident strains. The resulting interaction networks were found to self-organize into a network with strongly interacting specialized predator-prey pairs, resembling that of the “kill the winner” model. Furthermore, the “kill the winner” dynamics is enforced with the occasional elimination of even the fastest-growing bacteria strains due to a phage infecting the fast and slow growers. The frequency of slower-growing strains was increased with the introduction of even a few non-diagonal interactions. Hence, phages capable of infecting multiple hosts play significant roles both in the evolution of the ecosystem by eliminating the winner and in supporting diversity by allowing slow growers to coexist with faster growers.

We demonstrate that in an open system of phages and bacteria with very limited resources, a bacterial strain that has a high growth rate can still be outcompeted by a slower-growing strain if they have a common phage. The impact of this on ecosystem structure is significant as soon as there is a small probability to have a common phage among bacterial strains. Furthermore, by analysing the structure of the interaction network we show that it self-organizes into a network with strongly interacting specialized predator-prey pairs, in order to reduce phages competition. Nevertheless, the presence of the remaining links is very important for the network dynamics since even a few of them significantly enhance the frequency of slower-growing strains.

Persistence of plasmids targeted by CRISPR interference in bacterial populations

RNA-guided CRISPR-Cas nucleases efficiently protect bacterial cells from phage infection and plasmid transformation. Yet, the efficiency of CRISPR-Cas defense is not absolute. Mutations in either CRISPR-Cas components of the host or mobile genetic elements regions targeted by CRISPR-Cas inactivate the defensive action. Here, we show that even at conditions of active CRISPR-Cas and unaltered targeted plasmids, a kinetic equilibrium between CRISPR-Cas nucleases action and plasmid replication processes allows for existence of a small subpopulation of plasmid-bearing cells on the background of cells that have been cured from the plasmid. In nature, the observed diversification of phenotypes may allow rapid changes in the population structure to meet the demands of the environment.

CRISPR-Cas systems provide prokaryotes with an RNA-guided defense against foreign mobile genetic elements (MGEs) such as plasmids and viruses. A common mechanism by which MGEs avoid interference by CRISPR consists of acquisition of escape mutations in regions targeted by CRISPR. Here, using microbiological, live microscopy and microfluidics analyses we demonstrate that plasmids can persist for multiple generations in some Escherichia coli cell lineages at conditions of continuous targeting by the type I-E CRISPR-Cas system. We used mathematical modeling to show how plasmid persistence in a subpopulation of cells mounting CRISPR interference is achieved due to the stochastic nature of CRISPR interference and plasmid replication events. We hypothesize that the observed complex dynamics provides bacterial populations with long-term benefits due to continuous maintenance of mobile genetic elements in some cells, which leads to diversification of phenotypes in the entire community and allows rapid changes in the population structure to meet the demands of a changing environment.

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.

Unveiling Human Non-Random Genome Editing Mechanisms Activated in Response to Chronic Environmental Changes: I. Where Might These Mechanisms Come from and What Might They Have Led To?

This article challenges the notion of the randomness of mutations in eukaryotic cells by unveiling stress-induced human non-random genome editing mechanisms. To account for the existence of such mechanisms, I have developed molecular concepts of the cell environment and cell environmental stressors and, making use of a large quantity of published data, hypothesised the origin of some crucial biological leaps along the evolutionary path of life on Earth under the pressure of natural selection, in particular, (1) virus-cell mating as a primordial form of sexual recombination and symbiosis; (2) Lamarckian CRISPR-Cas systems; (3) eukaryotic gene development; (4) antiviral activity of retrotransposon-guided mutagenic enzymes; and finally, (5) the exaptation of antiviral mutagenic mechanisms to stress-induced genome editing mechanisms directed at "hyper-transcribed" endogenous genes. Genes transcribed at their maximum rate (hyper-transcribed), yet still unable to meet new chronic environmental demands generated by "pollution", are inadequate and generate more and more intronic retrotransposon transcripts. In this scenario, RNA-guided mutagenic enzymes (e.g., Apolipoprotein B mRNA editing catalytic polypeptide-like enzymes, APOBECs), which have been shown to bind to retrotransposon RNA-repetitive sequences, would be surgically targeted by intronic retrotransposons on opened chromatin regions of the same "hyper-transcribed" genes. RNA-guided mutagenic enzymes may therefore "Lamarkianly" generate single nucleotide polymorphisms (SNP) and gene copy number variations (CNV), as well as transposon transposition and chromosomal translocations in the restricted areas of hyper-functional and inadequate genes, leaving intact the rest of the genome. CNV and SNP of hyper-transcribed genes may allow cells to surgically explore a new fitness scenario, which increases their adaptability to stressful environmental conditions. Like the mechanisms of immunoglobulin somatic hypermutation, non-random genome editing mechanisms may generate several cell mutants, and those codifying for the most environmentally adequate proteins would have a survival advantage and would therefore be Darwinianly selected. Non-random genome editing mechanisms represent tools of evolvability leading to organismal adaptation including transgenerational non-Mendelian gene transmission or to death of environmentally inadequate genomes. They are a link between environmental changes and biological novelty and plasticity, finally providing a molecular basis to reconcile gene-centred and "ecological" views of evolution.

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’.

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.

Evolution of the U.S. Biological Select Agent Rathayibacter toxicus

Rathayibacter toxicus is a toxin-producing species found in Australia and is often fatal to grazing animals. The threat of introduction of the species into the United States led to its inclusion in the Federal Select Agent Program, which makes R. toxicus a highly regulated species. This work provides novel insights into the evolution of R. toxicus. R. toxicus is the only species in the genus to have acquired a CRISPR adaptive immune system to protect against bacteriophages. Results suggest that coexistence with the bacteriophage NCPPB3778 led to the massive shrinkage of the R. toxicus genome, species divergence, and the maintenance of low genetic diversity in extant bacterial groups. This work contributes to an understanding of the evolution and ecology of an agriculturally important species of bacteria.

Rathayibacter toxicus is a species of Gram-positive, corynetoxin-producing bacteria that causes annual ryegrass toxicity, a disease often fatal to grazing animals. A phylogenomic approach was employed to model the evolution of R. toxicus to explain the low genetic diversity observed among isolates collected during a 30-year period of sampling in three regions of Australia, gain insight into the taxonomy of Rathayibacter, and provide a framework for studying these bacteria. Analyses of a data set of more than 100 sequenced Rathayibacter genomes indicated that Rathayibacter forms nine species-level groups. R. toxicus is the most genetically distant, and evidence suggested that this species experienced a dramatic event in its evolution. Its genome is significantly reduced in size but is colinear to those of sister species. Moreover, R. toxicus has low intergroup genomic diversity and almost no intragroup genomic diversity between ecologically separated isolates. R. toxicus is the only species of the genus that encodes a clustered regularly interspaced short palindromic repeat (CRISPR) locus and that is known to host a bacteriophage parasite. The spacers, which represent a chronological history of infections, were characterized for information on past events. We propose a three-stage process that emphasizes the importance of the bacteriophage and CRISPR in the genome reduction and low genetic diversity of the R. toxicus species.

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.

Theory of invasion extinction dynamics in minimal food webs

When food webs are exposed to species invasion, secondary extinction cascades may be set off. Although much work has gone into characterizing the structure of food webs, systematic predictions on their evolutionary dynamics are still scarce. Here we present a theoretical framework that predicts extinctions in terms of an alternating sequence of two basic processes: resource depletion by or competitive exclusion between consumers. We first propose a conceptual invasion extinction model (IEM) involving random fitness coefficients. We bolster this IEM by an analytical, recursive procedure for calculating idealized extinction cascades after any species addition and simulate the long-time evolution. Our procedure describes minimal food webs where each species interacts with only a single resource through the generalized Lotka-Volterra equations. For such food webs ex- tinction cascades are determined uniquely and the system always relaxes to a stable steady state. The dynamics and scale invariant species life time resemble the behavior of the IEM, and correctly predict an upper limit for trophic levels as observed in the field.

Existence and construction of large stable food webs

Ecological diversity is ubiquitous despite the restrictions imposed by competitive exclusion and apparent competition. To explain the observed richness of species in a given habitat, food-web theory has explored nonlinear functional responses, self-interaction, or spatial structure and dispersal-model ingredients that have proven to promote stability and diversity. We return instead here to classical Lotka-Volterra equations, where species-species interaction is characterized by a simple product and spatial restrictions are ignored. We quantify how this idealization imposes constraints on coexistence and diversity for many species. To this end, we introduce the concept of free and controlled species and use this to demonstrate how stable food webs can be constructed by the sequential addition of species. The resulting food webs can reach dozens of species and generally yield nonrandom degree distributions in accordance with the constraints imposed through the assembly process. Our model thus serves as a formal starting point for the study of sustainable interaction patterns between species.

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.

Models of life: epigenetics, diversity and cycles

This review emphasizes aspects of biology that can be understood through repeated applications of simple causal rules. The selected topics include perspectives on gene regulation, phage lambda development, epigenetics, microbial ecology, as well as model approaches to diversity and to punctuated equilibrium in evolution. Two outstanding features are repeatedly described. One is the minimal number of rules to sustain specific states of complex systems for a long time. The other is the collapse of such states and the subsequent dynamical cycle of situations that restitute the system to a potentially new metastable state.

Population cycles and species diversity in dynamic Kill-the-Winner model of microbial ecosystems

Determinants of species diversity in microbial ecosystems remain poorly understood. Bacteriophages are believed to increase the diversity by the virtue of Kill-the-Winner infection bias preventing the fastest growing organism from taking over the community. Phage-bacterial ecosystems are traditionally described in terms of the static equilibrium state of Lotka-Volterra equations in which bacterial growth is exactly balanced by losses due to phage predation. Here we consider a more dynamic scenario in which phage infections give rise to abrupt and severe collapses of bacterial populations whenever they become sufficiently large. As a consequence, each bacterial population in our model follows cyclic dynamics of exponential growth interrupted by sudden declines. The total population of all species fluctuates around the carrying capacity of the environment, making these cycles cryptic. While a subset of the slowest growing species in our model is always driven towards extinction, in general the overall ecosystem diversity remains high. The number of surviving species is inversely proportional to the variation in their growth rates but increases with the frequency and severity of phage-induced collapses. Thus counter-intuitively we predict that microbial communities exposed to more violent perturbations should have higher diversity.

Lamarckian Illusions

In recent years the term 'Lamarckian evolution' has become a household name for processes that do not follow classical Mendelian pattern of inheritance, and it is seen as a relevant complement to Darwinism. In this article I argue that bringing back Lamarck is unjustified and misleading.

Resistance and tolerance to foreign elements by prokaryotic immune systems - curating the genome

To engage in adaptive symbioses or genetic exchange, organisms must interact with foreign, non-self elements despite the risks of predation and parasitism. By surveying the interface between self and non-self, immune systems can help ensure the benevolence of these interactions without isolating their hosts altogether. In this Essay, we examine prokaryotic restriction-modification and CRISPR-Cas (clustered, regularly interspaced palindromic repeat-CRISPR-associated proteins) activities and discuss their analogy to mammalian immune pathways. We further explain how their capacities for resistance and tolerance are optimized to reduce parasitism and immunopathology during encounters with non-self.

Restriction modification systems as engines of diversity

Restriction modification (RM) systems provide protection against a broad spectrum of phages. However, the likelihood of a phage permanently bypassing this can be as high as 0.1 per infection (Korona et al., 1993) which makes for a relatively weak defense. Here we argue that, apart from providing such transient defenses, RM systems can facilitate long-term coexistence of many bacterial strains. We show that this diversity can be as large as the burst size of the phage but no larger-a curious correspondence between a number at the level of species and another number at the level of individuals. Such a highly diverse and stably coexisting ecosystem is robust to substantial variation in both bacterial growth rates and strength of their RM systems, which might be one reason why quite weak RM systems exist in the wild.

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.

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.

CRISPR-Cas systems: Prokaryotes upgrade to adaptive immunity

Clustered regularly interspaced short palindromic repeats (CRISPR), and associated proteins (Cas) comprise the CRISPR-Cas system, which confers adaptive immunity against exogenic elements in many bacteria and most archaea. CRISPR-mediated immunization occurs through the uptake of DNA from invasive genetic elements such as plasmids and viruses, followed by its integration into CRISPR loci. These loci are subsequently transcribed and processed into small interfering RNAs that guide nucleases for specific cleavage of complementary sequences. Conceptually, CRISPR-Cas shares functional features with the mammalian adaptive immune system, while also exhibiting characteristics of Lamarckian evolution. Because immune markers spliced from exogenous agents are integrated iteratively in CRISPR loci, they constitute a genetic record of vaccination events and reflect environmental conditions and changes over time. Cas endonucleases, which can be reprogrammed by small guide RNAs have shown unprecedented potential and flexibility for genome editing and can be repurposed for numerous DNA targeting applications including transcriptional control.

Phage and bacteria support mutual diversity in a narrowing staircase of coexistence

The competitive exclusion principle states that phage diversity M should not exceed bacterial diversity N. By analyzing the steady-state solutions of multistrain equations, we find a new constraint: the diversity N of bacteria living on the same resources is constrained to be M or M+1 in terms of the diversity of their phage predators. We quantify how the parameter space of coexistence exponentially decreases with diversity. For diversity to grow, an open or evolving ecosystem needs to climb a narrowing 'diversity staircase' by alternatingly adding new bacteria and phages. The unfolding coevolutionary arms race will typically favor high growth rate, but a phage that infects two bacterial strains differently can occasionally eliminate the fastest growing bacteria. This context-dependent fitness allows abrupt resetting of the 'Red-Queen's race' and constrains the local diversity.

Speciation, diversification, and coexistence of sessile species that compete for space

Speciation, diversification, and competition between species challenge the stability of complex ecosystems. Laboratory experiments often focus on one or two species competing under conditions where they may grow exponentially. Field studies, in contrast, emphasize multi-species communities characterized by many types of ecological interactions. A general problem is to understand conditions that support a dynamically maintained coexistence of many species in an ecosystem over a long time span. In the present paper we propose a lattice model of multiple competing and evolving sessile species. When allowing the interspecies interactions to mutate, we obtain coexistence of many species in a complex ecosystem, provided that there is a cost for each interaction. The diversity reached by the model incorporating speciation is found to be substantially higher than in the case when entirely new species appear due to immigration from outside of the considered ecosystem. The species self-organize their spatial distribution through competitive interactions to create many patches, implicitly protecting each other from competitively superior species, and speciation in each patch leads the system to high diversity. We also show that species that exist a long time tend to have a relatively small population, as this allows them to avoid encounter with competitive invaders.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

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.

Mechanisms of multi-strain coexistence in host-phage systems with nested infection networks

Bacteria and their viruses (bacteriophages) coexist in natural environments forming complex infection networks. Recent empirical findings suggest that phage-bacteria infection networks often possess a nested structure such that there is a hierarchical relationship among who can infect whom. Here we consider how nested infection networks may affect phage and bacteria dynamics using a multi-type Lotka-Volterra framework with cross-infection. Analysis of similar models has, in the past, assumed simpler interaction structures as a first step towards tractability. We solve the proposed model, finding trade-off conditions on the life-history traits of both bacteria and viruses that allow coexistence in communities with nested infection networks. First, we find that bacterial growth rate should decrease with increasing defense against infection. Second, we find that the efficiency of viral infection should decrease with host range. Next, we establish a relationship between relative densities and the curvature of life history trade-offs. We compare and contrast the current findings to the "Kill-the-Winner" model of multi-species phage-bacteria communities. Finally, we discuss a suite of testable hypotheses stemming from the current model concerning relationships between infection range, life history traits and coexistence in complex phage-bacteria communities.

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.

Probabilistic models for CRISPR spacer content evolution

BACKGROUND

The CRISPR/Cas system is known to act as an adaptive and heritable immune system in Eubacteria and Archaea. Immunity is encoded in an array of spacer sequences. Each spacer can provide specific immunity to invasive elements that carry the same or a similar sequence. Even in closely related strains, spacer content is very dynamic and evolves quickly. Standard models of nucleotide evolution cannot be applied to quantify its rate of change since processes other than single nucleotide changes determine its evolution.

METHODS

We present probabilistic models that are specific for spacer content evolution. They account for the different processes of insertion and deletion. Insertions can be constrained to occur on one end only or are allowed to occur throughout the array. One deletion event can affect one spacer or a whole fragment of adjacent spacers. Parameters of the underlying models are estimated for a pair of arrays by maximum likelihood using explicit ancestor enumeration.

RESULTS

Simulations show that parameters are well estimated on average under the models presented here. There is a bias in the rate estimation when including fragment deletions. The models also estimate times between pairs of strains. But with increasing time, spacer overlap goes to zero, and thus there is an upper bound on the distance that can be estimated. Spacer content similarities are displayed in a distance based phylogeny using the estimated times.We use the presented models to analyze different Yersinia pestis data sets and find that the results among them are largely congruent. The models also capture the variation in diversity of spacers among the data sets. A comparison of spacer-based phylogenies and Cas gene phylogenies shows that they resolve very different time scales for this data set.

CONCLUSIONS

The simulations and data analyses show that the presented models are useful for quantifying spacer content evolution and for displaying spacer content similarities of closely related strains in a phylogeny. This allows for comparisons of different CRISPR arrays or for comparisons between CRISPR arrays and nucleotide substitution rates.

Adding a cost of resistance description extends the ability of virus-host model to explain observed patterns in structure and function of pelagic microbial communities

By adding a generic description of cost of resistance (COR) to the existing 'killing-the-winner' model, we show how this expands the model's explanatory power to include rank-abundance relationships in the host population. The model can predict a counter-intuitive relationship previously suggested in the literature, where abundant viruses are associated with rare hosts and vice versa. The model explains the observed dominance of slow-growing prokaryotes as the result of successful defence strategies, rather than as dormancy of hosts lacking essential substrates. In addition to these important conceptual aspects, the model is able to reproduce realistic values for virus : host ratios and partitioning of bacterial production between predatory loss and viral lysis. A high COR is also shown to increase the community's richness and Shannon diversity index. This model thus not only couples life strategies at the cellular level with system properties, but it also links the two system level properties of biogeochemical flows and diversity to each other. The model operates with host groups, and consequences for biodiversity when interpreting these groups in terms of species and strains are discussed.

Spreading of multiple epidemics with cross immunization

Pathogen-host relationships are the result of an ongoing coevolutionary race where the immune system of the host attempts to eliminate the pathogen, while the successful pathogen mutates to become invisible for the host's immune system. We here propose a minimal pathogen-host evolution model that takes into account cross immunization and allows for evolution of a spatially heterogeneous immune status of a population of hosts. With only the mutation rate as a determining parameter, the model allows us to produce an evolutionary tree of diseases which is highly branched, but hardly ever splits into separate long-lived trunks. Side branches remain short lived and seldom diverge to the extent of losing all cross immunizations.