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

Do-or-die life cycles and diverse post-infection resistance mechanisms limit the evolution of parasite host ranges

In light of the dynamic nature of parasite host ranges and documented potential for rapid host shifts, the observed high host specificity of most parasites remains an ecological paradox. Different variants of host-use trade-offs have become a mainstay of theoretical explanations of the prevalence of host specialism, but empirical evidence for such trade-offs is rare. We propose an alternative theory based on basic features of the parasite life cycle: host selection and subsequent intrahost replication. We introduce a new concept of effective burst size that accounts for the fact that successful host selection does not guarantee intrahost replication. Our theory makes a general prediction that a parasite will expand its host range if its effective burst size is positive. An in silico model of bacteria-phage coevolution verifies our predictions and demonstrates that the tendency for relatively narrow host ranges in parasites can be explained even in the absence of trade-offs.

The evolution of resistance against good and bad infections

Opportunities for genetic exchange are abundant between bacteria and foreign genetic elements (FGEs) such as conjugative plasmids, transposable elements and bacteriophages. The genetic novelty that may arise from these forms of genetic exchange is potentially beneficial to bacterial hosts, but there are also potential costs, which may be considerable in the case of phage infection. Some bacterial resistance mechanisms target both beneficial and deleterious forms of genetic exchange. Using a general epidemiological model, we explored under which conditions such resistance mechanisms may evolve. We considered a population of hosts that may be infected by FGEs that either confer a benefit or are deleterious to host fitness, and we analysed the epidemiological and evolutionary outcomes of resistance evolving under different cost/benefit scenarios. We show that the degree of co-infection between these two types of infection is particularly important in determining the evolutionarily stable level of host resistance. We explore these results using the example of CRISPR-Cas, a form of bacterial immunity that targets a variety of FGEs, and we show the potential role of bacteriophage infection in selecting for resistance mechanisms that in turn limit the acquisition of plasmid-borne antibiotic resistance. Finally, beyond microbes, we discuss how endosymbiotic associations may have shaped the evolution of host immune responses to pathogens.

Dealing with the evolutionary downside of CRISPR immunity: bacteria and beneficial plasmids

The immune systems that protect organisms from infectious agents invariably have a cost for the host. In bacteria and archaea CRISPR-Cas loci can serve as adaptive immune systems that protect these microbes from infectiously transmitted DNAs. When those DNAs are borne by lytic viruses (phages), this protection can provide a considerable advantage. CRISPR-Cas immunity can also prevent cells from acquiring plasmids and free DNA bearing genes that increase their fitness. Here, we use a combination of experiments and mathematical-computer simulation models to explore this downside of CRISPR-Cas immunity and its implications for the maintenance of CRISPR-Cas loci in microbial populations. We analyzed the conjugational transfer of the staphylococcal plasmid pG0400 into Staphylococcus epidermidis RP62a recipients that bear a CRISPR-Cas locus targeting this plasmid. Contrary to what is anticipated for lytic phages, which evade CRISPR by mutations in the target region, the evasion of CRISPR immunity by plasmids occurs at the level of the host through loss of functional CRISPR-Cas immunity. The results of our experiments and models indicate that more than 10⁻⁴ of the cells in CRISPR-Cas positive populations are defective or deleted for the CRISPR-Cas region and thereby able to receive and carry the plasmid. Most intriguingly, the loss of CRISPR function even by large deletions can have little or no fitness cost in vitro. These theoretical and experimental results can account for the considerable variation in the existence, number and function of CRISPR-Cas loci within and between bacterial species. We postulate that as a consequence of the opposing positive and negative selection for immunity, CRISPR-Cas systems are in a continuous state of flux. They are lost when they bear immunity to laterally transferred beneficial genes, re-acquired by horizontal gene transfer, and ascend in environments where phage are a major source of mortality.

In addition to the virtue of protecting archaea and bacteria from the ravages of lethal viruses (phage), the immunity generated by the CRISPR-Cas systems have an evolutionary downside; they can prevent the acquisition of genes and genetic elements required for the adaptation and even the survival of these microbes. Using mathematical models and experiments with Staphylococcus epidermidis and the staphylococcal conjugative plasmid pG0400, we explore how bacteria deal with this evolutionary downside of CRISPR-Cas immunity. Although there are mechanisms by which immune populations of bacteria can acquire essential plasmids without the loss of CRISPR-Cas immunity, the results of our conjugation and fitness cost experiments suggest the most likely mechanism is the deactivation and deletion of this region. These results provide an explanation for the considerable variation in the existence, number and function of CRISPR-Cas within and between species of microbes. Along with other observations our work also suggests that the CRISPR-Cas loci are in a continuous state of flux: acquired by horizontal gene transfer, ascend when populations are confronted with phage and are rapidly lost when infectiously transmitted genes and genetic elements are required for the adaptation and survival of the population.

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.

Holding a grudge: persisting anti-phage CRISPR immunity in multiple human gut microbiomes

The CRISPR (clustered regularly interspaced short palindromic repeats)/Cas (CRISPR-associated) system of bacteria and archaea constitutes a mechanism of acquired adaptive immunity against phages, which is based on genome-encoded markers of previously infecting phage sequences ("spacers"). As a repository of phage sequences, these spacers make the system particularly suitable for elucidating phage-bacteria interactions in metagenomic studies. Recent metagenomic analyses of CRISPRs associated with the human microbiome intriguingly revealed conserved "memory spacers" shared by bacteria in multiple unrelated, geographically separated individuals. Here, we discuss possible avenues for explaining this phenomenon by integrating insights from CRISPR biology and phage-bacteria ecology, with a special focus on the human gut. We further explore the growing body of evidence for the role of CRISPR/Cas in regulating the interplay between bacteria and lysogenic phages, which may be intimately related to the presence of memory spacers and sheds new light on the multifaceted biological and ecological modes of action of CRISPR/Cas.

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.

The population and evolutionary dynamics of phage and bacteria with CRISPR-mediated immunity

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), together with associated genes (cas), form the CRISPR–cas adaptive immune system, which can provide resistance to viruses and plasmids in bacteria and archaea. Here, we use mathematical models, population dynamic experiments, and DNA sequence analyses to investigate the host–phage interactions in a model CRISPR–cas system, Streptococcus thermophilus DGCC7710 and its virulent phage 2972. At the molecular level, the bacteriophage-immune mutant bacteria (BIMs) and CRISPR–escape mutant phage (CEMs) obtained in this study are consistent with those anticipated from an iterative model of this adaptive immune system: resistance by the addition of novel spacers and phage evasion of resistance by mutation in matching sequences or flanking motifs. While CRISPR BIMs were readily isolated and CEMs generated at high rates (frequencies in excess of 10⁻⁶), our population studies indicate that there is more to the dynamics of phage–host interactions and the establishment of a BIM–CEM arms race than predicted from existing assumptions about phage infection and CRISPR–cas immunity. Among the unanticipated observations are: (i) the invasion of phage into populations of BIMs resistant by the acquisition of one (but not two) spacers, (ii) the survival of sensitive bacteria despite the presence of high densities of phage, and (iii) the maintenance of phage-limited communities due to the failure of even two-spacer BIMs to become established in populations with wild-type bacteria and phage. We attribute (i) to incomplete resistance of single-spacer BIMs. Based on the results of additional modeling and experiments, we postulate that (ii) and (iii) can be attributed to the phage infection-associated production of enzymes or other compounds that induce phenotypic phage resistance in sensitive bacteria and kill resistant BIMs. We present evidence in support of these hypotheses and discuss the implications of these results for the ecology and (co)evolution of bacteria and phage.

The evidence that the CRISPR regions of the genomes of archaea and bacteria play a role in the ecology and (co)evolution of these microbes and their viruses is overwhelming: (i) the spacers (variable sequences of 26–72 bp of DNA between the repeats of this region) of these prokaryotes are homologous to the DNA of viruses in their communities; (ii) experimentally, the acquisition and incorporation of spacers of viral DNA can protect these organisms from subsequent infection by these viruses; (iii) experimentally, viruses evade this immunity by mutation in homologous protospacers or protospacer-adjacent motifs (PAMs). Not so clear are the nature and magnitude of the role CRISPR plays in this ecology and evolution. Here, we use mathematical models, experiments with Streptococcus thermophilus and the phage 2972, and DNA sequence analyses to explore the contribution of CRISPR–cas immunity to the ecology and (co)evolution of bacteria and their viruses. The results of this study suggest that the contribution of CRISPR to the ecology of bacteria and phage is more modest and limited, and the conditions for a CRISPR–mediated coevolutionary arms race between these organisms more restrictive, than anticipated from models based on the canonical view of phage infection and CRISPR–cas immunity.

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.

Understanding bacteriophage specificity in natural microbial communities

Studying the coevolutionary dynamics between bacteria and the bacteriophage viruses that infect them is critical to understanding both microbial diversity and ecosystem functioning. Phages can play a key role in shaping bacterial population dynamics and can significantly alter both intra- and inter-specific competition among bacterial hosts. Predicting how phages might influence community stability and apparent competition, however, requires an understanding of how bacteria-phage interaction networks evolve as a function of host diversity and community dynamics. Here, we first review the progress that has been made in understanding phage specificity, including the use of experimental evolution, we then introduce a new dataset on natural bacteriophages collected from the phyllosphere of horse chestnut trees, and finally we highlight that bacterial sensitivity to phage is rarely a binary trait and that this variation should be taken into account and reported. We emphasize that there is currently insufficient evidence to make broad generalizations about phage host range in natural populations, the limits of phage adaptation to novel hosts, or the implications of phage specificity in shaping microbial communities. However, the combination of experimental and genomic approaches with the study of natural communities will allow new insight to the evolution and impact of phage specificity within complex bacterial communities.

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.

Phage-bacteria infection networks

Phage and their bacterial hosts are the most abundant and genetically diverse group of organisms on the planet. Given their dominance, it is no wonder that many recent studies have found that phage-bacteria interactions strongly influence global biogeochemical cycles, incidence of human diseases, productivity of industrial microbial commodities, and patterns of microbial genome diversity. Unfortunately, given the extreme diversity and complexity of microbial communities, traditional analyses fail to characterize interaction patterns and underlying processes. Here, we review emerging systems approaches that combine empirical data with rigorous theoretical analysis to study phage-bacterial interactions as networks rather than as coupled interactions in isolation.

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.

Multi-scale structure and geographic drivers of cross-infection within marine bacteria and phages

Bacteriophages are the most abundant biological life forms on Earth. However, relatively little is known regarding which bacteriophages infect and exploit which bacteria. A recent meta-analysis showed that empirically measured phage-bacteria infection networks are often significantly nested, on average, and not modular. A perfectly nested network is one in which phages can be ordered from specialist to generalist such that the host range of a given phage is a subset of the host range of the subsequent phage in the ordering. The same meta-analysis hypothesized that modularity, in which groups of phages specialize on distinct groups of hosts, should emerge at larger geographic and/or taxonomic scales. In this paper, we evaluate the largest known phage-bacteria interaction data set, representing the interaction of 215 phage types with 286 host types sampled from geographically separated sites in the Atlantic Ocean. We find that this interaction network is highly modular. In addition, some of the modules identified in this data set are nested or contain submodules, indicating the presence of multi-scale structure, as hypothesized in the earlier meta-analysis. We examine the role of geography in driving these patterns and find evidence that the host range of phages and the phage permissibility of bacteria is driven, in part, by geographic separation. We conclude by discussing approaches to disentangle the roles of ecology and evolution in driving complex patterns of interaction between phages and bacteria.

Ocean viruses and their effects on microbial communities and biogeochemical cycles

Viruses are the most abundant life forms on Earth, with an estimated 10³¹ total viruses globally. The majority of these viruses infect microbes, whether bacteria, archaea or microeukaryotes. Given the importance of microbes in driving global biogeochemical cycles, it would seem, based on numerical abundances alone, that viruses also play an important role in the global cycling of carbon and nutrients. However, the importance of viruses in controlling host populations and ecosystem functions, such as the regeneration, storage and export of carbon and other nutrients, remains unresolved. Here, we report on advances in the study of ecological effects of viruses of microbes. In doing so, we focus on an area of increasing importance: the role that ocean viruses play in shaping microbial population sizes as well as in regenerating carbon and other nutrients.

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.

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.

Intricate interactions between the bloom-forming cyanobacterium Microcystis aeruginosa and foreign genetic elements, revealed by diversified clustered regularly interspaced short palindromic repeat (CRISPR) signatures

Clustered regularly interspaced short palindromic repeats (CRISPR) confer sequence-dependent, adaptive resistance in prokaryotes against viruses and plasmids via incorporation of short sequences, called spacers, derived from foreign genetic elements. CRISPR loci are thus considered to provide records of past infections. To describe the host-parasite (i.e., cyanophages and plasmids) interactions involving the bloom-forming freshwater cyanobacterium Microcystis aeruginosa, we investigated CRISPR in four M. aeruginosa strains and in two previously sequenced genomes. The number of spacers in each locus was larger than the average among prokaryotes. All spacers were strain specific, except for a string of 11 spacers shared in two closely related strains, suggesting diversification of the loci. Using CRISPR repeat-based PCR, 24 CRISPR genotypes were identified in a natural cyanobacterial community. Among 995 unique spacers obtained, only 10 sequences showed similarity to M. aeruginosa phage Ma-LMM01. Of these, six spacers showed only silent or conservative nucleotide mutations compared to Ma-LMM01 sequences, suggesting a strategy by the cyanophage to avert CRISPR immunity dependent on nucleotide identity. These results imply that host-phage interactions can be divided into M. aeruginosa-cyanophage combinations rather than pandemics of population-wide infectious cyanophages. Spacer similarity also showed frequent exposure of M. aeruginosa to small cryptic plasmids that were observed only in a few strains. Thus, the diversification of CRISPR implies that M. aeruginosa has been challenged by diverse communities (almost entirely uncharacterized) of cyanophages and plasmids.