Exceptional convergent evolution in a virus

Genomics of adaptation during experimental evolution of the opportunistic pathogen Pseudomonas aeruginosa

Adaptation is likely to be an important determinant of the success of many pathogens, for example when colonizing a new host species, when challenged by antibiotic treatment, or in governing the establishment and progress of long-term chronic infection. Yet, the genomic basis of adaptation is poorly understood in general, and for pathogens in particular. We investigated the genetics of adaptation to cystic fibrosis-like culture conditions in the presence and absence of fluoroquinolone antibiotics using the opportunistic pathogen Pseudomonas aeruginosa. Whole-genome sequencing of experimentally evolved isolates revealed parallel evolution at a handful of known antibiotic resistance genes. While the level of antibiotic resistance was largely determined by these known resistance genes, the costs of resistance were instead attributable to a number of mutations that were specific to individual experimental isolates. Notably, stereotypical quinolone resistance mutations in DNA gyrase often co-occurred with other mutations that, together, conferred high levels of resistance but no consistent cost of resistance. This result may explain why these mutations are so prevalent in clinical quinolone-resistant isolates. In addition, genes involved in cyclic-di-GMP signalling were repeatedly mutated in populations evolved in viscous culture media, suggesting a shared mechanism of adaptation to this CF–like growth environment. Experimental evolutionary approaches to understanding pathogen adaptation should provide an important complement to studies of the evolution of clinical isolates.

Pathogens face a hostile and often novel environment when infecting a new host, and adaptation to this environment can be critical to a pathogen's survival. The genetic basis of pathogen adaptation is in turn important for treatment, since the consistency with which therapies succeed may depend on the extent to which a pathogen adapts via the same routes in different patients. In this study, we investigate adaptation of the bacterium Pseudomonas aeruginosa to laboratory conditions that resemble the lungs of cystic fibrosis patients and to quinolone antibiotics. We find that a handful of genes and genetic pathways are repeatedly involved in adaptation to each condition. Nonetheless, other, less common mutations can play important roles in determining fitness, complicating strategies aimed at reducing the prevalence of antibiotic resistance.

Cancer in light of experimental evolution

Cancer initiation, progression, and the emergence of therapeutic resistance are evolutionary phenomena of clonal somatic cell populations. Studies in microbial experimental evolution and the theoretical work inspired by such studies are yielding deep insights into the evolutionary dynamics of clonal populations, yet there has been little explicit consideration of the relevance of this rapidly growing field to cancer biology. Here, we examine how the understanding of mutation, selection, and spatial structure in clonal populations that is emerging from experimental evolution may be applicable to cancer. Along the way, we discuss some significant ways in which cancer differs from the model systems used in experimental evolution. Despite these differences, we argue that enhanced prediction and control of cancer may be possible using ideas developed in the context of experimental evolution, and we point out some prospects for future research at the interface between these traditionally separate areas.

The consistency of beneficial fitness effects of mutations across diverse genetic backgrounds

Parallel and convergent evolution have been remarkably common observations in molecular adaptation but primarily in the context of the same genotype adapting to the same conditions. These phenomena therefore tell us about the stochasticity and limitations of adaptation. The limited data on convergence and parallelism in the adaptation of different genotypes conflict as to the importance of such events. If the effects of beneficial mutations are highly context dependent (i.e., if they are epistatic), different genotypes should adapt through different mutations. Epistasis for beneficial mutations has been investigated but mainly through measurement of interactions between individually beneficial mutations for the same genotype. We examine epistasis for beneficial mutations at a broader genetic scale by measuring the fitness effects of two mutations beneficial for the ssDNA bacteriophage ID11 in eight different, related genotypes showing 0.3-3.7% nucleotide divergence from ID11. We found no evidence for sign epistasis, but the mutations tended to have much smaller or no effects on fitness in the new genotypes. We found evidence for diminishing-returns epistasis; the effects were more beneficial for lower-fitness genotypes. The patterns of epistasis were not determined by phylogenetic relationships to the original genotype. To improve our understanding of the patterns of epistasis, we fit the data to a model in which each mutation had a constant, nonepistatic phenotypic effect across genotypes and the phenotype-fitness map had a single optimum. This model fit the data well, suggesting that epistasis for these mutations was due to nonlinearity in the phenotype-fitness mapping and that the likelihood of parallel evolution depends more on phenotype than on genotype.

Effects of adaptation, chance, and history on the evolution of the toxic dinoflagellate Alexandrium minutum under selection of increased temperature and acidification

The roles of adaptation, chance, and history on evolution of the toxic dinoflagellate Alexandrium minutum Halim, under selective conditions simulating global change, have been addressed. Two toxic strains (AL1V and AL2V), previously acclimated for two years at pH 8.0 and 20°C, were transferred to selective conditions: pH 7.5 to simulate acidification and 25°C. Cultures under selective conditions were propagated until growth rate and toxin cell quota achieved an invariant mean value at 720 days (ca. 250 and ca. 180 generations for strains AL1V and AL2V, respectively). Historical contingencies strongly constrained the evolution of growth rate and toxin cell quota, but the forces involved in the evolution were not the same for both traits. Growth rate was 1.5–1.6 times higher than the one measured in ancestral conditions. Genetic adaptation explained two-thirds of total adaptation while one-third was a consequence of physiological adaptation. On the other hand, the evolution of toxin cell quota showed a pattern attributable to neutral mutations because the final variances were significantly higher than those measured at the start of the experiment. It has been hypothesized that harmful algal blooms will increase under the future scenario of global change. Although this study might be considered an oversimplification of the reality, it can be hypothesized that toxic blooms will increase but no predictions can be advanced about toxicity.

Evolutionary genomics of host-use in bifurcating demes of RNA virus phi-6

BACKGROUND

Viruses are exceedingly diverse in their evolved strategies to manipulate hosts for viral replication. However, despite these differences, most virus populations will occasionally experience two commonly-encountered challenges: growth in variable host environments, and growth under fluctuating population sizes. We used the segmented RNA bacteriophage ϕ6 as a model for studying the evolutionary genomics of virus adaptation in the face of host switches and parametrically varying population sizes. To do so, we created a bifurcating deme structure that reflected lineage splitting in natural populations, allowing us to test whether phylogenetic algorithms could accurately resolve this 'known phylogeny'. The resulting tree yielded 32 clones at the tips and internal nodes; these strains were fully sequenced and measured for phenotypic changes in selected traits (fitness on original and novel hosts).

RESULTS

We observed that RNA segment size was negatively correlated with the extent of molecular change in the imposed treatments; molecular substitutions tended to cluster on the Small and Medium RNA chromosomes of the virus, and not on the Large segment. Our study yielded a very large molecular and phenotypic dataset, fostering possible inferences on genotype-phenotype associations. Using further experimental evolution, we confirmed an inference on the unanticipated role of an allelic switch in a viral assembly protein, which governed viral performance across host environments.

CONCLUSIONS

Our study demonstrated that varying complexities can be simultaneously incorporated into experimental evolution, to examine the combined effects of population size, and adaptation in novel environments. The imposed bifurcating structure revealed that some methods for phylogenetic reconstruction failed to resolve the true phylogeny, owing to a paucity of molecular substitutions separating the RNA viruses that evolved in our study.

Experimental evolution

Experimental evolution is the study of evolutionary processes occurring in experimental populations in response to conditions imposed by the experimenter. This research approach is increasingly used to study adaptation, estimate evolutionary parameters, and test diverse evolutionary hypotheses. Long applied in vaccine development, experimental evolution also finds new applications in biotechnology. Recent technological developments provide a path towards detailed understanding of the genomic and molecular basis of experimental evolutionary change, while new findings raise new questions that can be addressed with this approach. However, experimental evolution has important limitations, and the interpretation of results is subject to caveats resulting from small population sizes, limited timescales, the simplified nature of laboratory environments, and, in some cases, the potential to misinterpret the selective forces and other processes at work.

Comparative analysis of two phenotypically-similar but genomically-distinct Burkholderia cenocepacia-specific bacteriophages

Background

Genomic analysis of bacteriophages infecting the Burkholderia cepacia complex (BCC) is an important preliminary step in the development of a phage therapy protocol for these opportunistic pathogens. The objective of this study was to characterize KL1 (vB_BceS_KL1) and AH2 (vB_BceS_AH2), two novel Burkholderia cenocepacia-specific siphoviruses isolated from environmental samples.

Results

KL1 and AH2 exhibit several unique phenotypic similarities: they infect the same B. cenocepacia strains, they require prolonged incubation at 30°C for the formation of plaques at low titres, and they do not form plaques at similar titres following incubation at 37°C. However, despite these similarities, we have determined using whole-genome pyrosequencing that these phages show minimal relatedness to one another. The KL1 genome is 42,832 base pairs (bp) in length and is most closely related to Pseudomonas phage 73 (PA73). In contrast, the AH2 genome is 58,065 bp in length and is most closely related to Burkholderia phage BcepNazgul. Using both BLASTP and HHpred analysis, we have identified and analyzed the putative virion morphogenesis, lysis, DNA binding, and MazG proteins of these two phages. Notably, MazG homologs identified in cyanophages have been predicted to facilitate infection of stationary phase cells and may contribute to the unique plaque phenotype of KL1 and AH2.

Conclusions

The nearly indistinguishable phenotypes but distinct genomes of KL1 and AH2 provide further evidence of both vast diversity and convergent evolution in the BCC-specific phage population.

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.

Differential transcription of bacteriophage φX174 genes at 37 °C and 42 °C

To investigate how high temperature affects viral transcription, the absolute amounts of mRNA for six bacteriophage φX174 genes were compared at 37°C and 42°C using Q-PCR. At 37°C, mRNA levels for all genes were consistent with previous studies, but at 42°C mRNA levels for four genes were significantly different from levels at 37°C. Transcript levels were higher for genes B and D; the promoter before gene B appears to be up-regulated at high temperature. Levels for genes F and G were reduced at high temperature, possibly due to increased efficiency of the transcription termination signal immediately upstream of gene F. These functional changes in φX174 gene regulation at high temperature have not been described previously. Studies of phage evolution at high temperatures indicate that this difference in transcript levels is subject to adaptation.

Real time forecasting of near-future evolution

A metaphor for adaptation that informs much evolutionary thinking today is that of mountain climbing, where horizontal displacement represents change in genotype, and vertical displacement represents change in fitness. If it were known a priori what the 'fitness landscape' looked like, that is, how the myriad possible genotypes mapped onto fitness, then the possible paths up the fitness mountain could each be assigned a probability, thus providing a dynamical theory with long-term predictive power. Such detailed genotype-fitness data, however, are rarely available and are subject to change with each change in the organism or in the environment. Here, we take a very different approach that depends only on fitness or phenotype-fitness data obtained in real time and requires no a priori information about the fitness landscape. Our general statistical model of adaptive evolution builds on classical theory and gives reasonable predictions of fitness and phenotype evolution many generations into the future.

The evolution of endogenous viral elements

Endogenous retroviruses are a common component of the eukaryotic genome, and their evolution and potential function have attracted considerable interest. More surprising was the recent discovery that eukaryotic genomes contain sequences from RNA viruses that have no DNA stage in their life cycle. Similarly, several single-stranded DNA viruses have left integrated copies in their host genomes. This review explores some major evolutionary aspects arising from the discovery of these endogenous viral elements (EVEs). In particular, the reasons for the bias toward EVEs derived from negative-sense RNA viruses are considered, as well as what they tell us about the long-term "arms races" between hosts and viruses, characterized by episodes of selection and counter-selection. Most dramatically, the presence of orthologous EVEs in divergent hosts demonstrates that some viral families have ancestries dating back almost 100 million years, and hence are far older than expected from the phylogenetic analysis of their exogenous relatives.

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.

Impact of host cell line adaptation on quasispecies composition and glycosylation of influenza A virus hemagglutinin

The genome of influenza A viruses is constantly changing (genetic drift) resulting in small, gradual changes in viral proteins. Alterations within antibody recognition sites of the viral membrane glycoproteins hemagglutinin (HA) and neuraminidase (NA) result in an antigenetic drift, which requires the seasonal update of human influenza virus vaccines. Generally, virus adaptation is necessary to obtain sufficiently high virus yields in cell culture-derived vaccine manufacturing. In this study detailed HA N-glycosylation pattern analysis was combined with in-depth pyrosequencing analysis of the virus genomic RNA. Forward and backward adaptation from Madin-Darby Canine Kidney (MDCK) cells to African green monkey kidney (Vero) cells was investigated for two closely related influenza A virus PR/8/34 (H1N1) strains: from the National Institute for Biological Standards and Control (NIBSC) or the Robert Koch Institute (RKI). Furthermore, stability of HA N-glycosylation patterns over ten consecutive passages and different harvest time points is demonstrated. Adaptation to Vero cells finally allowed efficient influenza A virus replication in Vero cells. In contrast, during back-adaptation the virus replicated well from the very beginning. HA N-glycosylation patterns were cell line dependent and stabilized fast within one (NIBSC-derived virus) or two (RKI-derived virus) successive passages during adaptation processes. However, during adaptation new virus variants were detected. These variants carried "rescue" mutations on the genomic level within the HA stem region, which result in amino acid substitutions. These substitutions finally allowed sufficient virus replication in the new host system. According to adaptation pressure the composition of the virus populations varied. In Vero cells a selection for "rescue" variants was characteristic. After back-adaptation to MDCK cells some variants persisted at indifferent frequencies, others slowly diminished and even dropped below the detection limit.

Stickbreaking: a novel fitness landscape model that harbors epistasis and is consistent with commonly observed patterns of adaptive evolution

In relating genotypes to fitness, models of adaptation need to both be computationally tractable and qualitatively match observed data. One reason that tractability is not a trivial problem comes from a combinatoric problem whereby no matter in what order a set of mutations occurs, it must yield the same fitness. We refer to this as the bookkeeping problem. Because of their commutative property, the simple additive and multiplicative models naturally solve the bookkeeping problem. However, the fitness trajectories and epistatic patterns they predict are inconsistent with the patterns commonly observed in experimental evolution. This motivates us to propose a new and equally simple model that we call stickbreaking. Under the stickbreaking model, the intrinsic fitness effects of mutations scale by the distance of the current background to a hypothesized boundary. We use simulations and theoretical analyses to explore the basic properties of the stickbreaking model such as fitness trajectories, the distribution of fitness achieved, and epistasis. Stickbreaking is compared to the additive and multiplicative models. We conclude that the stickbreaking model is qualitatively consistent with several commonly observed patterns of adaptive evolution.

Tempo and mode of plant RNA virus escape from RNA interference-mediated resistance

A biotechnological application of artificial microRNAs (amiRs) is the generation of plants that are resistant to virus infection. This resistance has proven to be highly effective and sequence specific. However, before these transgenic plants can be deployed in the field, it is important to evaluate the likelihood of the emergence of resistance-breaking mutants. Two issues are of particular interest: (i) whether such mutants can arise in nontransgenic plants that may act as reservoirs and (ii) whether a suboptimal expression level of the transgene, resulting in subinhibitory concentrations of the amiR, would favor the emergence of escape mutants. To address the first issue, we experimentally evolved independent lineages of Turnip mosaic virus (TuMV) (family Potyviridae) in fully susceptible wild-type Arabidopsis thaliana plants and then simulated the spillover of the evolving virus to fully resistant A. thaliana transgenic plants. To address the second issue, the evolution phase took place with transgenic plants that expressed the amiR at subinhibitory concentrations. Our results show that TuMV populations replicating in susceptible hosts accumulated resistance-breaking alleles that resulted in the overcoming of the resistance of fully resistant plants. The rate at which resistance was broken was 7 times higher for TuMV populations that experienced subinhibitory concentrations of the antiviral amiR. A molecular characterization of escape alleles showed that they all contained at least one nucleotide substitution in the target sequence, generally a transition of the G-to-A and C-to-U types, with many instances of convergent molecular evolution. To better understand the viral population dynamics taking place within each host, as well as to evaluate relevant population genetic parameters, we performed in silico simulations of the experiments. Together, our results contribute to the rational management of amiR-based antiviral resistance in plants.

Evolutionary changes in growth rate and toxin production in the cyanobacterium Microcystis aeruginosa under a scenario of eutrophication and temperature increase

Toxic blooms of the cyanobacterium Microcystis aeruginosa affect humans and animals in inland water systems worldwide, and it has been hypothesized that the development of these blooms will increase under the future scenario of global change, considering eutrophication and temperature increase as two important consequences. The importance of genetic adaptation, chance and history on evolution of growth rate, and toxin production of M. aeruginosa was studied under these new conditions. The experiment followed the idea of "replaying life's tape" by means of the simultaneous propagation of 15 independent isolates of three M. aeruginosa strains, which were grown under doubled nutrient concentration and temperature during c. 87 generations. Adaptation by new mutations that resulted in the enhancement of growth rate arose during propagation of derived cultures under the new environmental conditions was the main component of evolution; however, chance also contributed in a lesser extension to evolution of growth rate. Mutations were selected, displacing the wild-type ancestral genotypes. In contrast, the effect of selection on mutations affecting microcystin production was neutral. Chance and history were the pacemakers in evolution of toxin production. Although this study might be considered an oversimplification of the reality, it suggest that a future scenario of global change might lead to an increase in M. aeruginosa bloom frequency, but no predictions about the frequency of toxicity can be made.

The fitness effects of synonymous mutations in DNA and RNA viruses

Despite being silent with respect to protein sequence, synonymous nucleotide substitutions can be targeted by natural selection directly at the DNA or RNA level. However, there has been no systematic assessment of how frequent this type of selection is. Here, we have constructed 53 single random synonymous substitution mutants of the bacteriophages Qβ and ΦX174 by site-directed mutagenesis and assayed their fitness. Analysis of this mutant collection and of previous studies undertaken with a variety of single-stranded (ss) viruses demonstrates that selection at synonymous sites is stronger in RNA viruses than in DNA viruses. We estimate that this type of selection contributes approximately 18% of the overall mutational fitness effects in ssRNA viruses under our assay conditions and that random synonymous substitutions have a 5% chance of being lethal to the virus, whereas in ssDNA viruses, these figures drop to 1.4% and 0%, respectively. In contrast, the effects of nonsynonymous substitutions appear to be similar in ssRNA and ssDNA viruses.

Genomic and protein structural maps of adaptive evolution of human influenza A virus to increased virulence in the mouse

Adaptive evolution is characterized by positive and parallel, or repeated selection of mutations. Mouse adaptation of influenza A virus (IAV) produces virulent mutants that demonstrate positive and parallel evolution of mutations in the hemagglutinin (HA) receptor and non-structural protein 1 (NS1) interferon antagonist genes. We now present a genomic analysis of all 11 genes of 39 mouse adapted IAV variants from 10 replicate adaptation experiments. Mutations were mapped on the primary and structural maps of each protein and specific mutations were validated with respect to virulence, replication, and RNA polymerase activity. Mouse adapted (MA) variants obtained after 12 or 20–21 serial infections acquired on average 5.8 and 7.9 nonsynonymous mutations per genome of 11 genes, respectively. Among a total of 115 nonsynonymous mutations, 51 demonstrated properties of natural selection including 27 parallel mutations. The greatest degree of parallel evolution occurred in the HA receptor and ribonucleocapsid components, polymerase subunits (PB1, PB2, PA) and NP. Mutations occurred in host nuclear trafficking factor binding sites as well as sites of virus-virus protein subunit interaction for NP, NS1, HA and NA proteins. Adaptive regions included cap binding and endonuclease domains in the PB2 and PA polymerase subunits. Four mutations in NS1 resulted in loss of binding to the host cleavage and polyadenylation specificity factor (CPSF30) suggesting that a reduction in inhibition of host gene expression was being selected. The most prevalent mutations in PB2 and NP were shown to increase virulence but differed in their ability to enhance replication and demonstrated epistatic effects. Several positively selected RNA polymerase mutations demonstrated increased virulence associated with >300% enhanced polymerase activity. Adaptive mutations that control host range and virulence were identified by their repeated selection to comprise a defined model for studying IAV evolution to increased virulence in the mouse.

Experimental evolution of RNA versus DNA viruses

Based on their extremely high mutation rates, RNA viruses have been traditionally considered as the fastest evolving entities in nature. However, recent work has revealed that, despite their greater replication fidelity, single-stranded (ss) DNA viruses can evolve fast in a similar way. To further investigate this issue, we have compared the rates of adaptation and molecular evolution of ssRNA and ssDNA viruses under highly controlled laboratory conditions using the bacteriophages ΦX174, G4, f1, Qβ, SP, and MS2 as model systems. Our results indicate that ssRNA phages evolve faster than ssDNA phages under strong selective pressure, and that their extremely high mutation rates appear to be optimal for this kind of scenario. However, their performance becomes similar to that of ssDNA phages over the longer term or when the population is moderately well-adapted. Interestingly, the roughly 100-fold difference between the mutation rates of ssRNA and ssDNA phages yields less than a fivefold difference in adaptation and nucleotide substitution rates. The results are therefore consistent with the observation that, despite their lower mutation rates, ssDNA viruses can sometimes match the evolvability of RNA viruses.

What is parallelism?

Although parallel and convergent evolution are discussed extensively in technical articles and textbooks, their meaning can be overlapping, imprecise, and contradictory. The meaning of parallel evolution in much of the evolutionary literature grapples with two separate hypotheses in relation to phenotype and genotype, but often these two hypotheses have been inferred from only one hypothesis, and a number of subsidiary but problematic criteria, in relation to the phenotype. However, examples of parallel evolution of genetic traits that underpin or are at least associated with convergent phenotypes are now emerging. Four criteria for distinguishing parallelism from convergence are reviewed. All are found to be incompatible with any single proposition of homoplasy. Therefore, all homoplasy is equivalent to a broad view of convergence. Based on this concept, all phenotypic homoplasy can be described as convergence and all genotypic homoplasy as parallelism, which can be viewed as the equivalent concept of convergence for molecular data. Parallel changes of molecular traits may or may not be associated with convergent phenotypes but if so describe homoplasy at two biological levels-genotype and phenotype. Parallelism is not an alternative to convergence, but rather it entails homoplastic genetics that can be associated with and potentially explain, at the molecular level, how convergent phenotypes evolve.

From resistance to stimulation: the evolution of a virus in the presence of a dominant lethal inhibitory scaffolding protein

By acquiring resistance to an inhibitor, viruses can become dependent on that inhibitor for optimal fitness. However, inhibitors rarely, if ever, stimulate resistant strain fitness to values that equal or exceed the uninhibited wild-type level. This would require an adaptive mechanism that converts the inhibitor into a beneficial replication factor. Using a plasmid-encoded inhibitory external scaffolding protein that blocks ϕX174 assembly, we previously demonstrated that such mechanisms are possible. The resistant strain, referred to as the evolved strain, contains four mutations contributing to the resistance phenotype. Three mutations confer substitutions in the coat protein, whereas the fourth mutation alters the virus-encoded external scaffolding protein. To determine whether stimulation by the inhibitory protein coevolved with resistance or whether it was acquired after resistance was firmly established, the strain temporally preceding the previously characterized mutant, referred to as the intermediary strain, was isolated and characterized. The results of the analysis indicated that the mutation in the virus-encoded external scaffolding protein was primarily responsible for stimulating strain fitness. When the mutation was placed in a wild-type background, it did not confer resistance. The mutation was also placed in cis with the plasmid-encoded dominant lethal mutation. In this configuration, the stimulating mutation exhibited no activity, regardless of the genotype (wild type, evolved, or intermediary) of the infecting virus. Thus, along with the coat protein mutations, stimulation required two external scaffolding protein genes: the once inhibitory gene and the mutant gene acquired during evolution.

Experimental evolution, loss-of-function mutations, and "the first rule of adaptive evolution"

Adaptive evolution can cause a species to gain, lose, or modify a function; therefore, it is of basic interest to determine whether any of these modes dominates the evolutionary process under particular circumstances. Because mutation occurs at the molecular level, it is necessary to examine the molecular changes produced by the underlying mutation in order to assess whether a given adaptation is best considered as a gain, loss, or modification of function. Although that was once impossible, the advance of molecular biology in the past half century has made it feasible. In this paper, I review molecular changes underlying some adaptations, with a particular emphasis on evolutionary experiments with microbes conducted over the past four decades. I show that by far the most common adaptive changes seen in those examples are due to the loss or modification of a pre-existing molecular function, and I discuss the possible reasons for the prominence of such mutations.

Experimental adaptation of an RNA virus mimics natural evolution

Identification of virulence determinants of viruses is of critical importance in virology. In search of such determinants, virologists traditionally utilize comparative genomics between a virulent and an avirulent virus strain and construct chimeras to map their locations. Subsequent comparison reveals sequence differences, and through analyses of site-directed mutants, key residues are identified. In the absence of a naturally occurring virulent strain, an avirulent strain can be functionally converted to a virulent variant via an experimental evolutionary approach. However, the concern remains whether experimentally evolved virulence determinants mimic those that have evolved naturally. To provide a direct comparison, we exploited a plant RNA virus, soybean mosaic virus (SMV), and its natural host, soybean. Through a serial in vivo passage experiment, the molecularly cloned genome of an avirulent SMV strain was converted to virulent variants on functionally immune soybean genotypes harboring resistance factor(s) from the complex Rsv1 locus. Several of the experimentally evolved virulence determinants were identical to those discovered through a comparative genomic approach with a naturally evolved virulent strain. Thus, our observations validate an experimental evolutionary approach to identify relevant virulence determinants of an RNA virus.

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.

Genetic basis of infectivity evolution in a bacteriophage

Antagonistic coevolution between hosts and parasites is probably ubiquitous. However, very little is known of the genetic changes associated with parasite infectivity evolution during adaptation to a coevolving host. We followed the phenotypic and genetic changes in a lytic virus population (bacteriophage; phage Φ2) that coevolved with its bacterial host, Pseudomonas fluorescens SBW25. First, we show the rapid evolution of numerous unique phage infectivity phenotypes, and that both phage host range and bacterial resistance to individual phage increased over coevolutionary time. Second, each of the distinct phage phenotypes in our study had a unique genotype, and molecular evolution did not act uniformly across the phage genome during coevolution. In particular, we detected numerous substitutions on the tail fibre gene, which is involved in the first step of the host-parasite interaction: host adsorption. None of the observed mutations could be directly linked with infection against a particular host, suggesting that the phenotypic effects of infectivity mutations are probably epistatic. However, phage genotypes with the broadest host ranges had the largest number of nonsynonymous amino acid changes on genes implicated in infectivity evolution. An understanding of the molecular genetics of phage infectivity has helped to explain the complex phenotypic coevolutionary dynamics in this system.

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.

Convergence in pigmentation at multiple levels: mutations, genes and function

Convergence--the independent evolution of the same trait by two or more taxa--has long been of interest to evolutionary biologists, but only recently has the molecular basis of phenotypic convergence been identified. Here, we highlight studies of rapid evolution of cryptic coloration in vertebrates to demonstrate that phenotypic convergence can occur at multiple levels: mutations, genes and gene function. We first show that different genes can be responsible for convergent phenotypes even among closely related populations, for example, in the pale beach mice inhabiting Florida's Gulf and Atlantic coasts. By contrast, the exact same mutation can create similar phenotypes in distantly related species such as mice and mammoths. Next, we show that different mutations in the same gene need not be functionally equivalent to produce similar phenotypes. For example, separate mutations produce divergent protein function but convergent pale coloration in two lizard species. Similarly, mutations that alter the expression of a gene in different ways can, nevertheless, result in similar phenotypes, as demonstrated by sister species of deer mice. Together these studies underscore the importance of identifying not only the genes, but also the precise mutations and their effects on protein function, that contribute to adaptation and highlight how convergence can occur at different genetic levels.

Evolutionary trajectory of white spot syndrome virus (WSSV) genome shrinkage during spread in Asia

BACKGROUND

White spot syndrome virus (WSSV) is the sole member of the novel Nimaviridae family, and the source of major economic problems in shrimp aquaculture. WSSV appears to have rapidly spread worldwide after the first reported outbreak in the early 1990s. Genomic deletions of various sizes occur at two loci in the WSSV genome, the ORF14/15 and ORF23/24 variable regions, and these have been used as molecular markers to study patterns of viral spread over space and time. We describe the dynamics underlying the process of WSSV genome shrinkage using empirical data and a simple mathematical model.

METHODOLOGY/PRINCIPAL FINDINGS

We genotyped new WSSV isolates from five Asian countries, and analyzed this information together with published data. Genome size appears to stabilize over time, and deletion size in the ORF23/24 variable region was significantly related to the time of the first WSSV outbreak in a particular country. Parameter estimates derived from fitting a simple mathematical model of genome shrinkage to the data support a geometric progression (k<1) of the genomic deletions, with k = 0.371 ± 0.150.

CONCLUSIONS/SIGNIFICANCE

The data suggest that the rate of genome shrinkage decreases over time before attenuating. Bioassay data provided support for a link between genome size and WSSV fitness in an aquaculture setting. Differences in genomic deletions between geographic WSSV isolates suggest that WSSV spread did not follow a smooth pattern of geographic radiation, suggesting spread of WSSV over long distances by commercial activities. We discuss two hypotheses for genome shrinkage, an adaptive and a neutral one. We argue in favor of the adaptive hypothesis, given that there is support for a link between WSSV genome size and fitness.

Diminishing returns from beneficial mutations and pervasive epistasis shape the fitness landscape for rifampicin resistance in Pseudomonas aeruginosa

Because adaptation depends upon the fixation of novel beneficial mutations, the fitness effects of beneficial mutations that are substituted by selection are key to our understanding of the process of adaptation. In this study, we experimentally investigated the fitness effects of beneficial mutations that are substituted when populations of the pathogenic bacterium Pseudomonas aeruginosa adapt to the antibiotic rifampicin. Specifically, we isolated the first beneficial mutation to be fixed by selection when 96 populations of three different genotypes of P. aeruginosa that vary considerably in fitness in the presence of rifampicin were challenged with adapting to a high dose of this antibiotic. The simple genetics of rifampicin resistance allowed us to determine the genetic basis of adaptation in the majority of our populations. We show that the average fitness effects of fixed beneficial mutations show a simple and clear pattern of diminishing returns, such that selection tends to fix mutations with progressively smaller effects as populations approach a peak on the adaptive landscape. The fitness effects of individual mutations, on the other hand, are highly idiosyncratic across genetic backgrounds, revealing pervasive epistasis. In spite of this complexity of genetic interactions in this system, there is an overall tendency toward diminishing-returns epistasis. We argue that a simple overall pattern of diminishing-returns adaptation emerges, despite pervasive epistasis between beneficial mutations, because many beneficial mutations are available, and while the fitness landscape is rugged at the fine scale, it is smooth and regular when we consider the average over possible routes to adaptation. In the context of antibiotic resistance, these results show that acquiring mutations that confer low levels of antibiotic resistance does not impose any constraint on the ability to evolve high levels of resistance.

Mutational fitness effects in RNA and single-stranded DNA viruses: common patterns revealed by site-directed mutagenesis studies

The fitness effects of mutations are central to evolution, yet have begun to be characterized in detail only recently. Site-directed mutagenesis is a powerful tool for achieving this goal, which is particularly suited for viruses because of their small genomes. Here, I discuss the evolutionary relevance of mutational fitness effects and critically review previous site-directed mutagenesis studies. The effects of single-nucleotide substitutions are standardized and compared for five RNA or single-stranded DNA viruses infecting bacteria, plants or animals. All viruses examined show very low tolerance to mutation when compared with cellular organisms. Moreover, for non-lethal mutations, the mean fitness reduction caused by single mutations is remarkably constant (0.10-0.13), whereas the fraction of lethals varies only modestly (0.20-0.41). Other summary statistics are provided. These generalizations about the distribution of mutational fitness effects can help us to better understand the evolution of RNA and single-stranded DNA viruses.

Selection for thermostability can lead to the emergence of mutational robustness in an RNA virus

Mutational robustness has important evolutionary implications, yet the mechanisms leading to its emergence remain poorly understood. One possibility is selection acting on a correlated trait, as for instance thermostability (plastogenetic congruence). Here, we examine the correlation between mutational robustness and thermostability in experimental populations of the RNA bacteriophage Qβ. Thermostable viruses evolved after only six serial passages in the presence of heat shocks, and genome sequencing suggested that thermostability can be conferred by several alternative mutations. To test whether thermostable viruses have increased mutational robustness, we performed additional passages in the presence of nitrous acid. Whereas in control lines this treatment produced the expected reduction in growth rate caused by the accumulation of deleterious mutations, thermostable viruses showed no such reduction, indicating that they are more resistant to mutagenesis. Our results suggest that selection for thermostability can lead to the emergence of mutational robustness driven by plastogenetic congruence. As temperature is a widespread selective pressure in nature, the mechanism described here may be relevant to the evolution of mutational robustness.

Some consequences of demographic stochasticity in population genetics

Much of population genetics is based on the diffusion limit of the Wright-Fisher model, which assumes a fixed population size. This assumption is violated in most natural populations, particularly for microbes. Here we study a more realistic model that decouples birth and death events and allows for a stochastically varying population size. Under this model, classical quantities such as the probability of and time before fixation of a mutant allele can differ dramatically from their Wright-Fisher expectations. Moreover, inferences about natural selection based on Wright-Fisher assumptions can yield erroneous and even contradictory conclusions: at small population densities one allele will appear superior, whereas at large densities the other allele will dominate. Consequently, competition assays in laboratory conditions may not reflect the outcome of long-term evolution in the field. These results highlight the importance of incorporating demographic stochasticity into basic models of population genetics.

Simple genomes, complex interactions: epistasis in RNA virus

Owed to their reduced size and low number of proteins encoded, RNA viruses and other subviral pathogens are often considered as being genetically too simple. However, this structural simplicity also creates the necessity for viral RNA sequences to encode for more than one protein and for proteins to carry out multiple functions, all together resulting in complex patterns of genetic interactions. In this work we will first review the experimental studies revealing that the architecture of viral genomes is dominated by antagonistic interactions among loci. Second, we will also review mathematical models and provide a description of computational tools for the study of RNA virus dynamics and evolution. As an application of these tools, we will finish this review article by analyzing a stochastic bit-string model of in silico virus replication. This model analyzes the interplay between epistasis and the mode of replication on determining the population load of deleterious mutations. The model suggests that, for a given mutation rate, the deleterious mutational load is always larger when epistasis is predominantly antagonistic than when synergism is the rule. However, the magnitude of this effect is larger if replication occurs geometrically than if it proceeds linearly.

Evolutionary trajectory of turnip mosaic virus populations adapting to a new host

Little is known about how some plant viruses establish successful cross-species transmission whilst others do not; the genetic basis for adaptation is largely unknown. This study investigated the genetic changes that occurred using the progeny of an infectious clone, p35Tunos, derived from the turnip mosaic virus (TuMV) UK 1 isolate, which has a Brassica host type, but rarely infects Raphanus systemically and then only asymptomatically. The genetic trajectory leading to viral adaptation was studied in a TuMV isolate passaged in Nicotiana benthamiana (parental), Brassica rapa, the old (susceptible) host and Raphanus sativus, the new (almost insusceptible) host. Almost-complete consensus genomic sequences were obtained by RT-PCR of viral populations passaged up to 35 times together with 59 full sequences of 578,200 nt. There were significant differences in the nucleotide and encoded amino acid changes in the consensus genomes from the old and new hosts. Furthermore, a 3264 nt region corresponding to nt 3222-6485 of the UK 1 genome was cloned, and 269 clones from 23 populations were sequenced; this region covered 33 % of the genome and represented a total of 878,016 nt. The results showed that the nucleotide diversity and the non-synonymous/synonymous ratio of the populations from the new host were higher than those from the old host. An analysis of molecular variance showed significant differences among the populations from the old and new hosts. As far as is known, this is the first report comparing the evolutionary trajectory dynamics of plant virus populations in old and new hosts.

Antagonistic coevolution accelerates molecular evolution

The Red Queen hypothesis proposes that coevolution of interacting species (such as hosts and parasites) should drive molecular evolution through continual natural selection for adaptation and counter-adaptation. Although the divergence observed at some host-resistance and parasite-infectivity genes is consistent with this, the long time periods typically required to study coevolution have so far prevented any direct empirical test. Here we show, using experimental populations of the bacterium Pseudomonas fluorescens SBW25 and its viral parasite, phage Phi2 (refs 10, 11), that the rate of molecular evolution in the phage was far higher when both bacterium and phage coevolved with each other than when phage evolved against a constant host genotype. Coevolution also resulted in far greater genetic divergence between replicate populations, which was correlated with the range of hosts that coevolved phage were able to infect. Consistent with this, the most rapidly evolving phage genes under coevolution were those involved in host infection. These results demonstrate, at both the genomic and phenotypic level, that antagonistic coevolution is a cause of rapid and divergent evolution, and is likely to be a major driver of evolutionary change within species.

Evolution: like any other science it is predictable

Evolutionary biology rejoices in the diversity of life, but this comes at a cost: other than working in the common framework of neo-Darwinian evolution, specialists in, for example, diatoms and mammals have little to say to each other. Accordingly, their research tends to track the particularities and peculiarities of a given group and seldom enquires whether there are any wider or deeper sets of explanations. Here, I present evidence in support of the heterodox idea that evolution might look to a general theory that does more than serve as a tautology ('evolution explains evolution'). Specifically, I argue that far from its myriad of products being fortuitous and accidental, evolution is remarkably predictable. Thus, I urge a move away from the continuing obsession with Darwinian mechanisms, which are entirely uncontroversial. Rather, I emphasize why we should seek explanations for ubiquitous evolutionary convergence, as well as the emergence of complex integrated systems. At present, evolutionary theory seems to be akin to nineteenth-century physics, blissfully unaware of the imminent arrival of quantum mechanics and general relativity. Physics had its Newton, biology its Darwin: evolutionary biology now awaits its Einstein.

Identification and dynamics of a beneficial mutation in a long-term evolution experiment with Escherichia coli

BACKGROUND

Twelve populations of E. coli were serially propagated for 20,000 generations in a glucose-supplemented minimal medium in order to study the dynamics of evolution. We sought to find and characterize one of the beneficial mutations responsible for the adaptation and other phenotypic changes, including increased cell size, in one of these populations.

RESULTS

We used transposon-tagging followed by P1-transduction into the ancestor, screening for increased cell size and fitness, co-transduction analysis, and DNA sequencing. We identified a 1-bp insertion in the BoxG1 region located upstream of glmUS, an operon involved in cell-wall biosynthesis. When transduced into the ancestor, this mutation increased competitive fitness by about 5%. This mutation spread through its population of origin between 500 and 1500 generations. Mutations in this region were not found in the other 11 evolving populations, even after 20,000 generations.

CONCLUSION

The 1-bp insertion in the BoxG1 region near glmUS was demonstrably beneficial in the environment in which it arose. The absence of similar mutations in the other evolved populations suggests that they substituted other mutations that rendered this particular mutation unimportant. These results show the unpredictability of adaptive evolution, whereas parallel substitutions at other loci in these same populations reveal the predictability.

Adaptation rates of lytic viruses depend critically on whether host cells survive the bottleneck

We use a branching process approach to estimate the substitution rate, the rate at which beneficial mutations occur and fix, in populations of lytic viruses whose growth is controlled by periodic population bottlenecks. Our model predicts that substitution rates, and by extension adaptation rates, are profoundly affected by the survival of infected host cells at the bottleneck. In particular, we find that direct transfer (or environmental) bottlenecks, in which some fraction of both free virus and host cells are preserved, are associated with relatively slow adaptation rates for the virus. In contrast, viruses can adapt much more quickly when only free virus is transferred to a new host population, as is typical in an epidemiological setting. Finally, when premature lysis of the host-cell population is induced at the bottleneck, we predict that adaptation rates for the virus will, in general, be faster still. These results hold irrespective of the life-history trait affected by the beneficial mutation. The substitution rates in the presence of environmental bottlenecks are predicted to be as much as an order of magnitude lower than in the other two cases.

The dynamics of adaptation on correlated fitness landscapes

Evolutionary theory predicts that a population in a new environment will accumulate adaptive substitutions, but precisely how they accumulate is poorly understood. The dynamics of adaptation depend on the underlying fitness landscape. Virtually nothing is known about fitness landscapes in nature, and few methods allow us to infer the landscape from empirical data. With a view toward this inference problem, we have developed a theory that, in the weak-mutation limit, predicts how a population's mean fitness and the number of accumulated substitutions are expected to increase over time, depending on the underlying fitness landscape. We find that fitness and substitution trajectories depend not on the full distribution of fitness effects of available mutations but rather on the expected fixation probability and the expected fitness increment of mutations. We introduce a scheme that classifies landscapes in terms of the qualitative evolutionary dynamics they produce. We show that linear substitution trajectories, long considered the hallmark of neutral evolution, can arise even when mutations are strongly selected. Our results provide a basis for understanding the dynamics of adaptation and for inferring properties of an organism's fitness landscape from temporal data. Applying these methods to data from a long-term experiment, we infer the sign and strength of epistasis among beneficial mutations in the Escherichia coli genome.

Viral adaptation to an antiviral protein enhances the fitness level to above that of the uninhibited wild type

Viruses often evolve resistance to antiviral agents. While resistant strains are able to replicate in the presence of the agent, they generally exhibit lower fitness than the wild-type strain in the absence of the inhibitor. In some cases, resistant strains become dependent on the antiviral agent. However, the agent rarely, if ever, elevates dependent strain fitness above the uninhibited wild-type level. This would require an adaptive mechanism to convert the antiviral agent into a beneficial growth factor. Using an inhibitory scaffolding protein that specifically blocks phiX174 capsid assembly, we demonstrate that such mechanisms are possible. To obtain the quintuple-mutant resistant strain, the wild-type virus was propagated for approximately 150 viral life cycles in the presence of increasing concentrations of the inhibitory protein. The expression of the inhibitory protein elevated the strain's fitness significantly above the uninhibited wild-type level. Thus, selecting for resistance coselected for dependency, which was characterized and found to operate on the level of capsid nucleation. To the best of our knowledge, this is the first report of a virus evolving a mechanism to productively utilize an antiviral agent to stimulate its fitness above the uninhibited wild-type level. The results of this study may be predictive of the types of resistant phenotypes that could be selected by antiviral agents that specifically target capsid assembly.

The genetics of adaptation for eight microvirid bacteriophages

Theories of adaptive molecular evolution have recently experienced significant expansion, and their predictions and assumptions have begun to be subjected to rigorous empirical testing. However, these theories focus largely on predicting the first event in adaptive evolution, the fixation of a single beneficial mutation. To address long-term adaptation it is necessary to include new assumptions, but empirical data are needed for guidance. To empirically characterize the general properties of adaptive walks, eight recently isolated relatives of the single-stranded DNA (ssDNA) bacteriophage φX174 (family Microviridae) were adapted to identical selective conditions. Three of the eight genotypes were adapted in replicate, for a total of 11 adaptive walks. We measured fitness improvement and identified the genetic changes underlying the observed adaptation. Nearly all phages were evolvable; nine of the 11 lineages showed a significant increase in fitness. However, fitness plateaued quickly, and adaptation was achieved through only three substitutions on average. Parallel evolution was rampant, both across replicates of the same genotype as well as across different genotypes, yet adaptation of replicates never proceeded through the exact same set of mutations. Despite this, final fitnesses did not vary significantly among replicates. Final fitnesses did vary significantly across genotypes but not across phylogenetic groupings of genotypes. A positive correlation was found between the number of substitutions in an adaptive walk and the magnitude of fitness improvement, but no correlation was found between starting and ending fitness. These results provide an empirical framework for future adaptation theory.

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.