|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: BRL. Performed the experiments: BRL. Analyzed the data: BRL. Wrote the paper: BRL. Did the programming: BRL.
Clustered, Regularly Interspaced Short Palindromic Repeats (CRISPR) abound in the genomes of almost all archaebacteria and nearly half the eubacteria sequenced. Through a genetic interference mechanism, bacteria with CRISPR regions carrying copies of the DNA of previously encountered phage and plasmids abort the replication of phage and plasmids with these sequences. Thus it would seem that protection against infecting phage and plasmids is the selection pressure responsible for establishing and maintaining CRISPR in bacterial populations. But is it? To address this question and provide a framework and hypotheses for the experimental study of the ecology and evolution of CRISPR, I use mathematical models of the population dynamics of CRISPR-encoding bacteria with lytic phage and conjugative plasmids. The results of the numerical (computer simulation) analysis of the properties of these models with parameters in the ranges estimated for Escherichia coli and its phage and conjugative plasmids indicate: (1) In the presence of lytic phage there are broad conditions where bacteria with CRISPR-mediated immunity will have an advantage in competition with non-CRISPR bacteria with otherwise higher Malthusian fitness. (2) These conditions for the existence of CRISPR are narrower when there is envelope resistance to the phage. (3) While there are situations where CRISPR-mediated immunity can provide bacteria an advantage in competition with higher Malthusian fitness bacteria bearing deleterious conjugative plasmids, the conditions for this to obtain are relatively narrow and the intensity of selection favoring CRISPR weak. The parameters of these models can be independently estimated, the assumption behind their construction validated, and the hypotheses generated from the analysis of their properties tested in experimental populations of bacteria with lytic phage and conjugative plasmids. I suggest protocols for estimating these parameters and outline the design of experiments to evaluate the validity of these models and test these hypotheses.
CRISPR is the acronym for the adaptive immune system that has been found in almost all archaebacteria and nearly half the eubacteria examined. Unlike the other defenses bacteria have for protection from phage and other deleterious DNAs, CRISPR has the virtues of specificity, memory, and the capacity to abort infections with a virtually indefinite diversity of deleterious DNAs. In this report, mathematical models of the population dynamics of bacteria, phage, and plasmids are used to determine the conditions under which CRISPR can become established and will be maintained in bacterial populations and the contribution of this adaptive immune system to the ecology and (co)evolution of bacteria and bacteriophage. The models predict realistic and broad conditions under which bacteria bearing CRISPR regions can invade and be maintained in populations of higher fitness bacteria confronted with bacteriophage and narrower conditions when the confrontation is with competitors carrying conjugative plasmids. The models predict that CRISPR can facilitate long-term co-evolutionary arms races between phage and bacteria and between phage- rather than resource-limited bacterial communities. The parameters of these models can be independently estimated, the assumptions behind their construction validated, and the hypotheses generated from the analysis of their properties tested with experimental populations of bacteria.
For many species of bacteria, adaptive evolution is through the expression of chromosomal and extrachromosomal (plasmid- and prophage - borne) genes or clusters of genes (pathogenicity and nicer islands) acquired by horizontal gene transfer (HGT) from the same or even quite distant species , . Thus, on first consideration it may seem that bacteria and their accessory genetic elements would have mechanism to promote the acquisition, incorporation and expression of genes from without. And, indeed there are mechanisms like integrons – that appear to have that function. On the other side, DNA acquired from external sources may be deleterious. This is certainly the case when that DNA is borne on lytic bacteriophage, but also for plasmids that engender fitness costs ,  or chromosomal DNA from the wrong source , . To deal with these contingencies, it would seem that bacteria would have mechanisms to protect themselves against infection by deleterious foreign DNA . And indeed there are systems like restriction-modification (restriction endonucleases) which appear to have that role , .
The most recently discovered mechanism postulated to provide bacteria immunity to infectious genetic elements are Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR). For recent reviews see , . CRISPR is particularly intriguing because of its ubiquity, appearing in ~90% and ~40% of archaeal and eubacterial sequenced genomes, respectively, and because of the adaptive mechanism by which it provides immunity to infections by a virtually indefinite diversity of bacteriophage and plasmids. DNA from infecting phage and plasmids is incorporated into the CRISPR array. Through a yet to be fully elucidated mechanism, bacteria abort the replication of infecting phage  or the establishment of conjugative plasmids  bearing copies of the DNA incorporated into their CRISPR arrays, also see . Further support for CRISPR being an adaptive immune system that is maintained because it protects bacteria from infection with phage comes from studies of the community ecology of bacteria and phage; DNA in the CRISPR regions of the bacteria from those communities corresponds to that in the co-existing phage –. For an intriguing perspective on CRISPR as a witness to the coevolutionary history of bacteria and phage, see .
CRISPR-mediated immunity has been likened to a Lamarckian mechanism , because the selection pressure, the infecting phage and plasmids, determine the genotype. This analogy however does not account for the evolution and maintenance of the machinery responsible for taking up the infecting phage and plasmid DNA and the mechanism employed to prevent the replication or establishment of infecting genetic elements with those sequences. Under what conditions will adaptive immunity to phage and plasmid infection be the selection pressure responsible for establishing and maintaining CRISPR-mediated immunity in populations of archeae and bacteria? What about other mechanisms of resistance, like structural modification blocking phage adsorption (envelope resistance) and restriction-modification? How do these mechanisms interact with CRISPR – acquired immunity and contribute to its establishment and maintenance?
To address these questions and provide a framework and hypotheses for their study experimentally, I use mathematical models of the population dynamics of bacteria, phage and plasmids to explore the conditions under which a CRISPR–like adaptive immune mechanism will provide bacteria a selective advantage in competition with bacteria without this immune system. The results of the numerical analysis of the properties of these models suggest that with bacterial replication and phage infection parameters in realistic ranges, there are broad but not universal conditions where a CRISPR–like adaptive immune system can be favored and will be maintained in populations of bacteria confronted with lytic phage. While this model predicts conditions where CRISPR-mediated immunity will be favored when bacteria compete with populations bearing conjugative plasmids, these conditions are relatively restrictive. The parameters of these models can be independently estimated, the validity of the assumptions behind their construction and the hypotheses generated from the analysis of the properties can be tested in experimental populations of bacteria with lytic phage and conjugative plasmids. Procedures for doing these experiments are outlined and their potential outcomes described and/or speculated upon. Also discussed are the broader implications of CRISR-mediated adaptive immunity to the population and evolutionary biology and ecology of bacteria and phage.
Both the lytic phage and conjugative plasmid models used here assume a chemostat-like habitat. The bacteria grow at a rate that is a monotonically increasing function of the concentration of a limiting resource, R µg/ml .
where Vi hr−1 is the maximum growth rate of the ith strain of bacteria and k the concentration of the resource when the growth rate is half its maximum value (the “Monod constant”). The populations are maintained in a vessel of unit volume, (1ml) into which medium containing the limiting resource from a reservoir where it is maintained at a concentration A µg/ml flows in at a rate w per hour. Excess resource and wastes are removed from the vessel at the same rate. As in , the rate of uptake of the resource by the bacteria is proportional to the density, the resource concentration-dependent growth rates of the different populations of bacteria and a conversion efficiency parameter, e µg/per cell.
The model developed here is an extension of that in . There are four populations of bacteria. Two are sensitive to the phage, N, non–CRISPR and C, CRISPR and two that are either fully resistant (envelope resistance), or immune because of CRISPR, NR and CR, respectively. The variables N, C, NR and CR are the both the densities (bacteria per ml) of these populations and used as their designations. There is one population of phage, with density and designation, P particles per ml.
The phage adsorb to the N and C and CR bacteria with rate constants, δN and δC (ml per phage per cell per hour) respectively. Phage do not adsorb to bacteria with envelope resistant, i.e. the NR cells. To account for a possible multiplicity of infection (MOI) effect on survival of phage-infected CR, the effective killing rate constant for phage adsorption to CRISPR can be an increasing function of the ratio of free phage and CR cells, M=P/CR.
where δMIN and δMAX are the minimum and maximum adsorption rates. The parameter x is a coefficient (0≤x≤1) that specifies the magnitude of the MOI effect, q is the MOI where the adsorption rate is half its maximum value and n is an exponent which contributes to the shape of the distribution. At low multiplicities, δCR (M) the CRISPR cells would be effectively immune (resistant) (Figure 1). At high multiplicities, however, immune CRISPR cells can be overburdened by phage, their immunity would be overridden, and the phage would replicate, killing the cells. On the other side, we assume that the phage are removed from the population by adsorption to immune CRISPR cells at the maximum adsorption rate, δMAX.
For convenience I neglect the latent periods of the phage infection but assume that the phage have potentially different burst sizes, βN, β C, and βCR particles per cell, for N, C and CR cells, respectively.
Phage-immune CRISPR cells, CR are produced from C at a rate proportional to the rate at which the phage adsorb to them and a constant m (0≤m≤1) which is the probability that a phage infection will be aborted and a CRISPR strain will be produced. At a rate v per cell per hour, CRISPR lose their immunity, CR→C. For the N and C populations the loss of the adsorbed phage is subsumed in the value of the burst size (which is one less than the number of phage produced). For the CR population, the loss of the phage due to adsorption is specifically considered because only a small fraction of the adsorbed phage replicate when the MOI is low.
The model developed here is an extension of that in . There are five bacterial populations. Two populations do not code for CRISPR, N and NP, and three populations code for CRISPR, C and CP and CX. The NP and CP populations bear the conjugative plasmid and CX, carries CRISPR and plasmid sequences that make it completely immune to the receipt of these plasmids. Plasmids are transferred by conjugation at rates proportional to the product of the densities of the plasmid-bearing and plasmid-free populations and rate constants, γNN, γNC, γCN and γCC (ml per cell per hour) respectively for the transfer of the plasmid from NP to N, NP to C, CP to N and CP to C., respectively. Plasmids are lost by vegetative segregation at rates τN and τC per cell per hour, with NP→N and Cp→C. C are converted to CX at a rate proportional to the rate at which C acquires the plasmid and a probability m (0≤m≤1). Cx lose the CRISPR plasmid immunity region and become C at rate ν per cell per hour. Each of the cell lines, have a maximum growth rate, VN, VNP, VC, and VCP, and VX per hour. In Figure 3, I illustrate the interactions between the different cell lines in this model, and, in Table 2, I separately define the parameters and variables. The equations for this model are:
For the numerical solutions to these equations (computer simulations) I use a differential equation-solving software package, Berkeley Madonna. For the phage simulations there is a refuge density, below which the phage are unable to adsorb to the bacteria. The purpose of this is to control the system from oscillating without limits, see . In these simulations, if the phage density falls below 10−1 particles per ml, the phage are considered to be lost. Copies of these simulations are available online, www.eclf.net/programs.
The bacterial growth, resource-uptake, phage adsorption parameters and burst sizes used in these simulations (Table 1) are in a range similar to that which we observed for E. coli and the phages T2 and T7 , .
In a chemostat with susceptible bacteria at an equilibrium density N*, a lytic phage can become established and will maintain a population with sensitive bacteria as long as the rate of phage production exceeds the rate of washout, δNβNN*>w . With the parameters used in these simulations, N*~108 (see ). As long as δNβNN*>2×10−9, the phage will become established and can maintain a population by replicating on sensitive bacteria (Figure 4A). The oscillations in the densities of bacteria and phage in these and the following simulations are those anticipated for the predator-prey nature of these dynamics.
To explore the conditions under which a CRISPR population will become established and be maintained in the presence of phage, I consider situations where the C and CR populations have an intrinsic selective disadvantage relative to N (VN>VC, VCR) and therefore cannot invade an established N population in the absence of these bacterial viruses. Because of the immunity of CR, with phage present and in the absence of a multiplicity effect, an initially rare CRISPR population will invade and ascend to dominance despite its lower intrinsic fitness (Figure 4B). With these parameters, the phage are maintained along with N and C, the latter being continually generated by the loss of immunity by the dominant CR population. The N population is maintained because of its higher intrinsic fitness (growth rate) relative to CR, and resources, rather than phage predation, limit the bacteria at large. The phage continue to be maintained by replicating on the N and C cells. Although the oscillations are damped and in time would no longer be noticed, that time would be considerably greater than would be feasible to study experimentally with chemostats. If we allow for a strong multiplicity effect (x=0.5), the CRISPR population becomes established, and both immune and non-immune CRISPR cells maintain their populations with sensitive non–CRISPR in a phage- rather than resource- limited community (Figure 4C). When the magnitude of the multiplicity effect is reduced (x=0.2), the phage continue to be maintained but immune CRISPR cells ascend to dominance and the community with three populations of bacteria, N, C and CR are maintained in a resource- rather than a phage-limited state (Figure 4D).
In addition to CRISPR immunity, when confronted with phage, bacteria may generate mutants to which phage are unable to adsorb or are resistant by other mechanisms . To explore how this envelope resistance will affect the conditions for the establishment and maintenance of CRISPR, we consider the invasion of an envelope resistant strain of N, NR, into a population of N and phage. In these simulations, the C and CR are less fit than N (VN>VC, VCR) and NR are less fit than C and CR, (VNR<VC, VCR). Were the NR cells more fit than C and CR, they would dominate and the CRISPR population would not invade an would not be established. Whether this fitness relationship will be seen with real bacteria and what those fitness will be is an empirical question.
As can be seen in Figure 5A, although the resistant, NR strain is the least intrinsically fit bacteria in the community (lowest maximum growth rate), in the presence of phage it ascends rapidly and achieves dominance. During this initial phase, as a consequence of the production of immune CR cells, the CRISPR population also increases in density, but remains a minority population relative to the resistant non–CRISPR NR. With these parameters, the phage density declines after the ascent of resistance and the densities of both the N and C populations increase. Shortly after the phage are eliminated the highest fitness N population ascends and lower fitness C, CR and NR decline. If the phage resistant population is substantially less fit than the other bacterial populations, the CR population ascends to dominance and continues to co-exist with the phage, N, and C populations (Figure 5B).
In accord with , conjugative plasmids will be maintained as long as the rate of infectious transfer exceeds the rates of loss of the plasmid due to selection against the cells carrying it, vegetative segregation, and the rate of flow through the chemostat. In terms of the above parameters, the plasmid will be maintained in an N-NP population as long as
where N* is the density of plasmid-free cells at the chemostat equilibrium. For example, if VN=1.0, VNP=0.95, w=0.2, τN=10−3, the plasmid will be maintained in a population of density N*=108 as long as γNN>1.1×10−10. If the plasmid augments the growth rate (which in this model is the sole parameter of cell fitness) of the bacteria that carry it, VNP>VN, as we would anticipate for antibiotic resistance encoding plasmids in the presence of the selecting antibiotic, bacteria bearing the plasmid will be able to invade even without transfer, as long as the segregation rate, τN, is sufficiently small.
The population dynamics of selection and plasmid transfer in an equilibrium chemostat in the absence of CRISPR are presented in Figure 6A. If the conditions specified in equation (2) are met, the plasmid- bearing cells become established and ascend to dominate the N-NP community, whether cells bearing the plasmid are favored or not. If the plasmid is maintained by transfer or selection for the genes it carries and τN>0, there will be a stable population of plasmid-free cells. When the rate constant of plasmid transfer is too low, the deleterious plasmid will be lost.
To consider the effects of CRISPR on the population dynamics of bacteria with conjugative plasmids and the conditions under which CRISPR immunity will provide an advantage to bacteria, I let the maximum growth rates of the CRISPR strains (the sole measure of intrinsic, phage-independent fitness) be somewhat lower than the corresponding non–CRISPR cells. In Figure 6B, the population is initially at equilibrium with a plasmid-free, non-immune CRISPR population and a low density of plasmid-bearing non–CRISPR bacteria are introduced. The plasmid spreads rapidly from NP to C producing a CP population which in turn generates immune CRISPR, CX. While the C and C P populations die out, CX ascends to dominance and minority populations of N and NP are maintained. Although the CX population has a lower growth rate than N, in the presence of a deleterious conjugative plasmid they have an advantage because they cannot be infected by that element. They do not eliminate the N and NP populations due to the loss of the CRISPR region and the conversion into N. As can be seen in Figure 6C, with these parameters and a lower growth rate, CX can invade an equilibrium N-NP population, but the rate of increase in the density of Cx is low. The invasion rate for CX would even be further reduced if, instead of CX, a plasmid-free C invaded an NP population, because it would be some time before the CX is produced and, in a finite population, may not be produced at all (“data” not shown). A very different situation obtains when the plasmid confers a growth rate advantage to the infected host (Figure 6D). Under these conditions, the C populations and its derivatives, CP and CX, are eliminated.
“All models are wrong, some are useful.” (George Box)
It has been less than eight years since the ubiquitous clusters of palindromic repeats now known as CRISPR first acquired this moniker . Although there had been compelling circumstantial evidence that CRISPR was part of an adaptive immune system that provides protection against infecting phage and plasmids, it has been less than four and three years respectively since the publication of the first direct (read experimental) evidence that CRISPR can provide immunity to infection by lytic phage  and conjugative plasmids .
In the course of this time a great deal has been learned about the molecular biology of CRISPR and the mechanisms by which it provides adaptive immunity to plasmid and phage infection. But there remain many unanswered questions about these processes. Most important for this consideration is a dearth of the quantitative information needed to understand the population dynamics of CRISPR-mediated adaptive immunity and thereby the conditions for the establishment and maintenance of CRISPR in bacterial populations. To my knowledge, this study is the first formal consideration of these dynamics.
The models developed in this report incorporate what has been learned about CRISPR-mediated adaptive immunity to phage and conjugative plasmids, primarily from the studies of Barrangou and colleagues  and Marraffini and Sontheimer , into models of the population dynamics of lytic phage  and conjugative plasmids . Although they may appear complex, at best they are simplistic caricatures of interactions between these infectious genetic elements and bacteria with CRISPR-mediated adaptive immunity. These models are not intended or anticipated to be numerically precise analogs of these processes and dynamics.
The role of these mathematical models is similar to that of the diagrammatic models (cartoons) used to illustrate the molecular basis and mode of action of CRISPR, i.e., to provide a framework for understanding these processes, designing experiments, and interpreting their results. In this case, these experiments are on population and evolutionary dynamics of bacteria with CRISPR-mediated immunity confronted with lytic phage and competing bacteria bearing conjugative plasmids. The purpose of these models for this experimental enterprise is: (i) to identify and, in a quantitative way, evaluate the role of the different factors (parameters) contributing to these dynamics and the conditions for the establishment and maintenance of CRISPR in bacterial populations, and (ii) to generate hypotheses about these dynamics and existence conditions that can be tested (and rejected) in experimental populations.
The results of the analysis of the properties of the phage - CRISPR model are consistent with the proposition that in the presence of lytic bateriophage there are broad conditions under which a CRISPR–like adaptive immune system can become established and will be maintained in bacterial populations. With population densities, growth rates, and phage infection parameters in realistic ranges, these models predict that despite a growth rate disadvantage, bacteria with CRISPR–like acquired immunity to infecting phage will increase in frequency when initially rare and will be maintained. The necessary condition for this is that the phage population continues to persist at a sufficiently high density for CRISPR-mediated adaptive immunity to overcome an intrinsic disadvantage associated with the costs of carrying and expressing these genes.
When will the phage maintain their populations at sufficient levels for this outcome? With the parameters used to address this question, the phage will be maintained under broad conditions, but may eventually be lost if a population with envelope or other resistance ascends to dominance. I emphasized the word may for two reasons. The first is theoretical, if the relative growth rate of the resistant population is adequately low, the phage and thereby CRISPR will be maintained. The second is empirical, even when resistant bacteria dominate experimental populations of bacteria and phage, in general the phage continue to be maintained , , .
The CRISPR plasmid model predicts that because of the immunity to infection with conjugative plasmids, a lower growth rate (Malthusian fitness) CRISPR population can become established and will be maintained when competing with bacteria with a greater Malthusian fitness but bearing deleterious (fitness-reducing) conjugative plasmids. Although these conditions are met with plasmid fitness costs in the range estimated for “laboratory” plasmids , , it is not clear that naturally occurring plasmids would be as burdensome as those maintained in the Lab. The greater the Malthusian fitness burden attributed to the plasmid, the greater the advantage of CRISPR-mediated immunity.
The rate constants of plasmid transfer used in these simulations are those for plasmids with permanently derepressed conjugative pili synthesis. Wild type conjugative plasmids are more likely to be repressed for the production of these transfer organelles and would have substantially lower rates of transmission than plasmids that are permanently derepressed for plasmid transfer , . Indeed, it is not clear whether in natural populations conjugative plasmids that engender fitness cost can be maintained by transfer alone. Their persistence may require periodic episodes where bacteria carrying them have an advantage , , but also see . If the rate of infectious transfer is not sufficient to maintain deleterious plasmid in a population and they persist by continually or periodically enhancing the cells Malthusian fitness, immunity to these plasmids would not be sufficient to maintain CRISPR-encoding cells that have an intrinsic fitness disadvantage.
It would be nearly impossible to determine whether the quantitative conditions predicted by these models for the establishment and maintenance of CRISPR-mediated immunity are met in natural populations. On the other hand, the values of the parameters of these models can be estimated and the validity of the assumptions behind their construction and hypotheses generated from the analysis of their properties can be tested in laboratory culture using CRISPR–positive and CRISPR–negative bacterial constructs, phage and plasmids of the types used respectively by Barrangou and colleagues  and Marraffini and Sontheimer,  in chemostat culture.
All of the parameters of these models (Table 1 and Table 2) can be independently estimated and procedures for doing so have been published for the majority of them: (1) for the bacterial growth and resource utilization parameters, the VS, k, and e, see , ; (2) for the phage latent periods, adsorption rates δs, and burst sizes, the βs, see , (3) for the rate constants of plasmid transfer, the γs, see , , and (4) for the mutation rate to envelope resistance, see , . Estimates of the plasmid segregation rate, τ, can be obtained by plating low-density cultures of plasmid-bearing cells, and testing colonies for the plasmid marker. However, unless τ is very high (τ>0.005 per cell per division), this procedure would be excessively labor intensive. However, if low, this parameter would have a negligible contribution to the dynamics of the plasmid and estimating its value would not be worthwhile.
Protocols for isolating bacteria with CRISPR-mediated resistance to phage and plasmid infection, can be found in  and , respectively. I am, however, unaware of published studies providing estimates of the fractions of phage and plasmid infected cells that become immune, the parameter m, or the rates of loss of these immunities, ν, in the models (Figure 2 and Figure 3). In Text S1, I outline potential ways to estimate these parameters. I emphasize the word potential because without actually doing these experiments, it is difficult to anticipate pitfalls and problems with the proposed procedures.
In developing the model, I made a series of assumptions about CRISPR – mediated immunity and the population dynamics of bacteria with lytic phage and conjugative plasmids. In the following, I list these assumptions and briefly describe what would be anticipated experimentally if these assumptions are correct.
One way to evaluate how well these models serve as analogs of the population dynamics of bacteria with CRISPR adaptive immunity to bacteria and phage is to compare the results of simulations with independently estimated parameters to that observed in chemostat populations. Although it would be gratifying to see quantitative agreement between the anticipated dynamics and those observed in experimental populations, populations with CRISPR constructs of bacteria, conjugative plasmids and phage, it would also be surprising. These models are far too simple to expect the predicted and observed dynamics to be numerically coincident. A more modest, realistic, and, I believe, more useful goal is test predictions made from the analysis of the properties of these models in a qualitative – semi-quantitative way and identify those elements of the model that have to be modified to make the models more realistic and accurate. In the following, I list these predictions.
(i) When mixtures of otherwise isogenic CRISPR positive and negative phage –sensitive constructs are introduced into chemostats in approximately equal frequencies:
(ii) When introduced at low frequencies into chemostats with sensitive non–CRISPR cells in the presence of phage, as long as immune CRISPR cells are produced, the CRISPR population will increase in frequency. This will not be the case in the absence of phage.
When mixtures of non–CRISPR cells bearing fitness reducing conjugative plasmids and plasmid-free CRISPR cells are introduced into chemostats:
In this report, I elected to restrict the model and its analysis to the simplest cases with lowest realistic number of states of bacteria, phage and plasmids. I have done so because at this time these minimum number of states models and the predictions generated from their analysis are more amenable to evaluating and testing experimentally than models with more states of bacteria, phage and plasmids. Moreover, these tests, and particularly the population dynamic experiments, should indicate the importance of the generation of additional population states by mutation, like host range phage and host range plasmids, are to these dynamics. Be that as it may, I also realize that this minimum number of states model will not account for what may turn out to be the most important contributions of CRISPR-mediated immunity to the ecology as well as the population and evolutionary biology of bacteria and phage.
Luciano Marraffini (personal communication) suggested one potentially important contribution of CRISPR to the population and evolutionary dynamics of bacteria and phage. Unlike envelope resistance, which is almost always restricted to phage that utilize single adsorption organelles, , CRISPR–immunity can be effective against multiple phages with different adsorption organelles (independent resistance). Moreover, envelope resistance is likely to engender a cost in Malthusian fitness, e.g. see , ,  and that cost will almost certainly be greater if this resistance is for multiple phages that employ different receptors for infection.
If these interpretations are correct, it would seem experimental populations with CRISPR-encoding bacteria with envelope resistance to all the phage will not evolve and CRISPR will prevail in competition with sensitive non–CRISPR cells. If, however, the results of a test of this multi-phage hypothesis Ryzard Koroana and did in a study of the conditions for the maintenance of restriction endonuclease (restriction-modification, R-M) immunity are general , this hypothesis may be rejected. E. coli bearing an R-M system conferring immunity to three phage with different organelles were challenged with a mixture of all three of these phages. As a consequence of a hierarchy of phage replication , there was sequential selection for the different resistant states and within a day of exposure, bacteria with envelope resistance to all three phages dominated the community .
A number of years ago, Richard Lenski and I postulated that the arms race between bacterial resistance and host range phage would be limited to few cycles and is likely to end with resistant bacteria to which phage would not be able to generate host range mutations . The empirical basis of our hypothesis was the results of experiments with E. coli and its phage and envelope resistance, , , , . While this interpretation was also supported by experiments with V. cholerae and its phage JSF4 , experiments with Pseudomonas fluorescens and its phage SBW25  suggest extended arms races are possible. Although, to my knowledge, the mechanisms responsible for the continuous changes in resistance and host-range reported in this study with this strain of Pseudomonas and phage have yet to be elucidated, CRISPR does provide a mechanism for long-term arms races between bacteria and phage , , . By single base changes in sequences of DNA into the spacer regions of CRISPR, a phage can infect and replicate on previously immune CRISPR cells. By incorporating the mutated or other region of that phage into another spacer, CRISPR cells can generate resistance to these host range phages. At this time, it is not at all clear how long or through how many cycles a CRISPR-mediated arms race can proceed. I would it certainly be interesting, tenable experimentally and fun to find out. Be it by CRISPR or by sequential resistance and host-range mutation ,  an extended arms race could provide a way for phage, rather than resources, to limit the densities of bacterial populations (see Text S2), which is an ecological outcome with practical as well as theoretical implications, e.g. see –.
Protocols to estimate the probability of formation, m, and rate of loss, ν, of CRISPR-mediated immunity to phage and conjugative plasmids.
(0.03 MB DOC)
Arms races and phage-limited bacterial populations.
(0.33 MB DOC)
I would like to thank Luciano Marraffini and Joshua Weitz for stimulating discussions and correspondence and useful comments and suggestions on an earlier draft of this report. I am grateful to Amy Kirby, Amoolya Singh, Klas Udekwu, and Jim Bull for helpful suggestions, comments, and corrections. I apologize to the members of the EcLF and others sipping coffee and eating lunch in the “Peoples room” and to the nice Grandmother sitting next to me on an airplane for my ranting on about the wonders of CRISPR.
The authors have declared that no competing interests exist.
This work was supported by US National Institutes of Health GM091875 (BRL). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.