Evolutionary Feedback Mediated through Population Density, Illustrated with Viruses in Chemostats

Warming drives feedback between plant phenotypes and ecosystem functioning in sub-Antarctic ponds

Ample evidence indicates that warming affects individuals in plant communities, ultimately threatening biodiversity. Individual plants in communities are also exposed to plant-plant interaction that may affect their performance. However, trait responses to these two constraints have usually been studied separately, while they may influence processes at the ecosystem level. In turn, these ecological modifications may impact the phenotypes of plants through nutrient availability and uptake. We developed an experimental approach based on the macrophyte communities in the ponds of the sub-Antarctic Iles Kerguelen. Individuals of the species Limosella australis were grown under different temperature × plant-plant interaction treatments to assess their trait responses and create litters with different characteristics. The litters were then decomposed in the presence of individual plants at different temperatures to examine effects on ecosystem functioning and potential feedback affecting plant trait values. Leaf resource-acquisition- and -conservation-related traits were altered in the context of temperature × plant-plant interaction. At 13 °C, SLA and leaf C:N were higher under interspecific and intraspecific interactions than without interaction, whereas at 23 °C, these traits increased under intraspecific interaction only. These effects only slightly improved the individual performance, suggesting that plant-plant interaction is an additional selective pressure on individuals in the context of climate warming. The decay rate of litter increased with the Leaf Carbon Content at 13 °C and 18 °C, but decreased at 23 °C. The highest decay rate was recorded at 18 °C. Besides, we observed evidence of positive feedback of the decay rate alone, and in interaction with the temperature, respectively on the leaf C:N and Leaf Dry Matter Content, suggesting that variations in ecological processes affect plant phenotypes. Our findings demonstrate that warming can directly and indirectly affect the evolutionary and ecological processes occurring in aquatic ecosystems through plants.

Unconstrained coevolution of bacterial size and the latent period of plastic phage

Viruses play critical roles in the dynamics of microbial communities. Lytic viruses, for example, kill significant fractions of autotrophic and heterotrophic microbes daily. The dynamic interplay between viruses and microbes results from an overlap of physiological, ecological, and evolutionary responses: environmental changes trigger host physiological changes, affecting the ecological interactions of host and virus and, ultimately, the evolutionary pressures influencing the two populations. Recent theoretical work studied how the dependence of viral traits on host physiology (viral plasticity) affects the evolutionarily stable host cell size and viral infection time emerging from coevolution. Here, we broaden the scope of the framework to consider any coevolutionary outcome, including potential evolutionary collapses of the system. We used the case study of Escherichia coli and T-like viruses under chemostat conditions, but the framework can be adapted to any microbe-virus system. Oligotrophic conditions led to smaller, lower-quality but more abundant hosts, and infections that were longer but produced a reduced viral offspring. Conversely, eutrophic conditions resulted in fewer but larger higher-quality hosts, and shorter but more productive infections. The virus influenced host evolution decreasing host size more noticeably for low than for high dilution rates, and for high than for low nutrient input concentration. For low dilution rates, the emergent infection time minimized host need/use, but higher dilution led to an opportunistic strategy that shortened the duration of infections. System collapses driven by evolution resulted from host failure to adapt quickly enough to the evolving virus. Our results contribute to understanding the eco-evolutionary dynamics of microbes and virus, and to improving the predictability of current models for host-virus interactions. The large quantitative and qualitative differences observed with respect to a classic description (in which viral traits are assumed to be constant) highlights the importance of including viral plasticity in theories describing short- and long-term host-virus dynamics.

Ecological Approach to Understanding Superinfection Inhibition in Bacteriophage

In microbial communities, viruses compete with each other for host cells to infect. As a consequence of competition for hosts, viruses evolve inhibitory mechanisms to suppress their competitors. One such mechanism is superinfection exclusion, in which a preexisting viral infection prevents a secondary infection. The bacteriophage ΦX174 exhibits a potential superinfection inhibition mechanism (in which secondary infections are either blocked or resisted) known as the reduction effect. In this auto-inhibitory phenomenon, a plasmid containing a fragment of the ΦX174 genome confers resistance to infection among cells that were once permissive to ΦX174. Taking advantage of this plasmid system, we examine the inhibitory properties of the ΦX174 reduction effect on a range of wild ΦX174-like phages. We then assess how closely the reduction effect in the plasmid system mimics natural superinfection inhibition by carrying out phage-phage competitions in continuous culture, and we evaluate whether the overall competitive advantage can be predicted by phage fitness or by a combination of fitness and reduction effect inhibition. Our results show that viral fitness often correctly predicts the winner. However, a phage's reduction sequence also provides an advantage to the phage in some cases, modulating phage-phage competition and allowing for persistence where competitive exclusion was expected. These findings provide strong evidence for more complex dynamics than were previously thought, in which the reduction effect may inhibit fast-growing viruses, thereby helping to facilitate coexistence.

Ancestral ecological regime shapes reaction to food limitation in the Least Killifish, Heterandria formosa

Populations with different densities often show genetically based differences in life histories. The divergent life histories could be driven by several agents of selection, one of which is variation in per-capita food levels. Its relationship with population density is complex, as it depends on overall food availability, individual metabolic demand, and food-independent factors potentially affecting density, such as predation intensity. Here, we present a case study of two populations of a small live-bearing freshwater fish, one characterized by high density, low predation risk, low overall food availability, and presumably low per-capita food levels, and the other by low density, high predation risk, high overall food availability, and presumably high per-capita food levels. Using a laboratory experiment, we examined whether fish from these populations respond differently to food limitation, and whether size at birth, a key trait with respect to density variation in this species, is associated with any such differential responses. While at the lower food level growth was slower, body size smaller, maturation delayed, and survival reduced in both populations, these fitness costs were smaller in fish from the high-density population. At low food, only 15% of high-density fish died, compared to 75% of low-density fish. This difference was much smaller at high food (0% vs. 15% mortality). The increased survival of high-density fish may, at least partly, be due to their larger size at birth. Moreover, being larger at birth enabled fish to mature relatively early even at the lower food level. We demonstrate that sensitivities to food limitation differ between study populations, consistent with selection for a greater ability to tolerate low per-capita food availability in the high-density population. While we cannot preclude other agents of selection from operating in these populations simultaneously, our results suggest that variation in per-capita food levels is one of those agents.

A High-Resolution View of Adaptive Event Dynamics in a Plasmid

Coadaptation between bacterial hosts and plasmids frequently results in adaptive changes restricted exclusively to host genome leaving plasmids unchanged. To better understand this remarkable stability, we transformed naïve Escherichia coli cells with a plasmid carrying an antibiotic-resistance gene and forced them to adapt in a turbidostat environment. We then drew population samples at regular intervals and subjected them to duplex sequencing—a technique specifically designed for identification of low-frequency mutations. Variants at ten sites implicated in plasmid copy number control emerged almost immediately, tracked consistently across the experiment’s time points, and faded below detectable frequencies toward the end. This variation crash coincided with the emergence of mutations on the host chromosome. Mathematical modeling of trajectories for adaptive changes affecting plasmid copy number showed that such mutations cannot readily fix or even reach appreciable frequencies. We conclude that there is a strong selection against alterations of copy number even if it can provide a degree of growth advantage. This incentive is likely rooted in the complex interplay between mutated and wild-type plasmids constrained within a single cell and underscores the importance of understanding of intracellular plasmid variability.

Bacteriophage Production Models: An Overview

The use of bacteriophages has been proposed as an alternative method to control pathogenic bacteria. During recent years several reports have been published about the successful use of bacteriophages in different fields such as food safety, agriculture, aquaculture, and even human health. Several companies are now commercializing bacteriophages or bacteriophage-based products for therapeutic purposes. However, this technology is still in development and there are challenges to overcome before bacteriophages can be widely used to control pathogenic bacteria. One big hurdle is the development of efficient methods for bacteriophage production. To date, several models for bacteriophage production have been reported, some of them evaluated experimentally. This mini-review offers an overview of different models and methods for bacteriophage production, contrasting their principal differences.

How Does the Evolution of Universal Ecological Traits Affect Population Size? Lessons from Simple Models

This article argues that adaptive evolutionary change in a consumer species should frequently decrease (and maladaptive change should increase) population size, producing adaptive decline. This conclusion is based on analysis of multiple consumer-resource models that examine evolutionary change in consumer traits affecting the universal ecological parameters of attack rate, conversion efficiency, and mortality. Two scenarios are investigated. In one, evolutionary equilibrium is initially maintained by opposing effects on the attack rate and other growth rate parameters; the environment or trait is perturbed, and the trait then evolves to a new (or back to a previous) equilibrium. Here evolution exhibits adaptive decline in up to one-half of all cases. The other scenario assumes a genetic perturbation having purely fitness-increasing effects. Here adaptive decline in the consumer requires that the resource be self-reproducing and overexploited and requires a sufficient increase in the attack rate. However, if the resource exhibits adaptive defense via behavior or evolution, adaptive decline may characterize consumer traits affecting all parameters. Favorable environmental change producing parameter shifts similar to those produced by adaptive evolution has similar counterintuitive effects on consumer population size. Many different food web models have already been shown to exhibit such counterintuitive changes in some species.

Ecological and Evolutionary Consequences of Viral Plasticity

Viruses use the host machinery to replicate, and their performance thus depends on the host's physiological state. For bacteriophages, this link between host and viral performance has been characterized empirically and with intracellular theories. Such theories are too detailed to be included in models that study host-phage interactions in the long term, which hinders our understanding of systems that range from pathogens infecting gut bacteria to marine phage shaping the oceans. Here, we combined data and models to study the short- and long-term consequences that host physiology has on bacteriophage performance. We compiled data showing the dependence of lytic-phage traits on host growth rate (referred to as viral phenotypic plasticity) to deduce simple expressions that represent such plasticity. Including these expressions in a standard host-phage model allowed us to understand mechanistically how viral plasticity affects emergent evolutionary strategies and the population dynamics associated with different environmental scenarios including, for example, nutrient pulses or host starvation. Moreover, we show that plasticity on the offspring number drives the phage ecological and evolutionary dynamics by reinforcing feedbacks between host, virus, and environment. Standard models neglect viral plasticity, which therefore handicaps their predictive ability in realistic scenarios. Our results highlight the importance of viral plasticity to unravel host-phage interactions and the need of laboratory and field experiments to characterize viral plastic responses across systems.

High Throughput Manufacturing of Bacteriophages Using Continuous Stirred Tank Bioreactors Connected in Series to Ensure Optimum Host Bacteria Physiology for Phage Production

Future industrial demand for large quantities of bacteriophages e.g., for phage therapy, necessitates the development of scalable Good Manufacturing Practice compliant (cGMP) production platforms. The continuous production of high titres of E coli T3 phages (10¹¹ PFU mL⁻¹) was achieved using two continuous stirred tank bioreactors connected in series, and a third bioreactor was used as a final holding tank operated in semi-batch mode to finish the infection process. The first bioreactor allowed the steady-state propagation of host bacteria using a fully synthetic medium with glucose as the limiting substrate. Host bacterial growth was decoupled from the phage production reactor downstream of it to suppress the production of phage-resistant mutants, thereby allowing stable operation over a period of several days. The novelty of this process is that the manipulation of the host reactor dilution rates (range 0.1–0.6 hr⁻¹) allows control over the physiological state of the bacterial population. This results in bacteria with considerably higher intracellular phage production capability whilst operating at high dilution rates yielding significantly higher overall phage process productivity. Using a pilot-scale chemostat system allowed optimisation of the upstream phage amplification conditions conducive for high intracellular phage production in the host bacteria. The effect of the host reactor dilution rates on the phage burst size, lag time, and adsorption rate were evaluated. The host bacterium physiology was found to influence phage burst size, thereby affecting the productivity of the overall process. Mathematical modelling of the dynamics of the process allowed parameter sensitivity evaluation and provided valuable insights into the factors affecting the phage production process. The approach presented here may be used at an industrial scale to significantly improve process control, increase productivity via process intensification, and reduce process manufacturing costs through process footprint reduction.

Evolutionarily stable communities: a framework for understanding the role of trait evolution in the maintenance of diversity

Biological diversity depends on the interplay between evolutionary diversification and ecological mechanisms allowing species to coexist. Current research increasingly integrates ecology and evolution over a range of timescales, but our common conceptual framework for understanding species coexistence requires better incorporation of evolutionary processes. Here, we focus on the idea of evolutionarily stable communities (ESCs), which are theoretical endpoints of evolution in a community context. We use ESCs as a unifying framework to highlight some important but under-appreciated theoretical results, and we review empirical research relevant to these theoretical predictions. We explain how, in addition to generating diversity, evolution can also limit diversity by reducing the effectiveness of coexistence mechanisms. The coevolving traits of competing species may either diverge or converge, depending on whether the number of species in the community is low (undersaturated) or high (oversaturated) relative to the ESC. Competition in oversaturated communities can lead to extinction or neutrally coexisting, ecologically equivalent species. It is critical to consider trait evolution when investigating fundamental ecological questions like the strength of different coexistence mechanisms, the feasibility of ecologically equivalent species, and the interpretation of different patterns of trait dispersion.

Bacteriophage Production in Bioreactors

The optimal conditions for the production of virulent bacteriophages in bioreactors can vary greatly depending on the host-bacteriophage system used. We present a general method for the production of virulent bacteriophages in bioreactors that can be adapted to many host-bacteriophage systems and various operating conditions (reactor volume, medium composition, temperature, etc.). The procedures detail how to establish optimal initial infection conditions (infection load and initial multiplicity of infection (MOI)), prepare the host pre-culture and bioreactor, operate the bioreactor, and harvest the bacteriophage product. Batch operation is detailed but a short discussion addresses other modes of operation, namely two-stage continuous bioreactors and two-stage cycling bioreactors.

Hot spots become cold spots: coevolution in variable temperature environments

Antagonistic coevolution between hosts and parasites is a key process in the genesis and maintenance of biological diversity. Whereas coevolutionary dynamics show distinct patterns under favourable environmental conditions, the effects of more realistic, variable conditions are largely unknown. We investigated the impact of a fluctuating environment on antagonistic coevolution in experimental microcosms of Pseudomonas fluorescens SBW25 and lytic phage SBWΦ2. High-frequency temperature fluctuations caused no deviations from typical coevolutionary arms race dynamics. However, coevolution was stalled during periods of high temperature under intermediate- and low-frequency fluctuations, generating temporary coevolutionary cold spots. Temperature variation affected population density, providing evidence that eco-evolutionary feedbacks act through variable bacteria-phage encounter rates. Our study shows that environmental fluctuations can drive antagonistic species interactions into and out of coevolutionary cold and hot spots. Whether coevolution persists or stalls depends on the frequency of change and the environmental optima of both interacting players.

Genetically Determined Variation in Lysis Time Variance in the Bacteriophage φX174

Researchers in evolutionary genetics recently have recognized an exciting opportunity in decomposing beneficial mutations into their proximal, mechanistic determinants. The application of methods and concepts from molecular biology and life history theory to studies of lytic bacteriophages (phages) has allowed them to understand how natural selection sees mutations influencing life history. This work motivated the research presented here, in which we explored whether, under consistent experimental conditions, small differences in the genome of bacteriophage φX174 could lead to altered life history phenotypes among a panel of eight genetically distinct clones. We assessed the clones' phenotypes by applying a novel statistical framework to the results of a serially sampled parallel infection assay, in which we simultaneously inoculated each of a large number of replicate host volumes with ∼1 phage particle. We sequentially plated the volumes over the course of infection and counted the plaques that formed after incubation. These counts served as a proxy for the number of phage particles in a single volume as a function of time. From repeated assays, we inferred significant, genetically determined heterogeneity in lysis time and burst size, including lysis time variance. These findings are interesting in light of the genetic and phenotypic constraints on the single-protein lysis mechanism of φX174. We speculate briefly on the mechanisms underlying our results, and we discuss the potential importance of lysis time variance in viral evolution.

Population and Evolutionary Dynamics based on Predator-Prey Relationships in a 3D Physical Simulation

Recent studies have reported that population dynamics and evolutionary dynamics, occurring at different time scales, can be affected by each other. Our purpose is to explore the interaction between population and evolutionary dynamics using an artificial life approach based on a 3D physically simulated environment in the context of predator-prey and morphology-behavior coevolution. The morphologies and behaviors of virtual prey creatures are evolved using a genetic algorithm based on the predation interactions between predators and prey. Both population sizes are also changed, depending on the fitness. We observe two types of cyclic behaviors, corresponding to short-term and long-term dynamics. The former can be interpreted as a simple population dynamics of Lotka-Volterra type. It is shown that the latter cycle is based on the interaction between the changes in the prey strategy against predators and the long-term change in both population sizes, resulting partly from a tradeoff between their defensive success and the cost of defense.

The habits of highly effective phages: population dynamics as a framework for identifying therapeutic phages

The use of bacteriophages as antibacterial agents is being actively researched on a global scale. Typically, the phages used are isolated from the wild by plating on the bacteria of interest, and a far larger set of candidate phages is often available than can be used in any application. When an excess of phages is available, how should the best phages be identified? Here we consider phage-bacterial population dynamics as a basis for evaluating and predicting phage success. A central question is whether the innate dynamical properties of phages are the determinants of success, or instead, whether extrinsic, indirect effects can be responsible. We address the dynamical perspective, motivated in part by the absence of dynamics in previously suggested principles of phage therapy. Current mathematical models of bacterial-phage dynamics do not capture the realities of in vivo dynamics, nor is this likely to change, but they do give insight to qualitative properties that may be generalizable. In particular, phage adsorption rate may be critical to treatment success, so understanding the effects of the in vivo environment on host availability may allow prediction of useful phages prior to in vivo experimentation. Principles for predicting efficacy may be derived by developing a greater understanding of the in vivo system, or such principles could be determined empirically by comparing phages with known differences in their dynamic properties. The comparative approach promises to be a powerful method of discovering the key to phage success. We offer five recommendations for future study: (i) compare phages differing in treatment efficacy to identify the phage properties associated with success, (ii) assay dynamics in vivo, (iii) understand mechanisms of bacterial escape from phages, (iv) test phages in model infections that are relevant to the intended clinical applications, and (v) develop new classes of models for phage growth in spatially heterogeneous environments.

A new Microviridae phage isolated from a failed biotechnological process driven by Escherichia coli

Bacteriophages are present in every environment that supports bacterial growth, including man made ecological niches. Virulent phages may even slow or, in more severe cases, interrupt bioprocesses driven by bacteria. Escherichia coli is one of the most widely used bacteria for large-scale bioprocesses; however, literature describing phage-host interactions in this industrial context is sparse. Here, we describe phage MED1 isolated from a failed industrial process. Phage MED1 (Microviridae family, with a single-stranded DNA [ssDNA] genome) is highly similar to the archetypal phage phiX174, sharing >95% identity between their genomic sequences. Whole-genome phylogenetic analysis of 52 microvirus genomes from public databases revealed three genotypes (alpha3, G4, and phiX174). Phage MED1 belongs to the phiX174 group. We analyzed the distribution of single nucleotide variants in MED1 and 18 other phiX174-like genomes and found that there are more missense mutations in genes G, B, and E than in the other genes of these genomes. Gene G encodes the spike protein, involved in host attachment. The evolution of this protein likely results from the selective pressure on phages to rapidly adapt to the molecular diversity found at the surface of their hosts.

Evolutionary comparison between viral lysis rate and latent period

Marine viruses shape the structure of the microbial community. They are, thus, a key determinant of the most important biogeochemical cycles in the planet. Therefore, a correct description of the ecological and evolutionary behavior of these viruses is essential to make reliable predictions about their role in marine ecosystems. The infection cycle, for example, is indistinctly modeled in two very different ways. In one representation, the process is described including explicitly a fixed delay between infection and offspring release. In the other, the offspring are released at exponentially distributed times according to a fixed release rate. By considering obvious quantitative differences pointed out in the past, the latter description is widely used as a simplification of the former. However, it is still unclear how the dichotomy "delay versus rate description" affects long-term predictions of host-virus interaction models. Here, we study the ecological and evolutionary implications of using one or the other approaches, applied to marine microbes. To this end, we use mathematical and eco-evolutionary computational analysis. We show that the rate model exhibits improved competitive abilities from both ecological and evolutionary perspectives in steady environments. However, rate-based descriptions can fail to describe properly long-term microbe-virus interactions. Moreover, additional information about trade-offs between life-history traits is needed in order to choose the most reliable representation for oceanic bacteriophage dynamics. This result affects deeply most of the marine ecosystem models that include viruses, especially when used to answer evolutionary questions.

Lethal mutagenesis failure may augment viral adaptation

Lethal mutagenesis, the attempt to extinguish a population by elevating its mutation rate, has been endorsed in the virology literature as a promising approach for treating viral infections. In support of the concept, in vitro studies have forced viral extinction with high doses of mutagenic drugs. However, the one known mutagenic drug used on patients commonly fails to cure infections, and in vitro studies typically find a wide range of mutagenic conditions permissive for viral growth. A key question becomes how subsequent evolution is affected if the viral population is mutated but avoids extinction--Is viral adaptation augmented rather than suppressed? Here we consider the evolution of highly mutated populations surviving mutagenesis, using the DNA phage T7. In assays using inhibitory hosts, whenever resistance mutants were observed, the mutagenized populations exhibited higher frequencies, but some inhibitors blocked plaque formation by even the mutagenized stock. Second, outgrowth of previously mutagenized populations led to rapid and potentially complete fitness recovery but polymorphism was slow to decay, and mutations exhibited inconsistent patterns of change. Third, the combination of population bottlenecks with mutagenesis did cause fitness declines, revealing a vulnerability that was not apparent from mutagenesis of large populations. The results show that a population surviving high mutagenesis may exhibit enhanced adaptation in some environments and experience little negative fitness consequences in many others.

Selection affects genes involved in replication during long-term evolution in experimental populations of the bacteriophage φX174

Observing organisms that evolve in response to strong selection over very short time scales allows the determination of the molecular mechanisms underlying adaptation. Although dissecting these molecular mechanisms is expensive and time-consuming, general patterns can be detected from repeated experiments, illuminating the biological processes involved in evolutionary adaptation. The bacteriophage φX174 was grown for 50 days in replicate chemostats under two culture conditions: Escherichia coli C as host growing at 37°C and Salmonella typhimurium as host growing at 43.5°C. After 50 days, greater than 20 substitutions per chemostat had risen to detectable levels. Of the 97 substitutions, four occurred in all four chemostats, five arose in both culture conditions, eight arose in only the high temperature S. typhimurium chemostats, and seven arose only in the E. coli chemostats. The remaining substitutions were detected only in a single chemostat, however, almost half of these have been seen in other similar experiments. Our findings support previous studies that host recognition and capsid stability are two biological processes that are modified during adaptation to novel hosts and high temperature. Based upon the substitutions shared across both environments, it is apparent that genome replication and packaging are also affected during adaptation to the chemostat environment, rather than to temperature or host per se. This environment is characterized by a large number of phage and very few hosts, leading to competition among phage within the host. We conclude from these results that adaptation to a high density environment selects for changes in genome replication at both protein and DNA sequence levels.

Evolution of virulence in emerging epidemics

Theory predicts that selection for pathogen virulence and horizontal transmission is highest at the onset of an epidemic but decreases thereafter, as the epidemic depletes the pool of susceptible hosts. We tested this prediction by tracking the competition between the latent bacteriophage λ and its virulent mutant λcI857 throughout experimental epidemics taking place in continuous cultures of Escherichia coli. As expected, the virulent λcI857 is strongly favored in the early stage of the epidemic, but loses competition with the latent virus as prevalence increases. We show that the observed transient selection for virulence and horizontal transmission can be fully explained within the framework of evolutionary epidemiology theory. This experimental validation of our predictions is a key step towards a predictive theory for the evolution of virulence in emerging infectious diseases.

Evolutionary stability of a refactored phage genome

Engineered genetic systems are commonly unstable; if propagated, they evolve to reverse or modify engineered elements because the elements impair fitness. A goal of synthetic biology is thus to anticipate and avoid detrimental engineering, but little is yet known about which types of elements cause problems in different contexts. In prior work, 30% of the genome of bacteriophage T7 was “refactored” by the insertion or modification of 65 short sequences that included a useful restriction enzyme site in order to, among other goals, separate genes and their translational initiation regions from each other and from other genetic elements. Although gene sequences and known important regions of regulatory elements were kept intact, the translational efficiency of some genes or element regulatory function might have been compromised. We adapted the refactored phage for rapid growth in two conditions, observing fitness and sequence evolution. As anticipated from the original work, refactoring had major fitness effects in both environments, but most of the fitness costs were recovered upon adaptation. The evolved phages retained 60–70% of the design elements, suggesting they had only minor fitness effects. Approximately half the elements that were lost lie within large deletions commonly observed during adaptation of the wild-type genome. Some elements were lost or modified in parallel between the adaptations without affecting T7 gene sequences, but no obvious correlates can be made. Nevertheless, experimental adaptations are useful for identifying specific synthetic design problems, and we suggest that experimental evolution in conjunction with alternative engineering may also be useful in overcoming those problems.

Competition-colonization trade-off promotes coexistence of low-virulence viral strains

RNA viruses exist as genetically diverse populations displaying a range of virulence degrees. The evolution of virulence in viral populations is, however, poorly understood. On the basis of the experimental observation of an RNA virus clone in cell culture diversifying into two subpopulations of different virulence, we study the dynamics of mutating virus populations with varying virulence. We introduce a competition-colonization trade-off into standard mathematical models of intra-host viral infection. Colonizers are fast-spreading virulent strains, whereas the competitors are less-virulent variants but more successful within co-infected cells. We observe a two-step dynamics of the population. Early in the infection, the population is dominated by colonizers, which later are outcompeted by competitors. Our simulations suggest the existence of steady state in which all virulence classes coexist but are dominated by the most competitive ones. This equilibrium implies collective virulence attenuation in the population, in contrast to previous models predicting evolution of the population towards increased virulence.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

Slow fitness recovery in a codon-modified viral genome

Extensive synonymous codon modification of viral genomes appears to be an effective way of attenuating strains for use as live vaccines. An assumption of this method is that codon changes have individually small effects, such that codon-attenuated viruses will be slow to evolve back to high fitness (and thus to high virulence). The major capsid gene of the bacterial virus T7 was modified to have varying levels of suboptimal synonymous codons in different constructs, and fitnesses declined linearly with the number of changes. Adaptation of the most extreme design, with 182 codon changes, resulted in a slow fitness recovery by standards of previous experimental evolution with this virus, although fitness effects of substitutions were higher than expected from the average effect of an engineered codon modification. Molecular evolution during recovery was modest, and changes evolved both within the modified gene and outside it. Some changes within the modified gene evolved in parallel across replicates, but with no obvious explanation. Overall, the study supports the premise that codon-modified viruses recover fitness slowly, although the evolution is substantially more rapid than expected from the design principle.

Multiple genetic pathways to similar fitness limits during viral adaptation to a new host

The gain in fitness during adaptation depends on the supply of beneficial mutations. Despite a good theoretical understanding of how evolution proceeds for a defined set of mutations, there is little understanding of constraints on net fitness-whether fitness will reach a limit despite ongoing selection and mutation, and if there is a limit, what determines it. Here, the dsDNA bacteriophage SP6, a virus of Salmonella, was adapted to Escherichia coli K-12. From an isolate capable of modest growth on E. coli, four lines were adapted for rapid growth by protocols differing in use of mutagen, propagation method, and duration, but using the same media, temperature, and a continual excess of the novel host. Nucleotide changes underlying those adaptations differed greatly in number and identity, but the four lines achieved similar absolute fitness at the end, an increase of more than 4000-fold phage descendants per hour. Thus, the fitness landscape allows multiple genetic paths to the same approximate fitness limit. The existence and causes of fitness limits have ramifications to genome engineering, vaccine design, and "lethal mutagenesis" treatments to cure viral infections.

Trophic network structure emerges through antagonistic coevolution in temporally varying environments

Understanding the mechanisms underlying ecological specialization is central to our understanding of community ecology and evolution. Although theoretical work has investigated how variable environments may affect specialization in single species, little is known about how such variation impacts bipartite network structure in antagonistically coevolving systems. Here, we develop and analyse a general model of victim-enemy coevolution that explicitly includes resource and population dynamics. We investigate how temporal environmental heterogeneity affects the evolution of specialization and associated community structure. Environmental productivity influences victim investment in resistance, which will shape patterns of specialization through its regulating effect on enemy investment in infectivity. We also investigate the epidemiological consequences of environmental variability and show that enemy population density is maximized for intermediate lengths of productive seasons, which corresponds to situations where enemies can evolve higher infectivity than victims can evolve defence. We discuss our results in the light of empirical studies, and further highlight ways in which our model applies to a range of natural systems.

Competition-colonization dynamics: An ecology approach to quasispecies dynamics and virulence evolution in RNA viruses

A single and purified clone of foot-and-mouth disease virus diversified in cell culture into two subpopulations that were genetically distinct. The subpopulation with higher virulence was a minority and was suppressed by the dominant but less virulent one. These two populations follow the competitioncolonization dynamics described in ecology. Virulent viruses can be regarded as colonizers because they killed the cells faster and they spread faster. The attenuated subpopulation resembles competitors because of its higher replication efficiency in coinfected cells. Our results suggest a new model for the evolution of virulence which is based on interactions between components of the quasispecies. Competition between viral mutants takes place at two levels, intracellular competition and competition for new cells. The two strategies are subjected to densitydependent selection.

Using artificial systems to explore the ecology and evolution of symbioses

The web of life is weaved from diverse symbiotic interactions between species. Symbioses vary from antagonistic interactions such as competition and predation to beneficial interactions such as mutualism. What are the bases for the origin and persistence of symbiosis? What affects the ecology and evolution of symbioses? How do symbiotic interactions generate ecological patterns? How do symbiotic partners evolve and coevolve? Many of these questions are difficult to address in natural systems. Artificial systems, from abstract to living, have been constructed to capture essential features of natural symbioses and to address these key questions. With reduced complexity and increased controllability, artificial systems can serve as useful models for natural systems. We review how artificial systems have contributed to our understanding of symbioses.

Positive selection at high temperature reduces gene transcription in the bacteriophage ϕX174

Background

Gene regulation plays a central role in the adaptation of organisms to their environments. There are many molecular components to gene regulation, and it is often difficult to determine both the genetic basis of adaptation and the evolutionary forces that influence regulation. In multiple evolution experiments with the bacteriophage ϕX174, adaptive substitutions in cis-acting regulatory sequences sweep through the phage population as the result of strong positive selection at high temperatures that are non-permissive for laboratory-adapted phage. For one cis-regulatory region, we investigate the individual effects of four adaptive substitutions on transcript levels and fitness for phage growing on three hosts at two temperatures.

Results

The effect of the four individual substitutions on transcript levels is to down-regulate gene expression, regardless of temperature or host. To ascertain the conditions under which these substitutions are adaptive, fitness was measured by a variety of methods for several bacterial hosts growing at two temperatures, the control temperature of 37°C and the selective temperature of 42°C. Time to lysis and doublings per hour indicate that the four substitutions individually improve fitness over the ancestral strain at high temperature independent of the bacterial host in which the fitness was measured. Competition assays between the ancestral strain and either of two mutant strains indicate that both mutants out-compete the ancestor at high temperature, but the relative frequencies of each phage remain the same at the control temperature.

Conclusions

Our results strongly suggest that gene transcription plays an important role in influencing fitness in the bacteriophage ϕX174, and different point mutations in a single cis-regulatory region provided the genetic basis for this role in adaptation to high temperature. We speculate that the adaptive nature of these substitutions is due to the physiology of the host at high temperature or the need to maintain particular ratios of phage proteins during capsid assembly. Our investigation of regulatory evolution contributes to interpreting genome-level assessments of regulatory variation, as well as to understanding the molecular basis of adaptation.

Models of RNA virus evolution and their roles in vaccine design

Viruses are fast evolving pathogens that continuously adapt to the highly variable environments they live and reproduce in. Strategies devoted to inhibit virus replication and to control their spread among hosts need to cope with these extremely heterogeneous populations and with their potential to avoid medical interventions. Computational techniques such as phylogenetic methods have broadened our picture of viral evolution both in time and space, and mathematical modeling has contributed substantially to our progress in unraveling the dynamics of virus replication, fitness, and virulence. Integration of multiple computational and mathematical approaches with experimental data can help to predict the behavior of viral pathogens and to anticipate their escape dynamics. This piece of information plays a critical role in some aspects of vaccine development, such as viral strain selection for vaccinations or rational attenuation of viruses. Here we review several aspects of viral evolution that can be addressed quantitatively, and we discuss computational methods that have the potential to improve vaccine design.

Experimental evolution of viruses: Microviridae as a model system

phiX174 was developed as a model system for experimental studies of evolution because of its small genome size and ease of cultivation. It has been used extensively to address statistical questions about the dynamics of adaptive evolution. Molecular changes seen during experimental evolution of phiX174 under a variety of conditions were compiled from 10 experiments comprising 58 lineages, where whole genomes were sequenced. A total of 667 substitutions was seen. Parallel evolution was rampant, with over 50 per cent of substitutions occurring at sites with three or more events. Comparisons of experimentally evolved sites to variation seen among wild phage suggest that at least some of the adaptive mechanisms seen in the laboratory are relevant to adaptation in nature. Elucidation of these mechanisms is aided by the availability of capsid and pro-capsid structures for phiX174 and builds on years of genetic studies of the phage life history.

Optimality models in the age of experimental evolution and genomics

Optimality models have been used to predict evolution of many properties of organisms. They typically neglect genetic details, whether by necessity or design. This omission is a common source of criticism, and although this limitation of optimality is widely acknowledged, it has mostly been defended rather than evaluated for its impact. Experimental adaptation of model organisms provides a new arena for testing optimality models and for simultaneously integrating genetics. First, an experimental context with a well-researched organism allows dissection of the evolutionary process to identify causes of model failure--whether the model is wrong about genetics or selection. Second, optimality models provide a meaningful context for the process and mechanics of evolution, and thus may be used to elicit realistic genetic bases of adaptation--an especially useful augmentation to well-researched genetic systems. A few studies of microbes have begun to pioneer this new direction. Incompatibility between the assumed and actual genetics has been demonstrated to be the cause of model failure in some cases. More interestingly, evolution at the phenotypic level has sometimes matched prediction even though the adaptive mutations defy mechanisms established by decades of classic genetic studies. Integration of experimental evolutionary tests with genetics heralds a new wave for optimality models and their extensions that does not merely emphasize the forces driving evolution.

Viral and microbial community dynamics in four aquatic environments

The species composition and metabolic potential of microbial and viral communities are predictable and stable for most ecosystems. This apparent stability contradicts theoretical models as well as the viral-microbial dynamics observed in simple ecosystems, both of which show Kill-the-Winner behavior causing cycling of the dominant taxa. Microbial and viral metagenomes were obtained from four human-controlled aquatic environments at various time points separated by one day to >1 year. These environments were maintained within narrow geochemical bounds and had characteristic species composition and metabolic potentials at all time points. However, underlying this stability were rapid changes at the fine-grained level of viral genotypes and microbial strains. These results suggest a model wherein functionally redundant microbial and viral taxa are cycling at the level of viral genotypes and virus-sensitive microbial strains. Microbial taxa, viral taxa, and metabolic function persist over time in stable ecosystems and both communities fluctuate in a Kill-the-Winner manner at the level of viral genotypes and microbial strains.

Competition-colonization dynamics in an RNA virus

During replication, RNA viruses rapidly generate diverse mutant progeny which differ in their ability to kill host cells. We report that the progeny of a single RNA viral genome diversified during hundreds of passages in cell culture and self-organized into two genetically distinct subpopulations that exhibited the competition-colonization dynamics previously recognized in many classical ecological systems. Viral colonizers alone were more efficient in killing cells than competitors in culture. In cells coinfected with both competitors and colonizers, viral interference resulted in reduced cell killing, and competitors replaced colonizers. Mathematical modeling of this coinfection dynamics predicted selection to be density dependent, which was confirmed experimentally. Thus, as is known for other ecological systems, biodiversity and even cell killing of virus populations can be shaped by a tradeoff between competition and colonization. Our results suggest a model for the evolution of virulence in viruses based on internal interactions within mutant spectra of viral quasispecies.

High-resolution mapping of evolutionary trajectories in a phage

Experimental evolution in rapidly reproducing viruses offers a robust means to infer substitution trajectories during evolution. But with conventional approaches, this inference is limited by how many individual genotypes can be sampled from the population at a time. Low-frequency changes are difficult to detect, potentially rendering early stages of adaptation unobservable. Here we circumvent this using short-read sequencing technology in a fine-grained analysis of polymorphism dynamics in the sentinel organism: a single-stranded DNA phage PhiX174. Nucleotide differences were educed from noise with binomial filtering methods that harnessed quality scores and separate data from brief phage amplifications. Remarkably, a significant degree of variation was observed in all samples including those grown in brief 2-h cultures. Sites previously reported as subject to high-frequency polymorphisms over a course of weeks exhibited monotonic increases in polymorphism frequency within hours in this study. Additionally, even with limitations imposed by the short length of sequencing reads, we were able to observe statistically significant linkage among polymorphic sites in evolved lineages. Additional parallels between replicate lineages were apparent in the sharing of polymorphic sites and in correlated polymorphism frequencies. Missense mutations were more likely to occur than silent mutations. This study offers the first glimpse into "real-time" substitution dynamics and offers a robust conceptual framework for future viral resequencing studies.

Space, time, and host evolution facilitate coexistence of competing bacteriophages: theory and experiment

We present a joint experimental/theoretical investigation into the roles of spatial structure and time in the competition between two pathogens for a single host. We suggest a natural mechanism by which competing pathogens can coexist when host evolution and competitive dynamics occur on similar timescales. Our experimental system consisted of a single bacterial host species and two competing bacteriophage strains grown on agar plates, with a serial transfer of samples of the bacteriophage population to fresh host populations after each incubation cycle. The experiments included two incubation times and two transfer protocols that either maintained or disrupted the spatial structure of the viruses at each transfer. The same bacteriophage acted as the dominant competitor under both transfer protocols. A striking difference between the treatments is that the weak competitor was able to persist in the long-incubation experiments but not in the short-incubation experiments. Mathematical and experimental evidence suggest that coexistence is due to the appearance of resistant mutant host cells that provide a transient "spatiotemporal refuge" for the weaker competitor. Our mathematical model is individual based, captures the stochastic spatial dynamics down to the level of individual cells, and helps to explain the differences in behavior under the various experimental conditions.

Bacteriophage plaques: theory and analysis

Laboratory characterization of bacteriophage growth traditionally is done either in broth cultures or in semisolid agar media. These two environments may be distinguished in terms of their spatial structure, i.e., the degree to which they limit diffusion, motility, and environmental mixing. Well-mixed broth, for example, represents the microbiological ideal of a non-spatially structured environment. Agar, by contrast, imposes significant limitations on phage and bacterial movement and therefore gives rise to spatial structure. The study of phage growth within spatially structured environments, such as that seen during phage plaque formation, is important for three reasons. First, a large fraction of environmental bacteria live within spatially structured environments such as within biofilms, within soil, or when growing in or on the tissues of plants and animals. Second, phage growth as plaques is a central technique to phage studies, yet appears to be under appreciated by phage workers in terms of its complexity. Third, selective pressures acting on phage during plaque growth differ from those seen during broth growth. In this chapter we will discuss just what a plaque is, how one forms, and what can affect plaque size. We will describe methods, both experimental and theoretical, that have been employed to study plaque growth. As caveats we will discuss why plaque formation failure is not necessarily equivalent to virion inviability (Note 1). We also will consider problems with inferring phage broth growth fitness as a function of plaque characteristics (Note 2).