Reproducible nonlinear population dynamics and critical points during replicative competitions of RNA virus quasispecies

Culicoides-Specific Fitness Increase of Vesicular Stomatitis Virus in Insect-to-Insect Infections

Vesicular stomatitis virus (VSV) is an arthropod-borne virus affecting livestock. In the United States, sporadic outbreaks result in significant economic losses. During epizootics, Culicoides biting midges are biological vectors and key to the geographic expansion of outbreaks. Additionally, Culicoides may play a role in VSV overwintering because females and males are capable of highly efficient venereal transmission, despite their relatively low virus titers. We hypothesized that VSV propagated within a midge has increased fitness for subsequent midge infections. To evaluate the potential host-specific fitness increase, we propagated three viral isolates of VSV in porcine skin fibroblasts and Culicoides cell lines. We then evaluated the viral infection dynamics of the different cell-source groups in Culicoides sonorensis. Our results indicate that both mammalian- and insect-derived VSV replicate well in midges inoculated via intrathoracic injection, thereby bypassing the midgut barriers. However, when the virus was required to infect and escape the midgut barrier to disseminate after oral acquisition, the insect-derived viruses had significantly higher titers, infection, and dissemination rates than mammalian-derived viruses. Our research suggests that VSV replication in Culicoides cells increases viral fitness, facilitating midge-to-midge transmission and subsequent replication, and further highlights the significance of Culicoides midges in VSV maintenance and transmission dynamics.

Phase transitions in evolutionary dynamics

Sharp changes in state, such as transitions from survival to extinction, are hallmarks of evolutionary dynamics in biological systems. These transitions can be explored using the techniques of statistical physics and the physics of nonlinear and complex systems. For example, a survival-to-extinction transition can be characterized as a non-equilibrium phase transition to an absorbing state. Here, we review the literature on phase transitions in evolutionary dynamics. We discuss directed percolation transitions in cellular automata and evolutionary models, and models that diverge from the directed percolation universality class. We explore in detail an example of an absorbing phase transition in an agent-based model of evolutionary dynamics, including previously unpublished data demonstrating similarity to, but also divergence from, directed percolation, as well as evidence for phase transition behavior at multiple levels of the model system's evolutionary structure. We discuss phase transition models of the error catastrophe in RNA virus dynamics and phase transition models for transition from chemistry to biochemistry, i.e., the origin of life. We conclude with a review of phase transition dynamics in models of natural selection, discuss the possible role of phase transitions in unraveling fundamental unresolved questions regarding multilevel selection and the major evolutionary transitions, and assess the future outlook for phase transitions in the investigation of evolutionary dynamics.

Next-Generation Sequencing for Confronting Virus Pandemics

Virus pandemics have happened, are happening and will happen again. In recent decades, the rate of zoonotic viral spillover into humans has accelerated, mirroring the expansion of our global footprint and travel network, including the expansion of viral vectors and the destruction of natural spaces, bringing humans closer to wild animals. Once viral cross-species transmission to humans occurs, transmission cannot be stopped by cement walls but by developing barriers based on knowledge that can prevent or reduce the effects of any pandemic. Controlling a local transmission affecting few individuals is more efficient that confronting a community outbreak in which infections cannot be traced. Genetic detection, identification, and characterization of infectious agents using next-generation sequencing (NGS) has been proven to be a powerful tool allowing for the development of fast PCR-based molecular assays, the rapid development of vaccines based on mRNA and DNA, the identification of outbreaks, transmission dynamics and spill-over events, the detection of new variants and treatment of vaccine resistance mutations, the development of direct-acting antiviral drugs, the discovery of relevant minority variants to improve knowledge of the viral life cycle, strengths and weaknesses, the potential for becoming dominant to take appropriate preventive measures, and the discovery of new routes of viral transmission.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

Genomic and immunogenic changes of Piscine novirhabdovirus (Viral Hemorrhagic Septicemia Virus) over its evolutionary history in the Laurentian Great Lakes

A unique and highly virulent subgenogroup (-IVb) of Piscine novirhabdovirus, also known as Viral Hemorrhagic Septicemia Virus (VHSV), suddenly appeared in the Laurentian Great Lakes, causing large mortality outbreaks in 2005 and 2006, and affecting >32 freshwater fish species. Periods of apparent dormancy have punctuated smaller and more geographically-restricted outbreaks in 2007, 2008, and 2017. In this study, we conduct the largest whole genome sequencing analysis of VHSV-IVb to date, evaluating its evolutionary changes from 48 isolates in relation to immunogenicity in cell culture. Our investigation compares genomic and genetic variation, selection, and rates of sequence changes in VHSV-IVb, in relation to other VHSV genogroups (VHSV-I, VHSV-II, VHSV-III, and VHSV-IVa) and with other Novirhabdoviruses. Results show that the VHSV-IVb isolates we sequenced contain 253 SNPs (2.3% of the total 11,158 nucleotides) across their entire genomes, with 85 (33.6%) of them being non-synonymous. The most substitutions occurred in the non-coding region (NCDS; 4.3%), followed by the Nv- (3.8%), and M- (2.8%) genes. Proportionally more M-gene substitutions encoded amino acid changes (52.9%), followed by the Nv- (50.0%), G- (48.6%), N- (35.7%) and L- (23.1%) genes. Among VHSV genogroups and subgenogroups, VHSV-IVa from the northeastern Pacific Ocean has shown the fastest substitution rate (2.01x10^-3), followed by VHSV-IVb (6.64x10^-5) and by the VHSV-I, -II and-III genogroups from Europe (4.09x10^-5). A 2016 gizzard shad (Dorosoma cepedianum) from Lake Erie possessed the most divergent VHSV-IVb sequence. The in vitro immunogenicity analysis of that sample displayed reduced virulence (as did the other samples from 2016), in comparison to the original VHSV-IVb isolate (which had been traced back to 2003, as an origin date). The 2016 isolates that we tested induced milder impacts on fish host cell innate antiviral responses, suggesting altered phenotypic effects. In conclusion, our overall findings indicate that VHSV-IVb has undergone continued sequence change and a trend to lower virulence over its evolutionary history (2003 through present-day), which may facilitate its long-term persistence in fish host populations.

Evolutionary trajectory of fish Piscine novirhabdovirus (=Viral Hemorrhagic Septicemia Virus) across its Laurentian Great Lakes history: Spatial and temporal diversification

Piscine novirhabdovirus = Viral Hemorrhagic Septicemia Virus (VHSV) first appeared in the Laurentian Great Lakes with large outbreaks from 2005 to 2006, as a new and novel RNA rhabdovirus subgenogroup (IVb) that killed >30 fish species. Interlude periods punctuated smaller more localized outbreaks in 2007, 2010, and 2017, although some fishes tested positive in the intervals. There have not been reports of outbreaks or positives from 2018, 2019, or 2020. Here, we employ a combined population genetics and phylogenetic approach to evaluate spatial and temporal evolutionary trajectory on its G-gene sequence variation, in comparison with whole-genome sequences (11,083 bp) from a subset of 44 individual isolates (including 40 newly sequenced ones). Our results show that IVb (N = 184 individual fish isolates) diversified into 36 G-gene haplotypes from 2003 to 2017, stemming from two originals ("a" and "b"). G-gene haplotypes "a" and "b" differed by just one synonymous single-nucleotide polymorphism (SNP) substitution, remained the most abundant until 2011, then disappeared. Group "a" descendants (14 haplotypes) remained most prevalent in the Upper and Central Great Lakes, with eight (51%) having nonsynonymous substitutions. Group "b" descendants primarily have occurred in the Lower Great Lakes, including 22 haplotypes, of which 15 (68%) contained nonsynonymous changes. Evolutionary patterns of the whole-genome sequences (which had 34 haplotypes among 44 isolates) appear congruent with those from the G-gene. Virus populations significantly diverged among the Upper, Central, and Lower Great Lakes, diversifying over time. Spatial divergence was apparent in the overall patterns of nucleotide substitutions, while amino acid changes increased temporally. VHSV-IVb thus significantly differentiated across its less than two decades in the Great Lakes, accompanied by declining outbreaks and virulence. Continuing diversification likely allowed the virus to persist at low levels in resident fish populations, and may facilitate its potential for further and future spread to new habitats and nonacclimated hosts.

Virus population dynamics examined with experimental model systems

Experimental evolution permits exploring the effect of controlled environmental variables in virus evolution. Several designs in cell culture and in vivo have established basic concepts that can assist in the interpretation of evolutionary events in the field. Important information has come from cytolytic and persistent infections in cell culture that have unveiled the power of virus-cell coevolution in virus and cell diversification. Equally informative are comparisons of the response of viral populations when subjected to different passage régimens. In particular, plaque-to-plaque transfers in cell culture have revealed unusual genotypes and phenotypes that populate minority layers of viral quasispecies. Some of these viruses display properties that contradict features established in virology textbooks. Several hypotheses and principles of population genetics have found experimental confirmation in experimental designs with viruses. The possibilities of using experimental evolution to understand virus behavior are still largely unexploited.

Vesicular Stomatitis Virus Transmission: A Comparison of Incriminated Vectors

Vesicular stomatitis (VS) is a viral disease of veterinary importance, enzootic in tropical and subtropical regions of the Americas. In the U.S., VS produces devastating economic losses, particularly in the southwestern states where the outbreaks display an occurrence pattern of 10-year intervals. To date, the mechanisms of the geographic spread and maintenance cycles during epizootics remain unclear. This is due, in part, to the fact that VS epidemiology has a complex of variables to consider, including a broad range of vertebrate hosts, multiple routes of transmission, and an extensive diversity of suspected vector species acting as both mechanical and biological vectors. Infection and viral progression within vector species are highly influenced by virus serotype, as well as environmental factors, including temperature and seasonality; however, the mechanisms of viral transmission, including non-conventional pathways, are yet to be fully studied. Here, we review VS epidemiology and transmission mechanisms, with comparisons of transmission evidence for the four most incriminated hematophagous dipteran taxa: Aedes mosquitoes, Lutzomyia sand flies, Simulium black flies, and Culicoides biting midges.

Quasispecies and virus

Quasispecies theory has been instrumental in the understanding of RNA virus population dynamics because it considered for the first time mutation as an integral part of the replication process. The key influences of quasispecies theory on experimental virology have been: (1) to disclose the mutant spectrum nature of viral populations and to evaluate its consequences; (2) to unveil collective properties of genome ensembles that can render a mutant spectrum a unit of selection; and (3) to identify new vulnerability points of pathogenic RNA viruses on three fronts: the need to apply multiple selective constraints (in the form of drug combinations) to minimize selection of treatment-escape variants, to translate the error threshold concept into antiviral designs, and to construct attenuated vaccine viruses through alterations of viral polymerase copying fidelity or through displacements of viral genomes towards unfavorable regions of sequence space. These three major influences on the understanding of viral pathogens preceded extensions of quasispecies to non-viral systems such as bacterial and tumor cell collectivities and prions. These developments are summarized here.

Differences in adaptive dynamics determine the success of virus variants that propagate together

Virus fitness is a complex parameter that results from the interaction of virus-specific characters (e.g. intracellular growth rate, adsorption rate, virion extracellular stability, and tolerance to mutations) with others that depend on the underlying fitness landscape and the internal structure of the whole population. Individual mutants usually have lower fitness values than the complex population from which they come from. When they are propagated and allowed to attain large population sizes for a sufficiently long time, they approach mutation-selection equilibrium with the concomitant fitness gains. The optimization process follows dynamics that vary among viruses, likely due to differences in any of the parameters that determine fitness values. As a consequence, when different mutants spread together, the number of generations experienced by each of them prior to co-propagation may determine its particular fate. In this work we attempt a clarification of the effect of different levels of population diversity in the outcome of competition dynamics. To this end, we analyze the behavior of two mutants of the RNA bacteriophage Qβ that co-propagate with the wild-type virus. When both competitor viruses are clonal, the mutants rapidly outcompete the wild type. However, the outcome in competitions performed with partially optimized virus populations depends on the distance of the competitors to their clonal origin. We also implement a theoretical population dynamics model that describes the evolution of a heterogeneous population of individuals, each characterized by a fitness value, subjected to subsequent cycles of replication and mutation. The experimental results are explained in the framework of our theoretical model under two non-excluding, likely complementary assumptions: (1) The relative advantage of both competitors changes as populations approach mutation-selection equilibrium, as a consequence of differences in their growth rates and (2) one of the competitors is more robust to mutations than the other. The main conclusion is that the nearness of an RNA virus population to mutation-selection equilibrium is a key factor determining the fate of particular mutants arising during replication.

Virus Population Dynamics Examined with Experimental Model Systems

Experimental evolution permits exploring the effect of controlled environmental variables in virus evolution. Several designs in cell culture and in vivo have established basic concepts that can assist in the interpretation of evolutionary events in the field. Important information has come from cytolytic and persistent infections in cell culture that have unveiled the power of virus-cell coevolution in virus and cell diversification. Equally informative are comparisons of the response of viral populations when subjected to different passage régimes. In particular, plaque-to-plaque transfers in cell culture have revealed unusual genotypes and phenotypes that populate minority layers of viral quasispecies. Some of these viruses display properties that contradict features established in virology textbooks. Several hypotheses and principles of population genetics have found experimental confirmation in experimental designs with viruses. The possibilities of using experimental evolution to understand virus behavior are still largely unexploited.

Quasispecies dynamics in hepatitis C liver transplant recipients receiving grafts from hepatitis C virus infected donors

The allocation of liver grafts from hepatitis C virus (HCV)-positive donors in HCV-infected liver transplant (LT) recipients leads to infection with two different viral populations. In a previous study, we examined quasispecies dynamics during reinfection by clonal sequencing, which did not allow an accurate characterization of coexistence and competition events. To overcome this limitation, here we used deep-sequencing analysis of a fragment of the HCV NS5B gene in six HCV-infected LT recipients who received HCV-infected grafts. Successive expansions and contractions of quasispecies complexity were observed, evolving in all cases towards a more homogeneous population. The population that became dominant was the one displaying the highest mutant spectrum complexity. In four patients, coexistence of minority mutants, derived from the donor or the recipient, were detected. In conclusion, our study shows that, during reinfection with a different HCV strain in LT recipients, the viral population with the highest diversity always becomes dominant.

Determination of the complete genome sequence of infectious hematopoietic necrosis virus (IHNV) Ch20101008 and viral molecular evolution in China

This study determined the complete genomic sequence of the infectious hematopoietic necrosis virus (IHNV) strain Ch20101008 isolated from farmed brook trout (Salvelinus fontinalis) that died from a disease caused by the virus in northeast China. The sequence was determined from 10 overlapping clones obtained through RT-PCR amplification. The whole genome length of Ch20101008 comprised 11,129 nucleotides (nt), and the overall organization was typical of that observed for all other IHNV strains. The phylogenetic analysis results of the 65 IHNV glycoprotein genes and 47 IHNV partial nucleoprotein genes presented five major genogroups (J, U, L, E and M). The J genogroup included the J Nagano and J Shizuoka subgroups. The IHNV Ch20101008 strain belonged to the J Nagano subgroup of the J genogroup and was significantly related to other Chinese IHNV strains. All Chinese IHNV isolates are monophyletic, with a recent common ancestor, except for the BjLL strain. The N, P, M, G, NV and L genes of Ch20101008 were compared with the available IHNV sequences in GenBank. The results indicated that 198 nt were substituted, 53 of which exhibited amino acid change in the Ch20101008 genome. An adenine nucleotide deletion was found at position 4959 of the 5' UTR of the L gene. In the G gene, specific nucleotides and amino acid variations of the Chinese IHNV strains were observed when compared with 23 isolates from other countries. Of the 15 nucleotide sites that changed, seven resulted in amino acid substitution. The data further demonstrated that the J genogroup IHNV was introduced to and evolved in salmon farm environments in China.

A replicase of Potato virus X acts as the resistance-breaking determinant for JAX1-mediated resistance

Lectin-mediated resistance (LMR) has been suggested to comprise an uncharacterized branch of antiviral plant innate immunity. To unveil the feature of resistance conferred by jacalin-type lectin required for potexvirus resistance 1 (JAX1), a recently isolated LMR gene against potexviruses, we analyzed the resistance-breaking variants to find the viral component involved in resistance. We employed grafting-mediated inoculation, a high-pressure virus inoculation method, to obtain Potato virus X (PVX) variants that can overcome JAX1-mediated resistance. Whole-genome sequencing of the variants suggested that a single amino acid in the methyl transferase domain of the replicase encoded by PVX is responsible for this resistance-breaking property. Reintroduction of the amino-acid substitution to avirulent wild-type PVX was sufficient to overcome the JAX1-mediated resistance. These results suggest that viral replicase is involved in JAX1-mediated resistance. The residue that determines the resistance-breaking properties was highly conserved among potexviruses, suggesting a general role of the residue in potexvirus-JAX1 interactions.

Viral quasispecies evolution

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

Specific and nonspecific host adaptation during arboviral experimental evolution

During the past decade or so, there has been a substantial body of work to dissect arboviral evolution and to develop models of adaptation during host switching. Regardless of what species serve as host or vectors, and of the geographic distribution and the mechanisms of replication, arboviruses tend to have slow evolutionary rates in nature. The hypothesis that this is the result of replication in the disparate environments provided by host and vector did not receive solid experimental support in any of the many viral species tested. Instead, it seems that from the virus's point of view, either the two environments are sufficiently similar or one of the environments so dominates viral evolution that there is tolerance for suboptimal adaptation to the other environment. Replication in alternating environments has an unexpected cost in that there is decreased genetic variance that translates into a compromised adaptability for bypassed environments. Arboviruses under strong and continuous positive selection may have unusual patterns of genomic changes, with few or no mutations accumulated in the consensus sequence or with dN/dS values typically consistent with random drift in DNA-based organisms.

Evolution and biogeography of an emerging quasispecies: diversity patterns of the fish Viral Hemorrhagic Septicemia virus (VHSv)

Viral Hemorrhagic Septicemia virus (VHSv) is an RNA rhabdovirus that causes one of the most important finfish diseases, affecting over 70 marine and freshwater species. It was discovered in European cultured fish in 1938 and since has been described across the Northern Hemisphere. Four strains and several substrains have been hypothesized, whose phylogenetic relationships and evolutionary radiation are evaluated here in light of a quasispecies model, including an in-depth analysis of the novel and especially virulent new substrain (IVb) that first appeared in the North American Laurentian Great Lakes in 2003. We analyze the evolutionary patterns, genetic diversity, and biogeography of VHSv using all available RNA sequences from the glycoprotein (G), nucleoprotein (N), and non-virion (Nv) genes, with Maximum Likelihood and bayesian approaches. Results indicate that the G gene evolves at an estimated rate of μ=2.58×10(-4) nucleotide substitutions per site per year, the N gene at μ=4.26×10(-4), and Nv fastest at μ=1.25×10(-3). Phylogenetic trees from the three genes largely are congruent, distinguishing strains I-IV as reciprocally monophyletic with high bootstrap and posterior probability support. VHSv appears to have originated from a marine ancestor in the North Atlantic Ocean, diverging into two primary clades: strain IV in North America (the Northwestern Atlantic Ocean), and strains I-III in the Northeastern Atlantic region (Europe). Strain II may comprise the basal group of the latter clade and diverged in Baltic Sea estuarine waters; strains I and III appear to be sister groups (according to the G and Nv genes), with the former mostly in European freshwaters and the latter in North Sea marine/estuarine waters. Strain IV is differentiated into three monophyletic substrains, with IVa infecting Northeastern Pacific salmonids and many marine fishes (with 44 unique G gene haplotypes), IVb endemic to the freshwater Great Lakes (11 haplotypes), and a newly-designated IVc in marine/estuarine North Atlantic waters (five haplotypes). Two separate substrains independently appeared in the Northwestern Pacific region (Asia) in 1996, with Ib originating from the west and IVa from the east. Our results depict an evolutionary history of relatively rapid population diversifications in star-like patterns, following a quasispecies model. This study provides a baseline for future tracking of VHSv spread and interpreting its evolutionary diversification pathways.

Mutant spectra in virus behavior

RNA viruses replicate as complex mutant spectra, also termed ‘mutant clouds’, known as viral quasispecies. While this is a widely observed viral population structure, it is less known that a number of biologically relevant features of this important group of viral pathogens depend on (or are strongly influenced by) the complexity and composition of mutant spectra. Among them, fitness increase or decrease depending on intrapopulation complementation or interference, selection triggered by memory genomes, pathogenic potential of viruses, disease evolution and the response to antiviral treatments. Quasispecies represent the recognition of complex behavior in viruses, and it is an oversimplification to equate such a population structure with the classic polymorphism of population biology. Darwinian principles acting on genome collectivities that replicate with high error rates provide a unique population structure prone to flexible and largely unpredictable behavior.

Mutational effects and population dynamics during viral adaptation challenge current models

Adaptation in haploid organisms has been extensively modeled but little tested. Using a microvirid bacteriophage (ID11), we conducted serial passage adaptations at two bottleneck sizes (10(4) and 10(6)), followed by fitness assays and whole-genome sequencing of 631 individual isolates. Extensive genetic variation was observed including 22 beneficial, several nearly neutral, and several deleterious mutations. In the three large bottleneck lines, up to eight different haplotypes were observed in samples of 23 genomes from the final time point. The small bottleneck lines were less diverse. The small bottleneck lines appeared to operate near the transition between isolated selective sweeps and conditions of complex dynamics (e.g., clonal interference). The large bottleneck lines exhibited extensive interference and less stochasticity, with multiple beneficial mutations establishing on a variety of backgrounds. Several leapfrog events occurred. The distribution of first-step adaptive mutations differed significantly from the distribution of second-steps, and a surprisingly large number of second-step beneficial mutations were observed on a highly fit first-step background. Furthermore, few first-step mutations appeared as second-steps and second-steps had substantially smaller selection coefficients. Collectively, the results indicate that the fitness landscape falls between the extremes of smooth and fully uncorrelated, violating the assumptions of many current mutational landscape models.

Genomic evolution of vesicular stomatitis virus strains with differences in adaptability

Virus strains with a history of repeated genetic bottlenecks frequently show a diminished ability to adapt compared to strains that do not have such a history. These differences in adaptability suggest differences in either the rate at which beneficial mutations are produced, the effects of beneficial mutations, or both. We tested these possibilities by subjecting four populations (two controls and two mutants with lower adaptabilities) to multiple replicas of a regimen of positive selection and then determining the fitnesses of the progeny through time and the changes in the consensus, full-length sequences of 56 genomes. We observed that at a given number of passages, the overall fitness gains observed for control populations were larger than fitness gains in mutant populations. However, these changes did not correlate with differences in the numbers of mutations accumulated in the two types of genomes. This result is consistent with beneficial mutations having a lower beneficial effect on mutant strains. Despite the overall fitness differences, some replicas of one mutant strain at passage 50 showed fitness increases similar to those observed for the wild type. We hypothesized that these evolved, high-fitness mutants may have a lower robustness than evolved, high-fitness controls. Robustness is the ability of a virus to avoid phenotypic changes in the face of mutation. We confirmed our hypothesis in mutation-accumulation experiments that showed a normalized fitness loss that was significantly larger in mutant bottlenecked populations than in control populations.

Unfinished stories on viral quasispecies and Darwinian views of evolution

Experimental evidence that RNA virus populations consist of distributions of mutant genomes, termed quasispecies, was first published 31 years ago. This work provided the earliest experimental support for a theory to explain a system that replicated with limited fidelity and to understand the self-organization and adaptability of early life forms on Earth. High mutation rates and quasispecies dynamics of RNA viruses are intimately related to both viral disease and antiviral treatment strategies. Moreover, the quasispecies concept is being applied to other biological systems such as cancer research in which cellular mutant spectra can be also detected. This review addresses some of the unanswered questions regarding viral and theoretical quasispecies concepts as well as more practical aspects concerning resistance to antiviral treatments and pathogenesis.

Incongruent fitness landscapes, not tradeoffs, dominate the adaptation of vesicular stomatitis virus to novel host types

Host radiation refers to the ability of parasites to adapt to new environments and expand or change their niches. Adaptation to one specific environment may involve a loss in adaptation to a second environment. Thus, fitness costs may impose limits to niche expansion and constitute the cost of specialization. Several reports have addressed the cost of host radiation in vesicular stomatitis virus (VSV), but in some cases the experimental setup may have resulted in the overestimation of fitness costs. To clarify this issue, experiments were carried out in which a reference strain of VSV was allowed to adapt to HeLa, MDCK and BHK-21 cells, and to a regime of alternation between HeLa and Madin–Darby canine kidney (MDCK) cells. Measurement of viral fitness on each cell type showed that most virus populations behaved as generalists, and increased in fitness in all environments. Tradeoffs, where a fitness increase in one environment led to a fitness decrease in another environment, were rare. These results highlight the importance of using appropriate methods to measure fitness in evolved virus populations, and provide further support to a model of evolutionary dynamics in which costs due to incongruent landscapes provided by different environments are more common than tradeoffs.

No evidence of selection for mutational robustness during lethal mutagenesis of lymphocytic choriomeningitis virus

Lethal mutagenesis is a transition towards virus extinction mediated by enhanced mutation rates during viral genome replication. Theoretical studies suggest that viruses can evolve towards regions of their fitness landscapes at which they display resistance to the deleterious effects of mutations. It has been suggested that such mutational robustness could jeopardize lethal mutagenesis. We have used the Arenavirus lymphocytic choriomeningitis virus (LCMV) to explore whether treatment with the mutagenic base analogue 5-fluorouracil (FU) selected for viral populations displaying resistance to lethal mutagenesis. Neither average LCMV populations with a history of FU mutagenesis, nor individual biological LCMV clones derived from those populations, displayed any resistance to lethal mutagenesis by FU. They were as sensitive to FU-induced extinction as LCMV populations and clones treated in parallel, but without a history of FU mutagenesis. Current evidence of the molecular events affecting quasispecies dynamics suggests that it is unlikely that a viral population can acquire mutational robustness under the increased mutation rates associated with mutagenic treatments. We consider mechanisms by which viruses could escape extinction by lethal mutagenesis, and provide evidence that mutational robustness is unlikely to be one of them.

Persistence of foot-and-mouth disease virus in cell culture revisited: implications for contingency in evolution

If we could rewind the tape of evolution and play it again, would it turn out to be similar to or different from what we know? Obviously, this key question can only be addressed by fragmentary experimental approaches. Twenty-two years ago, we described the establishment of BHK-21 cells persistently infected with foot-and-mouth disease virus (FMDV), a system that displayed as its major biological feature a coevolution of the cells and the resident virus in the course of persistence. Now we report the establishment of two persistently infected cell lines in parallel, starting with the same clones of FMDV and BHK-21 cells used 22 years ago. We have asked whether the evolution of the two newly established cell lines and of the earlier cell line would be similar or different. The main conclusions of the study are: (i) the basic behaviour characterized by virus–cell coevolution is similar in the three carrier cell lines, despite differences in some genetic alterations of FMDV; (ii) a strikingly parallel behaviour has been observed with the two newly established cell lines passaged in parallel, unveiling a deterministic virus behaviour during persistence; and (iii) selective RT-PCR amplifications have detected imbalances in the proportion of positive- versus negative-strand viral RNA, mediated by both viral and cellular factors. The results confirm coevolution of cells and virus as a major and reproducible feature of FMDV persistence in cell culture, and suggest that rapidly evolving viruses may constitute adequate test systems to probe the influence of historical contingency on evolutionary events.

Analysis of the genetic and the corresponding antigenic variability of the VP1 3' end of ECHO virus type 11 and ECHO virus type 30

The enteroviruses (EV), RNA viruses belonging to the Picornaviridae family, have a high genetic variability due to the absence of the efficient proofreading and post replicative repair activities associated with the RNA polymerase. In the present work, we studied the genetic and the antigenic variability of ECHO virus types 11 (E11) and 30 (E30), which are the most isolated echoviruses serotypes in clinical and environmental samples. We established on the 3' end of the VP1 gene, consensus sequences of E11 and E30 by alignment of 67 E11 and 247 E30 published sequences in GenBank. Our results of sequences comparison showed that the majority of the mutational sites are situated on the third nucleotide of the codon. These mutations were without consequence on the antigenic sequences of the VP1 protein. Thus, E11 and E30 have a high genetic variability (1/3 of the nucleotides are variable), but a relative antigenic conservation. The analysis of the intertypic antigenic variability between E11 and E30 was obtained by the alignment of the corresponding amino acids sequences relative to the N-terminal part of the VP1 protein. Two discriminating parts were highlighted, probably representing antigenic sites for neutralisation antibodies.

Population bottlenecks in quasispecies dynamics

The characteristics of natural populations result from different stochastic and deterministic processes that include reproduction with error, selection, and genetic drift. In particular, population fluctuations constitute a stochastic process that may play a very relevant role in shaping the structure of populations. For example, it is expected that small asexual populations will accumulate mutations at a higher rate than larger ones. As a consequence, in any population the fixation of mutations is accelerated when environmental conditions cause population bottlenecks. Bottlenecks have been relatively frequent in the history of life and it is generally accepted that they are highly relevant for speciation. Although population bottlenecks can occur in any species, their effects are more noticeable in organisms that form large and heterogeneous populations, such as RNA viral quasispecies. Bottlenecks can also positively select and isolate particles that still keep the ability to infect cells from a disorganized population created by crossing the error threshold.

Effect of bottlenecking on evolution of the nonstructural protein 3 gene of hepatitis C virus during sexually transmitted acute resolving infection

Sexual partners of patients infected with the hepatitis C virus (HCV) often have detectable HCV-specific T-cell responses in the absence of seroconversion, suggesting unapparent, spontaneously resolving infection. To determine whether differences in the evolutionary potential of bottlenecked inoculum may explain the low rate of HCV persistence after sexual exposure, we have investigated changes in the entire HCV nonstructural 3 (NS3) gene over time in a chronic carrier and compared his viral quasispecies with that of the acute-phase isolate of his sexual partner, who developed acute resolving hepatitis C. The overall rate of accumulation of mutations, estimated by regression analysis of six consecutive consensus NS3 sequences over 8 years, was 1.5 x 10(-3) mutations per site per year, with small intersample fluctuations related to changes in environmental conditions. Comparison of quasispecies parameters in one isolate of the chronic carrier with those of the acute-phase isolate of the infected partner revealed a higher heterogeneity and lower proportion of nonsynonymous mutations in the former. All NS3 sequences from the acute-phase isolate clustered with a single sequence from the chronic isolate, despite complete HLA mismatch between the patients, suggesting bottlenecking during transmission. The low risk of viral persistence after sexual exposure to HCV may be related to the selection of a limited number of viral particles carrying a particular combination of mutations which may further limit the potential of a relatively homogeneous quasispecies to rapidly diversify and overcome the immune response of the exposed host.

Transitions in understanding of RNA viruses: a historical perspective

This chapter documents that RNA viruses have been known for over a century to be genetically variable. In recent decades, genetic and molecular analyses demonstrate that they form RNA quasispecies populations; the most rapidly mutating, highly variable and genetically versatile life forms on earth. Their enormous populations, rapid replication and extreme genetic plasticity can allow rates of evolution that exceed those of their eukaryotic host populations by millions-fold.

RNA viruses as complex adaptive systems

RNA viruses have high mutation rates and so their populations exist as dynamic and complex mutant distributions. It has been consistently observed that when challenged with a new environment, viral populations adapt following hyperbolic-like kinetics: adaptation is initially very rapid, but then slows down as fitness reaches an asymptotic value. These adaptive dynamics have been explained in terms of populations moving towards the top of peaks on rugged fitness landscapes. Fitness fluctuations of varying magnitude are observed during adaptation. Often the presence of fluctuations in the evolution of physical systems indicates some form of self-organization, or where many components of the system are simultaneously involved. Here we analyze data from several in vitro evolution experiments carried out with vesicular stomatitis virus (VSV) looking for the signature of criticality and scaling. Long-range fitness correlations have been detected during the adaptive process. We also found that the magnitude of fitness fluctuations, far from being trivial, conform to a Weibull probability distribution function, suggesting that viral adaptation belongs to a broad category of phenomena previously documented in other fields and related with emergence.

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.

Evolution of the fish rhabdovirus viral haemorrhagic septicaemia virus

Viral haemorrhagic septicaemia (VHS) caused by the rhabdovirus VHSV is economically the most important viral disease in European rainbow trout farming. Until 1989, this virus was mainly isolated from freshwater salmonids but in the last decade, it has also been isolated from an increasing number of free-living marine fish species. To study the genetic evolution of VHSV, the entire G gene from 74 isolates was analysed. VHSV from wild marine species caught in the Baltic Sea, Skagerrak, Kattegat, North Sea, and English Channel and European freshwater isolates, appeared to share a recent common ancestor. Based on the estimated nucleotide substitution rate, the ancestor of the European fresh water isolates was dated some 50 years ago. This finding fits with the initial reports in the 1950s on clinical observations of VHS in Danish freshwater rainbow trout farms. The study also indicates that European marine VHSV and the North American marine line separated approx. 500 years ago. The codon substitution rate among the freshwater VHSV isolates was found to be 2.5 times faster than among marine isolates. The data support the hypothesis of the marine environment being the original reservoir of VHSV and that the change in host range (to include rainbow trout) may have occurred several times. Virus from the marine environment will therefore continue to represent a threat to the trout aquaculture industry.

Negative effect of genetic bottlenecks on the adaptability of vesicular stomatitis virus

Muller's ratchet is a principle of evolutionary genetics describing mutant accumulation in populations that are repeatedly subjected to genetic bottleneck. The immediate effect of Muller's ratchet, overall loss of fitness, has been confirmed in several viral systems belonging to different groups. This report shows that in addition to fitness loss, genetic bottlenecks also have longer-term effects, namely changes in the capacity of viral populations to adapt. Thus, vesicular stomatitis virus strains with a history of genetic bottleneck have lower adaptability than strains maintained at relatively large population sizes. This lower adaptability is illustrated by their reduced ability to regain fitness and by their inability to outcompete wild-type populations in situations where the initial fitness of the bottlenecked mutant is the same or even higher than the initial fitness of the wild-type.

Synchronous loss of quasispecies memory in parallel viral lineages: a deterministic feature of viral quasispecies

Viral quasispecies are endowed with a memory of their past evolutionary history in the form of minority genomes of their mutant spectra. To determine the fate of memory genomes in evolving viral quasispecies, we have measured memory levels of antigenic variant of foot-and-mouth disease virus (FMDV) RED, which includes an Arg-Glu-Asp (RED) at a surface antigenic loop of the viral capsid. The RED reverted to the standard Arg-Gly-Asp (RGD), and the RED remained as memory in the evolving quasispecies. In four parallel evolutionary lineages, memory reduction followed a strikingly similar pattern, and at passage 60 memory levels were indistinguishable from those of control populations (devoid of memory). Nucleotide sequence analyses indicated that memory loss occurred synchronously despite its ultimate molecular basis being the stochastic occurrence of mutations in the evolving quasispecies. These results on the kinetics of memory levels have unveiled a deterministic feature of viral quasispecies. Molecular mechanisms that may underlie synchronous memory loss are the averaging of noise signals derived from mutational input, and constraints to genome diversification imposed by a nucleotide sequence context in the viral genome. Possible implications of the behaviour of complex, adaptive viral systems as experimental models to address primary mechanisms of neurological memory are discussed.

Contributions of vesicular stomatitis virus to the understanding of RNA virus evolution

Vesicular stomatitis virus has been a preferred system to study evolution for several decades. New approaches to antiviral treatment, such as lethal mutagenesis, stem from investigations done with VSV. Recent work has shed new light in the way we view neutrality, a fundamental concept to understand the evolutionary history of RNA viruses.

Coinfection and superinfection in RNA virus populations: a selection-mutation model

In this paper, we present a general selection-mutation model of evolution on a one-dimensional continuous fitness space. The formulation of our model includes both the classical diffusion approach to mutation process as well as an alternative approach based on an integral operator with a mutation kernel. We show that both approaches produce fundamentally equivalent results. To illustrate the suitability of our model, we focus its analytical study into its application to recent experimental studies of in vitro viral evolution. More specifically, these experiments were designed to test previous theoretical predictions regarding the effects of multiple infection dynamics (i.e., coinfection and superinfection) on the virulence of evolving viral populations. The results of these experiments, however, did not match with previous theory. By contrast, the model we present here helps to understand the underlying viral dynamics on these experiments and makes new testable predictions about the role of parameters such the time between successive infections and the growth rates of resident and invading populations.

Quasispecies Theory in Virology

Holmes and Moya claim that quasispecies is an unnecessary and misleading description of RNA virus evolution, that virologists refer to quasispecies inappropriately, and that there is little evidence of quasispecies in RNA virus evolution. They wish to look for other ideas in evolutionary biology and to set down an agenda for future research. I argue here that real virus quasispecies often differ from the theoretical quasispecies as initially formulated and that this difference does not invalidate quasispecies as a suitable theoretical framework to understand viruses at the population level.

Longitudinal evaluation of the structure of replicating and circulating hepatitis C virus quasispecies in nonprogressive chronic hepatitis C patients

In previous cross-sectional studies, we demonstrated that, in most patients with chronic hepatitis C, the composition and complexity of the circulating hepatitis C virus (HCV) population do not coincide with those of the virus replicating in the liver. In the subgroup of patients with similar complexities in both compartments, the ratio of quasispecies complexity in the liver to that in serum (liver/serum complexity ratio) of paired samples correlated with disease stage. In the present study we investigated the dynamic behavior of viral population parameters in consecutive paired liver and serum samples, obtained 3 to 6 years apart, from four chronic hepatitis C patients with persistently normal transaminases and stable liver histology. We sequenced 359 clones of a genomic fragment encompassing the E2(p7)-NS2 junction, in two consecutive liver-serum sample pairs from the four patients and in four intermediate serum samples from one of the patients. The results show that the liver/serum complexity ratio is not stable but rather fluctuates widely over time. Hence, the liver/serum complexity ratio does not identify a particular group of patients but a particular state of the infecting quasispecies. Phylogenetic analysis and signature mutation patterns showed that virtually all circulating sequences originated from sequences present in the liver specimens. The overall behavior of the circulating viral quasispecies appears to originate from changes in the relative replication kinetics of the large mutant spectrum present in the infected liver.

Contingent neutrality in competing viral populations

The replicative fitness of a genetically marked (MARM-C) population of vesicular stomatitis virus was examined in competition assays in BHK-21 cells. In standard fitness assays involving up to eight competition passages of the mixed populations, MARM-C competes equally with the wild type (wt), but very prolonged competitions always led to the wt gaining dominance over MARM-C in a very slowed, nonlinear manner (J. Quer et al., J. Mol. Biol. 264:465-471, 1996). In the present study we show that a number of quite unrelated environmental perturbations, which decreased virus replication during competitions, all led to an accelerated dominance of the wt over MARM-C. These perturbations were (i) the presence of added (or endogenously generated) defective interfering particles, (ii) the presence of the chemical mutagen 5-fluorouracil (5-FU), or (iii) an increase in temperature to 40.5 degrees C. Thus, the "neutral fitness" of the MARM-C population is contingent. We have determined the entire genomic consensus sequence of MARM-C and have identified only six mutations. Clearly, some or all of these mutations allowed the MARM-C quasispecies population to compete equally with wt in a defined constant host environment, but the period of neutrality was shortened when the environment was perturbed during competitions. Interestingly, when four passages of each population were carried out independently in the presence of 5-FU (but in the absence of competition), no significant differences were detected in the fitness changes of wt and MARM-C, nor was there a difference in their subsequent abilities to compete with each other in a standard fitness assay. We propose a model for this contingent neutrality. The conditions employed to generate the MARM-C quasispecies population selected a small number of mutations in the consensus sequence. It appears that the MARM-C quasispecies population has moved into a segment of sequence space in which the average fitness value is neutral but, under environmental stress, beneficial mutations cannot be generated rapidly enough to compete with those being generated concurrently by competing wt virus quasispecies populations.

Dynamics of dominance of a dipeptide insertion in reverse transcriptase of HIV-1 from patients subjected to prolonged therapy

A small proportion (0.8%) of individuals of a cohort of HIV-1 infected patients subjected to prolonged therapy with nucleoside analogues included a recently recognised dipeptide insertion in their RT (Ser-Ser or Ser-Gly between RT codons 69 and 70). To study the dynamics of dominance of genomes with this genetic change, sequential HIV-1 isolates from two patients were analyzed with regard to consensus sequences and complexity of mutant spectra. The two patients displayed completely different, complex evolutionary patterns leading to temporary dominance of dipeptide insertions. In one patient, a virus very closely related to an ancestor virus from the same patient overtook the population at late times, displacing genomes encoding a Ser-Ser insertion. In another patient the sequential dominance of genomes with Ser-Ser insertion-->no insertion-->Ser-Gly insertion was observed. These three types of genomes coexisted in the mutant spectrum of one HIV-1 isolate. Complexity was also reflected in the shape of phylogenetic trees derived with genomes from the mutant spectrum at each time point. The results suggest that HIV-1 genomes encoding a dipeptide insertion between RT codons 69 and 70 do not show a clear selective advantage over other genomes lacking the insertion. Such an absence of a clear selective advantage will favor that such genomes encoding this RT insertion become dominant only in a transient fashion, and following disparate kinetics in different patients.

Criticality and scaling in evolutionary ecology

Fluctuations in ecological systems are known to involve a wide range of spatial and temporal scales, often displaying self-similar (fractal) properties. Recent theoretical approaches are trying to shed light on the nature of these complex dynamics. The results suggest that complexity in ecology and evolution comes from the network-like structure of multispecies communities that are close to instability. If true, these ideas might change our understanding of how complexity emerges in the biosphere and how macroevolutionary events could be decoupled from microevolutionary ones.

Exponential fitness gains of RNA virus populations are limited by bottleneck effects

Fitness is a parameter that quantitatively measures adaptation of a virus to a given environment. We have previously reported exponential fitness gains of large populations of vesicular stomatitis virus replicating in a constant environment (I. S. Novella et al., Proc. Natl. Acad. Sci. USA 92:5841-5844, 1995). In this paper, we report that during long-term passage of such large viral populations, fitness values reached a high-fitness plateau during which stochastic fitness variations were observed. This effect appears likely to be due to bottleneck effects on very high fitness populations.

Rabies virus quasispecies: implications for pathogenesis

Passage of the mouse-adapted rabies virus strain CVS-24 (where CVS is challenge virus standard) in BHK cells results in the rapid selection of a dominant variant designated CVS-B2c that differs genotypically and phenotypically from the dominant variant CVS-N2c present in mouse-brain- or neuroblastoma-cell-passaged CVS-24. The glycoprotein of CVS-B2c has 10 amino acid substitutions compared with that of CVS-N2c. Because CVS-B2c can be reproducibly selected in BHK cells, it is likely to be a conserved minor subpopulation of CVS-24. CVS-N2c is more neurotropic in vitro and in vivo than CVS-B2c, which replicates more readily in nonneuronal cells in vitro and in vivo. These characteristics appear to be relevant to the pathogenicity of the two variants. CVS-N2c is more pathogenic for adult mice than CVS-B2c. In contrast, CVS-B2c is more pathogenic for neonatal mice. These differences in pathogenicity are reflected in the selection pattern when mixtures of CVS-N2c and CVS-B2c were used to infect neonatal and adult mice. Although CVS-N2c was highly selected in adult mice, no selection for either variant was seen in neonates, suggesting that certain aspects of development, such as maturation of the nervous and immune systems, may contribute to the selection process. We speculate that the existence of different variants within a rabies virus strain may facilitate the virus in overcoming barriers to its spread, both within the host and between species.

RNA virus mutations and fitness for survival

RNA viruses exploit all known mechanisms of genetic variation to ensure their survival. Distinctive features of RNA virus replication include high mutation rates, high yields, and short replication times. As a consequence, RNA viruses replicate as complex and dynamic mutant swarms, called viral quasispecies. Mutation rates at defined genomic sites are affected by the nucleotide sequence context on the template molecule as well as by environmental factors. In vitro hypermutation reactions offer a means to explore the functional sequence space of nucleic acids and proteins. The evolution of a viral quasispecies is extremely dependent on the population size of the virus that is involved in the infections. Repeated bottleneck events lead to average fitness losses, with viruses that harbor unusual, deleterious mutations. In contrast, large population passages result in rapid fitness gains, much larger than those so far scored for cellular organisms. Fitness gains in one environment often lead to fitness losses in an alternative environment. An important challenge in RNA virus evolution research is the assignment of phenotypic traits to specific mutations. Different constellations of mutations may be associated with a similar biological behavior. In addition, recent evidence suggests the existence of critical thresholds for the expression of phenotypic traits. Epidemiological as well as functional and structural studies suggest that RNA viruses can tolerate restricted types and numbers of mutations during any specific time point during their evolution. Viruses occupy only a tiny portion of their potential sequence space. Such limited tolerance to mutations may open new avenues for combating viral infections.