Why put up with immunity when there is resistance: an excursion into the population and evolutionary dynamics of restriction-modification and CRISPR-Cas

Multi-layered genome defences in bacteria

Bacteria have evolved a variety of defence mechanisms to protect against mobile genetic elements, including restriction-modification systems and CRISPR-Cas. In recent years, dozens of previously unknown defence systems (DSs) have been discovered. Notably, diverse DSs often coexist within the same genome, and some co-occur at frequencies significantly higher than would be expected by chance, implying potential synergistic interactions. Recent studies have provided evidence of defence mechanisms that enhance or complement one another. Here, we review the interactions between DSs at the mechanistic, regulatory, ecological and evolutionary levels.

Conditions for the spread of CRISPR-Cas immune systems into bacterial populations

Bacteria contain a wide variety of innate and adaptive immune systems which provide protection to the host against invading genetic material, including bacteriophages (phages). It is becoming increasingly clear that bacterial immune systems are frequently lost and gained through horizontal gene transfer. However, how and when new immune systems can become established in a bacterial population have remained largely unstudied. We developed a joint epidemiological and evolutionary model that predicts the conditions necessary for the spread of a CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated) immune system into a bacterial population lacking this system. We found that whether bacteria carrying CRISPR-Cas will spread (increase in frequency) into a bacterial population depends on the abundance of phages and the difference in the frequency of phage resistance mechanisms between bacteria carrying a CRISPR-Cas immune system and those not (denoted as ${f}_{\Delta }$). Specifically, the abundance of cells carrying CRISPR-Cas will increase if there is a higher proportion of phage resistance (either via CRISPR-Cas immunity or surface modification) in the CRISPR-Cas-possessing population than in the cells lacking CRISPR-Cas. We experimentally validated these predictions in a model using Pseudomonas aeruginosa PA14 and phage DMS3vir. Specifically, by varying the initial ratios of different strains of bacteria that carry alternative forms of phage resistance, we confirmed that the spread of cells carrying CRISPR-Cas through a population can be predicted based on phage density and the relative frequency of resistance phenotypes. Understanding which conditions promote the spread of CRISPR-Cas systems helps to predict when and where these defences can become established in bacterial populations after a horizontal gene transfer event, both in ecological and clinical contexts.

Microbial evolution through horizontal gene transfer by mobile genetic elements

Mobile genetic elements (MGEs) are crucial for horizontal gene transfer (HGT) in bacteria and facilitate their rapid evolution and adaptation. MGEs include plasmids, integrative and conjugative elements, transposons, insertion sequences and bacteriophages. Notably, the spread of antimicrobial resistance genes (ARGs), which poses a serious threat to public health, is primarily attributable to HGT through MGEs. This mini-review aims to provide an overview of the mechanisms by which MGEs mediate HGT in microbes. Specifically, the behaviour of conjugative plasmids in different environments and conditions was discussed, and recent methodologies for tracing the dynamics of MGEs were summarised. A comprehensive understanding of the mechanisms underlying HGT and the role of MGEs in bacterial evolution and adaptation is important to develop strategies to combat the spread of ARGs.

Development and characterization of a bacteriophage cocktail with high lytic efficacy against field-isolated Salmonella enterica

Salmonella spp. is a prevalent pathogen that causes great public health concern worldwide. Bacteriophage-based cocktails have arisen as an alternative to antibiotics to inhibit the growth of Salmonella. However, the bactericidal effect of bacteriophage cocktails in vivo largely differs from their observed effect in vitro. This is partly because in vitro developments of cocktails do not always consider the bacterial diversity nor the environmental conditions where bacteriophages will have to replicate. Here, we isolated and sequenced 47 bacteriophages that showed variable degrees of lytic activity against 258 Salmonella isolates from a commercial broiler company in Brazil. Three of these bacteriophages were characterized and selected to assemble a cocktail. In vitro quantitative assays determined the cocktail to be highly effective against multiple serovars of Salmonella, including Minnesota and Heidelberg. Remarkably, the in vitro lytic activity of the cocktail was retained or improved in conditions that more closely resembled the chicken gut, such as anaerobiosis, 42°C, and Salmonella mono-strain biofilms. Analysis of bacterial cross-resistance between the 3 bacteriophages composing the cocktail revealed limited or no generation of cross-resistance. Our results highlight the relevance of an optimized flux of work to develop bacteriophage cocktails against Salmonella with high lytic efficacy and strong potential to be applied in vivo in commercial broiler farms.

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.

Engineering bacteriophages for enhanced host range and efficacy: insights from bacteriophage-bacteria interactions

Bacteriophages, the most abundant organisms on earth, have the potential to address the rise of multidrug-resistant bacteria resulting from the overuse of antibiotics. However, their high specificity and limited host range can hinder their effectiveness. Phage engineering, through the use of gene editing techniques, offers a means to enhance the host range of bacteria, improve phage efficacy, and facilitate efficient cell-free production of phage drugs. To engineer phages effectively, it is necessary to understand the interaction between phages and host bacteria. Understanding the interaction between the receptor recognition protein of bacteriophages and host receptors can serve as a valuable guide for modifying or replacing these proteins, thereby altering the receptor range of the bacteriophage. Research and development focused on the CRISPR-Cas bacterial immune system against bacteriophage nucleic acids can provide the necessary tools to promote recombination and counter-selection in engineered bacteriophage programs. Additionally, studying the transcription and assembly functions of bacteriophages in host bacteria can facilitate the engineered assembly of bacteriophage genomes in non-host environments. This review highlights a comprehensive summary of phage engineering methods, including in-host and out-of-host engineering, and the use of high-throughput methods to understand their role. The main aim of these techniques is to harness the intricate interactions between bacteriophages and hosts to inform and guide the engineering of bacteriophages, particularly in the context of studying and manipulating the host range of bacteriophages. By employing advanced high-throughput methods to identify specific bacteriophage receptor recognition genes, and subsequently introducing modifications or performing gene swapping through in-host recombination or out-of-host synthesis, it becomes possible to strategically alter the host range of bacteriophages. This capability holds immense significance for leveraging bacteriophages as a promising therapeutic approach against antibiotic-resistant bacteria.

The coordination of anti-phage immunity mechanisms in bacterial cells

Bacterial cells are equipped with a variety of immune strategies to fight bacteriophage infections. Such strategies include unspecific mechanisms directed against any phage infecting the cell, ranging from the identification and cleavage of the viral DNA by restriction nucleases (restriction-modification systems) to the suicidal death of infected host cells (abortive infection, Abi). In addition, CRISPR-Cas systems generate an immune memory that targets specific phages in case of reinfection. However, the timing and coordination of different antiviral systems in bacterial cells are poorly understood. Here, we use simple mathematical models of immune responses in individual bacterial cells to propose that the intracellular dynamics of phage infections are key to addressing these questions. Our models suggest that the rates of viral DNA replication and cleavage inside host cells define functional categories of phages that differ in their susceptibility to bacterial anti-phage mechanisms, which could give raise to alternative phage strategies to escape bacterial immunity. From this viewpoint, the combined action of diverse bacterial defenses would be necessary to reduce the chances of phage immune evasion. The decision of individual infected cells to undergo suicidal cell death or to incorporate new phage sequences into their immune memory would be determined by dynamic interactions between the host's immune mechanisms and the phage DNA. Our work highlights the importance of within-cell dynamics to understand bacterial immunity, and formulates hypotheses that may inspire future research in this area.

Fitness costs of CRISPR-Cas systems in bacteria

CRISPR-Cas systems provide bacteria with both specificity and adaptability in defence against invading genetic elements. From a theoretical perspective, CRISPR-Cas systems confer many benefits. However, they are observed at an unexpectedly low prevalence across the bacterial domain. While these defence systems can be gained horizontally, fitness costs may lead to selection against their carriage. Understanding the source of CRISPR-related fitness costs will help us to understand the evolutionary dynamics of CRISPR-Cas systems and their role in shaping bacterial genome evolution. Here, we review our current understanding of the potential fitness costs associated with CRISPR-Cas systems. In addition to potentially restricting the acquisition of genetic material that could confer fitness benefits, we explore five alternative biological factors that from a theoretical perspective may influence the fitness costs associated with CRISPR-Cas system carriage: (1) the repertoire of defence mechanisms a bacterium has available to it, (2) the potential for a metabolic burden, (3) larger-scale population and environmental factors, (4) the phenomenon of self-targeting spacers, and (5) alternative non-defence roles for CRISPR-Cas.

High viral abundance and low diversity are associated with increased CRISPR-Cas prevalence across microbial ecosystems

CRISPR-Cas are adaptive immune systems that protect their hosts against viruses and other parasitic mobile genetic elements.¹ Although widely distributed among prokaryotic taxa, CRISPR-Cas systems are not ubiquitous.^2-4 Like most defense-system genes, CRISPR-Cas are frequently lost and gained, suggesting advantages are specific to particular environmental conditions.⁵ Selection from viruses is assumed to drive the acquisition and maintenance of these immune systems in nature, and both theory^6-8 and experiments have identified phage density and diversity as key fitness determinants.^9,10 However, these approaches lack the biological complexity inherent in nature. Here, we exploit metagenomic data from 324 samples across diverse ecosystems to analyze CRISPR abundance in natural environments. For each metagenome, we quantified viral abundance and diversity to test whether these contribute to CRISPR-Cas abundance across ecosystems. We find a strong positive association between CRISPR-Cas abundance and viral abundance. In addition, when controlling for differences in viral abundance, CRISPR-Cas systems are more abundant when viral diversity is low, suggesting that such adaptive immune systems may offer limited protection when required to target a diverse viral community. CRISPR-Cas abundance also differed among environments, with environmental classification explaining roughly a quarter of the variation in CRISPR-Cas relative abundance. The relationships between CRISPR-Cas abundance, viral abundance, and viral diversity are broadly consistent across environments, providing robust evidence from natural ecosystems that supports predictions of when CRISPR is beneficial. These results indicate that viral abundance and diversity are major ecological factors that drive the selection and maintenance of CRISPR-Cas in microbial ecosystems.

Ecological memory preserves phage resistance mechanisms in bacteria

Bacterial defenses against phage, which include CRISPR-mediated immunity and other mechanisms, can carry substantial growth rate costs and can be rapidly lost when pathogens are eliminated. How bacteria preserve their molecular defenses despite their costs, in the face of variable pathogen levels and inter-strain competition, remains a major unsolved problem in evolutionary biology. Here, we present a multilevel model that incorporates biophysics of molecular binding, host-pathogen population dynamics, and ecological dynamics across a large number of independent territories. Using techniques of game theory and non-linear dynamical systems, we show that by maintaining a non-zero failure rate of defenses, hosts sustain sufficient levels of pathogen within an ecology to select against loss of the defense. This resistance switching strategy is evolutionarily stable, and provides a powerful evolutionary mechanism that maintains host-pathogen interactions, selects against cheater strains that avoid the costs of immunity, and enables co-evolutionary dynamics in a wide range of systems.

High viral abundance and low diversity are associated with increased CRISPR-Cas prevalence across microbial ecosystems

CRISPR-Cas are adaptive immune systems that protect their hosts against viruses and other parasitic mobile genetic elements.¹ Although widely distributed among prokaryotic taxa, CRISPR-Cas systems are not ubiquitous.2, 3, 4 Like most defense-system genes, CRISPR-Cas are frequently lost and gained, suggesting advantages are specific to particular environmental conditions.⁵ Selection from viruses is assumed to drive the acquisition and maintenance of these immune systems in nature, and both theory6, 7, 8 and experiments have identified phage density and diversity as key fitness determinants.^9,10 However, these approaches lack the biological complexity inherent in nature. Here, we exploit metagenomic data from 324 samples across diverse ecosystems to analyze CRISPR abundance in natural environments. For each metagenome, we quantified viral abundance and diversity to test whether these contribute to CRISPR-Cas abundance across ecosystems. We find a strong positive association between CRISPR-Cas abundance and viral abundance. In addition, when controlling for differences in viral abundance, CRISPR-Cas systems are more abundant when viral diversity is low, suggesting that such adaptive immune systems may offer limited protection when required to target a diverse viral community. CRISPR-Cas abundance also differed among environments, with environmental classification explaining roughly a quarter of the variation in CRISPR-Cas relative abundance. The relationships between CRISPR-Cas abundance, viral abundance, and viral diversity are broadly consistent across environments, providing robust evidence from natural ecosystems that supports predictions of when CRISPR is beneficial. These results indicate that viral abundance and diversity are major ecological factors that drive the selection and maintenance of CRISPR-Cas in microbial ecosystems.

•

Metagenomic data from diverse ecosystems are used to analyze CRISPR prevalence

•

Environment type explains ∼a quarter of the variation in CRISPR-Cas abundance

•

There is a positive association between CRISPR-Cas abundance and viral abundance

•

CRISPR-Cas is more abundant when viral diversity is comparatively lower

Evolutionary Ecology and Interplay of Prokaryotic Innate and Adaptive Immune Systems

Like many organisms, bacteria and archaea have both innate and adaptive immune systems to defend against infection by viruses and other parasites. Innate immunity most commonly relies on the endonuclease-mediated cleavage of any incoming DNA that lacks a specific epigenetic modification, through a system known as restriction-modification. CRISPR-Cas-mediated adaptive immunity relies on the insertion of short DNA sequences from parasite genomes into CRISPR arrays on the host genome to provide sequence-specific protection. The discovery of each of these systems has revolutionised our ability to carry out genetic manipulations, and, as a consequence, the enzymes involved have been characterised in exquisite detail. In comparison, much less is known about the importance of these two arms of the defence for the ecology and evolution of prokaryotes and their parasites. Here, we review our current ecological and evolutionary understanding of these systems in isolation, and discuss the need to study how innate and adaptive immune responses are integrated when they coexist in the same cell.

Nutrient Loading and Viral Memory Drive Accumulation of Restriction Modification Systems in Bloom-Forming Cyanobacteria

The mechanisms driving cyanobacterial harmful algal blooms (HABs) like those caused by Microcystis aeruginosa remain elusive, but improved defense against viral predation has been implicated for success in eutrophic environments. Our genus-level analyses of 139,023 genomes revealed that HAB-forming cyanobacteria carry vastly more restriction modification systems per genome (RMPG) than nearly all other prokaryotic genera, suggesting that viral defense is a cornerstone of their ecological success. In contrast, picocyanobacteria that numerically dominate nutrient-poor systems have the fewest RMPG within the phylum Cyanobacteria. We used classic resource competition models to explore the hypothesis that nutrient enrichments drive ecological selection for high RMPG due to increased host-phage contact rate. These classic models, agnostic to the mechanism of defense, explain how nutrient loading can select for increased RMPG but, importantly, fail to explain the extreme accumulation of these defense systems. However, extreme accumulation of RMPG can be achieved in a novel “memory” model that accounts for a unique activity of restriction modification systems: the accidental methylation of viral DNA by the methyltransferase. The methylated virus “remembers” the RM defenses of its former host and can evade these defenses if they are present in the next host. This viral memory leads to continual RM system devaluation; RMs accumulate extensively because the benefit of each addition is diminished. Our modeling leads to the hypothesis that nutrient loading and virion methylation drive the extreme accumulation of RMPG in HAB-forming cyanobacteria. Finally, our models suggest that hosts with different RMPG values can coexist when hosts have unique sets of RM systems.

Coevolution between bacterial CRISPR-Cas systems and their bacteriophages

CRISPR-Cas systems provide bacteria and archaea with adaptive, heritable immunity against their viruses (bacteriophages and phages) and other parasitic genetic elements. CRISPR-Cas systems are highly diverse, and we are only beginning to understand their relative importance in phage defense. In this review, we will discuss when and why CRISPR-Cas immunity against phages evolves, and how this, in turn, selects for the evolution of immune evasion by phages. Finally, we will discuss our current understanding of if, and when, we observe coevolution between CRISPR-Cas systems and phages, and how this may be influenced by the mechanism of CRISPR-Cas immunity.

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.

It is unclear how important CRISPR-Cas systems are for protecting natural populations of bacteria against infections by mobile genetic elements

Articles on CRISPR commonly open with some variant of the phrase "these short palindromic repeats and their associated endonucleases (Cas) are an adaptive immune system that exists to protect bacteria and archaea from viruses and infections with other mobile genetic elements." There is an abundance of genomic data consistent with the hypothesis that CRISPR plays this role in natural populations of bacteria and archaea, and experimental demonstrations with a few species of bacteria and their phage and plasmids show that CRISPR-Cas systems can play this role in vitro. Not at all clear are the ubiquity, magnitude, and nature of the contribution of CRISPR-Cas systems to the ecology and evolution of natural populations of microbes and the strength of selection mediated by different types of phage and plasmids to the evolution and maintenance of CRISPR-Cas systems. In this perspective, with the aid of heuristic mathematical-computer simulation models, we explore the a priori conditions under which exposure to lytic and temperate phage and conjugative plasmids will select for and maintain CRISPR-Cas systems in populations of bacteria and archaea. We review the existing literature addressing these ecological and evolutionary questions and highlight the experimental and other evidence needed to fully understand the conditions responsible for the evolution and maintenance of CRISPR-Cas systems and the contribution of these systems to the ecology and evolution of bacteria, archaea, and the mobile genetic elements that infect them.

The network structure and eco-evolutionary dynamics of CRISPR-induced immune diversification

As a heritable sequence-specific adaptive immune system, CRISPR-Cas is a powerful force shaping strain diversity in host-virus systems. While the diversity of CRISPR alleles has been explored, the associated structure and dynamics of host-virus interactions have not. We explore the role of CRISPR in mediating the interplay between host-virus interaction structure and eco-evolutionary dynamics in a computational model and compare the results with three empirical datasets from natural systems. We show that the structure of the networks describing who infects whom and the degree to which strains are immune, are respectively modular (containing groups of hosts and viruses that interact strongly) and weighted-nested (specialist hosts are more susceptible to subsets of viruses that in turn also infect the more generalist hosts with many spacers matching many viruses). The dynamic interplay between these networks influences transitions between dynamical regimes of virus diversification and host control. The three empirical systems exhibit weighted-nested immunity networks, a pattern our theory shows is indicative of hosts able to suppress virus diversification. Previously missing from studies of microbial host-pathogen systems, the immunity network plays a key role in the coevolutionary dynamics.

Phage steering of antibiotic-resistance evolution in the bacterial pathogen, Pseudomonas aeruginosa

Background and objectives

Antimicrobial resistance is a growing global concern and has spurred increasing efforts to find alternative therapeutics. Bacteriophage therapy has seen near constant use in Eastern Europe since its discovery over a century ago. One promising approach is to use phages that not only reduce bacterial pathogen loads but also select for phage resistance mechanisms that trade-off with antibiotic resistance-so called 'phage steering'.

Methodology

Recent work has shown that the phage OMKO1 can interact with efflux pumps and in so doing select for both phage resistance and antibiotic sensitivity of the pathogenic bacterium Pseudomonas aeruginosa. We tested the robustness of this approach to three different antibiotics in vitro (tetracycline, erythromycin and ciprofloxacin) and one in vivo (erythromycin).

Results

We show that in vitro OMKO1 can reduce antibiotic resistance of P. aeruginosa (Washington PAO1) even in the presence of antibiotics, an effect still detectable after ca.70 bacterial generations in continuous culture with phage. Our in vivo experiment showed that phage both increased the survival times of wax moth larvae (Galleria mellonella) and increased bacterial sensitivity to erythromycin. This increased antibiotic sensitivity occurred both in lines with and without the antibiotic.

Conclusions and implications

Our study supports a trade-off between antibiotic resistance and phage sensitivity. This trade-off was maintained over co-evolutionary time scales even under combined phage and antibiotic pressure. Similarly, OMKO1 maintained this trade-off in vivo, again under dual phage/antibiotic pressure. Our findings have implications for the future clinical use of steering in phage therapies. Lay Summary: Given the rise of antibiotic-resistant bacterial infection, new approaches to treatment are urgently needed. Bacteriophages (phages) are bacterial viruses. The use of such viruses to treat infections has been in near-continuous use in several countries since the early 1900s. Recent developments have shown that these viruses are not only effective against routine infections but can also target antibiotic resistant bacteria in a novel, unexpected way. Similar to other lytic phages, these so-called 'steering phages' kill the majority of bacteria directly. However, steering phages also leave behind bacterial variants that resist the phages, but are now sensitive to antibiotics. Treatment combinations of these phages and antibiotics can now be used to greater effect than either one independently. We evaluated the impact of steering using phage OMKO1 and a panel of three antibiotics on Pseudomonas aeruginosa, an important pathogen in hospital settings and in people with cystic fibrosis. Our findings indicate that OMKO1, either alone or in combination with antibiotics, maintains antibiotic sensitivity both in vitro and in vivo, giving hope that phage steering will be an effective treatment option against antibiotic-resistant bacteria.

Promises and Pitfalls of In Vivo Evolution to Improve Phage Therapy

Phage therapy is the use of bacterial viruses (phages) to treat bacterial infections, a medical intervention long abandoned in the West but now experiencing a revival. Currently, therapeutic phages are often chosen based on limited criteria, sometimes merely an ability to plate on the pathogenic bacterium. Better treatment might result from an informed choice of phages. Here we consider whether phages used to treat the bacterial infection in a patient may specifically evolve to improve treatment on that patient or benefit subsequent patients. With mathematical and computational models, we explore in vivo evolution for four phage properties expected to influence therapeutic success: generalized phage growth, phage decay rate, excreted enzymes to degrade protective bacterial layers, and growth on resistant bacteria. Within-host phage evolution is strongly aligned with treatment success for phage decay rate but only partially aligned for phage growth rate and growth on resistant bacteria. Excreted enzymes are mostly not selected for treatment success. Even when evolution and treatment success are aligned, evolution may not be rapid enough to keep pace with bacterial evolution for maximum benefit. An informed use of phages is invariably superior to naive reliance on within-host evolution.

Molecular noise of innate immunity shapes bacteria-phage ecologies

Mathematical models have been used successfully at diverse scales of biological organization, ranging from ecology and population dynamics to stochastic reaction events occurring between individual molecules in single cells. Generally, many biological processes unfold across multiple scales, with mutations being the best studied example of how stochasticity at the molecular scale can influence outcomes at the population scale. In many other contexts, however, an analogous link between micro- and macro-scale remains elusive, primarily due to the challenges involved in setting up and analyzing multi-scale models. Here, we employ such a model to investigate how stochasticity propagates from individual biochemical reaction events in the bacterial innate immune system to the ecology of bacteria and bacterial viruses. We show analytically how the dynamics of bacterial populations are shaped by the activities of immunity-conferring enzymes in single cells and how the ecological consequences imply optimal bacterial defense strategies against viruses. Our results suggest that bacterial populations in the presence of viruses can either optimize their initial growth rate or their population size, with the first strategy favoring simple immunity featuring a single restriction modification system and the second strategy favoring complex bacterial innate immunity featuring several simultaneously active restriction modification systems.