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

In-Host Flat-like Quasispecies: Characterization Methods and Clinical Implications

The repeated failure to treat patients chronically infected with hepatitis E (HEV) and C (HCV) viruses, despite the absence of resistance-associated substitutions (RAS), particularly in response to prolonged treatments with the mutagenic agents of HEV, suggests that quasispecies structure may play a crucial role beyond single point mutations. Quasispecies structured in a flat-like manner (referred to as flat-like) are considered to possess high average fitness, occupy a significant fraction of the functional genetic space of the virus, and exhibit a high capacity to evade specific or mutagenic treatments. In this paper, we studied HEV and HCV samples using high-depth next-generation sequencing (NGS), with indices scoring the different properties describing flat-like quasispecies. The significance of these indices was demonstrated by comparing the values obtained from these samples with those from acute infections caused by respiratory viruses (betacoronaviruses, enterovirus, respiratory syncytial viruses, and metapneumovirus). Our results revealed that flat-like quasispecies in HEV and HCV chronic infections without RAS are characterized by numerous low-frequency haplotypes with no dominant one. Surprisingly, these low-frequency haplotypes (at the nucleotide level) exhibited a high level of synonymity, resulting in much lower diversity at the phenotypic level. Currently, clinical approaches for managing flat-like quasispecies are lacking. Here, we propose methods to identifying flat-like quasispecies, which represents an essential initial step towards exploring alternative treatment protocols for viruses resistant to conventional therapies.

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.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

Comparison of Endemic and Epidemic Vesicular Stomatitis Virus Lineages in Culicoides sonorensis Midges

Vesicular stomatitis virus (VSV) primarily infects livestock and is transmitted by direct contact and vectored by Culicoides midges (Diptera: Ceratopogonidae). Endemic to Central and South America, specific VSV lineages spread northward out of endemic regions of Mexico and into the U.S. sporadically every five to ten years. In 2012, a monophyletic epidemic lineage 1.1 successfully spread northward into the U.S. In contrast, the closest endemic ancestor, lineage 1.2, remained circulating exclusively in endemic regions in Mexico. It is not clear what roles virus-animal interactions and/or virus-vector interactions play in the ability of specific viral lineages to escape endemic regions in Mexico and successfully cause outbreaks in the U.S., nor the genetic basis for such incursions. Whole-genome sequencing of epidemic VSV 1.1 and endemic VSV 1.2 revealed significant differences in just seven amino acids. Previous studies in swine showed that VSV 1.1 was more virulent than VSV 1.2. Here, we compared the efficiency of these two viral lineages to infect the vector Culicoides sonorensis (Wirth and Jones) and disseminate to salivary glands for subsequent transmission. Our results showed that midges orally infected with the epidemic VSV 1.1 lineage had significantly higher infection dissemination rates compared to those infected with the endemic VSV 1.2 lineage. Thus, in addition to affecting virus-animal interactions, as seen with higher virulence in pigs, small genetic changes may also affect virus-vector interactions, contributing to the ability of specific viral lineages to escape endemic regions via vector-borne transmission.

A high-throughput drug screening strategy against coronaviruses

The emergence and re-emergence of coronaviruses (CoV) continually cause circulating epidemics and pandemics worldwide, such as the on-going outbreak of the novel coronavirus SARS-CoV-2. The resultant disease, coronavirus disease 2019 (COVID-19), has rapidly developed into a worldwide pandemic, leading to severe health and economic burdens. Although the recently announced vaccines against COVID-19 has rekindled hope, there is still a major challenge to urgently meet the global need for rapid treatment of the pandemic. Given the urgency of the CoV outbreak, we propose a strategy to screen potential broad-spectrum drugs against CoV in a high-throughput manner, particularly against SARS-CoV-2. Since the essential functional domains of CoV are extensively homologous, the availability of two types of mild CoV, HCoV-OC43 and MHV, should provide a valuable tool for the rapid identification of promising drugs against CoV without the drawbacks of level three biological confinements. The luciferase reporter gene is introduced into HCoV-OC43 and MHV to indicate viral activity, and hence the antiviral efficiency of screened drugs can be quantified by luciferase activity. Compounds with antiviral activity against both HCoV-OC43 and MHV are further evaluated in SARS-CoV-2 after structural optimizations. This system allows large-scale compounds to be screened to search for broad-spectrum drugs against CoV in a high-throughput manner, providing potential alternatives for clinical management of SARS-CoV-2 or other CoV.

Experimental Evolution Generates Novel Oncolytic Vesicular Stomatitis Viruses with Improved Replication in Virus-Resistant Pancreatic Cancer Cells

Vesicular stomatitis virus (VSV)-based oncolytic viruses are promising agents against pancreatic ductal adenocarcinoma (PDAC). However, some PDAC cell lines are resistant to VSV. Here, using a directed viral evolution approach, we generated novel oncolytic VSVs with an improved ability to replicate in virus-resistant PDAC cell lines, while remaining highly attenuated in nonmalignant cells. Two independently evolved VSVs obtained 2 identical VSV glycoprotein mutations, K174E and E238K. Additional experiments indicated that these acquired G mutations improved VSV replication, at least in part due to improved virus attachment to SUIT-2 cells. Importantly, no deletions or mutations were found in the virus-carried transgenes in any of the passaged viruses. Our findings demonstrate long-term genomic stability of complex VSV recombinants carrying large transgenes and support further clinical development of oncolytic VSV recombinants as safe therapeutics for cancer.

Vesicular stomatitis virus (VSV) based oncolytic viruses are promising agents against various cancers. We have shown that pancreatic ductal adenocarcinoma (PDAC) cell lines exhibit great diversity in susceptibility and permissibility to VSV. Here, using a directed evolution approach with our two previously described oncolytic VSV recombinants, VSV-p53wt and VSV-p53-CC, we generated novel oncolytic VSVs with an improved ability to replicate in virus-resistant PDAC cell lines. VSV-p53wt and VSV-p53-CC encode a VSV matrix protein (M) with a ΔM51 mutation (M-ΔM51) and one of two versions of a functional human tumor suppressor, p53, fused to a far-red fluorescent protein, eqFP650. Each virus was serially passaged 32 times (which accounts for more than 60 viral replication cycles) on either the SUIT-2 (moderately resistant to VSV) or MIA PaCa-2 (highly permissive to VSV) human PDAC cell lines. While no phenotypic changes were observed for MIA PaCa-2-passaged viruses, both SUIT-2-passaged VSV-p53wt and VSV-p53-CC showed improved replication in SUIT-2 and AsPC-1, another human PDAC cell line also moderately resistant to VSV, while remaining highly attenuated in nonmalignant cells. Surprisingly, two identical VSV glycoprotein (VSV-G) mutations, K174E and E238K, were identified in both SUIT-2-passaged viruses. Additional experiments indicated that the acquired G mutations improved VSV replication, at least in part due to improved virus attachment to SUIT-2 cells. Importantly, no mutations were found in the M-ΔM51 protein, and no deletions or mutations were found in the p53 or eqFP650 portions of virus-carried transgenes in any of the passaged viruses, demonstrating long-term genomic stability of complex VSV recombinants carrying large transgenes.

IMPORTANCE Vesicular stomatitis virus (VSV)-based oncolytic viruses are promising agents against pancreatic ductal adenocarcinoma (PDAC). However, some PDAC cell lines are resistant to VSV. Here, using a directed viral evolution approach, we generated novel oncolytic VSVs with an improved ability to replicate in virus-resistant PDAC cell lines, while remaining highly attenuated in nonmalignant cells. Two independently evolved VSVs obtained 2 identical VSV glycoprotein mutations, K174E and E238K. Additional experiments indicated that these acquired G mutations improved VSV replication, at least in part due to improved virus attachment to SUIT-2 cells. Importantly, no deletions or mutations were found in the virus-carried transgenes in any of the passaged viruses. Our findings demonstrate long-term genomic stability of complex VSV recombinants carrying large transgenes and support further clinical development of oncolytic VSV recombinants as safe therapeutics for cancer.

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.

The role of co-infection and swarm dynamics in arbovirus transmission

Arthropod-borne viruses (arboviruses) are transmitted by hematophagous insects, primarily mosquitoes. The geographic range and prevalence of mosquito-borne viruses and their vectors has dramatically increased over the last 50 years. As a result, the most medically important arboviurses now co-exist in many regions, resulting in an increased frequency of co-infections in hosts and vectors. In addition to concurrent infections with human pathogens, mosquito-only viruses and/or enzootic viruses not associated with human disease are ubiquitous in mosquito populations. Moreover, mosquito-borne viruses are largely RNA viruses that exist within individual hosts as a diverse and dynamic swarm of closely related genotypes. Interactions among co-infecting viruses and genotypes can have profound effects on virulence, fitness and evolution. Here, we review our understanding of how these complex interactions influence transmission of mosquito-borne viruses, focusing on the often-neglected virus interactions in the mosquito vector, and identify gaps in our knowledge that should guide future studies.

Novel Synthesis and Phenotypic Analysis of Mutant Clouds for Hepatitis E Virus Genotype 1

Many RNA viruses exist as an ensemble of genetically diverse, replicating populations known as a mutant cloud. The genetic diversity (cloud size) and composition of this mutant cloud may influence several important phenotypic features of the virus, including its replication capacity. We applied a straightforward, bacterium-free approach using error-prone PCR coupled with reverse genetics to generate infectious mutant RNA clouds with various levels of genetic diversity from a genotype 1 strain of hepatitis E virus (HEV). Cloning and sequencing of a genomic fragment encompassing 70% of open reading frame 1 (ORF1) or of the full genome from variants in the resultant clouds showed the occurrence of nucleotide mutations at a frequency on the order of 10^-3 per nucleotide copied and the existence of marked genetic diversity, with a high normalized Shannon entropy value. The mutant clouds showed transient replication in cell culture, while wild-type HEV did not. Cross-sectional data from these cell cultures supported the existence of differential effects of clouds of various sizes and compositions on phenotypic characteristics, such as the replication level of (+)-RNA progeny, the amounts of double-stranded RNA (a surrogate for the rate of viral replication) and ORF1 protein, and the expression of interferon-stimulated genes. Since mutant cloud size and composition influenced the viral phenotypic properties, a better understanding of this relationship may help to provide further insights into virus evolution and prediction of emerging viral diseases.IMPORTANCE Several biological or practical limitations currently prevent the study of phenotypic behavior of a mutant cloud in vitro We developed a simple and rapid method for synthesizing mutant clouds of hepatitis E virus (HEV), a single-stranded (+)-RNA [ss(+) RNA] virus, with various and controllable levels of genetic diversity, which could then be used in a cell culture system to study the effects of cloud size and composition on viral phenotype. In a cross-sectional analysis, we demonstrated that a particular mutant cloud which had an extremely high genetic diversity had a replication rate exceeding that of wild-type HEV. This method should thus provide a useful model for understanding the phenotypic behavior of ss(+) RNA viruses.

Sequence polymorphism in an insect RNA virus field population: A snapshot from a single point in space and time reveals stochastic differences among and within individual hosts

Population structure of Homalodisca coagulata Virus-1 (HoCV-1) among and within field-collected insects sampled from a single point in space and time was examined. Polymorphism in complete consensus sequences among single-insect isolates was dominated by synonymous substitutions. The mutant spectrum of the C2 helicase region within each single-insect isolate was unique and dominated by nonsynonymous singletons. Bootstrapping was used to correct the within-isolate nonsynonymous:synonymous arithmetic ratio (N:S) for RT-PCR error, yielding an N:S value ~one log-unit greater than that of consensus sequences. Probability of all possible single-base substitutions for the C2 region predicted N:S values within 95% confidence limits of the corrected within-isolate N:S when the only constraint imposed was viral polymerase error bias for transitions over transversions. These results indicate that bottlenecks coupled with strong negative/purifying selection drive consensus sequences toward neutral sequence space, and that most polymorphism within single-insect isolates is composed of newly-minted mutations sampled prior to selection.

Antigenic diversification is correlated with increased thermostability in a mammalian virus

The theory of plastogenetic congruence posits that ultimately, the pressure to maintain function in the face of biomolecular destabilization produces robustness. As temperature goes up so does destabilization. Thus, genetic robustness, defined as phenotypic constancy despite mutation, should correlate with survival during thermal challenge. We tested this hypothesis using vesicular stomatitis virus (VSV). We produced two sets of evolved strains after selection for higher thermostability by either preincubation at 37°C or by incubation at 40°C during infection. These VSV populations became more thermostable and also more fit in the absence of thermal selection, demonstrating an absence of tradeoffs. Eleven out of 12 evolved populations had a fixed, nonsynonymous substitution in the nucleocapsid (N) open reading frame. There was a partial correlation between thermostability and mutational robustness that was observed when the former was measured at 42°C, but not at 37°C. These results are consistent with our earlier work and suggest that the relationship between robustness and thermostability is complex. Surprisingly, many of the thermostable strains also showed increased resistance to monoclonal antibody and polyclonal sera, including sera from natural hosts. These data suggest that evolved thermostability may lead to antigenic diversification and an increased ability to escape immune surveillance in febrile hosts, and potentially to an improved robustness. These relationships have important implications not only in terms of viral pathogenesis, but also for the development of vaccine vectors and oncolytic agents.

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.

In Vivo Replication and Pathogenesis of Vesicular Stomatitis Virus Recombinant M40 Containing Ebola Virus L-Domain Sequences

The M40 VSV recombinant was engineered to contain overlapping PTAP and PPxY L-domain motifs and flanking residues from the VP40 protein of Ebola virus. Replication of M40 in cell culture is virtually indistinguishable from that of control viruses. However, the presence of the Ebola PTAP motif in the M40 recombinant enabled this virus to interact with and recruit host Tsg101, which was packaged into M40 virions. In this brief report, we compared replication and the pathogenic profiles of M40 and the parental virus M51R in mice to determine whether the presence of the Ebola L-domains and flanking residues altered in vivo characteristics of the virus. Overall, the in vivo characteristics of M40 were similar to those of the parental M51R virus, indicating that the Ebola sequences did not alter pathogenesis of VSV in this small animal model of infection.

Relationship between within-host fitness and virulence in the vesicular stomatitis virus: correlation with partial decoupling

Given the parasitic nature of viruses, it is sometimes assumed that rates of viral replication and dissemination within hosts (within-host fitness) correlate with virulence. However, there is currently little empirical evidence supporting this principle. To test this, we quantified the fitness and virulence of 21 single- or double-nucleotide mutants of the vesicular stomatitis virus in baby hamster kidney cells (BHK-21). We found that, overall, these two traits correlated positively, but significant outliers were identified. Particularly, a single mutation in the conserved C terminus of the N nucleocapsid (U1323A) had a strongly deleterious fitness effect but did not alter or even slightly increased virulence. We also found a double mutant of the M matrix protein and G glycoprotein (U2617G/A3802G mutant) with high fitness yet low virulence. We further characterized these mutants in primary cultures from mouse brain cells and in vivo and found that their relative fitness values were similar to those observed in BHK-21 cells. The mutations had weak effects on the virus-induced death rate of total brain cells, although they specifically reduced neuron death rates. Furthermore, increased apoptosis levels were detected in neurons infected with the U2617G/A3802G mutant, consistent with its known inability to block interferon secretion. In vivo, this mutant had reduced virulence and, despite its low brain titer, it retained a relatively high fitness value owing to its ability to suppress competitor viruses. Overall, our results are in broad agreement with the notion that viral fitness and virulence should be positively correlated but show that certain mutations can break this association and that the fitness-virulence relationship can depend on complex virus-host and virus-virus interactions.

A single amino acid change resulting in loss of fluorescence of eGFP in a viral fusion protein confers fitness and growth advantage to the recombinant vesicular stomatitis virus

Using a recombinant vesicular stomatitis virus encoding eGFP fused in-frame with an essential viral replication protein, the phosphoprotein P, we show that during passage in culture, the virus mutates the nucleotide C289 within eGFP of the fusion protein PeGFP to A or T, resulting in R97S/C amino acid substitution and loss of fluorescence. The resultant non-fluorescent virus exhibits increased fitness and growth advantage over its fluorescent counterpart. The growth advantage of the non-fluorescent virus appears to be due to increased transcription and replication activities of the PeGFP protein carrying the R97S/C substitution. Further, our results show that the R97S/C mutation occurs prior to accumulation of mutations that can result in loss of expression of the gene inserted at the G-L gene junction. These results suggest that fitness gain is more important for the recombinant virus than elimination of expression of the heterologous gene.

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.

Quasispecies as a matter of fact: viruses and beyond

We review the origins of the quasispecies concept and its relevance for RNA virus evolution, viral pathogenesis and antiviral treatment strategies. We emphasize a critical point of quasispecies that refers to genome collectivities as the unit of selection, and establish parallels between RNA viruses and some cellular systems such as bacteria and tumor cells. We refer also to tantalizing new observations that suggest quasispecies behavior in prions, perhaps as a result of the same quantum-mechanical indeterminations that underlie protein conformation and error-prone replication in genetic systems. If substantiated, these observations with prions could lead to new research on the structure-function relationship of non-nucleic acid biological molecules.

Dynamics of in vitro fitness recovery of HIV-1

The study on the evolutionary consequences of an RNA viral population's fluctuations can be approached by in vitro experiments. This work describes the fitness recovery of HIV-1 after 20 large-population passages in 10 debilitated clones. The serial passages promoted an increase in viral fitness. In addition, we detected a significant number of mutations fixed in the complete genome consensus sequence of the final viral populations. Among the mutations, events of convergent evolution with important phenotypic characteristics occurred in several independent clones. One common change, V35I, in the nuclear localization signal of the p17 protein appeared in four viruses of three different lineages. Other common alterations mapped in position E196K of the reverse transcriptase or in position S316K of the V3 loop of the gp120 residue that is associated with the X4/R5 phenotype. Together with this mutational analysis, we studied the quasispecies heterogeneity of the initial and final viruses, revealing that fitness increase correlated with an augmentation in the genetic heterogeneity of viral quasispecies. However, while heterogeneity was mostly composed of synonymous (dS) mutations in the first 10 passages performed, at passage 21 it switched to nonsynonymous (dN) substitutions, with significant differences in dN - dS values between passages 11 and 21. In summary, the HIV-1 in vitro fitness recovery depicts a multiphase process occurring first by generation of mutations followed by fixation of the beneficial ones, depicting a classical Darwinian process.

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.

Mechanisms of viral emergence

A number of virologic and environmental factors are involved in the emergence and re-emergence of viral disease. Viruses do not conservatively occupy a single and permanent ecological niche. Rather, due to their intrinsic capacity for genetic change, and to the evolvability of fitness levels, viruses display a potential to parasitize alternative host species. Mutation, recombination and genome segment reassortment, and combination of these molecular events, produce complex and phenotypically diverse populations of viruses, which constitute the raw material on which selection acts. The majority of emerging viral diseases of humans have a zoonotic origin. Sociologic and ecologic factors produce diverse and changing environments in which viral subpopulations have ample opportunities to be selected from intrinsically heterogeneous viral populations, particularly in the case of RNA viruses. In this manner, new human, animal and plant viruses have emerged periodically and, from all evidence, will continue to emerge. This article reviews some of the mechanisms that have been identified in viral emergence, with a focus on the importance of genetic variation of viruses, and on the general concept of biological complexity.

Combining mathematics and empirical data to predict emergence of RNA viruses that differ in reservoir use

RNA viruses may be particularly capable of contributing to the increasing biomedical problem of infectious disease emergence. Empirical studies and epidemiological models are informative for the understanding of evolutionary processes that promote pathogen emergence, but rarely are these approaches combined in the same study. Here, we used an epidemiology model containing observations of pathogen productivity in reservoirs, as a means to predict which pathogens should be most prone to emerge in a primary host such as humans. We employed as a model system a collection of vesicular stomatitis virus populations that had previously diverged in host-use strategy: specialists, directly selected generalists and indirectly selected (fortuitous) generalists. Using data from experiments where these viral strategists were challenged to grow on unencountered novel hosts in vitro, logistic growth models determined that the directly selected generalist viruses tended to grow best on model reservoirs. Furthermore, when we used the growth data to estimate average reproductive rate across secondary reservoirs, we showed that the combined approach could be used to estimate relative success of the differing virus strategists when encountering a primary host. Our study suggests that synergistic approaches combining epidemiological modelling with empirical data from experimental evolution may be useful for developing efforts to predict which types of pathogens pose the greatest probability of emerging in the future.

Inhibition of encephalomyocarditis virus and poliovirus replication by quinacrine: implications for the design and discovery of novel antiviral drugs

The 9-aminoacridine (9AA) derivative quinacrine (QC) has a long history of safe human use as an antiprotozoal and antirheumatic agent. QC intercalates into DNA and RNA and can inhibit DNA replication, RNA transcription, and protein synthesis. The extent of QC intercalation into RNA depends on the complexity of its secondary and tertiary structure. Internal ribosome entry sites (IRESs) that are required for initiation of translation of some viral and cellular mRNAs typically have complex structures. Recent work has shown that some intercalating drugs, including QC, are capable of inhibiting hepatitis C virus IRES-mediated translation in a cell-free system. Here, we show that QC suppresses translation directed by the encephalomyocarditis virus (EMCV) and poliovirus IRESs in a cell-free system and in virus-infected HeLa cells. In contrast, IRESs present in the mammalian p53 transcript that are predicted to have less-complex structures were not sensitive to QC. Inhibition of IRES-mediated translation by QC correlated with the affinity of binding between QC and the particular IRES. Expression of viral capsid proteins, replication of viral RNAs, and production of virus were all strongly inhibited by QC (and 9AA). These results suggest that QC and similar intercalating drugs could potentially be used for treatment of viral infections.

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.

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.

Enhanced adaptation of vesicular stomatitis virus in cells infected with vaccinia virus

Infections involving different viruses (multiple infections) are common in nature and can take place between different strains of the same virus or between different virus species, including DNA and RNA viruses. The influence of multiple infections on viral evolution has been previously studied using different populations of the same virus. Here, we took a step forward by studying the evolution of an RNA virus (vesicular stomatitis virus, VSV) in the presence of a resident DNA virus (vaccinia virus, VV). Cell cultures were infected with a constant amount of VV, and VSV was added at four different post-VV-inoculation times and four different population sizes. The results showed that the presence of VV accelerates the adaptation of VSV to a cellular environment, especially at high population sizes. The effect of VV on VSV evolution was stronger when cells were incubated for longer times with VV prior to the addition of VSV. Our results suggest that cooperation between the two viruses rather than competition might be responsible for the enhanced rate of adaptation of VSV. Further studies are needed to discern whether infections involving different viruses could have an increased ability to escape antiviral strategies.

Antagonistic pleiotropy involving promoter sequences in a virus

Selection of specialist genotypes, that is, populations with limited niche width, promotes the maintenance of diversity. Specialization to a particular environment may have a cost in other environments, including fitness tradeoffs. When the tradeoffs are the result of mutations that have a beneficial effect in the selective environment but a deleterious effect in other environments, we have antagonistic pleiotropy. Alternatively, tradeoffs can result from the fixation of mutations that are neutral in the selective environment but have a negative effect in other environments, and thus the tradeoff is due to mutation accumulation. We tested the mechanisms underlying the fitness tradeoffs observed during adaptation to persistent infection of vesicular stomatitis virus in insect cells by sequencing the full-length genomes of 12 strains with a history of replication in a single niche (acute mammalian infection or persistent insect infection) or in temporally heterogeneous niches and correlated genetic and fitness changes. Ecological theory predicts a correlation between the selective environment and the niche width of the evolved populations, such that adaptation to single niches should lead to the selection of specialists and niche cycling should result in the selection of generalists. Contrary to this expectation, adaptation to one of the single niches resulted in a generalist and adaptation to a heterogeneous environment led to the selection of a specialist. Only one-third of the mutations that accumulated during persistent infection had a fitness cost that could be explained in all cases by antagonistic pleiotropy. Mutations involved in fitness tradeoffs included changes in regulatory sequences, particularly at the 3' termini of the genomes, which contain the single promoter that controls viral transcription and replication.

Replicating viral vectors as HIV vaccines Summary Report from IAVI Sponsored Satellite Symposium, International AIDS Society Conference, July 22, 2007

At the International AIDS Society Conference on Pathogenesis, Treatment and Prevention held in Sydney, Australia, in July 2007, the International AIDS Vaccine Initiative (IAVI) convened a satellite symposium entitled 'Accelerating the Development of Replicating Viral Vectors for AIDS Vaccines.' Its purpose was to highlight the rationale for accelerating the development of replicating viral vectors for use as vaccines against HIV-1, and to bring together vaccine scientists, regulatory officials, and public health specialists from industrialized and developing nations to discuss the major issues facing the development and testing of replicating viral vector-based vaccines.

Quasispecies and its impact on viral hepatitis

Quasispecies dynamics mediates adaptability of RNA viruses through a number of mechanisms reviewed in the present article, with emphasis on the medical implications for the hepatitis viruses. We discuss replicative and non-replicative molecular mechanisms of genome variation, modulating effects of mutant spectra, and several modes of viral evolution that can affect viral pathogenesis. Relevant evolutionary events include the generation of minority virus variants with altered functional properties, and alterations of mutant spectrum complexity that can affect disease progression or response to treatment. The widespread occurrence of resistance to antiviral drugs encourages new strategies to control hepatic viral disease such as combination therapies and lethal mutagenesis. In particular, ribavirin may be exerting in some cases its antiviral activity with participation of its mutagenic action. Despite many unanswered questions, here we document that quasispecies dynamics has provided an interpretation of the adaptability of the hepatitis viruses, with features conceptually similar to those observed with other RNA viruses, a reflection of the common underlying Darwinian principles.

A real-time RT-PCR assay for quantifying the fitness of tobacco etch virus in competition experiments

Relative fitness determination has become a standard tool in experimental virus evolution studies. In this type of studies, the tested strain is mixed with a reference strain, which differs in an easy-to-score and genetically stable marker, and allowed to compete for a limited common pool of resources during a given number of generations. In this report, a TaqMan real-time PCR methodology is proposed for quantifying the relative fitness of tobacco etch potyvirus strains (TEV) in in planta mixed infections with a reference TEV strain. Two different forward primers along with a common reverse one are used into separated reactions mixes from the same RNA preparation. The reference strain, named TEV-PC1, was genetically engineered to carry a neutral marker in a highly conserved region of the RNA polymerase NIb gene. This marker allows tracking the frequency of both competitors during competition experiments by real-time quantitative PCR using specific primers. Both the reproducibility and sensitivity of the method have been explored. Reproducibility was assessed by running multiple competition experiments for the same genotype. Sensitivity was assessed by comparing the results of competition experiments against TEV-PC1 of 24 single-nucleotide substitutions mutants.

Quantifying antagonistic epistasis in a multifunctional RNA secondary structure of the Rous sarcoma virus

Recent studies have suggested that antagonistic epistasis (i.e. mutations having smaller effects in combination than alone) may be common among RNA viruses, in contrast to other biological systems. Here, by re-analysing previously published data from a random viral library, selection and epistasis coefficients were estimated in the U5-IR stem and loop of theRous sarcoma virus, a region that adopts a conserved secondary structure and is involved in various essential steps of viral infection. The estimated mutational fitness effects are extremely high and genetic interactions are antagonistic on average. This pattern might be representative of RNA virus genomes, which show high compaction and frequent secondary structures. The implications for RNA virus adaptability are explored.

Viruses as quasispecies: biological implications

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

Symbiosis versus competition in plant virus evolution

Darwin's theory of evolution by natural selection has been supported by molecular evidence and by experimental evolution of viruses. However, it might not account for the evolution of all life, and an alternative model of evolution through symbiotic relationships also has gained support. In this review, the evolution of plant viruses has been reinterpreted in light of these two seemingly opposing theories by using evidence from the earliest days of plant virology to the present. Both models of evolution probably apply in different circumstances, but evolution by symbiotic association (symbiogenesis) is the most likely model for many evolutionary events that have resulted in rapid changes or the formation of new species. In viruses, symbiogenesis results in genomic reassortment or recombination events among disparate species. These are most noticeable by phylogenetic comparisons of extant viruses from different taxonomic groups.

Selection of and recombination between minor variants lead to the adaptation of an avian coronavirus to primate cells

An interesting question posed by the current evidence that severe acute respiratory syndrome coronavirus may be originated from an animal coronavirus is how such an animal coronavirus breaks the host species barrier and becomes zoonotic. In this report, we study the chronological order of genotypic changes in the spike protein of avian coronavirus infectious bronchitis virus (IBV) during its adaptation to a primate cell line. Adaptation of the Beaudette strain of IBV from chicken embryo to Vero cells showed the accumulation of 49 amino acid mutations. Among them, 26 (53.06%) substitutions were located in the S protein. Sequencing analysis and comparison of the S gene demonstrated that the majority of the mutations were accumulated and fixed at passage 7 on Vero cells and minor variants were isolated in several passages. Evidence present suggests that the dominant Vero cell-adapted IBV strain may be derived from the chicken embryo passages by selection of and potential recombination between the minor variants. This may explain why adaptation is a rapid process and the dominant strain, once adapted to a new host cell, becomes relatively stable.

Population structure within lineages of Wheat streak mosaic virus derived from a common founding event exhibits stochastic variation inconsistent with the deterministic quasi-species model

Structure of Wheat streak mosaic virus (WSMV) populations derived from a common founding event and subjected to serial passage at high multiplicity of infection (MOI) was evaluated. The founding population was generated by limiting dilution inoculation. Lineages of known pedigree were sampled at passage 9 (two populations) and at passage 15, with (three populations) or without mixing (four populations) of lineages at passage 10. Polymorphism within each population was assessed by sequencing 17-21 clones containing a 1371 nt region (WSMV-Sidney 81 nts 8001-9371) encompassing the entire coat protein cistron and flanking regions. Mutation frequency averaged approximately 5.0 x 10(-4)/nt across all populations and ranged from 2.4 to 11.6 x 10(-4)/nt within populations, but did not consistently increase or decrease with the number of passages removed from the founding population. Shared substitutions (19 nonsynonymous, 10 synonymous, and 3 noncoding) occurred at 32 sites among 44 haplotypes. Only four substitutions became fixed (frequency = 100%) within a population and nearly one third (10/32) never achieved a frequency of 10% or greater in any sampled population. Shared substitutions were randomly distributed with respect to genome position, with transitions outnumbering transversions 5.4:1 and a clear bias for A to G and U to C substitutions. Haplotype composition of each population was unique with complexity of each population varying unpredictably, in that the number and frequency of haplotypes within a lineage were not correlated with number of passages removed from the founding population or whether the population was derived from a single or mixed lineage. The simplest explanation is that plant virus lineages, even those propagated at high MOI, are subject to frequent, narrow genetic bottlenecks during systemic movement that result in low effective population size and stochastic changes in population structure upon serial passage.

Epistasis and the adaptability of an RNA virus

We have explored the patterns of fitness recovery in the vesicular stomatitis RNA virus. We show that, in our experimental setting, reversions to the wild-type genotype were rare and fitness recovery was at least partially driven by compensatory mutations. We compared compensatory adaptation for genotypes carrying (1) mutations with varying deleterious fitness effects, (2) one or two deleterious mutations, and (3) pairs of mutations showing differences in the strength and sign of epistasis. In all cases, we found that the rate of fitness recovery and the proportion of reversions were positively affected by population size. Additionally, we observed that mutations with large fitness effect were always compensated faster than mutations with small fitness effect. Similarly, compensatory evolution was faster for genotypes carrying a single deleterious mutation than for those carrying pairs of mutations. Finally, for genotypes carrying two deleterious mutations, we found evidence of a negative correlation between the epistastic effect and the rate of compensatory evolution.

Adaptability costs in immune escape variants of vesicular stomatitis virus

We have used vesicular stomatitis virus (VSV) to determine the cost of antiserum resistance during escape from a polyclonal immune response. Replication of VSV in the presence of polyclonal antiserum resulted in the selection of antibody-escape mutants, as shown by increased fitness in the presence of antiserum and by increased resistance to neutralization. However, resistance came at a cost of overall fitness loss in the BHK-21 host cells. Sequencing of the surface G glycoprotein showed that two to four mutations were fixed in each population, most of which mapped in the A1 and A2 antigenic sites. Selected resistant populations were passaged as large populations in BHK-21 cells under constant conditions, which would normally lead to fitness increases. Nevertheless, many of the populations showed little or no sign of recovery, although the resistant phenotype was maintained. These results suggest that while antiserum resistance can develop, it may come at a cost in fitness and further limitations in the adaptability of the populations.

Vesicular stomatitis virus evolution during alternation between persistent infection in insect cells and acute infection in mammalian cells is dominated by the persistence phase

Vesicular stomatitis virus has the potential for very rapid evolution in the laboratory, but like many other arboviruses, it evolves at a relatively slow rate in the natural environment. Previous work showed that alternating replication in different cell types does not promote stasis. In order to determine whether other factors promote stasis, we compared the fitness trajectories of populations evolving during acute infections in mammalian cells, populations evolving during persistent infections in insect cells, and populations evolving during alternating acute and persistent infection cycles. Populations evolving under constant conditions increased in fitness in the environment in which they replicated. An asymmetric trade-off was observed such that acute infection had no cost for persistence but persistent replication had a dramatic cost for acute infection in mammalian cells. After an initial period of increase, fitness remained approximately constant in all the populations that included persistent replication, but fitness continuously increased in populations evolving during acute infections. Determination of the consensus sequence of the genes encoding the N, P, M, and G proteins showed that the pattern of mutation accumulation was coherent with fitness changes during persistence so that once fitness reached a maximum, the rate of mutation accumulation dropped. Persistent replication dominated both the genetic and the phenotypic evolution of the populations that alternated between acute infection of mammalian cells and persistence in insect cells, and fitness loss was observed in the mammalian environment despite periodic replication in mammalian cells. These results show that stasis can be achieved without good levels of adaptation to both the mammalian and the insect environments.

Human respiratory coronavirus OC43: genetic stability and neuroinvasion

The complete genome sequences of the human coronavirus OC43 (HCoV-OC43) laboratory strain from the American Type Culture Collection (ATCC), and a HCoV-OC43 clinical isolate, designated Paris, were obtained. Both genomes are 30,713 nucleotides long, excluding the poly(A) tail, and only differ by 6 nucleotides. These six mutations are scattered throughout the genome and give rise to only two amino acid substitutions: one in the spike protein gene (I958F) and the other in the nucleocapsid protein gene (V81A). Furthermore, the two variants were shown to reach the central nervous system (CNS) after intranasal inoculation in BALB/c mice, demonstrating neuroinvasive properties. Even though the ATCC strain could penetrate the CNS more effectively than the Paris 2001 isolate, these results suggest that intrinsic neuroinvasive properties already existed for the HCoV-OC43 ATCC human respiratory isolate from the 1960s before it was propagated in newborn mouse brains. It also demonstrates that the molecular structure of HCoV-OC43 is very stable in the environment (the two variants were isolated ca. 40 years apart) despite virus shedding and chances of persistence in the host. The genomes of the two HCoV-OC43 variants display 71, 53.1, and 51.2% identity with those of mouse hepatitis virus A59, severe acute respiratory syndrome human coronavirus Tor2 strain (SARS-HCoV Tor2), and human coronavirus 229E (HCoV-229E), respectively. HCoV-OC43 also possesses well-conserved motifs with regard to the genome sequence of the SARS-HCoV Tor2, especially in open reading frame 1b. These results suggest that HCoV-OC43 and SARS-HCoV may share several important functional properties and that HCoV-OC43 may be used as a model to study the biology of SARS-HCoV without the need for level three biological facilities.

Fitness analyses of vesicular stomatitis strains with rearranged genomes reveal replicative disadvantages

Gene expression of the nonsegmented negative-strand RNA viruses is determined by the position of each gene relative to that the single 3' promoter. The general order of genes among all of the viruses of the order Mononegavirales is highly conserved. In previous work we generated recombinant viruses in which the order of the three central genes of the prototypical rhabdovirus, vesicular stomatitis virus, was rearranged to all six possible permutations. While some of these viruses replicated less well than the wild type when assayed by single-step growth analyses in BSC-1 cells, others replicated as well or slightly better. In the work reported here, we used competition assays to compare the fitness of the viruses with alternative gene orders to that of the wild-type (wt) virus. We found that the relative fitness of these recombinant viruses depended on the multiplicity of infection (MOI) but not on the population size. However, during competitions at low MOI, when complementation cannot compensate for the defects of the populations with rearranged genomes, the virus with the wt gene order was always the most fit.

The distribution of fitness effects caused by single-nucleotide substitutions in an RNA virus

Little is known about the mutational fitness effects associated with single-nucleotide substitutions on RNA viral genomes. Here, we used site-directed mutagenesis to create 91 single mutant clones of vesicular stomatitis virus derived from a common ancestral cDNA and performed competition experiments to measure the relative fitness of each mutant. The distribution of nonlethal deleterious effects was highly skewed and had a long, flat tail. As expected, fitness effects depended on whether mutations were chosen at random or reproduced previously described ones. The effect of random deleterious mutations was well described by a log-normal distribution, with -19% reduction of average fitness; the effects distribution of preobserved deleterious mutations was better explained by a beta model. The fit of both models was improved when combined with a uniform distribution. Up to 40% of random mutations were lethal. The proportion of beneficial mutations was unexpectedly high. Beneficial effects followed a gamma distribution, with expected fitness increases of 1% for random mutations and 5% for preobserved mutations.

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.

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.