Well-temperate phage: optimal bet-hedging against local environmental collapses

Lytic and temperate phage naturally coexist in a dynamic population model

When phage infect their bacterial hosts, they may either lyse the cell and generate a burst of new phage, or lysogenize the bacterium, incorporating the phage genome into it. Phage lysis/lysogeny strategies are assumed to be highly optimized, with the optimal tradeoff depending on environmental conditions. However, in nature, phage of radically different lysis/lysogeny strategies coexist in the same environment, preying on the same bacteria. How can phage preying on the same bacteria coexist if one is more optimal than the other? Here, we address this conundrum within a modeling framework, simulating the population dynamics of communities of phage and their lysogens. We find that coexistence between phage of different lysis/lysogeny strategies is a natural outcome of chaotic population dynamics that arise within sufficiently diverse communities, which ensure no phage is able to absolutely dominate its competitors. Our results further suggest a bet-hedging mechanism at the level of the phage pan-genome, wherein obligate lytic (virulent) strains typically outcompete temperate strains, but also more readily fluctuate to extinction within a local community.

Optimal phenotypic adaptation in fluctuating environments

Phenotypic adaptation is a universal feature of biological systems navigating highly variable environments. Recent empirical data support the role of memory-driven decision making in cellular systems navigating uncertain future nutrient landscapes, wherein a distinct growth phenotype emerges in fluctuating conditions. We develop a simple stochastic mathematical model to describe memory-driven cellular adaptation required for systems to optimally navigate such uncertainty. In this framework, adaptive populations traverse dynamic environments by inferring future variation from a memory of prior states, and memory capacity imposes a fundamental trade-off between the speed and accuracy of adaptation to new fluctuating environments. Our results suggest that the observed growth reductions that occur in fluctuating environments are a direct consequence of optimal decision making and result from bet hedging and occasional phenotypic-environmental mismatch. We anticipate that this modeling framework will be useful for studying the role of memory in phenotypic adaptation, including in the design of temporally varying therapies against adaptive systems.

The Unique Genome of the Virus and Alternative Strategies for its Realization

Dedicated to the 130th anniversary of Dmitry Ivanovsky's discovery of the virus kingdom as a new form of biological life. The genome of some RNA-containing viruses comprises ambipolar genes that are arranged in stacks (one above the other) encoding proteins in opposite directions. Ambipolar genes provide a new approach for developing viral diversity when virions possessing an identical genome may differ in its expression scheme (strategy) and have distinct types of progeny virions varying in the genomic RNA polarity and the composition of proteins expressed by positive- or negative-sense genes, the so-called ambipolar virions. So far, this pathway of viral genome expression remains hypothetical and hidden from us, like the dark side of the Moon, and deserves a detailed study.

Phenotypic flux: The role of physiology in explaining the conundrum of bacterial persistence amid phage attack

Bacteriophages, the viruses of bacteria, have been studied for over a century. They were not only instrumental in laying the foundations of molecular biology, but they are also likely to play crucial roles in shaping our biosphere and may offer a solution to the control of drug-resistant bacterial infections. However, it remains challenging to predict the conditions for bacterial eradication by phage predation, sometimes even under well-defined laboratory conditions, and, most curiously, if the majority of surviving cells are genetically phage-susceptible. Here, I propose that even clonal phage and bacterial populations are generally in a state of continuous 'phenotypic flux', which is caused by transient and nongenetic variation in phage and bacterial physiology. Phenotypic flux can shape phage infection dynamics by reducing the force of infection to an extent that allows for coexistence between phages and susceptible bacteria. Understanding the mechanisms and impact of phenotypic flux may be key to providing a complete picture of phage-bacteria coexistence. I review the empirical evidence for phenotypic variation in phage and bacterial physiology together with the ways they have been modeled and discuss the potential implications of phenotypic flux for ecological and evolutionary dynamics between phages and bacteria, as well as for phage therapy.

Mutualistic interplay between bacteriophages and bacteria in the human gut

Bacteriophages (phages) are often described as obligate predators of their bacterial hosts, and phage predation is one of the leading forces controlling the density and distribution of bacterial populations. Every 48 h half of all bacteria on Earth are killed by phages. Efficient killing also forms the basis of phage therapy in humans and animals and the use of phages as food preservatives. In turn, bacteria have a plethora of resistance systems against phage attack, but very few bacterial species, if any, have entirely escaped phage predation. However, in complex communities and environments such as the human gut, this antagonistic model of attack and counter-defence does not fully describe the scope of phage-bacterium interactions. In this Review, we explore some of the more mutualistic aspects of phage-bacterium interactions in the human gut, and we suggest that the relationship between phages and their bacterial hosts in the gut is best characterized not as a fight to the death between enemies but rather as a mutualistic relationship between partners.

Quiescence Generates Moving Average in a Stochastic Epidemiological Model with One Host and Two Parasites

Mathematical modelling of epidemiological and coevolutionary dynamics is widely being used to improve disease management strategies of infectious diseases. Many diseases present some form of intra-host quiescent stage, also known as covert infection, while others exhibit dormant stages in the environment. As quiescent/dormant stages can be resistant to drug, antibiotics, fungicide treatments, it is of practical relevance to study the influence of these two life-history traits on the coevolutionary dynamics. We develop first a deterministic coevolutionary model with two parasite types infecting one host type and study analytically the stability of the dynamical system. We specifically derive a stability condition for a five-by-five system of equations with quiescence. Second, we develop a stochastic version of the model to study the influence of quiescence on stochasticity of the system dynamics. We compute the steady state distribution of the parasite types which follows a multivariate normal distribution. Furthermore, we obtain numerical solutions for the covariance matrix of the system under symmetric and asymmetric quiescence rates between parasite types. When parasite strains are identical, quiescence increases the variance of the number of infected individuals at high transmission rate and vice versa when the transmission rate is low. However, when there is competition between parasite strains with different quiescent rates, quiescence generates a moving average behaviour which dampen off stochasticity and decreases the variance of the number of infected hosts. The strain with the highest rate of entering quiescence determines the strength of the moving average and the magnitude of reduction of stochasticity. Thus, it is worth investigating simple models of multi-strain parasite under quiescence/dormancy to improve disease management strategies.

Alternating lysis and lysogeny is a winning strategy in bacteriophages due to Parrondo's paradox

Bacteriophages, the most widespread reproducing biological entity on Earth, employ two strategies of virus–host interaction: lysis of the host cell and lysogeny whereby the virus genome integrates into the host genome and propagates vertically with it. We present a population model that reveals an effect known as Parrondo’s paradox in game theory: Alternating between lysis and lysogeny is a winning strategy for a bacteriophage, even when each strategy individually is at a disadvantage compared with a competing bacteriophage. Thus, evolution of bacteriophages appears to optimize the ratio between the lysis and lysogeny propensities rather than the phage burst size in any individual phase. This phenomenon is likely to be relevant for understanding evolution of other host–parasites systems.

Temperate bacteriophages lyse or lysogenize host cells depending on various parameters of infection, a key one being the ratio of the number of free viruses to the number of host cells. However, the effect of different propensities of phages for lysis and lysogeny on phage fitness remains an open problem. We explore a nonlinear dynamic evolution model of competition between two phages, one of which is disadvantaged in both the lytic and lysogenic phases. We show that the disadvantaged phage can win the competition by alternating between the lytic and lysogenic phases, each of which individually is a “loser.” This counterintuitive result is analogous to Parrondo’s paradox in game theory, whereby individually losing strategies combine to produce a winning outcome. The results suggest that evolution of phages optimizes the ratio between the lysis and lysogeny propensities rather than the phage burst size in any individual phase. These findings are likely to broadly apply to the evolution of host–parasite interactions.

OUP accepted manuscript

Temperate phage can initiate lysis or lysogeny after infecting a bacterial host. The genetic switch between lysis and lysogeny is mediated by phage regulatory genes as well as host and environmental factors. Recently, a new class of decision switches was identified in phage of the SPbeta group, mediated by the extracellular release of small, phage-encoded peptides termed arbitrium. Arbitrium peptides can be taken up by bacteria prior to infection, modulating the decision switch in the event of a subsequent phage infection. Increasing the concentration of arbitrium increases the chance that a phage infection will lead to lysogeny, rather than lysis. Although prior work has centered on the molecular mechanisms of arbitrium-induced switching, here we focus on how selective pressures impact the benefits of plasticity in switching responses. In this work, we examine the possible advantages of near-term adaptation of communication-based decision switches used by the SPbeta-like group. We combine a nonlinear population model with a control-theoretic approach to evaluate the relationship between a putative phage reaction norm (i.e. the probability of lysogeny as a function of arbitrium) and the extent of phage reproduction at a near-term time horizon. We measure phage reproduction in terms of a cellular-level metric previously shown to enable comparisons of near-term phage fitness across a continuum from lysis to latency. We show the adaptive potential of communication-based lysis–lysogeny responses and find that optimal switching between lysis and lysogeny increases the near-term phage reproduction compared to fixed responses, further supporting both molecular- and model-based analyses of the putative benefits of this class of decision switches. We further find that plastic responses are robust to the inclusion of cellular-level stochasticity, variation in life history traits, and variation in resource availability. These findings provide further support to explore the long-term evolution of plastic decision systems mediated by extracellular decision-signaling molecules and the feedback between phage reaction norms and ecological context.

Principles of seed banks and the emergence of complexity from dormancy

Across the tree of life, populations have evolved the capacity to contend with suboptimal conditions by engaging in dormancy, whereby individuals enter a reversible state of reduced metabolic activity. The resulting seed banks are complex, storing information and imparting memory that gives rise to multi-scale structures and networks spanning collections of cells to entire ecosystems. We outline the fundamental attributes and emergent phenomena associated with dormancy and seed banks, with the vision for a unifying and mathematically based framework that can address problems in the life sciences, ranging from global change to cancer biology.

The Baltimore Classification of Viruses 50 Years Later: How Does It Stand in the Light of Virus Evolution?

Fifty years ago, David Baltimore published a brief conceptual paper delineating the classification of viruses by the routes of genome expression. The six "Baltimore classes" of viruses, with a subsequently added 7th class, became the conceptual framework for the development of virology during the next five decades. During this time, it became clear that the Baltimore classes, with relatively minor additions, indeed cover the diversity of virus genome expression schemes that also define the replication cycles. Here, we examine the status of the Baltimore classes 50 years after their advent and explore their links with the global ecology and biology of the respective viruses. We discuss an extension of the Baltimore scheme and why many logically admissible expression-replication schemes do not appear to be realized in nature. Recent phylogenomic analyses allow tracing the complex connections between the Baltimore classes and the monophyletic realms of viruses. The five classes of RNA viruses and reverse-transcribing viruses share an origin, whereas both the single-stranded DNA viruses and double-stranded DNA (dsDNA) viruses evolved on multiple independent occasions. Most of the Baltimore classes of viruses probably emerged during the earliest era of life evolution, at the stage of the primordial pool of diverse replicators, and before the advent of modern-like cells with large dsDNA genomes. The Baltimore classes remain an integral part of the conceptual foundation of biology, providing the essential structure for the logical space of information transfer processes, which is nontrivially connected with the routes of evolution of viruses and other replicators.

Enhanced mutualistic symbiosis between soil phages and bacteria with elevated chromium-induced environmental stress

BACKGROUND

Microbe-virus interactions have broad implications on the composition, function, and evolution of microbiomes. Elucidating the effects of environmental stresses on these interactions is critical to identify the ecological function of viral communities and understand microbiome environmental adaptation. Heavy metal-contaminated soils represent a relevant ecosystem to study the interplay between microbes, viruses, and environmental stressors.

RESULTS

Metagenomic analysis revealed that Cr pollution adversely altered the abundance, diversity, and composition of viral and bacterial communities. Host-phage linkage based on CRISPR indicated that, in soils with high Cr contamination, the abundance of phages associated with heavy metal-tolerant hosts increased, as did the relative abundance of phages with broad host ranges (identified as host-phage linkages across genera), which would facilitate transfection and broader distribution of heavy metal resistance genes in the bacterial community. Examining variations along the pollutant gradient, enhanced mutualistic phage-bacterium interactions were observed in the face of greater environmental stresses. Specifically, the fractions of lysogens in bacterial communities (identified by integrase genes within bacterial genomes and prophage induction assay by mitomycin-C) were positively correlated with Cr contamination levels. Furthermore, viral genomic analysis demonstrated that lysogenic phages under higher Cr-induced stresses carried more auxiliary metabolic genes regulating microbial heavy metal detoxification.

CONCLUSION

With the intensification of Cr-induced environmental stresses, the composition, replication strategy, and ecological function of the phage community all evolve alongside the bacterial community to adapt to extreme habitats. These result in a transformation of the phage-bacterium interaction from parasitism to mutualism in extreme environments and underscore the influential role of phages in bacterial adaptation to pollution-related stress and in related biogeochemical processes. Video Abstract.

From Bench to Keyboard and Back Again: A Brief History of Lambda Phage Modeling

Cellular decision making is the process whereby cells choose one developmental pathway from multiple possible ones, either spontaneously or due to environmental stimuli. Examples in various cell types suggest an almost inexhaustible plethora of underlying molecular mechanisms. In general, cellular decisions rely on the gene regulatory network, which integrates external signals to drive cell fate choice. The search for general principles of such a process benefits from appropriate biological model systems that reveal how and why certain gene regulatory mechanisms drive specific cellular decisions according to ecological context and evolutionary outcomes. In this article, we review the historical and ongoing development of the phage lambda lysis-lysogeny decision as a model system to investigate all aspects of cellular decision making. The unique generality, simplicity, and richness of phage lambda decision making render it a constant source ofmathematical modeling-aided inspiration across all of biology. We discuss the origins and progress of quantitative phage lambda modeling from the 1950s until today, as well as its possible future directions. We provide examples of how modeling enabled methods and theory development, leading to new biological insights by revealing gaps in the theory and pinpointing areas requiring further experimental investigation. Overall, we highlight the utility of theoretical approaches both as predictive tools, to forecast the outcome of novel experiments, and as explanatory tools, to elucidate the natural processes underlying experimental data.

When to wake up? The optimal waking-up strategies for starvation-induced persistence

Prolonged lag time can be induced by starvation contributing to the antibiotic tolerance of bacteria. We analyze the optimal lag time to survive and grow the iterative and stochastic application of antibiotics. A simple model shows that the optimal lag time can exhibit a discontinuous transition when the severeness of the antibiotic application, such as the probability to be exposed the antibiotic, the death rate under the exposure, and the duration of the exposure, is increased. This suggests the possibility of reducing tolerant bacteria by controlled usage of antibiotics application. When the bacterial populations are able to have two phenotypes with different lag times, the fraction of the second phenotype that has different lag time shows a continuous transition. We then present a generic framework to investigate the optimal lag time distribution for total population fitness for a given distribution of the antibiotic application duration. The obtained optimal distributions have multiple peaks for a wide range of the antibiotic application duration distributions, including the case where the latter is monotonically decreasing. The analysis supports the advantage in evolving multiple, possibly discrete phenotypes in lag time for bacterial long-term fitness.

Repeated outbreaks drive the evolution of bacteriophage communication

Recently, a small-molecule communication mechanism was discovered in a range of Bacillus-infecting bacteriophages, which these temperate phages use to inform their lysis-lysogeny decision. We present a mathematical model of the ecological and evolutionary dynamics of such viral communication and show that a communication strategy in which phages use the lytic cycle early in an outbreak (when susceptible host cells are abundant) but switch to the lysogenic cycle later (when susceptible cells become scarce) is favoured over a bet-hedging strategy in which cells are lysogenised with constant probability. However, such phage communication can evolve only if phage-bacteria populations are regularly perturbed away from their equilibrium state, so that acute outbreaks of phage infections in pools of susceptible cells continue to occur. Our model then predicts the selection of phages that switch infection strategy when half of the available susceptible cells have been infected.

Bacteriophages, or phages for short, are viruses that need to infect bacteria to multiply. Once inside a cell, phages follow one of two strategies. They either start to replicate quickly, killing the host in the process; or they lay dormant, their genetic material slowly duplicating as the bacterium divides. These two strategies are respectively known as a ‘lytic’ or a ‘lysogenic’ infection.

In 2017, scientists discovered that, during infection, some phages produce a signalling molecule that influences the strategy other phages will use. Generally, a high concentration of the signal triggers lysogenic infection, while a low level prompts the lytic type. However, it is still unclear what advantages this communication system brings to the viruses, and how it has evolved.

Here, Doekes et al. used a mathematical model to explore how communication changes as phages infect a population of bacteria, rigorously testing earlier theories. The simulations showed that early in an outbreak, when only a few cells have yet been infected, the signalling molecule levels are low: lytic infections are therefore triggered and the phages quickly multiply, killing their hosts in the process. This is an advantageous strategy since many bacteria are available for the viruses to prey on. Later on, as more phages are being produced and available bacteria become few and far between, the levels of the signalling molecule increase. The viruses then switch to lysogenic infections, which allows them to survive dormant, inside their host.

Doekes et al. also discovered that this communication system only evolves if phages regularly cause large outbreaks in new, uninfected bacterial populations. From there, the model was able to predict that phages switch from lytic to lysogenic infections when about half the available bacteria have been infected.

As antibiotic resistance rises around the globe, phages are increasingly considered as a new way to fight off harmful bacteria. Deciphering the way these viruses communicate could help to understand how they could be harnessed to control the spread of bacteria.

The evolutionary dynamics of viruses: virion release strategies, time delays and fitness minima

Viruses exhibit a diverse array of strategies for infecting host cells and for virion release after replication. Cell exit strategies generally involve either budding from the cell membrane or killing the host cell. The conditions under which either is at a selective advantage is a key question in the evolutionary theory of viruses, with the outcome having potentially important impacts on the course of infection and pathogenicity. Although a plethora of external factors will influence the fitness of either strategy; here, we focus just on the effects of the physical properties of the system. We develop theoretical approaches to assess the effects of the time delays between initial infection and virion release. We show that the length of the delay before apoptosis is an important trait in virus evolutionary dynamics. Our results show that for a fixed time to apoptosis, intermediate delays lead to virus fitness that is lower than short times to apoptosis - leading to an apoptotic strategy - and long times to apoptosis - leading to a budding strategy at the between-cell level. At fitness minima, selection is expected to be disruptive and the potential for adaptive radiation in virus strategies is feasible. Hence, the physical properties of the system are sufficient to explain the existence of both budding and virus-induced apoptosis. The fitness functions presented here provide a formal basis for further work focusing on the evolutionary implications of trade-offs between time delays, intracellular replication and resulting mutation rates.

Making sense of virus size and the tradeoffs shaping viral fitness

Viruses span an impressive size range, with genome length varying a thousandfold and virion volume nearly a millionfold. For cellular organisms the scaling of traits with size is a pervasive influence on ecological processes, but whether size plays a central role in viral ecology is unknown. Here, we focus on viruses of aquatic unicellular organisms, which exhibit the greatest known range of virus size. We outline hypotheses within a quantitative framework, and analyse data where available, to consider how size affects the primary components of viral fitness. We argue that larger viruses have fewer offspring per infection and slower contact rates with host cells, but a larger genome tends to increase infection efficiency, broaden host range, and potentially increase attachment success and decrease decay rate. These countervailing selective pressures may explain why a breadth of sizes exist and even coexist when infecting the same host populations. Oligotrophic ecosystems may be enriched in "giant" viruses, because environments with resource-limited phagotrophs at low concentrations may select for broader host range, better control of host metabolism, lower decay rate and a physical size that mimics bacterial prey. Finally, we describe where further research is needed to understand the ecology and evolution of viral size diversity.

Characterization of the Vaginal DNA Virome in Health and Dysbiosis

Bacterial vaginosis (BV) is characterized by a reduction in Lactobacillus (L.) spp. abundance and increased abundance of facultative anaerobes, such as Gardnerella spp. BV aetiology is not fully understood; however, bacteriophages could play a pivotal role in the perturbation of the vaginal bacterial community. We investigated the vaginal viral community, including bacteriophages and the association to the bacterial community and BV-status. Vaginal samples from 48 patients undergoing IVF treatment for non-female factor infertility were subjected to metagenomic sequencing of purified virus-like particles. The vaginal viral community was characterized and correlated with the BV-status by Nugent score, bacterial community, structure, and the presence of key vaginal bacterial species. The majority of identified vaginal viruses belonged to the class of double-stranded DNA bacteriophages, with eukaryotic viruses constituting 4% of the total reads. Clear links between the viral community composition and BV (q = 0.006, R = 0.26) as well as the presence of L. crispatus (q = 0.001, R = 0.43), L. iners, Gardnerella spp., and Atopobium vaginae were found (q < 0.002, R > 0.15). The eukaryotic viral community also correlated with BV-status (q = 0.018, R = 0.20). In conclusion, the vaginal virome was clearly linked with bacterial community structure and BV-status.

Diversity of Vaginal Microbiome in Pregnancy: Deciphering the Obscurity

Human microbiota plays an indispensable role in physiology, nutrition and most significantly, in imparting immunity. The role of microbiota has remained cryptic for years, until recently meticulous studies revealed the interaction and dynamics of these microbial communities. This diversified state is governed by hormonal, behavioral and physio-chemical changes in the genital tract. Many inclusive studies have revealed "Lactobacillus" to be the most dominant member of vaginal flora in most of the healthy, reproductive age group and pregnant females. A total of five community state types have been described, out of which four are dominated by Lactobacillus while the fifth one by facultative or strict anaerobic species. A variation between species stability and gestational age has also been revealed. Studies have divulged a significant higher stability of vaginal microbiota in early stages of pregnancy and the same increased subsequently. Inter-species and racial variation has shown women belonging to White, Asian, and Caucasian race to harbor more of the anaerobic flora. The vaginal microbiome in pregnancy play a significant role in preterm and spontaneous labor. This Lactobacillus-rich microbiome falls tremendously, becoming more diverse in post-partum period. Apart from these known bacterial communities in human vagina, other microbial communities have also been traced. The major fragment is constituted by vaginal viral virome and very little information exists in relation to vaginal mycobiome. Studies have revealed the abundance of ds DNA viruses in vaginal microbiome, followed by ssDNA, and few unidentified viruses. The eukaryotic viruses detected were very few, with Herpesvirales, and Papillomaviridae being the only pathogenic ones. This flora is transmitted to infants either via maternal gut, vagina or breast milk. Recent studies have given an insight for vaginal microbiome, dissociating the old concept of "healthy" and "diseased." However, more extensive studies are required to study evolution of virome and mycobiome in relation to their association with bacterial communities; to establish and decode full array of vaginal virome under the influence of genotypic and environmental factors, using novel bioinformatic, multi-omic, statistical model, and CRISPR/Cas approaches.

Evolutionary Stabilization of Cooperative Toxin Production through a Bacterium-Plasmid-Phage Interplay

Colicins are toxins produced and released by Enterobacteriaceae to kill competitors in the gut. While group A colicins employ a division of labor strategy to liberate the toxin into the environment via colicin-specific lysis, group B colicin systems lack cognate lysis genes. In Salmonella enterica serovar Typhimurium (S. Tm), the group B colicin Ib (ColIb) is released by temperate phage-mediated bacteriolysis. Phage-mediated ColIb release promotes S. Tm fitness against competing Escherichia coli It remained unclear how prophage-mediated lysis is realized in a clonal population of ColIb producers and if prophages contribute to evolutionary stability of toxin release in S. Tm. Here, we show that prophage-mediated lysis occurs in an S. Tm subpopulation only, thereby introducing phenotypic heterogeneity to the system. We established a mathematical model to study the dynamic interplay of S. Tm, ColIb, and a temperate phage in the presence of a competing species. Using this model, we studied long-term evolution of phage lysis rates in a fluctuating infection scenario. This revealed that phage lysis evolves as bet-hedging strategy that maximizes phage spread, regardless of whether colicin is present or not. We conclude that the ColIb system, lacking its own lysis gene, is making use of the evolutionary stable phage strategy to be released. Prophage lysis genes are highly prevalent in nontyphoidal Salmonella genomes. This suggests that the release of ColIb by temperate phages is widespread. In conclusion, our findings shed new light on the evolution and ecology of group B colicin systems.IMPORTANCE Bacteria are excellent model organisms to study mechanisms of social evolution. The production of public goods, e.g., toxin release by cell lysis in clonal bacterial populations, is a frequently studied example of cooperative behavior. Here, we analyze evolutionary stabilization of toxin release by the enteric pathogen Salmonella The release of colicin Ib (ColIb), which is used by Salmonella to gain an edge against competing microbiota following infection, is coupled to bacterial lysis mediated by temperate phages. Here, we show that phage-dependent lysis and subsequent release of colicin and phage particles occurs only in part of the ColIb-expressing Salmonella population. This phenotypic heterogeneity in lysis, which represents an essential step in the temperate phage life cycle, has evolved as a bet-hedging strategy under fluctuating environments such as the gastrointestinal tract. Our findings suggest that prophages can thereby evolutionarily stabilize costly toxin release in bacterial populations.

When to be temperate: on the fitness benefits of lysis vs. lysogeny

Bacterial viruses, that is ‘bacteriophage’ or ‘phage’, can infect and lyse their bacterial hosts, releasing new viral progeny. In addition to the lytic pathway, certain bacteriophage (i.e. ‘temperate’ bacteriophage) can also initiate lysogeny, a latent mode of infection in which the viral genome is integrated into and replicated with the bacterial chromosome. Subsequently, the integrated viral genome, that is the ‘prophage’, can induce and restart the lytic pathway. Here, we explore the relationship among infection mode, ecological context, and viral fitness, in essence asking: when should viruses be temperate? To do so, we use network loop analysis to quantify fitness in terms of network paths through the life history of an infectious pathogen that start and end with infected cells. This analysis reveals that temperate strategies, particularly those with direct benefits to cellular fitness, should be favored at low host abundances. This finding applies to a spectrum of mechanistic models of phage–bacteria dynamics spanning both explicit and implicit representations of intra-cellular infection dynamics. However, the same analysis reveals that temperate strategies, in and of themselves, do not provide an advantage when infection imposes a cost to cellular fitness. Hence, we use evolutionary invasion analysis to explore when temperate phage can invade microbial communities with circulating lytic phage. We find that lytic phage can drive down niche competition amongst microbial cells, facilitating the subsequent invasion of latent strategies that increase cellular resistance and/or immunity to infection by lytic viruses—notably this finding holds even when the prophage comes at a direct fitness cost to cellular reproduction. Altogether, our analysis identifies broad ecological conditions that favor latency and provide a principled framework for exploring the impacts of ecological context on both the short- and long-term benefits of being temperate.

Coming-of-Age Characterization of Soil Viruses: A User’s Guide to Virus Isolation, Detection within Metagenomes, and Viromics

The study of soil viruses, though not new, has languished relative to the study of marine viruses. This is particularly due to challenges associated with separating virions from harboring soils. Generally, three approaches to analyzing soil viruses have been employed: (1) Isolation, to characterize virus genotypes and phenotypes, the primary method used prior to the start of the 21st century. (2) Metagenomics, which has revealed a vast diversity of viruses while also allowing insights into viral community ecology, although with limitations due to DNA from cellular organisms obscuring viral DNA. (3) Viromics (targeted metagenomics of virus-like-particles), which has provided a more focused development of ‘virus-sequence-to-ecology’ pipelines, a result of separation of presumptive virions from cellular organisms prior to DNA extraction. This separation permits greater sequencing emphasis on virus DNA and thereby more targeted molecular and ecological characterization of viruses. Employing viromics to characterize soil systems presents new challenges, however. Ones that only recently are being addressed. Here we provide a guide to implementing these three approaches to studying environmental viruses, highlighting benefits, difficulties, and potential contamination, all toward fostering greater focus on viruses in the study of soil ecology.

The enemy from within: a prophage of Roseburia intestinalis systematically turns lytic in the mouse gut, driving bacterial adaptation by CRISPR spacer acquisition

Despite an overall temporal stability in time of the human gut microbiota at the phylum level, strong variations in species abundance have been observed. We are far from a clear understanding of what promotes or disrupts the stability of microbiome communities. Environmental factors, like food or antibiotic use, modify the gut microbiota composition, but their overall impacts remain relatively low. Phages, the viruses that infect bacteria, might constitute important factors explaining temporal variations in species abundance. Gut bacteria harbour numerous prophages, or dormant viruses, which can evolve to become ultravirulent phage mutants, potentially leading to important bacterial death. Whether such phenomenon occurs in the mammal's microbiota has been largely unexplored. Here we studied temperate phage-bacteria coevolution in gnotoxenic mice colonised with Roseburia intestinalis, a dominant symbiont of the human gut microbiota, and Escherichia coli, a sub-dominant member of the same microbiota. We show that R. intestinalis L1-82 harbours two active prophages, Jekyll and Shimadzu. We observed the systematic evolution in mice of ultravirulent Shimadzu phage mutants, which led to a collapse of R. intestinalis population. In a second step, phage infection drove the fast counter-evolution of host phage resistance mainly through phage-derived spacer acquisition in a clustered regularly interspaced short palindromic repeats array. Alternatively, phage resistance was conferred by a prophage originating from an ultravirulent phage with a restored ability to lysogenize. Our results demonstrate that prophages are a potential source of ultravirulent phages that can successfully infect most of the susceptible bacteria. This suggests that prophages can play important roles in the short-term temporal variations observed in the composition of the gut microbiota.

How pirate phage interferes with helper phage: Comparison of the two distinct strategies

Pirate phages use the structural proteins encoded by unrelated helper phages to propagate. The best-studied example is the pirate P4 and helper P2 of coliphages, and it has been known that the Staphylococcus aureus pathogenicity islands (SaPIs) that can encode virulence factors act as pirate phages, too. When alone in the host, the pirate phages act as a prophage, but when the helper phage gene is also in the same host cell, the pirate phage has ability to exploit the helper phages structural proteins to produce pirate phage particles and spread, interfering with the helper phage production. The known helper phages in these systems are temperate phages. Interestingly, the interference of the pirate phage to the helper phage occurs in a different manner between the SaPI-helper system and the P4-P2 system. SaPIs cannot lyse a helper lysogen upon infection, while when a helper phage lyse a SaPI lysogen, most of the phage particles produced are the SaPI particles. On the contrary, in the P4-P2 system, a pirate phage P4 can lyse a helper P2 lysogen to produce mostly the P4 particles, while when P2 phage lyses a P4 lysogen, most of the produced phages are the P2 particles. Here, the consequences of these different strategies in the pirate and helper phage spreading among uninfected host is analyzed by using mathematical models. It is found that SaPI's strategy interferes with the helper phage spreading significantly more than the P4's strategy, because SaPI interferes with the helper phage's main reproduction step, while P4 interferes only by forcing the helper lysogens to lyse. However, the interference is found to be weaker in the spatially structured environment than in the well-mixed environment. This is because, in the spatial setting, the system tends to self-organize so that the helper phages take over the front of propagation due to the need of helper phage for the pirate phage spreading.

Optimality of the spontaneous prophage induction rate

Lysogens are bacterial cells that have survived after genomically incorporating the DNA of temperate bacteriophages infecting them. If an infection results in lysogeny, the lysogen continues to grow and divide normally, seemingly unaffected by the integrated viral genome known as a prophage. However, the prophage can still have an impact on the host's phenotype and overall fitness in certain environments. Additionally, the prophage within the lysogen can activate the lytic pathway via spontaneous prophage induction (SPI), killing the lysogen and releasing new progeny phages. These new phages can then lyse or lysogenize other susceptible nonlysogens, thereby impacting the competition between lysogens and nonlysogens. In a scenario with differing growth rates, it is not clear whether SPI would be beneficial or detrimental to the lysogens since it kills the host cell but also attacks nonlysogenic competitors, either lysing or lysogenizing them. Here we study the evolutionary dynamics of a mixture of lysogens and nonlysogens and derive general conditions on SPI rates for lysogens to displace nonlysogens. We show that there exists an optimal SPI rate for bacteriophage λ and explain why it is so low. We also investigate the impact of stochasticity and conclude that even at low cell numbers SPI can still provide an advantage to the lysogens. These results corroborate recent experimental studies showing that lower SPI rates are advantageous for phage-phage competition, and establish theoretical bounds on the SPI rate in terms of ecological and environmental variables associated with lysogens having a competitive advantage over their nonlysogenic counterparts.

Evolutionary Stability of Salmonella Competition with the Gut Microbiota: How the Environment Fosters Heterogeneity in Exploitative and Interference Competition

Following ingestion, gastrointestinal pathogens compete against the gastrointestinal microbiota and overcome host immune defenses in order to cause infections. Besides employing direct killing mechanisms, the commensal microbiota occupies metabolic niches to outcompete invading pathogens. Salmonella enterica serovar Typhimurium (S. Typhimurium) uses several strategies to successfully colonize the gut and establish infection, of which an increasing number is based on phenotypic heterogeneity within the S. Typhimurium population. The utilization of myo-inositol (MI) and the production of colicin confer a selective advantage over the microbiota in terms of exploitative and interference competition, respectively. In this review, we summarize the genetic basis underlying bistability of MI catabolism and colicin production. As demonstrated by single-cell analyses, a stochastic switch in the expression of the genes responsible for colicin production and MI degradation constitutes the heterogeneity of the two phenotypes. Both genetic systems are tightly regulated to avoid their expression under non-appropriate conditions and possible detrimental effects on bacterial fitness. Moreover, evolutionary mechanisms underlying formation and stability of these phenotypes in S. Typhimurium are discussed. We propose that both MI catabolism and colicin production create a bet-hedging strategy, which provides an adaptive benefit for S. Typhimurium in the fluctuating environment of the mammalian gut.

Viral invasion fitness across a continuum from lysis to latency

The prevailing paradigm in ecological studies of viruses and their microbial hosts is that the reproductive success of viruses depends on the proliferation of the 'predator', that is, the virus particle. Yet, viruses are obligate intracellular parasites, and the virus genome-the actual unit of selection-can persist and proliferate from one cell generation to the next without lysis or the production of new virus particles. Here, we propose a theoretical framework to quantify the invasion fitness of viruses using an epidemiological cell-centric metric that focuses on the proliferation of viral genomes inside cells instead of virus particles outside cells. This cell-centric metric enables direct comparison of viral strategies characterized by obligate killing of hosts (e.g. via lysis), persistence of viral genomes inside hosts (e.g. via lysogeny), and strategies along a continuum between these extremes (e.g. via chronic infections). As a result, we can identify environmental drivers, life history traits, and key feedbacks that govern variation in viral propagation in nonlinear population models. For example, we identify threshold conditions given relatively low densities of susceptible cells and relatively high growth rates of infected cells in which lysogenic and other chronic strategies have higher potential viral reproduction than lytic strategies. Altogether, the theoretical framework helps unify the ongoing study of eco-evolutionary drivers of viral strategies in natural environments.

Evolutionary stability of the lysis-lysogeny decision: Why be virulent?

Lytic viruses infect and kill host cells, producing a large number of viral copies. Temperate viruses, in contrast, are able to integrate viral genetic material into the host cell DNA, leaving a viable host cell. The evolutionary advantage of this strategy, lysogeny, has been demonstrated in complex environments that include spatial structure, oscillating population dynamics, or periodic environmental collapse. Here, we examine the evolutionary stability of the lysis-lysogeny decision, that is, we predict the long-term outcome of the evolution of lysogeny rates. We demonstrate that viruses with high rates of lysogeny are stable against invasion by more virulent viral strains even in simple environments, as long as the pool of susceptible hosts is not unlimited. This mirrors well-known results in both r-K selection theory and virulence evolution: although virulent viruses have a faster potential growth rate, temperate strains are able to maintain positive growth on a lower density of the limiting resource, susceptible hosts. We then outline scenarios in which the rate of lysogeny is predicted to evolve either toward full lysogeny or full lysis. Finally, we demonstrate conditions under which intermediate rates of lysogeny, as observed in temperate viruses in nature, can be sustained long-term. In general, intermediate lysogeny rates persist when the coupling between susceptible host density and virus density is relaxed.

Optimized bacteria are environmental prediction engines

Experimentalists observe phenotypic variability even in isogenic bacteria populations. We explore the hypothesis that in fluctuating environments this variability is tuned to maximize a bacterium's expected log-growth rate, potentially aided by epigenetic (all inheritable nongenetic) markers that store information about past environments. Crucially, we assume a time delay between sensing and action, so that a past epigenetic marker is used to generate the present phenotypic variability. We show that, in a complex, memoryful environment, the maximal expected log-growth rate is linear in the instantaneous predictive information-the mutual information between a bacterium's epigenetic markers and future environmental states. Hence, under resource constraints, optimal epigenetic markers are causal states-the minimal sufficient statistics for prediction-or lossy approximations thereof. We propose new theoretical investigations into and new experiments on bacteria phenotypic bet-hedging in fluctuating complex environments.

Regulated Stochasticity in a Bacterial Signaling Network Permits Tolerance to a Rapid Environmental Change

Microbial populations can maximize fitness in dynamic environments through bet hedging, a process wherein a subpopulation assumes a phenotype not optimally adapted to the present environment but well adapted to an environment likely to be encountered. Here, we show that oxygen induces fluctuating expression of the trimethylamine oxide (TMAO) respiratory system of Escherichia coli, diversifying the cell population and enabling a bet-hedging strategy that permits growth following oxygen loss. This regulation by oxygen affects the variance in gene expression but leaves the mean unchanged. We show that the oxygen-sensitive transcription factor IscR is the key regulator of variability. Oxygen causes IscR to repress expression of a TMAO-responsive signaling system, allowing stochastic effects to have a strong effect on the output of the system and resulting in heterogeneous expression of the TMAO reduction machinery. This work reveals a mechanism through which cells regulate molecular noise to enhance fitness.

A growing microcolony can survive and support persistent propagation of virulent phages

Bacteria form colonies and secrete extracellular polymeric substances that surround the individual cells. These spatial structures are often associated with collaboration and quorum sensing between the bacteria. Here we investigate the mutual protection provided by spherical growth of a monoclonal colony during exposure to phages that proliferate on its surface. As a proof of concept we exposed growing colonies of Escherichia coli to a virulent mutant of phage P1. When the colony consists of less than [Formula: see text]50,000 members it is eliminated, while larger initial colonies allow long-term survival of both phage-resistant mutants and, importantly, colonies of mostly phage-sensitive members. A mathematical model predicts that colonies formed solely by phage-sensitive bacteria can survive because the growth of bacteria throughout the colony exceeds the killing of bacteria on the surface and pinpoints how the critical colony size depends on key parameters in the phage infection cycle.

Viruses and mobile elements as drivers of evolutionary transitions

The history of life is punctuated by evolutionary transitions which engender emergence of new levels of biological organization that involves selection acting at increasingly complex ensembles of biological entities. Major evolutionary transitions include the origin of prokaryotic and then eukaryotic cells, multicellular organisms and eusocial animals. All or nearly all cellular life forms are hosts to diverse selfish genetic elements with various levels of autonomy including plasmids, transposons and viruses. I present evidence that, at least up to and including the origin of multicellularity, evolutionary transitions are driven by the coevolution of hosts with these genetic parasites along with sharing of ‘public goods’. Selfish elements drive evolutionary transitions at two distinct levels. First, mathematical modelling of evolutionary processes, such as evolution of primitive replicator populations or unicellular organisms, indicates that only increasing organizational complexity, e.g. emergence of multicellular aggregates, can prevent the collapse of the host–parasite system under the pressure of parasites. Second, comparative genomic analysis reveals numerous cases of recruitment of genes with essential functions in cellular life forms, including those that enable evolutionary transitions.

This article is part of the themed issue ‘The major synthetic evolutionary transitions’.

In silico Evolution of Lysis-Lysogeny Strategies Reproduces Observed Lysogeny Propensities in Temperate Bacteriophages

Bacteriophages are the most abundant organisms on the planet and both lytic and temperate phages play key roles as shapers of ecosystems and drivers of bacterial evolution. Temperate phages can choose between (i) lysis: exploiting their bacterial hosts by producing multiple phage particles and releasing them by lysing the host cell, and (ii) lysogeny: establishing a potentially mutually beneficial relationship with the host by integrating their chromosome into the host cell's genome. Temperate phages exhibit lysogeny propensities in the curiously narrow range of 5–15%. For some temperate phages, the propensity is further regulated by the multiplicity of infection, such that single infections go predominantly lytic while multiple infections go predominantly lysogenic. We ask whether these observations can be explained by selection pressures in environments where multiple phage variants compete for the same host. Our models of pairwise competition, between phage variants that differ only in their propensity to lysogenize, predict the optimal lysogeny propensity to fall within the experimentally observed range. This prediction is robust to large variation in parameters such as the phage infection rate, burst size, decision rate, as well as bacterial growth rate, and initial phage to bacteria ratio. When we compete phage variants whose lysogeny strategies are allowed to depend upon multiplicity of infection, we find that the optimal strategy is one which switches from full lysis for single infections to full lysogeny for multiple infections. Previous attempts to explain lysogeny propensity have argued for bet-hedging that optimizes the response to fluctuating environmental conditions. Our results suggest that there is an additional selection pressure for lysogeny propensity within phage populations infecting a bacterial host, independent of environmental conditions.

Models of life: epigenetics, diversity and cycles

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

Inevitability of Genetic Parasites

Almost all cellular life forms are hosts to diverse genetic parasites with various levels of autonomy including plasmids, transposons and viruses. Theoretical modeling of the evolution of primordial replicators indicates that parasites (cheaters) necessarily evolve in such systems and can be kept at bay primarily via compartmentalization. Given the (near) ubiquity, abundance and diversity of genetic parasites, the question becomes pertinent: are such parasites intrinsic to life? At least in prokaryotes, the persistence of parasites is linked to the rate of horizontal gene transfer (HGT). We mathematically derive the threshold value of the minimal transfer rate required for selfish element persistence, depending on the element duplication and loss rates as well as the cost to the host. Estimation of the characteristic gene duplication, loss and transfer rates for transposons, plasmids and virus-related elements in multiple groups of diverse bacteria and archaea indicates that most of these rates are compatible with the long term persistence of parasites. Notably, a small but non-zero rate of HGT is also required for the persistence of non-parasitic genes. We hypothesize that cells cannot tune their horizontal transfer rates to be below the threshold required for parasite persistence without experiencing highly detrimental side-effects. As a lower boundary to the minimum DNA transfer rate that a cell can withstand, we consider the process of genome degradation and mutational meltdown of populations through Muller's ratchet. A numerical assessment of this hypothesis suggests that microbial populations cannot purge parasites while escaping Muller's ratchet. Thus, genetic parasites appear to be virtually inevitable in cellular organisms.

Population Dynamics of Phage and Bacteria in Spatially Structured Habitats Using Phage λ and Escherichia coli

Bacteria living in physically structured habitats are exposed heterogeneously to both resources and different types of phages. While there have been numerous experimental approaches to examine spatially distributed bacteria exposed to phages, there is little theory to guide the design of these experiments, interpret their results, or expand the inferences drawn to a broader ecological and evolutionary context. Plaque formation provides a window into understanding phage-bacterium interactions in physically structured populations, including surfaces, semisolids, and biofilms. We develop models to address the plaque dynamics for a temperate phage and its virulent mutants. The models are compared with phage λ-Escherichia coli system to quantify their applicability. We found that temperate phages gave an increasing number of gradually smaller colonies as the distance increased from the plaque center. For low-lysogen frequency this resulted in plaques with most of the visible colonies at an intermediate distance between the center and periphery. Using spot inoculation, where phages in excess of bacteria were inoculated in a circular area, we measured the frequency and spatial distribution of lysogens. The spot morphology of cII-negative (cII(-)) and cIII(-) mutants of phage λ displays concentric rings of high-density lysogenic colonies. The simplest of these ring morphologies was reproduced by including multiplicity of infection (MOI) sensitivity in lysis-lysogeny decisions, but its failure to explain the occasional observation of multiple rings in cIII(-) mutant phages highlights unknown features of this phage. Our findings demonstrated advantages of temperate phages over virulent phages in exploiting limited resources in spatially distributed microbial populations.

IMPORTANCE

Phages are the most abundant organisms on earth, and yet little is known about how phages and bacterial hosts are influencing each other in density and evolution. Phages can be either virulent or temperate, a difference that is highlighted when a spatially structured bacterial population is infected. Phage λ is a temperate phage, with a capacity for dormancy that can be modified by single gene knockouts. The stochastic bias in the lysis-lysogeny decision's probability is reflected in plaque morphologies on bacterial lawns. We present a model for plaque morphology of both virulent and temperate phages, taking into account the underlying survival of bacterial microcolonies. It reproduces known plaque morphologies and speaks to advantages of temperate phages in a spatially structured environment.

Carriage of λ Latent Virus Is Costly for Its Bacterial Host due to Frequent Reactivation in Monoxenic Mouse Intestine

Temperate phages, the bacterial viruses able to enter in a dormant prophage state in bacterial genomes, are present in the majority of bacterial strains for which the genome sequence is available. Although these prophages are generally considered to increase their hosts’ fitness by bringing beneficial genes, studies demonstrating such effects in ecologically relevant environments are relatively limited to few bacterial species. Here, we investigated the impact of prophage carriage in the gastrointestinal tract of monoxenic mice. Combined with mathematical modelling, these experimental results provided a quantitative estimation of key parameters governing phage-bacteria interactions within this model ecosystem. We used wild-type and mutant strains of the best known host/phage pair, Escherichia coli and phage λ. Unexpectedly, λ prophage caused a significant fitness cost for its carrier, due to an induction rate 50-fold higher than in vitro, with 1 to 2% of the prophage being induced. However, when prophage carriers were in competition with isogenic phage susceptible bacteria, the prophage indirectly benefited its carrier by killing competitors: infection of susceptible bacteria led to phage lytic development in about 80% of cases. The remaining infected bacteria were lysogenized, resulting overall in the rapid lysogenization of the susceptible lineage. Moreover, our setup enabled to demonstrate that rare events of phage gene capture by homologous recombination occurred in the intestine of monoxenic mice. To our knowledge, this study constitutes the first quantitative characterization of temperate phage-bacteria interactions in a simplified gut environment. The high prophage induction rate detected reveals DNA damage-mediated SOS response in monoxenic mouse intestine. We propose that the mammalian gut, the most densely populated bacterial ecosystem on earth, might foster bacterial evolution through high temperate phage activity.

Dormant bacterial viruses, or prophages, are found in the genomes of almost all bacteria, but their impact on bacterial host fitness is largely unknown. Through experiments in mice, supported by a mathematical model, we quantified the activity of Escherichia coli prophage λ in monoxenic mouse gut, as well as its impact on its carrier bacteria. λ carriage negatively impacted its hosts due to frequent reactivation, but indirectly benefited its host by killing susceptible bacterial competitors. The high prophage activity unraveled in this study reflects a constant rate of SOS response, resulting from DNA damage in monoxenic mouse intestine. Our results should motivate researchers to take the presence of prophages into account when studying the action of specific bacteria in the gastrointestinal tract of mammals.

Diversity Waves in Collapse-Driven Population Dynamics

Populations of species in ecosystems are often constrained by availability of resources within their environment. In effect this means that a growth of one population, needs to be balanced by comparable reduction in populations of others. In neutral models of biodiversity all populations are assumed to change incrementally due to stochastic births and deaths of individuals. Here we propose and model another redistribution mechanism driven by abrupt and severe reduction in size of the population of a single species freeing up resources for the remaining ones. This mechanism may be relevant e.g. for communities of bacteria, with strain-specific collapses caused e.g. by invading bacteriophages, or for other ecosystems where infectious diseases play an important role. The emergent dynamics of our system is characterized by cyclic ''diversity waves'' triggered by collapses of globally dominating populations. The population diversity peaks at the beginning of each wave and exponentially decreases afterwards. Species abundances have bimodal time-aggregated distribution with the lower peak formed by populations of recently collapsed or newly introduced species while the upper peak--species that has not yet collapsed in the current wave. In most waves both upper and lower peaks are composed of several smaller peaks. This self-organized hierarchical peak structure has a long-term memory transmitted across several waves. It gives rise to a scale-free tail of the time-aggregated population distribution with a universal exponent of 1.7. We show that diversity wave dynamics is robust with respect to variations in the rules of our model such as diffusion between multiple environments, species-specific growth and extinction rates, and bet-hedging strategies.