Phage and bacteria support mutual diversity in a narrowing staircase of coexistence

Diverse phage communities are maintained stably on a clonal bacterial host

Bacteriophages are the most abundant and phylogenetically diverse biological entities on Earth, yet the ecological mechanisms that sustain this extraordinary diversity remain unclear. In this study, we discovered that phage diversity consistently outstripped the diversity of their bacterial hosts under simple experimental conditions. We assembled and passaged dozens of diverse phage communities on a single, nonevolving strain of Escherichia coli until the phage communities reached equilibrium. In all cases, we found that two or more phage species coexisted stably, despite competition for a single, clonal host population. Phage coexistence was supported through host phenotypic heterogeneity, whereby bacterial cells adopting different growth phenotypes served as niches for different phage species. Our experiments reveal that a rich community ecology of bacteriophages can emerge on a single bacterial host.

Phage predation accelerates the spread of plasmid-encoded antibiotic resistance

Phage predation is generally assumed to reduce microbial proliferation while not contributing to the spread of antibiotic resistance. However, this assumption does not consider the effect of phage predation on the spatial organization of different microbial populations. Here, we show that phage predation can increase the spread of plasmid-encoded antibiotic resistance during surface-associated microbial growth by reshaping spatial organization. Using two strains of the bacterium Escherichia coli, we demonstrate that phage predation slows the spatial segregation of the strains during growth. This increases the number of cell-cell contacts and the extent of conjugation-mediated plasmid transfer between them. The underlying mechanism is that phage predation shifts the location of fastest growth from the biomass periphery to the interior where cells are densely packed and aligned closer to parallel with each other. This creates straighter interfaces between the strains that are less likely to merge together during growth, consequently slowing the spatial segregation of the strains and enhancing plasmid transfer between them. Our results have implications for the design and application of phage therapy and reveal a mechanism for how microbial functions that are deleterious to human and environmental health can proliferate in the absence of positive selection.

Stochastic microbial dispersal drives local extinction and global diversity

Airborne dispersal of microorganisms is a ubiquitous migration mechanism, allowing otherwise independent microbial habitats to interact via biomass exchange. Here, we study the ecological implications of such advective transport using a simple spatial model for bacteria-phage interactions: the population dynamics at each habitat are described by classical Lotka-Volterra equations; however, species populations are taken as integer, that is, a discrete, positive extinction threshold exists. Spatially, species can spread from habitat to habitat by stochastic airborne dispersal. In any given habitat, the spatial biomass exchange causes incessant population density oscillations, which, as a consequence, occasionally drive species to extinction. The balance between local extinction events and dispersal-induced migration allows species to persist globally, even though diversity would be depleted by competitive exclusion, locally. The disruptive effect of biomass dispersal thus acts to increase microbial diversity, allowing system-scale coexistence of multiple species that would not coexist locally.

The dynamics of phage predation on a microcolony

Phage predation is an important factor for controlling the bacterial biomass. At face value, dense microbial habitats are expected to be vulnerable to phage epidemics due to the abundance of fresh hosts immediately next to any infected bacteria. Despite this, the bacterial microcolony is a common habitat for bacteria in nature. Here, we experimentally quantify the fate of microcolonies of Escherichia coli exposed to virulent phage T4. It has been proposed that the outer bacterial layers of the colony will shield the inner layers from the phage invasion and thereby constrain the phage to the colony's surface. We develop a dynamical model that incorporates this shielding mechanism and fit the results with experimental measurements to extract important phage-bacteria interaction parameters. The analysis suggests that, while the shielding mechanism delays phage attack, T4 phage are able to diffuse so deep into the dense bacterial environment that colony-level survival of the bacterial community is challenged.

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.

Plausible pathway for a host-parasite molecular replication network to increase its complexity through Darwinian evolution

How the complexity of primitive self-replication molecules develops through Darwinian evolution remains a mystery with regards to the origin of life. Theoretical studies have proposed that coevolution with parasitic replicators increases network complexity by inducing inter-dependent replication. Particularly, Takeuchi and Hogeweg proposed a complexification process of replicator networks by successive appearance of a parasitic replicator followed by the addition of a new host replicator that is resistant to the parasitic replicator. However, the feasibility of such complexification with biologically relevant molecules is still unknown owing to the lack of an experimental model. Here, we investigated the plausible complexification pathway of host-parasite replicators using both an experimental host-parasite RNA replication system and a theoretical model based on the experimental system. We first analyzed the parameter space that allows for sustainable replication in various replication networks ranging from a single molecule to three-member networks using computer simulation. The analysis shows that the most plausible complexification pathway from a single host replicator is the addition of a parasitic replicator, followed by the addition of a new host replicator that is resistant to the parasite, consistent with the previous study by Takeuchi and Hogeweg. We also provide evidence that the pathway actually occurred in our previous evolutionary experiment. These results provide experimental evidence that a population of a single replicator spontaneously evolves into multi-replicator networks through coevolution with parasitic replicators.

From kill the winner to eliminate the winner in open phage-bacteria systems

Phages and bacteria manage to coexist and sustain ecosystems with a high diversity of strains, despite limited resources and heavy predation. This diversity can be explained by the “kill the winner” model where virulent phages predominantly prey on fast-growing bacteria and thereby suppress the competitive exclusion of slower-growing bacteria. Here we computationally investigate the robustness of these systems against invasions, where new phages or bacteria may interact with more than one of the resident strains. The resulting interaction networks were found to self-organize into a network with strongly interacting specialized predator-prey pairs, resembling that of the “kill the winner” model. Furthermore, the “kill the winner” dynamics is enforced with the occasional elimination of even the fastest-growing bacteria strains due to a phage infecting the fast and slow growers. The frequency of slower-growing strains was increased with the introduction of even a few non-diagonal interactions. Hence, phages capable of infecting multiple hosts play significant roles both in the evolution of the ecosystem by eliminating the winner and in supporting diversity by allowing slow growers to coexist with faster growers.

We demonstrate that in an open system of phages and bacteria with very limited resources, a bacterial strain that has a high growth rate can still be outcompeted by a slower-growing strain if they have a common phage. The impact of this on ecosystem structure is significant as soon as there is a small probability to have a common phage among bacterial strains. Furthermore, by analysing the structure of the interaction network we show that it self-organizes into a network with strongly interacting specialized predator-prey pairs, in order to reduce phages competition. Nevertheless, the presence of the remaining links is very important for the network dynamics since even a few of them significantly enhance the frequency of slower-growing strains.

Fate of enteric viruses during leafy greens (romaine lettuce) production using treated municipal wastewater and AP205 bacteriophage as a surrogate

Water reuse programs are being explored to close the gap between supply and demand for irrigation in agriculture. However, these sources could contain hazardous microbial contaminants, and pose risks to public health. This study aimed to grow and irrigate romaine lettuce with inoculated wastewater effluent to track AP205 bacteriophage prevalence through cultivation and post-harvest storage. AP205 is a bacteriophage and was used as a surrogate for enteric viruses. Low and high dosages (mean ± standard deviation) of AP205 at 4.8 ± 0.4 log PFU/mL and 6.6 ± 0.2 log PFU/mL; respectively, were prepared to examine viral load influence on contamination levels. Foliage, leachate, and soil contamination levels were directly related to AP205 concentrations in the effluent. AP205 concentrations increased throughout cultivation for foliage and leachate, suggesting bacteriophage accumulation. During post-harvest storage (14 day at 4 °C), there was a significant decrease in AP205 concentration on the foliage. Results show that wastewater effluents usage for leafy greens cultivation can pose risks to humans and additional steps are required to safely apply wastewater effluents to soils and crops.

Diverse Viruses Carrying Genes for Microbial Extremotolerance in the Atacama Desert Hyperarid Soil

Viruses play an essential role in shaping microbial community structures and serve as reservoirs for genetic diversity in many ecosystems. In hyperarid desert environments, where life itself becomes scarce and loses diversity, the interactions between viruses and host populations have remained elusive. Here, we resolved host-virus interactions in the soil metagenomes of the Atacama Desert hyperarid core, one of the harshest terrestrial environments on Earth. We show evidence of diverse viruses infecting a wide range of hosts found in sites up to 205 km apart. Viral genomes carried putative extremotolerance features (i.e., spore formation proteins) and auxiliary metabolic genes, indicating that viruses could mediate the spread of microbial resilience against environmental stress across the desert. We propose a mutualistic model of host-virus interactions in the hyperarid core where viruses seek protection in microbial cells as lysogens or pseudolysogens, while viral extremotolerance genes aid survival of their hosts. Our results suggest that the host-virus interactions in the Atacama Desert soils are dynamic and complex, shaping uniquely adapted microbiomes in this highly selective and hostile environment.IMPORTANCE Deserts are one of the largest and rapidly expanding terrestrial ecosystems characterized by low biodiversity and biomass. The hyperarid core of the Atacama Desert, previously thought to be devoid of life, is one of the harshest environments, supporting only scant biomass of highly adapted microbes. While there is growing evidence that viruses play essential roles in shaping the diversity and structure of nearly every ecosystem, very little is known about the role of viruses in desert soils, especially where viral contact with viable hosts is significantly reduced. Our results demonstrate that diverse viruses are widely dispersed across the desert, potentially spreading key stress resilience and metabolic genes to ensure host survival. The desertification accelerated by climate change expands both the ecosystem cover and the ecological significance of the desert virome. This study sheds light on the complex virus-host interplay that shapes the unique microbiome in desert soils.

High genomic diversity of novel phages infecting the plant pathogen Ralstonia solanacearum, isolated in Mauritius and Reunion islands

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is among the most important plant diseases worldwide, severely affecting a high number of crops and ornamental plants in tropical regions. Only a limited number of phages infecting R. solanacearum have been isolated over the years, despite the importance of this bacterium and the associated plant disease. The antibacterial effect or morphological traits of these R. solanacearum viruses have been well studied, but not their genomic features, which need deeper consideration. This study reports the full genome of 23 new phages infecting RSSC isolated from agricultural samples collected in Mauritius and Reunion islands, particularly affected by this plant bacterial pathogen and considered biodiversity hotspots in the Southwest Indian Ocean. The complete genomic information and phylogenetic classification is provided, revealing high genetic diversity between them and weak similarities with previous related phages. The results support our proposal of 13 new species and seven new genera of R. solanacearum phages. Our findings highlight the wide prevalence of phages of RSSC in infected agricultural settings and the underlying genetic diversity. Discoveries of this kind lead more insight into the diversity of phages in general and to optimizing their use as biocontrol agents of bacterial diseases of plants in agriculture.

Bacteriophages presence in nature and their role in the natural selection of bacterial populations

Phages are the obligate parasite of bacteria and have complex interactions with their hosts. Phages can live in, modify, and shape bacterial communities by bringing about changes in their abundance, diversity, physiology, and virulence. In addition, phages mediate lateral gene transfer, modify host metabolism and reallocate bacterially-derived biochemical compounds through cell lysis, thus playing an important role in ecosystem. Phages coexist and coevolve with bacteria and have developed several antidefense mechanisms in response to bacterial defense strategies against them. Phages owe their existence to their bacterial hosts, therefore they bring about alterations in their host genomes by transferring resistance genes and genes encoding toxins in order to improve the fitness of the hosts. Application of phages in biotechnology, environment, agriculture and medicines demands a deep insight into the myriad of phage-bacteria interactions. However, to understand their complex interactions, we need to know how unique phages are to their bacterial hosts and how they exert a selective pressure on the microbial communities in nature. Consequently, the present review focuses on phage biology with respect to natural selection of bacterial populations.

Defence versus growth in a hostile world: lessons from phage and bacteria

Bacterial communities are often highly diverse with several closely related species (or strains) coexisting together. These bacteria compete for resources and the competitive exclusion principle predicts that all but the fastest-growing bacteria will go extinct. When exposed to phage, it is predicted that bacterial strains with restriction-modification (RM) systems can circumvent the competitive exclusion principle and reach diversity of the order of the phage burst size. We show that with a trade-off between bacterial growth rates and the strength of their RM systems, the diversity of such an ecosystem can further increase several fold beyond the burst size limit. Moreover, we find that the ratio of the growth rate of a bacterial strain to the imperfection of its RM system is an excellent predictor of (i) whether the strain will go extinct or not, and (ii) the biomass of the strain if it survives. In contrast, the growth rate alone is not a determinant of either of these properties. Our work provides a quantitative example of a model ecosystem where the fitness of a species is determined not by growth rate, but by a trade-off between growth and defence against predators.

Formation of phage lysis patterns and implications on co-propagation of phages and motile host bacteria

Coexistence of bacteriophages, or phages, and their host bacteria plays an important role in maintaining the microbial communities. In natural environments with limited nutrients, motile bacteria can actively migrate towards locations of richer resources. Although phages are not motile themselves, they can infect motile bacterial hosts and spread in space via the hosts. Therefore, in a migrating microbial community coexistence of bacteria and phages implies their co-propagation in space. Here, we combine an experimental approach and mathematical modeling to explore how phages and their motile host bacteria coexist and co-propagate. When lytic phages encountered motile host bacteria in our experimental set up, a sector-shaped lysis zone formed. Our mathematical model indicates that local nutrient depletion and the resulting inhibition of proliferation and motility of bacteria and phages are the key to formation of the observed lysis pattern. The model further reveals the straight radial boundaries in the lysis pattern as a telltale sign for coexistence and co-propagation of bacteria and phages. Emergence of such a pattern, albeit insensitive to extrinsic factors, requires a balance between intrinsic biological properties of phages and bacteria, which likely results from coevolution of phages and bacteria.

Sustainability of spatially distributed bacteria-phage systems

Virulent phages can expose their bacterial hosts to devastating epidemics, in principle leading to complete elimination of their hosts. Although experiments indeed confirm a large reduction of susceptible bacteria, there are no reports of complete extinctions. We here address this phenomenon from the perspective of spatial organization of bacteria and how this can influence the final survival of them. By modelling the transient dynamics of bacteria and phages when they are introduced into an environment with finite resources, we quantify how time delayed lysis, the spatial separation of initial bacterial positions, and the self-protection of bacteria growing in spherical colonies favour bacterial survival. Our results suggest that spatial structures on the millimetre and submillimetre scale play an important role in maintaining microbial diversity.

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.

Regime Shifts in a Phage-Bacterium Ecosystem and Strategies for Its Control

The competition between bacteria often involves both nutrients and phage predators and may give rise to abrupt regime shifts between the alternative stable states characterized by different species compositions. While such transitions have been previously studied in the context of competition for nutrients, the case of phage-induced bistability between competing bacterial species has not been considered yet. Here we demonstrate a possibility of regime shifts in well-mixed phage-bacterium ecosystems. In one of the bistable states, the fast-growing bacteria competitively exclude the slow-growing ones by depleting their common nutrient. Conversely, in the second state, the slow-growing bacteria with a large burst size generate such a large phage population that the other species cannot survive. This type of bistability can be realized as the competition between a strain of bacteria protected from phage by abortive infection and another strain with partial resistance to phage. It is often desirable to reliably control the state of microbial ecosystems, yet bistability significantly complicates this task. We discuss successes and limitations of one control strategy in which one adds short pulses to populations of individual species. Our study proposes a new type of phage therapy, where introduction of the phage is supplemented by the addition of a partially resistant host bacteria.IMPORTANCE Phage-microbe communities play an important role in human health as well as natural and industrial environments. Here we show that these communities can assume several alternative species compositions separated by abrupt regime shifts. Our model predicts these regime shifts in the competition between bacterial strains protected by two different phage defense mechanisms: abortive infection/CRISPR and partial resistance. The history dependence caused by regime shifts greatly complicates the task of manipulation and control of a community. We propose and study a successful control strategy via short population pulses aimed at inducing the desired regime shifts. In particular, we predict that a fast-growing pathogen could be eliminated by a combination of its phage and a slower-growing susceptible host.

Prophage LambdaSo uses replication interference to suppress reproduction of coexisting temperate phage MuSo2 in Shewanella oneidensis MR-1

Many bacterial genomes carry multiple prophages that compete with each other, potentially affecting the physiology, fitness, and pathogenicity of their hosts. However, molecular mechanisms of such prophage-prophage conflicts remain poorly understood. The genome of Shewanella oneidensis MR-1, a Gammaproteobacterium residing in aquatic environments and notable for its ability to reduce metal ions, harbours four prophages, two of which (LambdaSo and MuSo2) form infectious virions during biofilm formation. Here, we constructed indicator strains of LambdaSo and MuSo2 by deleting the corresponding prophages from the MR-1 chromosome and investigated their reproduction. Interestingly, the fitness of MuSo2 increased in the absence of LambdaSo, suggesting that prophage LambdaSo repressed MuSo2 reproduction. Partial deletion of LambdaSo from the MR-1 chromosome revealed that gene cluster R of LambdaSo, which was responsible for the switch to the lytic cycle and LambdaSo genome replication initiation, was necessary and sufficient to repress MuSo2. Furthermore, activation of cluster R genes facilitated replication of cluster R-encoding DNA and inhibited host and MuSo2 DNA replication. These findings suggest that LambdaSo represses MuSo2 propagation by inhibiting DNA replication during simultaneous induction. We predict that such a mechanism of inter-prophage interference is more widespread in bacteria than currently appreciated.

Simple models combining competition, defence and resource availability have broad implications in pelagic microbial food webs

In food webs, interactions between competition and defence control the partitioning of limiting resources. As a result, simple models of these interactions contain links between biogeochemistry, diversity, food web structure and ecosystem function. Working at hierarchical levels, these mechanisms also produce self-similarity and therefore suggest how complexity can be generated from repeated application of simple underlying principles. Reviewing theoretical and experimental literature relevant to the marine photic zone, we argue that there is a wide spectrum of phenomena, including single cell activity of prokaryotes, microbial biodiversity at different levels of resolution, ecosystem functioning, regional biogeochemical features and evolution at different timescales; that all can be understood as variations over a common principle, summarised in what has been termed the 'Killing-the-Winner' (KtW) motif. Considering food webs as assemblages of such motifs may thus allow for a more integrated approach to aquatic microbial ecology.

Biogeography and environmental conditions shape bacteriophage-bacteria networks across the human microbiome

Viruses and bacteria are critical components of the human microbiome and play important roles in health and disease. Most previous work has relied on studying bacteria and viruses independently, thereby reducing them to two separate communities. Such approaches are unable to capture how these microbial communities interact, such as through processes that maintain community robustness or allow phage-host populations to co-evolve. We implemented a network-based analytical approach to describe phage-bacteria network diversity throughout the human body. We built these community networks using a machine learning algorithm to predict which phages could infect which bacteria in a given microbiome. Our algorithm was applied to paired viral and bacterial metagenomic sequence sets from three previously published human cohorts. We organized the predicted interactions into networks that allowed us to evaluate phage-bacteria connectedness across the human body. We observed evidence that gut and skin network structures were person-specific and not conserved among cohabitating family members. High-fat diets appeared to be associated with less connected networks. Network structure differed between skin sites, with those exposed to the external environment being less connected and likely more susceptible to network degradation by microbial extinction events. This study quantified and contrasted the diversity of virome-microbiome networks across the human body and illustrated how environmental factors may influence phage-bacteria interactive dynamics. This work provides a baseline for future studies to better understand system perturbations, such as disease states, through ecological networks.

The human microbiome, the collection of microbial communities that colonize the human body, is a crucial component to health and disease. Two major components of the human microbiome are the bacterial and viral communities. These communities have primarily been studied separately using metrics of community composition and diversity. These approaches have failed to capture the complex dynamics of interacting bacteria and phage communities, which frequently share genetic information and work together to maintain ecosystem homestatsis (e.g. kill-the-winner dynamics). Removal of bacteria or phage can disrupt or even collapse those ecosystems. Relationship-based network approaches allow us to capture this interaction information. Using this network-based approach with three independent human cohorts, we were able to present an initial understanding of how phage-bacteria networks differ throughout the human body, so as to provide a baseline for future studies of how and why microbiome networks differ in disease states.

Theory of invasion extinction dynamics in minimal food webs

When food webs are exposed to species invasion, secondary extinction cascades may be set off. Although much work has gone into characterizing the structure of food webs, systematic predictions on their evolutionary dynamics are still scarce. Here we present a theoretical framework that predicts extinctions in terms of an alternating sequence of two basic processes: resource depletion by or competitive exclusion between consumers. We first propose a conceptual invasion extinction model (IEM) involving random fitness coefficients. We bolster this IEM by an analytical, recursive procedure for calculating idealized extinction cascades after any species addition and simulate the long-time evolution. Our procedure describes minimal food webs where each species interacts with only a single resource through the generalized Lotka-Volterra equations. For such food webs ex- tinction cascades are determined uniquely and the system always relaxes to a stable steady state. The dynamics and scale invariant species life time resemble the behavior of the IEM, and correctly predict an upper limit for trophic levels as observed in the field.

Existence and construction of large stable food webs

Ecological diversity is ubiquitous despite the restrictions imposed by competitive exclusion and apparent competition. To explain the observed richness of species in a given habitat, food-web theory has explored nonlinear functional responses, self-interaction, or spatial structure and dispersal-model ingredients that have proven to promote stability and diversity. We return instead here to classical Lotka-Volterra equations, where species-species interaction is characterized by a simple product and spatial restrictions are ignored. We quantify how this idealization imposes constraints on coexistence and diversity for many species. To this end, we introduce the concept of free and controlled species and use this to demonstrate how stable food webs can be constructed by the sequential addition of species. The resulting food webs can reach dozens of species and generally yield nonrandom degree distributions in accordance with the constraints imposed through the assembly process. Our model thus serves as a formal starting point for the study of sustainable interaction patterns between species.

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.

Apparent Competition

Most species have one or more natural enemies, e.g., predators, parasites, pathogens, and herbivores, among others. These species in turn typically attack multiple victim species. This leads to the possibility of indirect interactions among those victims, both positive and negative. The term apparent competition commonly denotes negative indirect interactions between victim species that arise because they share a natural enemy. This indirect interaction, which in principle can be reflected in many facets of the distribution and abundance of individual species and more broadly govern the structure of ecological communities in time and space, pervades many natural ecosystems. It also is a central theme in many applied ecological problems, including the control of agricultural pests, harvesting, the conservation of endangered species, and the dynamics of emerging diseases. At one end of the scale of life, apparent competition characterizes intriguing aspects of dynamics within individual organisms—for example, the immune system is akin in many ways to a predator that can induce negative indirect interactions among different pathogens. At intermediate scales of biological organization, the existence and strength of apparent competition depend upon many contingent details of individual behavior and life history, as well as the community and spatial context within which indirect interactions play out. At the broadest scale of macroecology and macroevolution, apparent competition may play a major, if poorly understood, role in the evolution of species’ geographical ranges and adaptive radiations.

Non-classical phase diagram for virus bacterial coevolution mediated by clustered regularly interspaced short palindromic repeats

CRISPR is a newly discovered prokaryotic immune system. Bacteria and archaea with this system incorporate genetic material from invading viruses into their genomes, providing protection against future infection by similar viruses. The condition for coexistence of prokaryots and viruses is an interesting problem in evolutionary biology. In this work, we show an intriguing phase diagram of the virus extinction probability, which is more complex than that of the classical predator-prey model. As the CRISPR incorporates genetic material, viruses are under pressure to evolve to escape recognition by CRISPR. When bacteria have a small rate of deleting spacers, a new parameter region in which bacteria and viruses can coexist arises, and it leads to a more complex coexistence patten for bacteria and viruses. For example, when the virus mutation rate is low, the virus extinction probability changes non-montonically with the bacterial exposure rate. The virus and bacteria coevolution not only alters the virus extinction probability, but also changes the bacterial population structure. Additionally, we show that recombination is a successful strategy for viruses to escape from CRISPR recognition when viruses have multiple proto-spacers, providing support for a recombination-mediated escape mechanism suggested experimentally. Finally, we suggest that the re-entrant phase diagram, in which phages can progress through three phases of extinction and two phases of abundance at low spacer deletion rates as a function of exposure rate to bacteria, is an experimentally testable phenomenon.

Spatially induced nestedness in a neutral model of phage-bacteria networks

Ecological networks, both displaying mutualistic or antagonistic interactions, seem to share common structural traits: the presence of nestedness and modularity. A variety of model approaches and hypothesis have been formulated concerning the significance and implications of these properties. In phage-bacteria bipartite infection networks, nestedness seems to be the rule in many different contexts. Modeling the coevolution of a diverse virus-host ensemble is a difficult task, given the dimensionality and multi parametric nature of a standard continuous approximation. Here, we take a different approach, by using a neutral, toy model of host-phage interactions on a spatial lattice. Each individual is represented by a bit string (a digital genome) but all strings in each class (i.e. hosts or phages) share the same sets of parameters. A matching allele model of phage-virus recognition rule is enough to generate a complex, diverse ecosystem with heterogeneous patterns of interaction and nestedness, provided that interactions take place under a spatially constrained setting. It is found that nestedness seems to be an emergent property of the co-evolutionary dynamics. Our results indicate that the enhanced diversity resulting from localized interactions strongly promotes the presence of nested infection matrices.

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.

Population cycles and species diversity in dynamic Kill-the-Winner model of microbial ecosystems

Determinants of species diversity in microbial ecosystems remain poorly understood. Bacteriophages are believed to increase the diversity by the virtue of Kill-the-Winner infection bias preventing the fastest growing organism from taking over the community. Phage-bacterial ecosystems are traditionally described in terms of the static equilibrium state of Lotka-Volterra equations in which bacterial growth is exactly balanced by losses due to phage predation. Here we consider a more dynamic scenario in which phage infections give rise to abrupt and severe collapses of bacterial populations whenever they become sufficiently large. As a consequence, each bacterial population in our model follows cyclic dynamics of exponential growth interrupted by sudden declines. The total population of all species fluctuates around the carrying capacity of the environment, making these cycles cryptic. While a subset of the slowest growing species in our model is always driven towards extinction, in general the overall ecosystem diversity remains high. The number of surviving species is inversely proportional to the variation in their growth rates but increases with the frequency and severity of phage-induced collapses. Thus counter-intuitively we predict that microbial communities exposed to more violent perturbations should have higher diversity.

Linking internal and external bacterial community control gives mechanistic framework for pelagic virus-to-bacteria ratios

For more than 25 years, virus-to-bacteria ratios (VBR) have been measured and interpreted as indicators of the importance of viruses in aquatic ecosystems, yet a generally accepted theory for understanding mechanisms controlling VBR is still lacking. Assuming that the denominator (total bacterial abundance) is primarily predator controlled, while viral lysis compensates for host growth rates exceeding this grazing loss, the numerator (viral abundance) reflects activity differences between prokaryotic hosts. VBR is then a ratio between mechanisms generating structure within the bacterial community and interactions between different plankton functional types controlling bacterial community size. We here show how these arguments can be formalized by combining a recently published model for co-evolutionary host-virus interactions, with a previously published "minimum" model for the microbial food web. The result is a framework where viral lysis links bacterial diversity to microbial food web structure and function, creating relationships between different levels of organization that are strongly modified by organism-level properties such as cost of resistance.

High-order species interactions shape ecosystem diversity

Classical theory shows that large communities are destabilized by random interactions among species pairs, creating an upper bound on ecosystem diversity. However, species interactions often occur in high-order combinations, whereby the interaction between two species is modulated by one or more other species. Here, by simulating the dynamics of communities with random interactions, we find that the classical relationship between diversity and stability is inverted for high-order interactions. More specifically, while a community becomes more sensitive to pairwise interactions as its number of species increases, its sensitivity to three-way interactions remains unchanged, and its sensitivity to four-way interactions actually decreases. Therefore, while pairwise interactions lead to sensitivity to the addition of species, four-way interactions lead to sensitivity to species removal, and their combination creates both a lower and an upper bound on the number of species. These findings highlight the importance of high-order species interactions in determining the diversity of natural ecosystems.

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.

Food Web Assembly Rules for Generalized Lotka-Volterra Equations

In food webs, many interacting species coexist despite the restrictions imposed by the competitive exclusion principle and apparent competition. For the generalized Lotka-Volterra equations, sustainable coexistence necessitates nonzero determinant of the interaction matrix. Here we show that this requirement is equivalent to demanding that each species be part of a non-overlapping pairing, which substantially constrains the food web structure. We demonstrate that a stable food web can always be obtained if a non-overlapping pairing exists. If it does not, the matrix rank can be used to quantify the lack of niches, corresponding to unpaired species. For the species richness at each trophic level, we derive the food web assembly rules, which specify sustainable combinations. In neighboring levels, these rules allow the higher level to avert competitive exclusion at the lower, thereby incorporating apparent competition. In agreement with data, the assembly rules predict high species numbers at intermediate levels and thinning at the top and bottom. Using comprehensive food web data, we demonstrate how omnivores or parasites with hosts at multiple trophic levels can loosen the constraints and help obtain coexistence in food webs. Hence, omnivory may be the glue that keeps communities intact even under extinction or ecological release of species.

Human impact currently induces rapid reductions in global biodiversity. Assessing the consequences of such modifications requires that ecological science better understand the conditions under which the species in a community can coexist and when not. Fundamentally, two species can not coexist indefinitely when they exclusively compete for the same prey—one must inevitably become extinct. This paradigm is known as the competitive exclusion principle. We consider communities of any number of species and multiple trophic levels, i.e. the average number of steps between a predator and basic nutrient, e.g. sunlight or sugars. We show that the extension of the competitive exclusion principle to such large systems means that each species must be part of a “non-overlapping pairing”. Such pairings are exclusive connections between two species, e.g. a predator and a prey. We demonstrate that a stable food web can always be obtained if a non-overlapping pairing exists. The food web assembly rules are explicit conditions that specify sustainable combinations of species at the different trophic levels. As also seen in field data, our rules imply high species numbers at intermediate levels and few at the top and bottom. We further show that omnivorous species—those with hosts at multiple trophic levels—may take a special role in stabilizing food webs, as they combine several trophic levels.

Variably lytic infection dynamics of large Bacteroidetes podovirus phi38:1 against two Cellulophaga baltica host strains

Bacterial viruses (phages) influence global biogeochemical cycles by modulating bacterial mortality, metabolic output and evolution. However, our understanding of phage infections is limited by few methods and environmentally relevant model systems. Prior work showed that Cellulophaga baltica phage ϕ38:1 infects its original host lytically, and an alternative host either delayed lytically or lysogenically. Here we investigate these infections through traditional and marker-based approaches, and introduce geneELISA for high-throughput examination of phage-host interactions. All methods confirmed the lytic, original host infection (70-80 min latent period; approximately eight phages produced per cell), but alternative host assays were more challenging. A 4.5 h experiment detected no phage production by plaque assay, whereas phageFISH and geneELISA revealed phage genome replication and a latent period ≥ 150 min. Longer experiments (26 h) suggested an 11 h latent period and a burst size of 871 by plaque assay, whereas phageFISH identified cell lysis starting at < 5 h and lasting to 11 h, but for only 7% to 21.5% of infected cells, respectively, and with ∼ 39 phages produced per cell. These findings help resolve the nature of the alternative host infection as delayed lytic and offer solutions to methodological challenges for studying inefficient phage-host interactions.

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.

Effects of Long Term Antibiotic Therapy on Human Oral and Fecal Viromes

Viruses are integral members of the human microbiome. Many of the viruses comprising the human virome have been identified as bacteriophage, and little is known about how they respond to perturbations within the human ecosystem. The intimate association of phage with their cellular hosts suggests their communities may change in response to shifts in bacterial community membership. Alterations to human bacterial biota can result in human disease including a reduction in the host's resilience to pathogens. Here we report the ecology of oral and fecal viral communities and their responses to long-term antibiotic therapy in a cohort of human subjects. We found significant differences between the viral communities of each body site with a more heterogeneous fecal virus community compared with viruses in saliva. We measured the relative diversity of viruses, and found that the oral viromes were significantly more diverse than fecal viromes. There were characteristic changes in the membership of oral and fecal bacterial communities in response to antibiotics, but changes in fecal viral communities were less distinguishing. In the oral cavity, an abundance of papillomaviruses found in subjects on antibiotics suggests an association between antibiotics and papillomavirus production. Despite the abundance of papillomaviruses identified, in neither the oral nor the fecal viromes did antibiotic therapy have any significant impact upon overall viral diversity. There was, however, an apparent expansion of the reservoir of genes putatively involved in resistance to numerous classes of antibiotics in fecal viromes that was not paralleled in oral viromes. The emergence of antibiotic resistance in fecal viromes in response to long-term antibiotic therapy in humans suggests that viruses play an important role in the resilience of human microbial communities to antibiotic disturbances.

Restriction modification systems as engines of diversity

Restriction modification (RM) systems provide protection against a broad spectrum of phages. However, the likelihood of a phage permanently bypassing this can be as high as 0.1 per infection (Korona et al., 1993) which makes for a relatively weak defense. Here we argue that, apart from providing such transient defenses, RM systems can facilitate long-term coexistence of many bacterial strains. We show that this diversity can be as large as the burst size of the phage but no larger-a curious correspondence between a number at the level of species and another number at the level of individuals. Such a highly diverse and stably coexisting ecosystem is robust to substantial variation in both bacterial growth rates and strength of their RM systems, which might be one reason why quite weak RM systems exist in the wild.

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

Upon infection of their bacterial hosts temperate phages must chose between lysogenic and lytic developmental strategies. Here we apply the game-theoretic bet-hedging strategy introduced by Kelly to derive the optimal lysogenic fraction of the total population of phages as a function of frequency and intensity of environmental downturns affecting the lytic subpopulation. "Well-temperate" phage from our title is characterized by the best long-term population growth rate. We show that it is realized when the lysogenization frequency is approximately equal to the probability of lytic population collapse. We further predict the existence of sharp boundaries in system's environmental, ecological, and biophysical parameters separating the regions where this temperate strategy is optimal from those dominated by purely virulent or dormant (purely lysogenic) strategies. We show that the virulent strategy works best for phages with large diversity of hosts, and access to multiple independent environments reachable by diffusion. Conversely, progressively more temperate or even dormant strategies are favored in the environments, that are subject to frequent and severe temporal downturns.

A continuum of specialists and generalists in empirical communities

Understanding the persistence of specialists and generalists within ecological communities is a topical research question, with far-reaching consequences for the maintenance of functional diversity. Although theoretical studies indicate that restricted conditions may be necessary to achieve co-occurrence of specialists and generalists, analyses of larger empirical (and species-rich) communities reveal the pervasiveness of coexistence. In this paper, we analyze 175 ecological bipartite networks of three interaction types (animal hosts-parasite, plant-herbivore and plant-pollinator), and measure the extent to which these communities are composed of species with different levels of specificity in their biotic interactions. We find a continuum from specialism to generalism. Furthermore, we demonstrate that diversity tends to be greatest in networks with intermediate connectance, and argue this is because of physical constraints in the filling of networks.

What difference does it make if viruses are strain-, rather than species-specific?

Theoretical work has suggested an important role of lytic viruses in controlling the diversity of their prokaryotic hosts. Yet, providing strong experimental or observational support (or refutation) for this has proven evasive. Such models have usually assumed "host groups" to correspond to the "species" level, typically delimited by 16S rRNA gene sequence data. Recent model developments take into account the resolution of species into strains with differences in their susceptibility to viral attack. With strains as the host groups, the models will have explicit viral control of abundance at strain level, combined with explicit predator or resource control at community level, but the direct viral control at species level then disappears. Abundance of a species therefore emerges as the combination of how many strains, and at what abundance, this species can establish in competition with other species from a seeding community. We here discuss how species diversification and strain diversification may introduce competitors and defenders, respectively, and that the balance between the two may be a factor in the control of species diversity in mature natural communities. These models can also give a dominance of individuals from strains with high cost of resistance; suggesting that the high proportion of "dormant" cells among pelagic heterotrophic prokaryotes may reflect their need for expensive defense rather than the lack of suitable growth substrates in their environment.