Effect of deleterious mutation-accumulation on the fitness of RNA bacteriophage MS2

The puzzling mutational landscape of the SARS-2-variant Omicron

The recently emerging SARS-CoV-2 variant omicron displays an unusual association of 30 mutations, 3 deletions, and 1 insertion. To analyze the impact of this atypic mutational landscape, we constructed a complete structure of the omicron spike protein. Compared with the delta variant, the receptor-binding domain (RBD) of omicron has an increased electrostatic surface potential, but a decreased affinity for the ACE-2 receptor. The N-terminal domain (NTD) has both a decreased surface potential and a lower affinity for lipid rafts. The omicron variant is predicted to be less fusogenic and thus less pathogenic than delta, due to a geometric reorganization of the S1-S2 cleavage site. Overall, these virological parameters suggest that omicron does not have a significant infectivity advantage over the delta variant. However, in omicron, neutralizing epitopes are greatly affected, suggesting that current vaccines will probably confer little protection against this variant. In conclusion, the puzzling mutational pattern of the omicron variant combines contradictory properties which may either decrease (virological properties) or increase (immunological escape/facilitation) the transmission of this variant in the human population. This Janus-like phenotype may explain some conflicting reports on the initial assessment of omicron and provide new insights about the molecular mechanisms controlling its dissemination and pathogenesis worldwide.

Direct sequencing of RNA with MinION Nanopore: detecting mutations based on associations

One of the key challenges in the field of genetics is the inference of haplotypes from next generation sequencing data. The MinION Oxford Nanopore sequencer allows sequencing long reads, with the potential of sequencing complete genes, and even complete genomes of viruses, in individual reads. However, MinION suffers from high error rates, rendering the detection of true variants difficult. Here, we propose a new statistical approach named AssociVar, which differentiates between true mutations and sequencing errors from direct RNA/DNA sequencing using MinION. Our strategy relies on the assumption that sequencing errors will be dispersed randomly along sequencing reads, and hence will not be associated with each other, whereas real mutations will display a non-random pattern of association with other mutations. We demonstrate our approach using direct RNA sequencing data from evolved populations of the MS2 bacteriophage, whose small genome makes it ideal for MinION sequencing. AssociVar inferred several mutations in the phage genome, which were corroborated using parallel Illumina sequencing. This allowed us to reconstruct full genome viral haplotypes constituting different strains that were present in the sample. Our approach is applicable to long read sequencing data from any organism for accurate detection of bona fide mutations and inter-strain polymorphisms.

Effects of life history and ecology on virus evolutionary potential

The life history traits of viruses pose many consequences for viral population structure. In turn, population structure may influence the evolutionary trajectory of a virus. Here we review factors that affect the evolutionary potential of viruses, including rates of mutation and recombination, bottlenecks, selection pressure, and ecological factors such as the requirement for hosts and vectors. Mutation, while supplying a pool of raw genetic material, also results in the generation of numerous unfit mutants. The infection of multiple host species may expand a virus' ecological niche, although it may come at a cost to genetic diversity. Vector-borne viruses often experience a diminished frequency of positive selection and exhibit little diversity, and resistance against vector-borne viruses may thus be more durable than against non-vectored viruses. Evidence indicates that adaptation to a vector is more evolutionarily difficult than adaptation to a host. Overall, a better understanding of how various factors influence viral dynamics in both plant and animal pathosystems will lead to more effective anti-viral treatments and countermeasures.

Phylogenetic and recombination analysis of Tobacco bushy top virus in China

BACKGROUND

During the past decade, tobacco bushy top disease, which is mainly caused by a combination of Tobacco bushy top virus (TBTV) and Tobacco vein-distorting virus (TVDV), underwent a sudden appearance, extreme virulence and degeneration of the epidemic in the Yunnan province of China. In addition to integrative control of its aphid vector, it is of interest to examine diversity and evolution among different TBTV isolates.

METHODS

5' and 3' RACE, combined with one step full-length RT-PCR, were used to clone the full-length genome of three new isolates of TBTV that exhibited mild pathogenicity in Chinese fields. Nucleotide and amino acid sequences for the TBTV isolates were analyzed by DNAMAN. MEGA 5.0 was used to construct phylogenetic trees. RDP4 was used to detect recombination events during evolution of these isolates.

RESULTS

The genomes of three isolates, termed TBTV-JC, TBTV-MD-I and TBTV-MD-II, were 4152 nt in length and included one distinctive difference from previously reported TBTV isolates: the first nucleotide of the genome was a guanylate instead of an adenylate. Diversity and phylogenetic analyses among these three new TBTV isolates and five other available isolates suggest that ORFs and 3'UTRs of TBTV may have evolved separately. Moreover, the RdRp-coding region was the most variable. Recombination analysis detected a total of 29 recombination events in the 8 TBTV isolates, in which 24 events are highly likely and 5 events have low-level likelihood based on their correlation with the phylogenetic trees. The three new TBTV isolates have individual recombination patterns with subtle divergences in parents and locations.

CONCLUSIONS

The genome sizes of TBTV isolates were constant while different ORF-coding regions and 3'UTRs may have evolved separately. The RdRp-coding region was the most variable. Frequent recombination occurred among TBTV isolates. Three new TBTV isolates have individual recombination patterns and may have different progenitors.

Experimentally guided models reveal replication principles that shape the mutation distribution of RNA viruses

Life history theory posits that the sequence and timing of events in an organism's lifespan are fine-tuned by evolution to maximize the production of viable offspring. In a virus, a life history strategy is largely manifested in its replication mode. Here, we develop a stochastic mathematical model to infer the replication mode shaping the structure and mutation distribution of a poliovirus population in an intact single infected cell. We measure production of RNA and poliovirus particles through the infection cycle, and use these data to infer the parameters of our model. We find that on average the viral progeny produced from each cell are approximately five generations removed from the infecting virus. Multiple generations within a single cell infection provide opportunities for significant accumulation of mutations per viral genome and for intracellular selection.

DOI:http://dx.doi.org/10.7554/eLife.03753.001

Viruses with genetic information made up of molecules of RNA can multiply quickly, but not very accurately. This means that many errors, or mutations, occur when the RNA is copied to create new viruses. The advantage of this rapid, but mistake-filled, RNA replication process is that some of the mutations will be beneficial to the virus. This allows viruses to rapidly evolve, for example, to develop resistance against drugs.

The poliovirus is an RNA virus that can cause paralysis and death in humans. To prevent such infections, scientists have extensively studied the poliovirus and have developed effective vaccines against it that have eliminated the virus from all but a few countries. Because so much is known about the poliovirus and because it has a very simple structure, scientists continue to use the poliovirus as a model to study virus behavior.

One unknown aspect of the poliovirus' behavior is how it replicates after invading a cell. Are all of its RNA copies made from the original viral RNA that first infected the cell, in what is known as a ‘stamping machine’ model? Or do the new copies of the RNA also get copied themselves in a ‘geometric replication mode’ that increases the likelihood of mutations and enables the virus to evolve more rapidly?

Viral RNA molecules are copied by one of the virus's own proteins and so before the viral RNA can be replicated, it must first be translated to form viral proteins. When and where replication begins depends on the concentration of translated proteins around the RNA and so replication tends to begin in particular areas of the cell at different times. Schulte, Draghi et al. used mathematical modeling to create computer simulations of the number of polioviruses in a cell that take into account these time and space constraints. By including random elements in the model, the simulated behavior more accurately follows experimentally recorded data than previously used models.

The results of the model led Schulte, Draghi et al. to conclude that the poliovirus replicates by the ‘geometric mode’; as new copies of the poliovirus RNA are made, each copy goes on to make more copies. This means that in a single infected cell there are multiple generations of RNA, and each generation may undergo distinct mutations that are passed on to the next set of RNA copies. In fact, Schulte, Draghi et al. found that the average virus released from an infected cell is the great-great-great-granddaughter of the original virus that infected the cell. With so many different generations of virus coexisting in a cell, there are a lot of opportunities for new genetic combinations to occur and for viruses to evolve new abilities.

DOI:http://dx.doi.org/10.7554/eLife.03753.002

Changes in protein domains outside the catalytic site of the bacteriophage Qβ replicase reduce the mutagenic effect of 5-azacytidine

The high genetic heterogeneity and great adaptability of RNA viruses are ultimately caused by the low replication fidelity of their polymerases. However, single amino acid substitutions that modify replication fidelity can evolve in response to mutagenic treatments with nucleoside analogues. Here, we investigated how two independent mutants of the bacteriophage Qβ replicase (Thr210Ala and Tyr410His) reduce sensitivity to the nucleoside analogue 5-azacytidine (AZC). Despite being located outside the catalytic site, both mutants reduced the mutation frequency in the presence of the drug. However, they did not modify the type of AZC-induced substitutions, which was mediated mainly by ambiguous base pairing of the analogue with purines. Furthermore, the Thr210Ala and Tyr410His substitutions had little or no effect on replication fidelity in untreated viruses. Also, both substitutions were costly in the absence of AZC or when the action of the drug was suppressed by adding an excess of natural pyrimidines (uridine or cytosine). Overall, the phenotypic properties of these two mutants were highly convergent, despite the mutations being located in different domains of the Qβ replicase. This suggests that treatment with a given nucleoside analogue tends to select for a unique functional response in the viral replicase.

IMPORTANCE

In the last years, artificial increase of the replication error rate has been proposed as an antiviral therapy. In this study, we investigated the mechanisms by which two substitutions in the Qβ replicase confer partial resistance to the mutagenic nucleoside analogue AZC. As opposed to previous work with animal viruses, where different mutations selected sequentially conferred nucleoside analogue resistance through different mechanisms, our results suggest that there are few or no alternative AZC resistance phenotypes in Qβ. Also, despite resistance mutations being highly costly in the absence of the drug, there was no sequential fixation of secondary mutations. Bacteriophage Qβ is the virus with the highest reported mutation rate, which should make it particularly sensitive to nucleoside analogue treatments, probably favoring resistance mutations even if they incur high costs. The results are also relevant for understanding the possible pathways by which fidelity of the replication machinery can be modified.

Variability in mutational fitness effects prevents full lethal transitions in large quasispecies populations

The distribution of mutational fitness effects (DMFE) is crucial to the evolutionary fate of quasispecies. In this article we analyze the effect of the DMFE on the dynamics of a large quasispecies by means of a phenotypic version of the classic Eigen's model that incorporates beneficial, neutral, deleterious, and lethal mutations. By parameterizing the model with available experimental data on the DMFE of Vesicular stomatitis virus (VSV) and Tobacco etch virus (TEV), we found that increasing mutation does not totally push the entire viral quasispecies towards deleterious or lethal regions of the phenotypic sequence space. The probability of finding regions in the parameter space of the general model that results in a quasispecies only composed by lethal phenotypes is extremely small at equilibrium and in transient times. The implications of our findings can be extended to other scenarios, such as lethal mutagenesis or genomically unstable cancer, where increased mutagenesis has been suggested as a potential therapy.

A new look at an old virus: patterns of mutation accumulation in the human H1N1 influenza virus since 1918

BACKGROUND

The H1N1 influenza A virus has been circulating in the human population for over 95 years, first manifesting itself in the pandemic of 1917-1918. Initial mortality was extremely high, but dropped exponentially over time. Influenza viruses have high mutation rates, and H1N1 has undergone significant genetic changes since 1918. The exact nature of H1N1 mutation accumulation over time has not been fully explored.

METHODS

We have made a comprehensive historical analysis of mutational changes within H1N1 by examining over 4100 fully-sequenced H1N1 genomes. This has allowed us to examine the genetic changes arising within H1N1 from 1918 to the present.

RESULTS

We document multiple extinction events, including the previously known extinction of the human H1N1 lineage in the 1950s, and an apparent second extinction of the human H1N1 lineage in 2009. These extinctions appear to be due to a continuous accumulation of mutations. At the time of its disappearance in 2009, the human H1N1 lineage had accumulated over 1400 point mutations (more than 10% of the genome), including approximately 330 non-synonymous changes (7.4% of all codons). The accumulation of both point mutations and non-synonymous amino acid changes occurred at constant rates (μ = 14.4 and 2.4 new mutations/year, respectively), and mutations accumulated uniformly across the entire influenza genome. We observed a continuous erosion over time of codon-specificity in H1N1, including a shift away from host (human, swine, and bird [duck]) codon preference patterns.

CONCLUSIONS

While there have been numerous adaptations within the H1N1 genome, most of the genetic changes we document here appear to be non-adaptive, and much of the change appears to be degenerative. We suggest H1N1 has been undergoing natural genetic attenuation, and that significant attenuation may even occur during a single pandemic. This process may play a role in natural pandemic cessation and has apparently contributed to the exponential decline in mortality rates over time, as seen in all major human influenza strains. These findings may be relevant to the development of strategies for managing influenza pandemics and strain evolution.

The causes of epistasis

Since Bateson's discovery that genes can suppress the phenotypic effects of other genes, gene interactions-called epistasis-have been the topic of a vast research effort. Systems and developmental biologists study epistasis to understand the genotype-phenotype map, whereas evolutionary biologists recognize the fundamental importance of epistasis for evolution. Depending on its form, epistasis may lead to divergence and speciation, provide evolutionary benefits to sex and affect the robustness and evolvability of organisms. That epistasis can itself be shaped by evolution has only recently been realized. Here, we review the empirical pattern of epistasis, and some of the factors that may affect the form and extent of epistasis. Based on their divergent consequences, we distinguish between interactions with or without mean effect, and those affecting the magnitude of fitness effects or their sign. Empirical work has begun to quantify epistasis in multiple dimensions in the context of metabolic and fitness landscape models. We discuss possible proximate causes (such as protein function and metabolic networks) and ultimate factors (including mutation, recombination, and the importance of natural selection and genetic drift). We conclude that, in general, pleiotropy is an important prerequisite for epistasis, and that epistasis may evolve as an adaptive or intrinsic consequence of changes in genetic robustness and evolvability.

Find the weakest link. A comparison between demographic, genetic and demo-genetic metapopulation extinction times

Background

While the ultimate causes of most species extinctions are environmental, environmental constraints have various secondary consequences on evolutionary and ecological processes. The roles of demographic, genetic mechanisms and their interactions in limiting the viabilities of species or populations have stirred much debate and remain difficult to evaluate in the absence of demography-genetics conceptual and technical framework. Here, I computed projected times to metapopulation extinction using (1) a model focusing on the effects of species properties, habitat quality, quantity and temporal variability on the time to demographic extinction; (2) a genetic model focusing on the dynamics of the drift and inbreeding loads under the same species and habitat constraints; (3) a demo-genetic model accounting for demographic-genetic processes and feedbacks.

Results

Results indicate that a given population may have a high demographic, but low genetic viability or vice versa; and whether genetic or demographic aspects will be the most limiting to overall viability depends on the constraints faced by the species (e.g., reduction of habitat quantity or quality). As a consequence, depending on metapopulation or species characteristics, incorporating genetic considerations to demographically-based viability assessments may either moderately or severely reduce the persistence time. On the other hand, purely genetically-based estimates of species viability may either underestimate (by neglecting demo-genetic interactions) or overestimate (by neglecting the demographic resilience) true viability.

Conclusion

Unbiased assessments of the viabilities of species may only be obtained by identifying and considering the most limiting processes (i.e., demography or genetics), or, preferentially, by integrating them.

Genetically engineered virulent phage banks in the detection and control of emergent pathogenic bacteria

Natural outbreaks of multidrug-resistant microorganisms can cause widespread devastation, and several can be used or engineered as agents of bioterrorism. From a biosecurity standpoint, the capacity to detect and then efficiently control, within hours, the spread and the potential pathological effects of an emergent outbreak, for which there may be no effective antibiotics or vaccines, become key challenges that must be met. We turned to phage engineering as a potentially highly flexible and effective means to both detect and eradicate threats originating from emergent (uncharacterized) bacterial strains. To this end, we developed technologies allowing us to (1) concurrently modify multiple regions within the coding sequence of a gene while conserving intact the remainder of the gene, (2) reversibly interrupt the lytic cycle of an obligate virulent phage (T4) within its host, (3) carry out efficient insertion, by homologous recombination, of any number of engineered genes into the deactivated genomes of a T4 wild-type phage population, and (4) reactivate the lytic cycle, leading to the production of engineered infective virulent recombinant progeny. This allows the production of very large, genetically engineered lytic phage banks containing, in an E. coli host, a very wide spectrum of variants for any chosen phage-associated function, including phage host-range. Screening of such a bank should allow the rapid isolation of recombinant T4 particles capable of detecting (ie, diagnosing), infecting, and destroying hosts belonging to gram-negative bacterial species far removed from the original E. coli host.

Computational models of HIV-1 resistance to gene therapy elucidate therapy design principles

Gene therapy is an emerging alternative to conventional anti-HIV-1 drugs, and can potentially control the virus while alleviating major limitations of current approaches. Yet, HIV-1's ability to rapidly acquire mutations and escape therapy presents a critical challenge to any novel treatment paradigm. Viral escape is thus a key consideration in the design of any gene-based technique. We develop a computational model of HIV's evolutionary dynamics in vivo in the presence of a genetic therapy to explore the impact of therapy parameters and strategies on the development of resistance. Our model is generic and captures the properties of a broad class of gene-based agents that inhibit early stages of the viral life cycle. We highlight the differences in viral resistance dynamics between gene and standard antiretroviral therapies, and identify key factors that impact long-term viral suppression. In particular, we underscore the importance of mutationally-induced viral fitness losses in cells that are not genetically modified, as these can severely constrain the replication of resistant virus. We also propose and investigate a novel treatment strategy that leverages upon gene therapy's unique capacity to deliver different genes to distinct cell populations, and we find that such a strategy can dramatically improve efficacy when used judiciously within a certain parametric regime. Finally, we revisit a previously-suggested idea of improving clinical outcomes by boosting the proliferation of the genetically-modified cells, but we find that such an approach has mixed effects on resistance dynamics. Our results provide insights into the short- and long-term effects of gene therapy and the role of its key properties in the evolution of resistance, which can serve as guidelines for the choice and optimization of effective therapeutic agents.

A primary obstacle to the success of any anti-HIV treatment is HIV's ability to rapidly resist it by generating new viral strains whose vulnerability to the treatment is reduced. Gene therapies represent a novel class of treatments for HIV infection that may supplement or replace present therapies, as they alleviate some of their major shortcomings. The design of gene therapeutic agents that effectively reduce viral resistance can be aided by a quantitative elucidation of the processes by which resistance is acquired following therapy initiation. We developed a computational model that describes a patient's response to therapy and used it to quantify the influence of therapy parameters and strategies on the development of viral resistance. We find that gene therapy induces different clinical conditions and a much slower viral response than present therapies. These dictate different design principles such as a greater significance to the virus' competence in the absence of therapy. We also show that one can effectively delay emergence of resistance by delivering distinct therapeutic genes into separate cell populations. Our results highlight the differences between traditional and gene therapies and provide a basic understanding of how key controllable parameters and strategies affect resistance development.

Estimation of the size of genetic bottlenecks in cell-to-cell movement of soil-borne wheat mosaic virus and the possible role of the bottlenecks in speeding up selection of variations in trans-acting genes or elements

Genetic bottlenecks facilitate the fixation and extinction of variants in populations, and viral populations are no exception to this theory. To examine the existence of genetic bottlenecks in cell-to-cell movement of plant RNA viruses, we prepared constructs for Soil-borne wheat mosaic virus RNA2 vectors carrying two different fluorescent proteins, yellow fluorescent protein (YFP) and cyan fluorescent protein (CFP). Coinoculation of host plant leaves with the two RNA2 vectors and the wild-type RNA1 showed separation of the two vector RNA2s, mostly within seven to nine cell-to-cell movements from individual initially coinfected cells. Our statistical analysis showed that the number of viral RNA genomes establishing infection in adjacent cells after the first cell-to-cell movement from an initially infected cell was 5.97 +/- 0.22 on average and 5.02 +/- 0.29 after the second cell-to-cell movement. These results indicate that plant RNA viruses may generally face narrow genetic bottlenecks in every cell-to-cell movement. Furthermore, our model suggests that, rather than suffering from fitness losses caused by the bottlenecks, the plant RNA viruses are utilizing the repeated genetic bottlenecks as an essential element of rapid selection of their adaptive variants in trans-acting genes or elements to respond to host shifting and changes in the growth conditions of the hosts.

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.

Optimal bacteriophage mutation rates for phage therapy

The mutability of bacteriophages offers a particular advantage in the treatment of bacterial infections not afforded by other antimicrobial therapies. When phage-resistant bacteria emerge, mutation may generate phage capable of exploiting and thus limiting population expansion among these emergent types. However, while mutation potentially generates beneficial variants, it also contributes to a genetic load of deleterious mutations. Here, we model the influence of varying phage mutation rate on the efficacy of phage therapy. All else being equal, phage types with historical mutation rates of approximately 0.1 deleterious mutations per genome per generation offer a reasonable balance between beneficial mutational diversity and deleterious mutational load. We determine that increasing phage inoculum density can undesirably increase the peak density of a mutant bacterial class by limiting the in situ production of mutant phage variants. For phage populations with minimal genetic load, engineering mutation rate increases beyond the mutation-selection balance optimum may provide even greater protection against emergent bacterial types, but only with very weak selective coefficients for de novo deleterious mutations (below approximately 0.01). Increases to the mutation rate beyond the optimal value at mutation-selection balance may therefore prove generally undesirable.

Genetic and phenotypic variation of foot-and-mouth disease virus during serial passages in a natural host

Foot-and-mouth disease virus (FMDV), like other RNA viruses, exhibits high mutation rates during replication that have been suggested to be of adaptive value. However, even though genetic variation in RNA viruses and, more specifically, FMDV has been extensively examined during virus replication in a wide variety of in vitro cell cultures, very little is known regarding the generation and effects of genetic variability of virus replication in the natural host under experimental conditions and no genetic data are available regarding the effects of serial passage in natural hosts. Here, we present the results of 20 serial contact transmissions of the highly pathogenic, pig-adapted O Taiwan 97 (O Tw97) isolate of FMDV in swine. We examined the virus genomic consensus sequences for a total of 37 full-length viral genomes recovered from 20 in vivo passages. The characteristics and distributions of changes in the sequences during the series of pig infections were analyzed in comparison to the O Tw97 genomes recovered from serially infected BHK-21 cell cultures. Unexpectedly, a significant reduction of virulence upon pig passages was observed, and finally, interruption of the viral transmission chain occurred after the14th pig passage (T14). Virus was, however, isolated from the tonsils and nasal swabs of the asymptomatic T15 pigs at 26 days postcontact, consistent with a natural establishment of the carrier state previously described only for ruminants. Surprisingly, the region encoding the capsid protein VP1 (1D) did not show amino acid changes during in vivo passages. These data demonstrate that contact transmission of FMDV O Tw97 in pigs mimics the fitness loss induced by the bottleneck effect, which was previously observed by others during plaque-to-plaque FMDV passage in vitro, suggesting that unknown mechanisms of virulence recovery might be necessary during the evolution and perpetuation of FMDV in nature.

Fitness declines in Tobacco etch virus upon serial bottleneck transfers

It has been well established that populations of RNA viruses transmitted throughout serial bottlenecks suffer from significant fitness declines as a consequence of the accumulation of deleterious mutations by the onset of Muller's ratchet. Bottlenecks are unavoidably linked to different steps of the infectious cycle of most plant RNA viruses, such as vector-mediated transmissions and systemic colonization of new leaves. Here we report evidence for fitness declines by the accumulation of deleterious mutations in the potyvirus Tobacco etch virus (TEV). TEV was inoculated into the nonsystemic host Chenopodium quinoa, and local lesions were isolated and used to initiate 20 independent mutation accumulation lineages. Weekly, a random lesion from each lineage was isolated and used to inoculate the next set of plants. At each transfer, the Malthusian growth rate was estimated. After 11 consecutive transfers, all lineages suffered significant fitness losses, and one even became extinct. The average rate of fitness decline was 5% per day. The average pattern of fitness decline was consistent with antagonistic epistasis between deleterious mutations, as postulated for antiredundant genomes. Temporal fitness fluctuations were not explained by random noise but reflected more complex underlying processes related to emergence and self-organization phenomena.

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.

The evolution of sex: empirical insights into the roles of epistasis and drift

Despite many years of theoretical and experimental work, the explanation for why sex is so common as a reproductive strategy continues to resist understanding. Recent empirical work has addressed key questions in this field, especially regarding rates of mutation accumulation in sexual and asexual organisms, and the roles of negative epistasis and drift as sources of adaptive constraint in asexually reproducing organisms. At the same time, new ideas about the evolution of sexual recombination are being tested, including intriguing suggestions of an important interplay between sex and genetic architecture, which indicate that sex and recombination could have affected their own evolution.

The Hill-Robertson effect and the evolution of recombination

In finite populations, genetic drift generates interference between selected loci, causing advantageous alleles to be found more often on different chromosomes than on the same chromosome, which reduces the rate of adaptation. This "Hill-Robertson effect" generates indirect selection to increase recombination rates. We present a new method to quantify the strength of this selection. Our model represents a new beneficial allele (A) entering a population as a single copy, while another beneficial allele (B) is sweeping at another locus. A third locus affects the recombination rate between selected loci. Using a branching process model, we calculate the probability distribution of the number of copies of A on the different genetic backgrounds, after it is established but while it is still rare. Then, we use a deterministic model to express the change in frequency of the recombination modifier, due to hitchhiking, as A goes to fixation. We show that this method can give good estimates of selection for recombination. Moreover, it shows that recombination is selected through two different effects: it increases the fixation probability of new alleles, and it accelerates selective sweeps. The relative importance of these two effects depends on the relative times of occurrence of the beneficial alleles.

Mechanisms of genetic robustness in RNA viruses

Two key features of RNA viruses are their compacted genomes and their high mutation rate. Accordingly, deleterious mutations are common and have an enormous impact on viral fitness. In their multicellular hosts, robustness can be achieved by genomic redundancy, including gene duplication, diploidy, alternative metabolic pathways and biochemical buffering mechanisms. However, here we review evidence suggesting that during RNA virus evolution, alternative robustness mechanisms may have been selected. After briefly describing how genetic robustness can be quantified, we discuss mechanisms of intrinsic robustness arising as consequences of RNA-genome architecture, replication peculiarities and quasi-species population dynamics. These intrinsic robustness mechanisms operate efficiently at the population level, despite the mutational sensitivity shown by individual genomes. Finally, we discuss the possibility that viruses might exploit cellular buffering mechanisms for their own benefit, producing a sort of extrinsic robustness.

Effect of varying epistasis on the evolution of recombination

Whether recombination decelerates or accelerates a population's response to selection depends, at least in part, on how fitness-determining loci interact. Realistically, all genomes likely contain fitness interactions both with positive and with negative epistasis. Therefore, it is crucial to determine the conditions under which the potential beneficial effects of recombination with negative epistasis prevail over the detrimental effects of recombination with positive epistasis. Here, we examine the simultaneous effects of diverse epistatic interactions with different strengths and signs in a simplified model system with independent pairs of interacting loci and selection acting only on the haploid phase. We find that the average form of epistasis does not predict the average amount of linkage disequilibrium generated or the impact on a recombination modifier when compared to results using the entire distribution of epistatic effects and associated single-mutant effects. Moreover, we show that epistatic interactions of a given strength can produce very different effects, having the greatest impact when selection is weak. In summary, we observe that the evolution of recombination at mutation-selection balance might be driven by a small number of interactions with weak selection rather than by the average epistasis of all interactions. We illustrate this effect with an analysis of published data of Saccharomyces cerevisiae. Thus to draw conclusions on the evolution of recombination from experimental data, it is necessary to consider the distribution of epistatic interactions together with the associated selection coefficients.

A measure of the within-chromosome synergistic epistasis for Drosophila viability

In order to detect possible synergistic epistasis for viability in Drosophila melanogaster we assayed the relative viability of chromosomes II in: (i) panmixia, (ii) forced total homozygosity, and (iii) homozygosity for, on the average, half of their loci. As these genotypes were constructed using exactly the same set of chromosomes in the three cases, the design allows us to estimate the inbreeding depression rate at two different inbreeding levels in the absence of purging natural selection. Overall, no consistent synergistic epistasis was found. However, there was a small fraction of chromosomes whose severely deleterious effect when homozygous was almost significantly larger than expected from their viability when homozygous for half of their loci. This suggests occasional but important synergistic epistasis, which might confer evolutionary advantage to recombination in tightly linked genomes. Nevertheless, such epistasis is unlikely to be an evolutionary advantage driving the evolution of sexual anisogamous reproduction, as its contribution to overall viability is small when compared with the two-fold cost of anisogamy.

Parasites and mutational load: an experimental test of a pluralistic theory for the evolution of sex

Ecological and mutational explanations for the evolution of sexual reproduction have usually been considered independently. Although many of these explanations have yielded promising theoretical results,experimental support for their ability to overcome a twofold cost of sex has been limited. For this reason, it has recently been argued that a pluralistic approach, combining effects from multiple models, may be necessary to explain the apparent advantage of sex. One such pluralistic model proposes that parasite load and synergistic epistasis between deleterious mutations might interact to create an advantage for recombination.Here, we test this proposal by comparing the fitness functions of parasitized and parasite-free genotypes of Escherichia coli bearing known numbers of transposon-insertion mutations. In both classes, we failed to detect any evidence for synergistic epistasis. However, the average effect of deleterious mutations was greater in parasitized than parasite-free genotypes. This effect might broaden the conditions under which another proposed model combining parasite-host coevolutionary dynamics and mutation accumulation can explain the maintenance of sex. These results suggest that, on average, deleterious mutations act multiplicatively with each other but in synergy with infection in determining fitness.

Self-fertilization and the evolution of recombination

In this article, we study the effect of self-fertilization on the evolution of a modifier allele that alters the recombination rate between two selected loci. We consider two different life cycles: under gametophytic selfing, a given proportion of fertilizations involves gametes produced by the same haploid individual, while under sporophytic selfing, a proportion of fertilizations involves gametes produced by the same diploid individual. Under both life cycles, we derive approximations for the change in frequency of the recombination modifier when selection is weak relative to recombination, so that the population reaches a state of quasi-linkage equilibrium. We find that gametophytic selfing increases the range of epistasis under which increased recombination is favored; however, this effect is substantial only for high selfing rates. Moreover, gametophytic selfing affects the relative influence of different components of epistasis (additive x additive, additive x dominance, dominance x dominance) on the evolution of the modifier. Sporophytic selfing has much stronger effects: even a small selfing rate greatly increases the parameter range under which recombination is favored, when there is negative dominance x dominance epistasis. This effect is due to the fact that selfing generates a correlation in homozygosity at linked loci, which is reduced by recombination.

High mutation rates, bottlenecks, and robustness of RNA viral quasispecies

Population bottlenecks are stochastic events that strongly condition the structure and evolution of natural populations. Their effects are readily observable in highly heterogeneous populations, such as RNA viruses, since bottlenecks cause a fast accumulation of mutations. Considering that most mutations are deleterious, it was predicted that the frequent application of bottlenecks would yield a population unable to replicate. However, in vitro as well as in vivo systems evolving through bottlenecks present a remarkable resistance to extinction. This observation reveals the robustness of RNA viruses and points to the existence of internal mechanisms which must confer a high degree of adaptability to fast mutating populations. In this contribution, we review experimental observations regarding the survival of RNA viruses, both in laboratory experiments and in natural populations. By means of a simple theoretical model of evolution which incorporates strong reductions of the population size, we explore the relationship between the number of replication rounds that a single founder particle undergoes before the next bottleneck is applied, and the mutation rate in a particular environment. Our numerical results reveal that the mutation rate has evolved in a concerted way with the degree of optimization achieved by the population originated from the founder particle. We hypothesize that this mechanism generates a mutation-selection equilibrium in natural populations that maximizes adaptability while maintaining their structure.

Molecular Basis of Fitness Loss and Fitness Recovery in Vesicular Stomatitis Virus

Viral populations subjected to repeated genetic bottleneck accumulate deleterious mutations in a process known as Muller's ratchet. Asexual viruses, such as vesicular stomatitis virus (VSV) can recover from Muller's ratchet by replication with large effective population sizes. However, mutants with a history of bottleneck transmissions often show decreased adaptability when compared to non-bottlenecked populations. We have generated a collection of bottlenecked mutants and allowed them to recover by large population passages. We have characterized fitness changes and the complete genomes of these strains. Mutations accumulated during the operation of Muller's ratchet led to the identification of two potential mutational hot spots in the VSV genome. As in other viral systems, transitions were more common than transversions. Both back mutation and compensatory mutations contributed to recovery, although a significant level of fitness increase was observed in nine of the 13 bottlenecked strains with no obvious changes in the consensus sequence. Additional replication of three strains resulted in the fixation of single point mutations. Only two mutations previously found in non-bottlenecked, high-fitness populations that had been adapting to the same environment were identified in the recovered strains.

Adaptation to the deleterious effects of antimicrobial drug resistance mutations by compensatory evolution

Compensatory mutations, due to their ability to mask the deleterious effects of another mutation, are important for the adaptation and evolution of most organisms. Resistance to antibiotics, antivirals, antifungals, herbicides and insecticides is usually associated with a fitness cost. As a result of compensatory evolution, the initial fitness costs conferred by resistance mutations (or other deleterious mutations) can often be rapidly and efficiently reduced. Such compensatory evolution is potentially of importance for (i) the long-term persistence of drug resistance, (ii) reducing the rate of fitness loss associated with the accumulation of deleterious mutations in small asexual populations, and (iii) the evolution of complexity of cellular processes.

Recombination in HIV and the evolution of drug resistance: for better or for worse?

The rapid evolution of drug resistance remains a major obstacle for HIV therapy. The capacity of the virus for recombination is widely believed to facilitate the evolution of drug resistance. Here, we challenge this intuitive view. We develop a population genetic model of HIV replication that incorporates the processes of mutation, cellular superinfection, and recombination. We show that cellular superinfection increases the abundance of low fitness viruses at the expense of the fittest strains due to the mixing of viral proteins during virion assembly. Moreover, we argue that whether recombination facilitates the evolution of drug resistance depends critically on how resistance mutations interact to determine viral fitness. Contrary to the commonly held belief, we find that, under the most plausible biological assumptions, recombination is expected to slow down the rate of evolution of multi-drug-resistant virus during therapy.

The population genetics and evolutionary epidemiology of RNA viruses

The authors discuss the main mechanisms of RNA virus evolution — mutation, recombination, natural selection, genetic drift and migration, and how these interact to shape the genetic structure of populations. The quasispecies model of RNA virus evolution is explained and the question of whether this model provides an accurate description of RNA virus evolution is discussed. Experiments that can be carried out to test the basic principles of evolutionary theory are briefly described. The authors review what such experiments have told us about virus evolution and, more widely, what these experiments have revealed in terms of general evolutionary principles. RNA viruses evolve quickly, so a detailed reconstruction of their epidemiological history can be undertaken. The authors show how epidemiological patterns of viruses result from their evolution at two different levels: within individual hosts (and vectors) and among hosts at the population level. Using several examples, including HIV and SARS, the authors describe how studying RNA virus evolution could be used to understand virus emergence. Finally, the important topics of the evolution of virulence and resistance to drugs are discussed.

RNA viruses are ubiquitous intracellular parasites that are responsible for many emerging diseases, including AIDS and SARS. Here, we discuss the principal mechanisms of RNA virus evolution and highlight areas where future research is required. The rapidity of sequence change in RNA viruses means that they are useful experimental models for the study of evolution in general and it enables us to watch them change in 'real time', and retrace the spread through populations with molecular phylogenies. An understanding of the mechanisms of RNA virus sequence change is also crucial to predicting important aspects of their emergence and long-term evolution.

Rates of deleterious mutation and the evolution of sex in Caenorhabditis

A variety of models propose that the accumulation of deleterious mutations plays an important role in the evolution of breeding systems. These models make predictions regarding the relative rates of protein evolution and deleterious mutation in taxa with contrasting modes of reproduction. Here we compare available coding sequences from one obligately outcrossing and two primarily selfing species of Caenorhabditis to explore the potential for mutational models to explain the evolution of breeding system in this clade. If deleterious mutations interact synergistically, the mutational deterministic hypothesis predicts that a high genomic deleterious mutation rate (U) will offset the reproductive disadvantage of outcrossing relative to asexual or selfing reproduction. Therefore, C. elegans and C. briggsae (both largely selfing) should both exhibit lower rates of deleterious mutation than the obligately outcrossing relative C. remanei. Using a comparative approach, we estimate U to be equivalent (and < 1) among all three related species. Stochastic mutational models, Muller's ratchet and Hill-Robertson interference, are expected to cause reductions in the effective population size in species that rarely outcross, thereby allowing deleterious mutations to accumulate at an elevated rate. We find only limited support for more rapid molecular evolution in selfing lineages. Overall, our analyses indicate that the evolution of breeding system in this group is unlikely to be explained solely by available mutational models.

Testing for epistasis between deleterious mutations in a parasitoid wasp

Determining the way in which deleterious mutations interact to effect fitness is crucial to numerous areas in evolutionary biology. For example, if each additional mutation leads to a greater decrease in log fitness than the last, termed synergistic epistasis, then sex and recombination provide an advantage because they enable deleterious mutations to be eliminated more efficiently. However, there is a severe shortage of relevant empirical data, especially of the form that can help test mutational explanations for the widespread occurrence of sex. Here, we test for epistasis in the parasitic wasp Nasonia vitripennis, examining the fitness consequences of chemically induced deleterious mutations. We examine two components of fitness, both of which are thought to be important in natural populations of parasitic wasps: longevity and egg production. Our results show synergistic epistasis for longevity, but not for egg production.

Evolution experiments with microorganisms: the dynamics and genetic bases of adaptation

Microorganisms have been mutating and evolving on Earth for billions of years. Now, a field of research has developed around the idea of using microorganisms to study evolution in action. Controlled and replicated experiments are using viruses, bacteria and yeast to investigate how their genomes and phenotypic properties evolve over hundreds and even thousands of generations. Here, we examine the dynamics of evolutionary adaptation, the genetic bases of adaptation, tradeoffs and the environmental specificity of adaptation, the origin and evolutionary consequences of mutators, and the process of drift decay in very small populations.

Compensatory mutations cause excess of antagonistic epistasis in RNA secondary structure folding

BACKGROUND

The rate at which fitness declines as an organism's genome accumulates random mutations is an important variable in several evolutionary theories. At an intuitive level, it might seem natural that random mutations should tend to interact synergistically, such that the rate of mean fitness decline accelerates as the number of random mutations is increased. However, in a number of recent studies, a prevalence of antagonistic epistasis (the tendency of multiple mutations to have a mitigating rather than reinforcing effect) has been observed.

RESULTS

We studied in silico the net amount and form of epistatic interactions in RNA secondary structure folding by measuring the fraction of neutral mutants as a function of mutational distance d. We found a clear prevalence of antagonistic epistasis in RNA secondary structure folding. By relating the fraction of neutral mutants at distance d to the average neutrality at distance d, we showed that this prevalence derives from the existence of many compensatory mutations at larger mutational distances.

CONCLUSIONS

Our findings imply that the average direction of epistasis in simple fitness landscapes is directly related to the density with which fitness peaks are distributed in these landscapes.

Experimental tests of the adaptive significance of sexual recombination

Numerous theories have been proposed to explain the advantages of sexual recombination the exchange of hereditary material between different genomes or homologous chromosomes. Many of these candidate benefits have been evaluated in controlled laboratory experiments, which, collectively, strongly indicate that sexual recombination provides important long-term advantages.

Distribution of spontaneous mutants and inferences about the replication mode of the RNA bacteriophage phi6

When a parent virus replicates inside its host, it must first use its own genome as the template for replication. However, once progeny genomes are produced, the progeny can in turn act as templates. Depending on whether the progeny genomes become templates, the distribution of mutants produced by an infection varies greatly. While information on the distribution is important for many population genetic models, it is also useful for inferring the replication mode of a virus. We have analyzed the distribution of mutants emerging from single bursts in the RNA bacteriophage phi6 and find that the distribution closely matches a Poisson distribution. The match suggests that replication in this bacteriophage is effectively by a stamping machine model in which the parental genome is the main template used for replication. However, because the distribution deviates slightly from a Poisson distribution, the stamping machine is not perfect and some progeny genomes must replicate. By fitting our data to a replication model in which the progeny genomes become replicative at a given rate or probability per round of replication, we estimated the rate to be very low and on the on the order of 10(-4). We discuss whether different replication modes may confer an adaptive advantage to viruses.

Soluble model for the accumulation of mutations in asexual populations

Finite asexual populations can accumulate an increasing number of deleterious mutations by a process known as Muller's ratchet, which consists of successive losses of the fittest or least-loaded classes of individuals in the population. We present here a simplified theoretical framework to describe the serial bottleneck passages setup used in experiments to demonstrate the decrease of the population mean fitness due to the operation of ratchet. In particular, we calculate the expected time between consecutive clicks of the ratchet and derive expressions relating the moments of the mean fitness distribution to the mutation and selection parameters.

Transmission bottlenecks and the evolution of fitness in rapidly evolving RNA viruses

We explored the evolutionary importance of two factors in the adaptation of RNA viruses to their cellular hosts, size of viral inoculum used to initiate a new infection, and mode of transmission (horizontal versus vertical). Transmission bottlenecks should occur in natural populations of viruses and their profound effects on viral adaptation have been previously documented. However, the role of transmission mode has not received the same attention. Here we used a factorial experimental design to test the combined effects of inoculum (bottleneck) size and mode of transmission in evolution of vesicular stomatitis virus (VSV) in tissue culture, and compared our results to the predictions of a recent theoretical model. Our data were in accord with basic genetic principles concerning the balance between mutation, selection and genetic drift. In particular, attenuation of vertically transmitted viruses was a consequence of the random accumulation of deleterious mutations, whereas horizontally transmitted viruses experiencing similar bottlenecks did not suffer the same fitness losses because effective bottleneck size was actually determined by the number of host individuals. In addition, high levels of viral fitness in horizontally transmitted populations were explained by competition among viral variants.

The evolution of RNA viruses: A population genetics view

RNA viruses are excellent experimental models for studying evolution under the theoretical framework of population genetics. For a proper justification of this thesis we have introduced some properties of RNA viruses that are relevant for studying evolution. On the other hand, population genetics is a reductionistic theory of evolution. It does not consider or make simplistic assumptions on the transformation laws within and between genotypic and phenotypic spaces. However, such laws are minimized in the case of RNA viruses because the phenotypic space maps onto the genotypic space in a much more linear way than on higher DNA-based organisms. Under experimental conditions, we have tested the role of deleterious and beneficial mutations in the degree of adaptation of vesicular stomatitis virus (VSV), a nonsegmented virus of negative strand. We also have studied how effective population size, initial genetic variability in populations, and environmental heterogeneity shapes the impact of mutations in the evolution of vesicular stomatitis virus. Finally, in an integrative attempt, we discuss pros and cons of the quasispecies theory compared with classic population genetics models for haploid organisms to explain the evolution of RNA viruses.