Clonal interference and the evolution of RNA viruses

Severe hindrance of viral infection propagation in spatially extended hosts

The production of large progeny numbers affected by high mutation rates is a ubiquitous strategy of viruses, as it promotes quick adaptation and survival to changing environments. However, this situation often ushers in an arms race between the virus and the host cells. In this paper we investigate in depth a model for the dynamics of a phenotypically heterogeneous population of viruses whose propagation is limited to two-dimensional geometries, and where host cells are able to develop defenses against infection. Our analytical and numerical analyses are developed in close connection to directed percolation models. In fact, we show that making the space explicit in the model, which in turn amounts to reducing viral mobility and hindering the infective ability of the virus, connects our work with similar dynamical models that lie in the universality class of directed percolation. In addition, we use the fact that our model is a multicomponent generalization of the Domany-Kinzel probabilistic cellular automaton to employ several techniques developed in the past in that context, such as the two-site approximation to the extinction transition line. Our aim is to better understand propagation of viral infections with mobility restrictions, e.g., in crops or in plant leaves, in order to inspire new strategies for effective viral control.

Frequent beneficial mutations during single-colony serial transfer of Streptococcus pneumoniae

The appearance of new mutations within a population provides the raw material for evolution. The consistent decline in fitness observed in classical mutation accumulation studies has provided support for the long-held view that deleterious mutations are more common than beneficial mutations. Here we present results of a study using a mutation accumulation design with the bacterium Streptococcus pneumoniae in which the fitness of the derived populations increased. This rise in fitness was associated specifically with adaptation to survival during brief stationary phase periods between single-colony population bottlenecks. To understand better the population dynamics behind this unanticipated adaptation, we developed a maximum likelihood model describing the processes of mutation and stationary-phase selection in the context of frequent population bottlenecks. Using this model, we estimate that the rate of beneficial mutations may be as high as 4.8×10⁻⁴ events per genome for each time interval corresponding to the pneumococcal generation time. This rate is several orders of magnitude higher than earlier estimates of beneficial mutation rates in bacteria but supports recent results obtained through the propagation of small populations of Escherichia coli. Our findings indicate that beneficial mutations may be relatively frequent in bacteria and suggest that in S. pneumoniae, which develops natural competence for transformation, a steady supply of such mutations may be available for sampling by recombination.

Beneficial mutations have long been considered extremely rare events and were thought to occur with a frequency of approximately one out of a billion times that a bacterium replicates its genome. Rare beneficial mutations would then be amplified by natural selection from the more frequent background of harmful mutations. Mutation accumulation experiments probe the nature of these spontaneous mutations by monitoring changes in fitness of model organisms propagated in the laboratory through numerous generations under conditions where the effects of selection are minimal. Previous mutation accumulation experiments have shown that organisms under study have declined in fitness as random mutations accrue in their genomes, consistent with a predominance of deleterious mutations. We conducted a mutation accumulation study with the bacterial pathogen S. pneumoniae in which a broad measure of fitness instead rose. We demonstrate that this unexpected adaptation was due to frequent beneficial mutations that were further amplified by selection in stationary-phase bacterial colonies. Together with recent work using E. coli, these results demonstrate that beneficial mutations can be common in bacteria and may contribute to our understanding of the evolution of traits such as antibiotic resistance and virulence.

Why do RNA viruses recombine?

Recombination occurs in many RNA viruses and can be of major evolutionary significance. However, rates of recombination vary dramatically among RNA viruses, which can range from clonal to highly recombinogenic. Here, we review the factors that might explain this variation in recombination frequency and show that there is little evidence that recombination is favoured by natural selection to create advantageous genotypes or purge deleterious mutations, as predicted if recombination functions as a form of sexual reproduction. Rather, recombination rates seemingly reflect larger-scale patterns of viral genome organization, such that recombination may be a mechanistic by-product of the evolutionary pressures acting on other aspects of virus biology.

Phylogenetic analysis of eastern equine encephalitis virus isolates from Florida

Florida has the highest degree of endemicity for eastern equine encephalitis virus (EEEV) of any state in the United States and is the only state with year-round transmission of EEEV. To further understand the viral population dynamics in Florida, the genome sequence of six EEEV isolates from central Florida were determined. These data were used to identify the most polymorphic regions of the EEEV genome from viruses isolated in Florida. The sequence of these polymorphic regions was then determined for 18 additional Florida isolates collected in four geographically distinct regions over a 20-year period. Phylogenetic analyses of these data suggested a rough temporal association of the Florida isolates, but no clustering by region or by source of the isolate. Some clustering of northeastern isolates with Florida isolates was seen, providing support for the hypothesis that Florida serves as a reservoir for the periodic introduction of EEEV into the northeastern United States.

Genetic variation and the fate of beneficial mutations in asexual populations

The fate of a newly arising beneficial mutation depends on many factors, such as the population size and the availability and fitness effects of other mutations that accumulate in the population. It has proved difficult to understand how these factors influence the trajectories of particular mutations, since experiments have primarily focused on characterizing successful clones emerging from a small number of evolving populations. Here, we present the results of a massively parallel experiment designed to measure the full spectrum of possible fates of new beneficial mutations in hundreds of experimental yeast populations, whether these mutations are ultimately successful or not. Using strains in which a particular class of beneficial mutation is detectable by fluorescence, we followed the trajectories of these beneficial mutations across 592 independent populations for 1000 generations. We find that the fitness advantage provided by individual mutations plays a surprisingly small role. Rather, underlying “background” genetic variation is quickly generated in our initially clonal populations and plays a crucial role in determining the fate of each individual beneficial mutation in the evolving population.

Role of evolved host breadth in the initial emergence of an RNA virus

Understanding how evolution promotes pathogen emergence would aid disease management, and prediction of future host shifts. Increased pathogen infectiousness of different hosts may occur through direct selection, or fortuitously via indirect selection. However, it is unclear which type of selection tends to produce host breadth promoting pathogen emergence. We predicted that direct selection for host breadth should foster emergence by causing higher population growth on new hosts, lower among-population variance in growth on new hosts, and lower population variance in growth across new hosts. We tested the predictions using experimentally evolved vesicular stomatitis virus populations, containing groups of host-use specialists, directly selected generalists, and indirectly selected generalists. In novel-host challenges, viruses directly selected for generalism showed relatively higher or equivalent host growth, lower among-population variance in host growth, and lower population variance in growth across hosts. Thus, two of three outcomes supported our prediction that directly selected host breadth should favor host colonization. Also, we observed that indirectly selected generalists were advantaged over specialist viruses, indicating that fortuitous changes in host breadth may also promote emergence. We discuss evolution of phenotypic plasticity versus environmental robustness in viruses, virus avoidance of extinction, and surveillance of pathogen niche breadth to predict future likelihood of emergence.

Tracking the evolution of multiple in vitro hepatitis C virus replicon variants under protease inhibitor selection pressure by 454 deep sequencing

Resistance to hepatitis C virus (HCV) inhibitors targeting viral enzymes has been observed in in vitro replicon studies and during clinical trials. The factors determining the emergence of resistance and the changes in the viral quasispecies population under selective pressure are not fully understood. To assess the dynamics of variants emerging in vitro under various selective pressures with TMC380765, a potent macrocyclic HCV NS3/4A protease inhibitor, HCV genotype 1b replicon-containing cells were cultured in the presence of a low, high, or stepwise-increasing TMC380765 concentration(s). HCV replicon RNA from representative samples thus obtained was analyzed using (i) population, (ii) clonal, and (iii) 454 deep sequencing technologies. Depending on the concentration of TMC380765, distinct mutational patterns emerged. In particular, culturing with low concentrations resulted in the selection of low-level resistance mutations (F43S and A156G), whereas high concentrations resulted in the selection of high-level resistance mutations (A156V, D168V, and D168A). Clonal and 454 deep sequencing analysis of the replicon RNA allowed the identification of low-frequency preexisting mutations possibly contributing to the mutational pattern that emerged. Stepwise-increasing TMC380765 concentrations resulted in the emergence and disappearance of multiple replicon variants in response to the changing selection pressure. Moreover, two different codons for the wild-type amino acids were observed at certain NS3 positions within one population of replicons, which may contribute to the emerging mutational patterns. Deep sequencing technologies enabled the study of minority variants present in the HCV quasispecies population present at baseline and during antiviral drug pressure, giving new insights into the dynamics of resistance acquisition by HCV.

Beneficial mutations and the dynamics of adaptation in asexual populations

We discuss the dynamics of adaptive evolution in asexual (clonal) populations. The classical 'periodic selection' model of clonal evolution assumed that beneficial mutations are very rare and therefore substitute unfettered into populations as occasional, isolated events. Newer models allow for the possibility that beneficial mutations are sufficiently common to coexist and compete for fixation within populations. Experimental evolution studies in microbes provide empirical support for stochastic models in which both selection and mutation are strong effects and clones compete for fixation; however, the relative importance of competition among clones bearing mutations of different selective effects versus competition among clones bearing multiple mutations remains unresolved. We provide some new theoretical results, moreover, suggesting that population dynamics consistent with the periodic selection model can arise even in a deterministic model that can accommodate a very high beneficial mutation rate.

Phenotypic differences in viral immune escape explained by linking within-host dynamics to host-population immunity

Viruses that do not cause life-long immunity persist by evolving rapidly in response to prevailing host immunity. The immune-escape mutants emerge frequently, displacing or co-circulating with native strains even though mutations conferring immune evasion are often detrimental to viral replication. The epidemiological dynamics of immune-escape in acute-infection viruses with high transmissibility have been interpreted mainly through immunity dynamics at the host population level, despite the fact that immune-escape evolution involves dynamical processes that feedback across the within- and between-host scales. To address this gap, we use a nested model of within- and between-host infection dynamics to examine how the interaction of viral replication rate and cross-immunity imprint host population immunity, which in turn determines viral immune escape. Our explicit consideration of direct and immune-mediated competitive interactions between strains within-hosts revealed three insights pertaining to risk and control of viral immune-escape: (1) replication rate and immune-stimulation deficiencies (i.e., original antigenic sin) act synergistically to increase immune escape, (2) immune-escape mutants with replication deficiencies relative to their wildtype progenitor are most successful under moderate cross-immunity and frequent re-infections, and (3) the immunity profile along short host-transmission chains (local host-network structure) is a key determinant of immune escape.

CRISPR-mediated phage resistance and the ghost of coevolution past

The past is never dead. It's not even past

William Faulkner (1951)

Bacteria can acquire heritable immunity to viral (phage) enemies by incorporating phage DNA into their own genome. This mechanism of anti-viral defence, known by the acronym CRISPR, simultaneously stores detailed information about current and past enemies and the evolved resistance to them. As a high-resolution genetic marker that is intimately tied with the host–pathogen interaction, the CRISPR system offers a unique, and relatively untapped, opportunity to study epidemiological and coevolutionary dynamics in microbial communities that were previously neglected because they could not be cultured in the laboratory. We briefly review the molecular mechanisms of CRISPR-mediated host–pathogen resistance, before assessing their potential importance for coevolution in nature, and their utility as a means of studying coevolutionary dynamics through metagenomics and laboratory experimentation.

Bacteriophages as model organisms for virus emergence research

Viruses fully emerge by gaining the ability to sustainably infect new host populations. When the hosts are humans, emerging viruses can present major public health issues, as exemplified by the AIDS pandemic. Therefore, heuristic approaches to identify nascent diseases before they become pandemic would be valuable. Unfortunately, the current patient-based and epidemiological approaches are ill-suited in this regard because they are largely responsive and not predictive. Alternative approaches based on virus evolutionary ecology might have greater potential to predict virus emergence. However, given the difficulties encountered when studying metazoan viruses in this context, the development of new model systems is greatly desirable. Here, I highlight studies that show that bacteriophages are appropriate model organisms for virus emergence research because of the ease in which important population parameters can be manipulated. Ideally this research will permit identifying major factors determining the persistence or extinction of emerging viruses. If such viruses could be recognized in advance, patient-based and epidemiological strategies could be better mobilized to deal with them.

Layers of evolvability in a bacteriophage life history trait

Functional redundancy in genomes arises from genes with overlapping functions, allowing phenotypes to persist after gene knockouts. Evolutionary redundancy or evolvability of a genome is one step removed, in that functional redundancy is absent but the genome has the potential to evolve to restore a lost phenotype. Exploring the extent to which this recovery alters gene networks can illuminate how functional gene interactions change through time. Here, the evolvability of lysis was studied in bacteriophage T7, revealing hidden functional interactions. Lysis is the destruction of host cell wall and membranes that releases progeny and is therefore essential for phage propagation. In most phages, lysis is mediated by a two-component genetic module: a muralytic enzyme that degrades the bacterial cell wall (endolysin) and a holin that permeabilizes the inner membrane to allow the endolysin access to the cell wall. T7 carries one known holin, one endolysin, and a second muralytic enzyme that plays little role in lysis by wild-type phage. If the primary endolysin is deleted, the second muralytic enzyme evolves to restore normal lysis after selection for faster growth. Here, a second level of evolutionary redundancy was revealed. When the second muralytic enzyme was prevented from adapting in a genome lacking the primary endolysin, the phage reevolved lysis de novo in the absence of any known muralytic enzymes by changes in multiple genes outside the original lysis module. This second level of redundancy proved to be evolutionarily inferior to the first, and both result in a lower fitness and slower lysis than wild-type T7. Deletion of the holin gene delayed lysis time modestly; fitness was restored by compensatory substitutions in genes that lack known roles in lysis of the wild type.

Growth of an RNA virus in single cells reveals a broad fitness distribution

Genetic and environmental factors will influence the growth of an RNA virus, but their relative contributions are challenging to resolve because standard culture methods mask how virus particles interact with individual host cells. Here, single particles of vesicular stomatitis virus, a prototype RNA virus, were used to infect individual BHK cells. Infected cells produced 50 to 8000 progeny virus particles, but these differences were lost upon subsequent culture, suggesting the diversity of yields reflected cell-to-cell differences rather than viral genetic variation. Cells infected at different phases of their cell cycle produced from 1400 (early S) to 8700 (G(2)M) infectious virus particles, coinciding with the middle-to-upper range of the observed distribution. Fluctuations in virus and cell compositions and noisy gene expression may also contribute to the broad distribution of virus yields. These findings take a step toward quantifying how environmental variation can impact the fitness distribution of an RNA virus.

The effects of population bottlenecks on clonal interference, and the adaptation effective population size

Clonal interference refers to the competition that arises in asexual populations when multiple beneficial mutations segregate simultaneously. A large body of theoretical and experimental work now addresses this issue. Although much of the experimental work is performed in populations that grow exponentially between periodic population bottlenecks, the theoretical work to date has addressed only populations of a constant size. We derive an analytical approximation for the rate of adaptation in the presence of both clonal interference and bottlenecks, and compare this prediction to the results of an individual-based simulation, showing excellent agreement in the parameter regime in which clonal interference prevails. We also derive an appropriate definition for the effective population size for adaptive evolution experiments in the presence of population bottlenecks. This "adaptation effective population size" allows for a good approximation of the expected rate of adaptation, either in the strong-selection weak-mutation regime, or when clonal interference comes into play. In the multiple mutation regime, when the product of the population size and mutation rate is extremely large, these results no longer hold.

Molecular characterization of clonal interference during adaptive evolution in asexual populations of Saccharomyces cerevisiae

The classical model of adaptive evolution in an asexual population postulates that each adaptive clone is derived from the one preceding it. However, experimental evidence has suggested more complex dynamics, with theory predicting the fixation probability of a beneficial mutation as dependent on the mutation rate, population size and the mutation's selection coefficient. Clonal interference has been demonstrated in viruses and bacteria but not in a eukaryote, and a detailed molecular characterization is lacking. Here we use three different fluorescent markers to visualize the dynamics of asexually evolving yeast populations. For each adaptive clone within one of our evolving populations, we identified the underlying mutations, monitored their population frequencies and used microarrays to characterize changes in the transcriptome. These results represent the most detailed molecular characterization of experimental evolution to date and provide direct experimental evidence supporting both the clonal interference and the multiple mutation models.

In vivo virus growth competition assays demonstrate equal fitness of fish rhabdovirus strains that co-circulate in aquaculture

A novel virus growth competition assay for determining relative fitness of RNA virus variants in vivo has been developed using the fish rhabdovirus, Infectious hematopoietic necrosis virus (IHNV), in juvenile rainbow trout (Oncorhynchus mykiss). We have conducted assays with IHNV isolates designated B, C, and D, representing the three most common genetic subtypes that co-circulate in Idaho trout farm aquaculture. In each assay, groups of 30 fish were immersed in a 1:1 mixture of two genotypes of IHNV, and then held in individual beakers for a 72h period of in vivo competitive virus replication. Progeny virus populations in each fish were analyzed for the presence and proportion of each viral genotype. In two independent assays of the B:C isolate pair, and two assays of the B:D pair, all fish were co-infected and there was a high level of fish-to-fish variation in the ratio of the two competing genotypes. However, in each assay the average ratio in the 30-fish group was not significantly different from the input ratio of 1:1, indicating equal or nearly equal viral fitness on a host population basis, under the conditions tested.

The evolution of foot-and-mouth disease virus: impacts of recombination and selection

Foot-and-mouth disease virus is an economically important animal virus that exhibits extensive genetic and antigenic heterogeneity. To examine the evolutionary forces that have influenced the population dynamics of foot-and-mouth disease virus, individual genes and the coding genomes for the Eurasian (Asia1, A, C, and O) serotypes were examined for phylogenetic relationships, recombination, genetic diversity and selection. Our analyses demonstrate that paraphyletic relationships among serotypes are not as prevalent as previously proposed and suggest that convergent evolution might be obscuring phylogenetic relationships. We provide evidence that identification of recombinant sequences and recombination breakpoint patterns among and within serotypes are heavily dependent on the level of genetic diversity and convergent characters present in a particular data set as well as the methods used to detect recombination. Here, we also investigate the impact of adaptive positive selection on the capsid proteins and the non-structural genes 2B, 2C, 3A, and 3Cpro to identify genome regions involved in genetic diversity and antigenic variation. Two different categories of positive selection at the amino acid level were examined; conservative (stabilizing) selection that maintains particular phenotypic properties of an amino acid residue and radical (destabilizing), and selection that dramatically alters the phenotype and potentially the functional and/or structural features of the protein. Approximately, 29% of residues in the capsid proteins were under positive selection. Of those, 64% were under the influence of destabilizing selection, 80% were under the influence of stabilizing selection, and 44% had phenotypic properties influenced by both selection types. The majority of residues under selection (74%) were located outside of known antigenic sites; suggestive of additional uncharacterized epitopes and genomic regions involved in antigenic drift.

Adaptation of Soybean mosaic virus avirulent chimeras containing P3 sequences from virulent strains to Rsv1-genotype soybeans is mediated by mutations in HC-Pro

In Rsv1-genotype soybean, Soybean mosaic virus (SMV)-N (an avirulent isolate of strain G2) elicits extreme resistance (ER) whereas strain SMV-G7 provokes a lethal systemic hypersensitive response (LSHR). SMV-G7d, an experimentally evolved variant of SMV-G7, induces systemic mosaic. Thus, for Rsv1-genotype soybean, SMV-N is avirulent whereas SMV-G7 and SMV-G7d are both virulent. Exploiting these differential interactions, we recently mapped the elicitor functions of SMV provoking Rsv1-mediated ER and LSHR to the N-terminal 271 amino acids of P3 from SMV-N and SMV-G7, respectively. The phenotype of both SMV-G7 and SMV-G7d were rendered avirulent on Rsv1-genotype soybean when the part of the genome encoding the N-terminus or the entire P3 cistron was replaced with that from SMV-N; however, reciprocal exchanges did not confer virulence to SMV-N-derived P3 chimeras. Here, we describe virulent SMV-N-derived P3 chimeras containing the full-length or the N-terminal P3 from SMV-G7 or SMV-G7d, with or without additional mutations in P3, that were selected on Rsv1-genotype soybean by sequential transfers on rsv1 and Rsv1-genotype soybean. Sequence analyses of the P3 and helper-component proteinase (HC-Pro) cistrons of progeny recovered from Rsv1-genotype soybean consistently revealed the presence of mutations in HC-Pro. Interestingly, the precise mutations in HC-Pro required for the adaptation varied among the chimeras. No mutation was detected in the HC-Pro of progeny passaged continuously in rsv1-genotype soybean, suggesting that selection is a consequence of pressure imposed by Rsv1. Mutations in HC-Pro alone failed to confer virulence to SMV-N; however, reconstruction of mutations in HC-Pro of the SMV-N-derived P3 chimeras resulted in virulence. Taken together, the data suggest that HC-Pro complementation of P3 is essential for SMV virulence on Rsv1-genotype soybean.

Evolutionary genomics of host adaptation in vesicular stomatitis virus

Populations experiencing similar selection pressures can sometimes diverge in the genetic architectures underlying evolved complex traits. We used RNA virus populations of large size and high mutation rate to study the impact of historical environment on genome evolution, thus increasing our ability to detect repeatable patterns in the evolution of genetic architecture. Experimental vesicular stomatitis virus populations were evolved on HeLa cells, on MDCK cells, or on alternating hosts. Turner and Elena (2000. Cost of host radiation in an RNA virus. Genetics. 156:1465-1470.) previously showed that virus populations evolved in single-host environments achieved high fitness on their selected hosts but failed to increase in fitness relative to their ancestor on the unselected host and that alternating-host-evolved populations had high fitness on both hosts. Here we determined the complete consensus sequence for each evolved population after 95 generations to gauge whether the parallel phenotypic changes were associated with parallel genomic changes. We also analyzed the patterns of allele substitutions to discern whether differences in fitness across hosts arose through true pleiotropy or the presence of not only a mutation that is beneficial in both hosts but also 1 or more mutations at other loci that are costly in the unselected environment (mutation accumulation [MA]). We found that ecological history may influence to what extent pleiotropy and MA contribute to fitness asymmetries across environments. We discuss the degree to which current genetic architecture is expected to constrain future evolution of complex traits, such as host use by RNA viruses.

Experimental evolution and genome sequencing reveal variation in levels of clonal interference in large populations of bacteriophage phiX174

Background

In large asexual populations where beneficial mutations may co-occur and recombination is absent, the fate of beneficial mutations can be significantly affected by competition (i.e., clonal interference). Theoretical models predict that clonal interference (CI) can slow adaptation, alter the distribution of fixed beneficial mutations, and affect disease progression by impacting within-host evolution of pathogens. While phenotypic data support that CI is a significant determinant of adaptive outcome, genetic data are needed to verify the patterns and to inform evolutionary models. We adapted replicate populations of the bacteriophage φX174 under two levels of CaCl₂ to create benign and harsh environments. Genomic sequences of multiple individuals from evolved populations were used to detect competing beneficial mutations.

Results

There were several competing genotypes in most of the populations where CaCl₂ was abundant, but no evidence of CI where CaCl₂ was scarce, even though rates of adaptation and population sizes among the treatments were similar. The sequence data revealed that observed mutations were limited to a single portion of one gene in the harsh treatment, but spanned a different and larger region of the genome under the benign treatments, suggesting that there were more adaptive solutions to the benign treatment.

Conclusion

Beneficial mutations with relatively large selection coefficients can be excluded by CI. CI may commonly determine the fate of beneficial mutations in large microbial populations, but its occurrence depends on selective conditions. CI was more frequent in benign selective conditions possibly due to a greater number of adaptive targets under this treatment. Additionally, the genomic sequence data showed that selection can target different types and numbers of phenotypes in environments that differ by only a single continuous variable.

Recombination speeds adaptation by reducing competition between beneficial mutations in populations of Escherichia coli

Identification of the selective forces contributing to the origin and maintenance of sex is a fundamental problem in biology. The Fisher–Muller model proposes that sex is advantageous because it allows beneficial mutations that arise in different lineages to recombine, thereby reducing clonal interference and speeding adaptation. I used the F plasmid to mediate recombination in the bacterium Escherichia coli and measured its effect on adaptation at high and low mutation rates. Recombination increased the rate of adaptation ∼3-fold more in the high mutation rate treatment, where beneficial mutations had to compete for fixation. Sequencing of candidate loci revealed the presence of a beneficial mutation in six high mutation rate lines. In the absence of recombination, this mutation took longer to fix and, over the course of its substitution, conferred a reduced competitive advantage, indicating interference between competing beneficial mutations. Together, these results provide experimental support for the Fisher–Muller model and demonstrate that plasmid-mediated gene transfer can accelerate bacterial adaptation.

Why have sex? One explanation is that sex is good because it allows beneficial mutations from different lineages to recombine. This reduces competition between mutations in a population and can increase the speed with which the population can adapt to environmental change. This explanation, known as the Fisher–Muller model, has extensive theoretical support; however, it is difficult to test experimentally. Using a simple microbial system I showed that recombination increased the rate of fitness improvement when beneficial mutations were common in the population and had to compete for fixation, but had little effect when mutations occurred rarely. Sequencing of candidate genes revealed the presence of the same beneficial mutation in a number of replicate populations. In the absence of recombination, this mutation took longer to spread and conferred a lower overall competitive advantage, indicating interference between competing beneficial mutations. Together, these results provide direct experimental support for the Fisher–Muller model.

Why there is sexual reproduction is an important question. Here the authors show that making bacteria have sex allows faster accumulation and fixation of beneficial mutations.

Evolutionary genetics: a piggyback ride to adaptation and diversity

Competition between adaptive mutations in different asexual lineages limits the rate of adaptation. But additional adaptive mutations can occur in lineages that already have one, altering the dynamics of evolving asexual populations.

Clonal interference in large populations

Clonal interference, the competition between lineages arising from different beneficial mutations in an asexually reproducing population, is an important factor determining the tempo and mode of microbial adaptation. The standard theory of this phenomenon neglects the occurrence of multiple mutations as well as the correlation between loss by genetic drift and clonal competition, which is questionable in large populations. Working within the Wright-Fisher model with multiplicative fitness (no epistasis), we determine the rate of adaptation asymptotically for very large population sizes and show that the standard theory fails in this regime. Our study also explains the success of the standard theory in predicting the rate of adaptation for moderately large populations. Furthermore, we show that the nature of the substitution process changes qualitatively when multiple mutations are allowed for, because several mutations can be fixed in a single fixation event. As a consequence, the index of dispersion for counts of the fixation process displays a minimum as a function of population size, whereas the origination process of fixed mutations becomes completely regular for very large populations. We find that the number of mutations fixed in a single event is geometrically distributed as in the neutral case. These conclusions are based on extensive simulations combined with analytic results for the limit of infinite population size.

Adaptive Mutations in Bacteria: High Rate and Small Effects

Evolution by natural selection is driven by the continuous generation of adaptive mutations. We measured the genomic mutation rate that generates beneficial mutations and their effects on fitness in Escherichia coli under conditions in which the effect of competition between lineages carrying different beneficial mutations is minimized. We found a rate on the order of 10(-5) per genome per generation, which is 1000 times as high as previous estimates, and a mean selective advantage of 1%. Such a high rate of adaptive evolution has implications for the evolution of antibiotic resistance and pathogenicity.

Beneficial mutation selection balance and the effect of linkage on positive selection

When beneficial mutations are rare, they accumulate by a series of selective sweeps. But when they are common, many beneficial mutations will occur before any can fix, so there will be many different mutant lineages in the population concurrently. In an asexual population, these different mutant lineages interfere and not all can fix simultaneously. In addition, further beneficial mutations can accumulate in mutant lineages while these are still a minority of the population. In this article, we analyze the dynamics of such multiple mutations and the interplay between multiple mutations and interference between clones. These result in substantial variation in fitness accumulating within a single asexual population. The amount of variation is determined by a balance between selection, which destroys variation, and beneficial mutations, which create more. The behavior depends in a subtle way on the population parameters: the population size, the beneficial mutation rate, and the distribution of the fitness increments of the potential beneficial mutations. The mutation-selection balance leads to a continually evolving population with a steady-state fitness variation. This variation increases logarithmically with both population size and mutation rate and sets the rate at which the population accumulates beneficial mutations, which thus also grows only logarithmically with population size and mutation rate. These results imply that mutator phenotypes are less effective in larger asexual populations. They also have consequences for the advantages (or disadvantages) of sex via the Fisher-Muller effect; these are discussed briefly.

Analysis of genetic systems using experimental evolution and whole-genome sequencing

Whole-genome sequencing of bacteria evolving in the laboratory promises to reveal the complex network of mutations that underlie adaptation.

The application of whole-genome sequencing to the study of microbial evolution promises to reveal the complex functional networks of mutations that underlie adaptation. A recent study of parallel evolution in populations of Escherichia coli shows how adaptation involves both functional changes to specific proteins as well as global changes in regulation.

Clonal interference is alleviated by high mutation rates in large populations

When a beneficial mutation is fixed in a population that lacks recombination, the genetic background linked to that mutation is fixed. As a result, beneficial mutations on different backgrounds experience competition, or "clonal interference," that can cause asexual populations to evolve more slowly than their sexual counterparts. Factors such as a large population size (N) and high mutation rates (mu) increase the number of competing beneficial mutations, and hence are expected to increase the intensity of clonal interference. However, recent theory suggests that, with very large values of Nmu, the severity of clonal interference may instead decline. The reason is that, with large Nmu, genomes including both beneficial mutations are rapidly created by recurrent mutation, obviating the need for recombination. Here, we analyze data from experimentally evolved asexual populations of a bacteriophage and find that, in these nonrecombining populations with very large Nmu, recurrent mutation does appear to ameliorate this cost of asexuality.

Adaptation of HIV-1 depends on the host-cell environment

Many viruses have the ability to rapidly develop resistance against antiviral drugs and escape from the host immune system. To which extent the host environment affects this adaptive potential of viruses is largely unknown. Here we show that for HIV-1, the host-cell environment is key to the adaptive potential of the virus. We performed a large-scale selection experiment with two HIV-1 strains in two different T-cell lines (MT4 and C8166). Over 110 days of culture, both virus strains adapted rapidly to the MT4 T-cell line. In contrast, when cultured on the C8166 T-cell line, the same strains did not show any increase in fitness. By sequence analyses and infections with viruses expressing either yellow or cyan fluorescent protein, we were able to show that the absence of adaptation was linked to a lower recombination rate in the C8166 T-cell line. Our findings suggest that if we can manipulate the host-cellular factors that mediate viral evolution, we may be able to significantly retard viral adaptability.

The speed of evolution and maintenance of variation in asexual populations

BACKGROUND

The rate at which beneficial mutations accumulate determines how fast asexual populations evolve, but this is only partially understood. Some recent clonal-interference models suggest that evolution in large asexual populations is limited because smaller beneficial mutations are outcompeted by larger beneficial mutations that occur in different lineages within the same population. This analysis assumes that the important mutations fix one at a time; it ignores multiple beneficial mutations that occur in the lineage of an earlier beneficial mutation, before the first mutation in the series can fix. We focus on the effects of such multiple mutations.

RESULTS

Our analysis predicts that the variation in fitness maintained by a continuously evolving population increases as the logarithm of the population size and logarithm of the mutation rate and thus yields a similar logarithmic increase in the speed of evolution. To test these predictions, we evolved asexual budding yeast in glucose-limited media at a range of population sizes and mutation rates.

CONCLUSIONS

We find that their evolution is dominated by the accumulation of multiple mutations of moderate effect. Our results agree with our theoretical predictions and are inconsistent with the one-by-one fixation of mutants assumed by recent clonal-interference analysis.

Differential impact of simultaneous migration on coevolving hosts and parasites

BACKGROUND

The dynamics of antagonistic host-parasite coevolution are believed to be crucially dependent on the rate of migration between populations. We addressed how the rate of simultaneous migration of host and parasite affected resistance and infectivity evolution of coevolving meta-populations of the bacterium Pseudomonas fluorescens and a viral parasite (bacteriophage). The increase in genetic variation resulting from small amounts of migration is expected to increase rates of adaptation of both host and parasite. However, previous studies suggest phages should benefit more from migration than bacteria; because in the absence of migration, phages are more genetically limited and have a lower evolutionary potential compared to the bacteria.

RESULTS

The results supported the hypothesis: migration increased the resistance of bacteria to their local (sympatric) hosts. Moreover, migration benefited phages more than hosts with respect to 'global' (measured with respect to the whole range of migration regimes) patterns of resistance and infectivity, because of the differential evolutionary responses of bacteria and phage to different migration regimes. Specifically, we found bacterial global resistance peaked at intermediate rates of migration, whereas phage global infectivity plateaued when migration rates were greater than zero.

CONCLUSION

These results suggest that simultaneous migration of hosts and parasites can dramatically affect the interaction of host and parasite. More specifically, the organism with the lower evolutionary potential may gain the greater evolutionary advantage from migration.

The effect of co- and superinfection on the adaptive dynamics of vesicular stomatitis virus

In many infectious diseases, hosts are often simultaneously infected with several genotypes of the same pathogen. Much theoretical work has been done on modelling multiple infection dynamics, but empirical evidences are relatively scarce. Previous studies have demonstrated that coinfection allows faster adaptation than single infection in RNA viruses. Here, we use experimental populations of the vesicular stomatitis Indiana virus derived from an infectious cDNA, to show that superinfection dynamics promotes faster adaptation than single infection. In addition, we have analysed two different periodicities of multiple infection, daily and separated 5 days in time. Daily multiple infections allow higher fitness increases than multiple infections taking place every 5 days. We propose that the effect of superinfection on fitness is mainly influenced by the time elapsed between the first and the second infection, since shorter time intervals offer more opportunities to competition between resident and invading populations.

The distribution of beneficial mutant effects under strong selection

For a general theory of adaptation, it is essential to know the distribution of fitness effects of beneficial mutations. Recent theoretical and empirical studies have made considerable progress in determining the characteristics of this distribution. To date, the experiments have largely verified the theoretical predictions. Despite the fact that the theoretical work has assumed small selection coefficients, strong selection has been observed in some experiments, especially those involving novel environments. Here, we derive the distribution of fitness effects among fixed beneficial mutants without the restriction of low selection coefficients. The fate of strongly favored alleles is less affected by stochastic drift while rare, causing the distribution of fitness effects among fixed beneficial mutations to reflect more closely the distribution among all newly arising beneficial mutations. We also find that when many alleles compete for fixation within an asexual population (clonal interference), the beneficial effects of a newly fixed mutant cannot be well estimated because of the high number of subsequent mutations that arise within the genome, regardless of whether selection is strong or weak.

The effect of population structure on the adaptive radiation of microbial populations evolving in spatially structured environments

Spatial structure is thought to be an important factor influencing the emergence and maintenance of genetic diversity. Previous studies have demonstrated that environmental heterogeneity, provided by spatial structure, leads to adaptive radiation of populations. In the present study, we investigate not only the impact of environmental heterogeneity on adaptive radiation, but also of population fragmentation and niche construction. Replicate populations founded by a single genotype of Escherichia coli were allowed to evolve for 900 generations by serial transfer in either a homogeneous environment, or a spatially structured environment that was either kept intact or destroyed with each daily transfer. Only populations evolving in the structured environment with intact population structure diversified: clones are significantly divergent in sugar catabolism, and show frequency-dependent fitness interactions indicative of stable coexistence. These findings demonstrate an important role for population fragmentation, a consequence of population structure in spatially structured environments, on the diversification of populations.

Effects of random mutations in the human immunodeficiency virus type 1 transcriptional promoter on viral fitness in different host cell environments

A mutation's effect on fitness or phenotype may in part depend on the interaction of the mutation with the environment. The resulting phenotype or fitness is important, since it determines the adaptive potential of a species. To date, most studies have focused on alterations to protein-coding regions of the genome and their consequential fitness effects. Non-protein-coding regulatory regions have been largely neglected, although they make up a large and important part of an organism's genome. Here, we use human immunodeficiency virus type 1 as a model system to investigate fitness effects of random mutations in noncoding DNA-binding sites of the transcriptional promoter. We determined 242 fitness values for 35 viral promoter mutants with one, two, or three mutations across seven distinct cellular environments and identified that (i) all mutants have an effect in at least one cellular environment; (ii) fitness effects are highly dependent on the cellular environment; (iii) disadvantageous and advantageous mutations occur at high and similar frequencies; and (iv) epistatic effects of multiple mutations are rare. Our results underline the evolutionary potential of regulatory regions and indicate that DNA-binding sites evolve under strong selection, while at the same time, they are very plastic to environmental change.

Virus fitness: concept, quantification, and application to HIV population dynamics

Viral fitness has been broadly studied during the past three decades, mainly to test evolutionary models and population theories difficult to analyze and interpret with more complex organisms. More recent studies, however, are focused in the role of fitness on viral transmission, pathogenesis, and drug resistance. Here, we used human immunodeficiency virus (HIV) as one of the most relevant models to evaluate the importance of viral quasispecies and fitness in HIV evolution, population dynamics, disease progression, and potential clinical implications.

Clonal interference and the periodic selection of new beneficial mutations in Escherichia coli

The conventional model of adaptation in asexual populations implies sequential fixation of new beneficial mutations via rare selective sweeps that purge all variation and preserve the clonal genotype. However, in large populations multiple beneficial mutations may co-occur, causing competition among them, a phenomenon called "clonal interference." Clonal interference is thus expected to lead to longer fixation times and larger fitness effects of mutations that ultimately become fixed, as well as to a genetically more diverse population. Here, we study the significance of clonal interference in populations consisting of mixtures of differently marked wild-type and mutator strains of Escherichia coli that adapt to a minimal-glucose environment for 400 generations. We monitored marker frequencies during evolution and measured the competitive fitness of random clones from each marker state after evolution. The results demonstrate the presence of multiple beneficial mutations in these populations and slower and more erratic invasion of mutants than expected by the conventional model, showing the signature of clonal interference. We found that a consequence of clonal interference is that fitness estimates derived from invasion trajectories were less than half the magnitude of direct estimates from competition experiments, thus revealing fundamental problems with this fitness measure. These results force a reevaluation of the conventional model of periodic selection for asexual microbes.

An equivalence principle for the incorporation of favorable mutations in asexual populations

Rapid evolution of asexual populations, such as that of cancer cells or of microorganisms developing drug resistance, can include the simultaneous spread of distinct beneficial mutations. We demonstrate that evolution in such cases is driven by the fitness effects and appearance times of only a small minority of favorable mutations. The complexity of the mutation-selection process is thereby greatly reduced, and much of the evolutionary dynamics can be encapsulated in two parameters-an effective selection coefficient and effective rate of beneficial mutations. We confirm this theoretical finding and estimate the effective parameters for evolving populations of fluorescently labeled Escherichia coli. The effective parameters constitute a simple description and provide a natural standard for comparing adaptation between species and across environments.

Future issues in RNA virus evolution

The study of RNA virus evolution has developed rapidly during the past 30 years. This review outlines some important recent findings, as well as a number of the remaining major challenges, particularly those that might explain why RNA viruses are the most important class of emerging diseases, yet often have difficulties adapting to sustained transmission cycles in new hosts. The author emphasizes the relevance of research on the underlying dynamics of mutation, fitness landscapes and the constraints to viral adaptation, as well as the evolution of recombination and reassortment. It is also suggested that a combination of theoretical, experimental and comparative approaches is essential for future studies of viral evolution, coupled with new genome sequence data on intrahost genetic variation.

Waiting times for the appearance of cytotoxic T-lymphocyte escape mutants in chronic HIV-1 infection

The failure of HIV-1 to escape at some cytotoxic T-lymphocyte (CTL) epitopes has generally been explained in terms of viral fitness costs or ineffective or attenuated CTL responses. Relatively little attention has been paid to the evolutionary time required for escape mutants to be detected. This time is significantly affected by selection, mutation rates, the presence of other advantageous mutations, and the effective population size of HIV-1 in vivo (typically estimated to be approximately 10(3) in chronically infected patients, though one study has estimated it to be approximately 10(5)). Here, we use a forward simulator with experimentally estimated HIV-1 parameters to show that these delays can be substantial. For an effective population size of 10(3), even highly advantageous mutants (s = 0.5) may not be detected for a couple of years in chronically infected patients, while moderately advantageous escape mutants (s = 0.1) may not be detected for up to 10 years. Even with an effective population size of 10(5), a moderately advantageous escape mutant (s = 0.1) may not be detected in the population within 2 years if it has to compete with other selectively advantageous mutants. Stochastic evolutionary forces, therefore, in addition to viral fitness costs and ineffective or attenuated CTL responses, must be taken into account when assessing the selection of CTL escape mutations.

Genetic drift of parvovirus B19 is found in AIDS patients with persistent B19 infection

It is generally thought that parvovirus B19 is stable genetically. Consistently, genetic drift has not been found in patients with persistent B19 infection. In this report, longitudinal genetic changes in NS1 and VP1 gene of B19 isolates from three AIDS patients with persistent B19 infection were studied. One of the three patients was not treated with highly active anti-retroviral therapy (HAART). B19 viral DNA from these patients was amplified by polymerase chain reaction (PCR) and then sequenced directly. A single genetic change was found in the B19 isolate obtained from the patient not treated with HAART on Day 10 after intravenous immunoglobulin (IVIG) treatment. The nucleotide sequences of B19 isolated from this patient, then remained unchanged over a period of 11 months. Analysis of NS1 clones derived from his longitudinal viral isolates showed the existence of quasi-species but genetic drift was not found. One of the other two patients treated with HAART experienced treatment failure; he was later treated with mega-HAART. In contrast to the genetic stability of B19 isolates from the patient not treated with HAART, multiple genetic changes were discovered in the viral isolates from the two other patients after HAART and mega-HAART, respectively. Through analysis of B19 clones, the frequency of clones containing these mutations confirmed the genetic drift. Nucleotide substitutions seen in VP2 gene of isolates with genetic drift from both patients were all non-conserved, suggesting that they are positively selected.

Viruses as quasispecies: biological implications

During viral infections, the complex and dynamic distributions of variants, termed viral quasispecies, play a key role in the adaptability of viruses to changing environments and the fate of the population as a whole. Mutant spectra are continuously and avoidably generated during RNA genome replication, and they are not just a by-product of error-prone replication, devoid of biological relevance. On the contrary, current evidence indicates that mutant spectra contribute to viral pathogenesis, can modulate the expression of phenotypic traits by subpopulations of viruses, can include memory genomes that reflect the past evolutionary history of the viral lineage, and, furthermore, can participate in viral extinction through lethal mutagenesis. Also, mutant spectra are the target on which selection and random drift act to shape the long-term evolution of viruses. The biological relevance of mutant spectra is the central topic of this chapter.

Limits to adaptation in asexual populations

In asexual populations, the rate of adaptation is basically limited by the frequency and properties of spontaneous beneficial mutations. Hence, knowledge of these mutational properties and how they are affected by particular evolutionary conditions is a precondition for understanding the process of adaptation. Here, we address how the rate of adaptation of asexual populations is limited by its population size and mutation rate, as well as by two factors affecting the fraction of mutations that confer a benefit, i.e. the initial adaptedness of the population and the variability of the environment. These factors both influence which mutations are likely to occur, as well as the probability that they will ultimately contribute to adaptation. We attempt to separate the consequences of these basic population features in terms of their effect on the rate of adaptation by using results from evolution experiments with microorganisms.

Experimental evolution of Pseudomonas fluorescens in simple and complex environments

In complex environments that contain several substitutable resources, lineages may become specialized to consume only one or a few of them. Here we investigate the importance of environmental complexity in determining the evolution of niche width over approximately 900 generations in a chemically defined experimental system. We propagated 120 replicate lines of the bacterium Pseudomonas fluorescens in environments of different complexity by using between one and eight carbon substrates in each environment. Genotypes from populations selected in complex environments evolved greater mean and variance in fitness than those from populations selected in simple environments. Thus, lineages were able to adapt to several substrates simultaneously without any appreciable loss of function with respect to other substrates present in the media. There was greater genetic and genotype-by-environment interaction variance for fitness within populations selected in complex environments. It is likely that genetic variance in populations grown on complex media was maintained because the identity of the fittest genotype varied among carbon substrates. Our results suggest that evolution in complex environments will result neither in narrow specialists nor in complete generalists but instead in overlapping imperfect generalists, each of which has become adapted to a certain range of substrates but not to all.

Adaptation in sexuals vs. asexuals: clonal interference and the Fisher-Muller model

Fisher and Muller's theory that recombination speeds adaptation by eliminating competition among beneficial mutations has proved a popular explanation for the advantage of sex. Recent theoretical studies have attempted to quantify the speed of adaptation under the Fisher-Muller model, partly in an attempt to understand the role of "clonal interference" in microbial experimental evolution. We reexamine adaptation in sexuals vs. asexuals, using a model of DNA sequence evolution. In this model, a modest number of sites can mutate to beneficial alleles and the fitness effects of these mutations are unequal. We study (1) transition probabilities to different beneficial mutations; (2) waiting times to the first and the last substitutions of beneficial mutations; and (3) trajectories of mean fitness through time. We find that some of these statistics are surprisingly similar between sexuals and asexuals. These results highlight the importance of the choice of substitution model in assessing the Fisher-Muller advantage of sex.

The speed of adaptation in large asexual populations

In large asexual populations, beneficial mutations have to compete with each other for fixation. Here, I derive explicit analytic expressions for the rate of substitution and the mean beneficial effect of fixed mutations, under the assumptions that the population size N is large, that the mean effect of new beneficial mutations is smaller than the mean effect of new deleterious mutations, and that new beneficial mutations are exponentially distributed. As N increases, the rate of substitution approaches a constant, which is equal to the mean effect of new beneficial mutations. The mean effect of fixed mutations continues to grow logarithmically with N. The speed of adaptation, measured as the change of log fitness over time, also grows logarithmically with N for moderately large N, and it grows double-logarithmically for extremely large N. Moreover, I derive a simple formula that determines whether at given N beneficial mutations are expected to compete with each other or go to fixation independently. Finally, I verify all results with numerical simulations.