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 (HGT). However, how and when new immune systems can become established in a bacterial population has remained largely unstudied. We developed a joint epidemiological and evolutionary model that predicts the conditions necessary for the spread of a CRISPR-Cas immune system into a bacterial population lacking this system. We find 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 using Pseudomonas aeruginosa PA14 and phage DMS3vir as a model. 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 establish in bacterial populations after a horizontal gene transfer event, both in ecological and clinical contexts.

Evolution and spread of multiadapted pathogens in a spatially heterogeneous environment.

Pathogen adaptation to multiple selective pressures challenges our ability to control their spread. Here we analyze the evolutionary dynamics of pathogens spreading in a heterogeneous host population where selection varies periodically in space. We study both the transient dynamics taking place at the front of the epidemic and the long-term evolution far behind the front. We identify five types of epidemic profiles arising for different levels of spatial heterogeneity and different costs of adaptation. In particular, we identify the conditions where a generalist pathogen carrying multiple adaptations can outrace a coalition of specialist pathogens. We also show that finite host populations promote the spread of generalist pathogens because demographic stochasticity enhances the extinction of locally maladapted pathogens. But higher mutation rates between genotypes can rescue the coalition of specialists and speed up the spread of epidemics for intermediate levels of spatial heterogeneity. Our work provides a comprehensive analysis of the interplay between migration, local selection, mutation, and genetic drift on the spread and on the evolution of pathogens in heterogeneous environments. This work extends our fundamental understanding of the outcome of the competition between two specialists and a generalist strategy (single- vs. multiadapted pathogens). These results have practical implications for the design of more durable control strategies against multiadapted pathogens in agriculture and in public health.

Building pyramids against the evolutionary emergence of pathogens.

Mutations allowing pathogens to escape host immunity promote the spread of infectious diseases in heterogeneous host populations and can lead to major epidemics. Understanding the conditions that slow down this evolution is key for the development of durable control strategies against pathogens. Here, we use theory and experiments to compare the efficacy of three strategies for the deployment of resistance: (i) a mixing strategy where the host population contains two single-resistant genotypes, (ii) a pyramiding strategy where the host carries a double-resistant genotype, (iii) a combining strategy where the host population is a mix of a single-resistant genotype and a double-resistant genotype. First, we use evolutionary epidemiology theory to clarify the interplay between demographic stochasticity and evolutionary dynamics to show that the pyramiding strategy always yields lower probability of evolutionary emergence. Second, we test experimentally these predictions with the introduction of bacteriophages into bacterial populations where we manipulated the diversity and the depth of immunity using a Clustered Regularly Interspaced Short Palindromic Repeats-CRISPR associated (CRISPR-Cas) system. These biological assays confirm that pyramiding multiple defences into the same host genotype and avoiding combination with single-defence genotypes is a robust way to reduce pathogen evolutionary emergence. The experimental validation of these theoretical recommendations has practical implications in various areas, including for the optimal deployment of resistance varieties in agriculture and for the design of durable vaccination strategies.

Effect of unequal vaccination coverage and migration on long-term pathogen evolution in a metapopulation.

Worldwide inequalities in vaccine availability are expected to affect the spread and spatial distribution of infectious diseases. It is unclear, however, how spatial variation in vaccination coverage can affect the long-term evolution of pathogens. Here we use an analytical model and numerical simulations to analyse the influence of different imperfect vaccines on the potential evolution of pathogen virulence in a two-population model where vaccination coverage varies between populations. We focus on four vaccines, with different modes of action on the life cycle of a pathogen infecting two host populations coupled by migration. We show that, for vaccines that reduce infection risk or transmissibility, spatial heterogeneity has little effect on pathogen prevalence and host mortality, and no effect on the evolution of pathogen virulence. In contrast, vaccines that reduce pathogen virulence can select for more virulent pathogens and may lead to the coexistence of different pathogen strains, depending on the degree of spatial heterogeneity in the metapopulation. This heterogeneity is driven by two parameters: pathogen migration and the difference in the vaccination rate between the two populations. We show that vaccines that only reduce pathogen virulence select mainly for a single pathogen strategy in the long term, while vaccines that reduce both transmission and virulence can favor the coexistence of two pathogen genotypes. We discuss the implications and potential extensions of our analysis.

Seasonality and the persistence of vector-borne pathogens.

Many vector-borne diseases are affected by the seasonality of the environment. Yet, seasonality can act on distinct steps of the life cycle of the pathogen and it is often difficult to predict the influence of the periodic fluctuations of the environment on the basic reproduction ratio R0 of vector-borne pathogens. Here, we analyse a general vector-borne disease model and we account for periodic fluctuations of different components of the pathogen's life cycle. We develop a perturbation analysis framework to obtain useful approximations to evaluate the overall consequences of seasonality on the R0 of the pathogen. This analysis reveals when seasonality is expected to increase or to decrease pathogen persistence. We show that seasonality in vector density or in the biting rate of the vector can have opposite effects on persistence and we provide a useful biological explanation for this result based on the covariance between key compartments of the epidemiological model. This framework could be readily extended to explore the influence of seasonality on other components of the life cycle of vector-borne pathogens.

Transient eco-evolutionary dynamics early in a phage epidemic have strong and lasting impact on the long-term evolution of bacterial defences.

Organisms have evolved a range of constitutive (always active) and inducible (elicited by parasites) defence mechanisms, but we have limited understanding of what drives the evolution of these orthogonal defence strategies. Bacteria and their phages offer a tractable system to study this: Bacteria can acquire constitutive resistance by mutation of the phage receptor (surface mutation, sm) or induced resistance through their CRISPR-Cas adaptive immune system. Using a combination of theory and experiments, we demonstrate that the mechanism that establishes first has a strong advantage because it weakens selection for the alternative resistance mechanism. As a consequence, ecological factors that alter the relative frequencies at which the different resistances are acquired have a strong and lasting impact: High growth conditions promote the evolution of sm resistance by increasing the influx of receptor mutation events during the early stages of the epidemic, whereas a high infection risk during this stage of the epidemic promotes the evolution of CRISPR immunity, since it fuels the (infection-dependent) acquisition of CRISPR immunity. This work highlights the strong and lasting impact of the transient evolutionary dynamics during the early stages of an epidemic on the long-term evolution of constitutive and induced defences, which may be leveraged to manipulate phage resistance evolution in clinical and applied settings.

Phage-antibiotic synergy: Cell filamentation is a key driver of successful phage predation.

Phages are promising tools to fight antibiotic-resistant bacteria, and as for now, phage therapy is essentially performed in combination with antibiotics. Interestingly, combined treatments including phages and a wide range of antibiotics lead to an increased bacterial killing, a phenomenon called phage-antibiotic synergy (PAS), suggesting that antibiotic-induced changes in bacterial physiology alter the dynamics of phage propagation. Using single-phage and single-cell techniques, each step of the lytic cycle of phage HK620 was studied in E. coli cultures treated with either ceftazidime, cephalexin or ciprofloxacin, three filamentation-inducing antibiotics. In the presence of sublethal doses of antibiotics, multiple stress tolerance and DNA repair pathways are triggered following activation of the SOS response. One of the most notable effects is the inhibition of bacterial division. As a result, a significant fraction of cells forms filaments that stop dividing but have higher rates of mutagenesis. Antibiotic-induced filaments become easy targets for phages due to their enlarged surface areas, as demonstrated by fluorescence microscopy and flow cytometry techniques. Adsorption, infection and lysis occur more often in filamentous cells compared to regular-sized bacteria. In addition, the reduction in bacterial numbers caused by impaired cell division may account for the faster elimination of bacteria during PAS. We developed a mathematical model to capture the interaction between sublethal doses of antibiotics and exposition to phages. This model shows that the induction of filamentation by sublethal doses of antibiotics can amplify the replication of phages and therefore yield PAS. We also use this model to study the consequences of PAS on the emergence of antibiotic resistance. A significant percentage of hyper-mutagenic filamentous bacteria are effectively killed by phages due to their increased susceptibility to infection. As a result, the addition of even a very low number of bacteriophages produced a strong reduction of the mutagenesis rate of the entire bacterial population. We confirm this prediction experimentally using reporters for bacterial DNA repair. Our work highlights the multiple benefits associated with the combination of sublethal doses of antibiotics with bacteriophages.

Phenotypic evolution of SARS-CoV-2: a statistical inference approach.

Since its emergence in late 2019, the SARS-CoV-2 virus has spread globally, causing the ongoing COVID-19 pandemic. In the fall of 2020, the Alpha variant (lineage B.1.1.7) was detected in England and spread rapidly, outcompeting the previous lineage. Yet, very little is known about the underlying modifications of the infection process that can explain this selective advantage. Here, we try to quantify how the Alpha variant differed from its predecessor on two phenotypic traits: The transmission rate and the duration of infectiousness. To this end, we analyzed the joint epidemiological and evolutionary dynamics as a function of the Stringency Index, a measure of the amount of Non-Pharmaceutical Interventions. Assuming that these control measures reduce contact rates and transmission, we developed a two-step approach based on ${{SEIR}}$ models and the analysis of a combination of epidemiological and evolutionary information. First, we quantify the link between the Stringency Index and the reduction in viral transmission. Second, based on a novel theoretical derivation of the selection gradient in an ${{SEIR}}$ model, we infer the phenotype of the Alpha variant from its frequency changes. We show that its selective advantage is more likely to result from a higher transmission than from a longer infectious period. Our work illustrates how the analysis of the joint epidemiological and evolutionary dynamics of infectious diseases can help understand the phenotypic evolution driving pathogen adaptation.

Bottlenecks can constrain and channel evolutionary paths.

Population bottlenecks are commonplace in experimental evolution, specifically in serial passaging experiments where microbial populations alternate between growth and dilution. Natural populations also experience such fluctuations caused by seasonality, resource limitation, or host-to-host transmission for pathogens. Yet, how unlimited growth with periodic bottlenecks influence the adaptation of populations is not fully understood. Here, we study theoretically the effects of bottlenecks on the accessibility of evolutionary paths and on the rate of evolution. We model an asexual population evolving on a minimal fitness landscape consisting of two types of beneficial mutations with the empirically supported trade-off between mutation rate and fitness advantage, in the regime where multiple beneficial mutations may segregate simultaneously. In the limit of large population sizes and small mutation rates, we show the existence of a unique most likely evolutionary scenario, determined by the size of the wild-type population at the beginning and at the end of each cycle. These two key demographic parameters determine which adaptive paths may be taken by the evolving population by controlling the supply of mutants during growth and the loss of mutants at the bottleneck. We do not only show that bottlenecks act as a deterministic control of evolutionary paths but also that each possible evolutionary scenario can be forced to occur by tuning demographic parameters. This work unveils the effects of demography on adaptation of periodically bottlenecked populations and can guide the design of evolution experiments.

Pathogen evolution during vaccination campaigns.

Following the initiation of the unprecedented global vaccination campaign against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), attention has now turned to the potential impact of this large-scale intervention on the evolution of the virus. In this Essay, we summarize what is currently known about pathogen evolution in the context of immune priming (including vaccination) from research on other pathogen species, with an eye towards the future evolution of SARS-CoV-2.

Competition and coevolution drive the evolution and the diversification of CRISPR immunity.

The diversity of resistance challenges the ability of pathogens to spread and to exploit host populations. Yet, how this host diversity evolves over time remains unclear because it depends on the interplay between intraspecific competition among host genotypes and coevolution with pathogens. Here we study experimentally the effect of coevolving phage populations on the diversification of bacterial CRISPR immunity across space and time. We demonstrate that the negative-frequency-dependent selection generated by coevolution is a powerful force that maintains host resistance diversity and selects for new resistance mutations in the host. We also find that host evolution is driven by asymmetries in competitive abilities among different host genotypes. Even if the fittest host genotypes are targeted preferentially by the evolving phages, they often escape extinctions through the acquisition of new CRISPR immunity. Together, these fluctuating selective pressures maintain diversity, but not by preserving the pre-existing host composition. Instead, we repeatedly observe the introduction of new resistance genotypes stemming from the fittest hosts in each population. These results highlight the importance of competition on the transient dynamics of host-pathogen coevolution.

Evolution of class-structured populations in periodic environments.

What is the influence of periodic environmental fluctuations on life-history evolution? We present a general theoretical framework to understand and predict the long-term evolution of life-history traits under a broad range of ecological scenarios. Specifically, we investigate how periodic fluctuations affect selection when the population is also structured in distinct classes. This analysis yields time-varying selection gradients that clarify the influence of the fluctuations of the environment on the competitive ability of a specific life-history mutation. We use this framework to analyse the evolution of key life-history traits of pathogens. We examine three different epidemiological scenarios and we show how periodic fluctuations of the environment can affect the evolution of virulence and transmission as well as the preference for different hosts. These examples yield new and testable predictions on pathogen evolution, and illustrate how our approach can provide a better understanding of the evolutionary consequences of time-varying environmental fluctuations in a broad range of scenarios.

Phage-host coevolution in natural populations.

Coevolution between bacteriophages (phages) and their bacterial hosts occurs through changes in resistance and counter-resistance mechanisms. To assess phage-host evolution in wild populations, we isolated 195 Vibrio crassostreae strains and 243 vibriophages during a 5-month time series from an oyster farm and combined these isolates with existing V. crassostreae and phage isolates. Cross-infection studies of 81,926 host-phage pairs delineated a modular network where phages are best at infecting co-occurring hosts, indicating local adaptation. Successful propagation of phage is restricted by the ability to adsorb to closely related bacteria and further constrained by strain-specific defence systems. These defences are highly diverse and predominantly located on mobile genetic elements, and multiple defences are active within a single genome. We further show that epigenetic and genomic modifications enable phage to adapt to bacterial defences and alter host range. Our findings reveal that the evolution of bacterial defences and phage counter-defences is underpinned by frequent genetic exchanges with, and between, mobile genetic elements.

Effects of epistasis and recombination between vaccine-escape and virulence alleles on the dynamics of pathogen adaptation.

Pathogen adaptation to public health interventions such as vaccination may take tortuous routes and involve multiple mutations at different locations in the pathogen genome, acting on distinct phenotypic traits. Yet how these multi-locus adaptations jointly evolve is poorly understood. Here we consider the joint evolution of two adaptations: pathogen escape from the vaccine-induced immune response and adjustments to pathogen virulence affecting transmission or clearance. We elucidate the role played by epistasis and recombination, with an emphasis on the different protective effects of vaccination. We show that vaccines blocking infection, reducing transmission and/or increasing clearance generate positive epistasis between the vaccine-escape and virulence alleles, favouring strains that carry both mutations, whereas vaccines reducing virulence mortality generate negative epistasis, favouring strains that carry either mutation but not both. High rates of recombination can affect these predictions. If epistasis is positive, frequent recombination can prevent the transient build-up of more virulent escape strains. If epistasis is negative, frequent recombination between loci can create an evolutionary bistability, favouring whichever adaptation is more accessible. Our work provides a timely alternative to the variant-centred perspective on pathogen adaptation and captures the effect of different types of vaccine on the interference between multiple adaptive mutations.

Targeted vaccination and the speed of SARS-CoV-2 adaptation.

The limited supply of vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) raises the question of targeted vaccination. Many countries have opted to vaccinate older and more sensitive hosts first to minimize the disease burden. However, what are the evolutionary consequences of targeted vaccination? We clarify the consequences of different vaccination strategies through the analysis of the speed of viral adaptation measured as the rate of change of the frequency of a vaccine-adapted variant. We show that such a variant is expected to spread faster if vaccination targets individuals who are likely to be involved in a higher number of contacts. We also discuss the pros and cons of dose-sparing strategies. Because delaying the second dose increases the proportion of the population vaccinated with a single dose, this strategy can both speed up the spread of the vaccine-adapted variant and reduce the cumulative number of deaths. Hence, strategies that are most effective at slowing viral adaptation may not always be epidemiologically optimal. A careful assessment of both the epidemiological and evolutionary consequences of alternative vaccination strategies is required to determine which individuals should be vaccinated first.

Regulation of prophage induction and lysogenization by phage communication systems.

Many viruses cause both lytic infections, where they release viral particles, and dormant infections, where they await future opportunities to reactivate.1 The benefits of each transmission mode depend on the density of susceptible hosts in the environment.2-4 Some viruses infecting bacteria use molecular signaling to respond plastically to changes in host availability.5 These viruses produce a signal during lytic infection and regulate, based on the signal concentration in the environment, the probability with which they switch to causing dormant infections.5,6 We present an analytical framework to examine the adaptive significance of plasticity in viral life-history traits in fluctuating environments. Our model generalizes and extends previous theory7 and predicts that host density fluctuations should select for plasticity in entering lysogeny as well as virus reactivation once signal concentrations decline. Using Bacillus subtilis and its phage phi3T, we experimentally confirm the prediction that phages use signal to make informed decisions over prophage induction. We also demonstrate that lysogens produce signaling molecules and that signal is degraded by hosts in a density-dependent manner. Declining signal concentrations therefore potentially indicate the presence of uninfected hosts and trigger prophage induction. Finally, we find that conflict over the responses of lysogenization and reactivation to signal is resolved through the evolution of different response thresholds for each trait. Collectively, these findings deepen our understanding of the ways viruses use molecular communication to regulate their infection strategies, which can be leveraged to manipulate host and phage population dynamics in natural environments.

Understanding the evolution of multiple drug resistance in structured populations.

The evolution of multidrug resistance (MDR) is a pressing public health concern. Yet many aspects, such as the role played by population structure, remain poorly understood. Here, we argue that studying MDR evolution by focusing upon the dynamical equations for linkage disequilibrium (LD) can greatly simplify the calculations, generate more insight, and provide a unified framework for understanding the role of population structure. We demonstrate how a general epidemiological model of MDR evolution can be recast in terms of the LD equations. These equations reveal how the different forces generating and propagating LD operate in a dynamical setting at both the population and metapopulation levels. We then apply these insights to show how the LD perspective: (i) explains equilibrium patterns of MDR, (ii) provides a simple interpretative framework for transient evolutionary dynamics, and (iii) can be used to assess the consequences of different drug prescription strategies for MDR evolution.

Gene regulation of the avian malaria parasite Plasmodium relictum, during the different stages within the mosquito vector.

The malaria parasite Plasmodium relictum is one of the most widespread species of avian malaria. As in the case of its human counterparts, bird Plasmodium undergoes a complex life cycle infecting two hosts: the arthropod vector and the vertebrate host. In this study, we examined transcriptomes of P. relictum (SGS1) during crucial timepoints within its vector, Culex pipiens quinquefasciatus. Differential gene-expression analyses identified genes linked to the parasites life-stages at: i) a few minutes after the blood meal is ingested, ii) during peak oocyst production phase, iii) during peak sporozoite phase and iv) during the late-stages of the infection. A large amount of genes coding for functions linked to host-immune invasion and multifunctional genes was active throughout the infection cycle. One gene associated with a conserved Plasmodium membrane protein with unknown function was upregulated throughout the parasite development in the vector, suggesting an important role in the successful completion of the sporogonic cycle. Gene expression analysis further identified genes, with unknown functions to be significantly differentially expressed during the infection in the vector as well as upregulation of reticulocyte-binding proteins, which raises the possibility of the multifunctionality of these RBPs. We establish the existence of highly stage-specific pathways being overexpressed during the infection. This first study of gene-expression of a non-human Plasmodium species in its vector provides a comprehensive insight into the molecular mechanisms of the common avian malaria parasite P. relictum and provides essential information on the evolutionary diversity in gene regulation of the Plasmodium's vector stages.

Phage gene expression and host responses lead to infection-dependent costs of CRISPR immunity.

CRISPR-Cas immune systems are widespread in bacteria and archaea, but not ubiquitous. Previous work has demonstrated that CRISPR immunity is associated with an infection-induced fitness cost, which may help explain the patchy distribution observed. However, the mechanistic basis of this cost has remained unclear. Using Pseudomonas aeruginosa PA14 and its phage DMS3vir as a model, we perform a 30-day evolution experiment under phage mediated selection. We demonstrate that although CRISPR is initially selected for, bacteria carrying mutations in the phage receptor rapidly invade the population following subsequent reinfections. We then test three potential mechanisms for the observed cost of CRISPR: (1) autoimmunity from the acquisition of self-targeting spacers, (2) immunopathology or energetic costs from increased cas gene expression and (3) toxicity caused by phage gene expression prior to CRISPR-mediated cleavage. We find that phages can express genes before the immune system clears the infection and that expression of these genes can have a negative effect on host fitness. While infection does not lead to increased expression of cas genes, it does cause differential expression of multiple other host processes that may further contribute to the cost of CRISPR immunity. In contrast, we found little support for infection-induced autoimmunological and immunopathological effects. Phage gene expression prior to cleavage of the genome by the CRISPR-Cas immune system is therefore the most parsimonious explanation for the observed phage-induced fitness cost.