High multiplicities of infection favor rapid and random evolution of vesicular stomatitis virus

Strategies for Assessing Arbovirus Genetic Variability in Vectors and/or Mammals

Animal arboviruses replicate in their invertebrate vectors and vertebrate hosts. They use several strategies to ensure replication/transmission. Their high mutation rates and propensity to generate recombinants and/or genome segment reassortments help them adapt to new hosts/emerge in new geographical areas. Studying arbovirus genetic variability has been used to identify indicators which predict their potential to adapt to new hosts and/or emergence and in particular quasi-species. Multiple studies conducted with insect-borne viruses laid the foundations for the "trade-off" hypothesis (alternation of host transmission cycle constrains arbovirus evolution). It was extrapolated to tick-borne viruses, where too few studies have been conducted, even though humans faced emergence of numerous tick-borne virus during the last decades. There is a paucity of information regarding genetic variability of these viruses. In addition, insects and ticks do not have similar lifecycles/lifestyles. Indeed, tick-borne viruses are longer associated with their vectors due to tick lifespan. The objectives of this review are: (i) to describe the state of the art for all strategies developed to study genetic variability of insect-borne viruses both in vitro and in vivo and potential applications to tick-borne viruses; and (ii) to highlight the specificities of arboviruses and vectors as a complex and diverse system.

Mutagenic Effects of Ribavirin on Hepatitis E Virus-Viral Extinction versus Selection of Fitness-Enhancing Mutations

Hepatitis E virus (HEV), an important agent of viral hepatitis worldwide, can cause severe courses of infection in pregnant women and immunosuppressed patients. To date, HEV infections can only be treated with ribavirin (RBV). Major drawbacks of this therapy are that RBV is not approved for administration to pregnant women and that the virus can acquire mutations, which render the intra-host population less sensitive or even resistant to RBV. One of the proposed modes of action of RBV is a direct mutagenic effect on viral genomes, inducing mismatches and subsequent nucleotide substitutions. These transition events can drive the already error-prone viral replication beyond an error threshold, causing viral population extinction. In contrast, the expanded heterogeneous viral population can facilitate selection of mutant viruses with enhanced replication fitness. Emergence of these mutant viruses can lead to therapeutic failure. Consequently, the onset of RBV treatment in chronically HEV-infected individuals can result in two divergent outcomes: viral extinction versus selection of fitness-enhanced viruses. Following an overview of RNA viruses treated with RBV in clinics and a summary of the different antiviral modes of action of this drug, we focus on the mutagenic effect of RBV on HEV intrahost populations, and how HEV is able to overcome lethal mutagenesis.

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.

Molecular intermediates of fitness gain of an RNA virus: characterization of a mutant spectrum by biological and molecular cloning

The mutant spectrum of a virus quasispecies in the process of fitness gain of a debilitated foot-and-mouth disease virus (FMDV) clone has been analysed. The mutant spectrum was characterized by nucleotide sequencing of three virus genomic regions (internal ribosome entry site; region between the two AUG initiation codons; VP1-coding region) from 70 biological clones (virus from individual plaques formed on BHK-21 cell monolayers) and 70 molecular clones (RT--PCR products cloned in E. coli). The biological and molecular clones provided statistically indistinguishable definitions of the mutant spectrum with regard to the distribution of mutations among the three genomic regions analysed and with regard to the types of mutations, mutational hot-spots and mutation frequencies. Therefore, the molecular cloning procedure employed provides a simple protocol for the characterization of mutant spectra of viruses that do not grow in cell culture. The number of mutations found repeated among the clones analysed was higher than expected from the mean mutation frequencies. Some components of the mutant spectrum reflected genomes that were dominant in the prior evolutionary history of the virus (previous passages), confirming the presence of memory genomes in virus quasispecies. Other components of the mutant spectrum were genomes that became dominant at a later stage of evolution, suggesting a predictive value of mutant spectrum analysis with regard to the outcome of virus evolution. The results underline the observation that greater insight into evolutionary processes of viruses may be gained from detailed clonal analyses of the mutant swarms at the sequence level.

Dynamics of rabies virus quasispecies during serial passages in heterologous hosts

To understand the mutations and genetic rearrangements that allow rabies virus infections of new hosts and adaptation in nature, the quasispecies structure of the nucleoprotein and glycoprotein genes as well as two noncoding sequences of a rabies virus genome were determined. Gene sequences were obtained from the brain and from the salivary glands of the original host, a naturally infected European fox, and after serial passages in mice, dogs, cats and cell culture. A relative genetic stasis of the consensus sequences confirmed previous results about the stability of rabies virus. At the quasispecies level, the mutation frequency varies, in the following order: glycoprotein region (21.9 x 10(-4) mutations per bp), noncoding sequence nucleoprotein-phosphoprotein region (7.2-7.9 x 10(-4) mutations per bp) and nucleoprotein gene region (2.9-3.7 x 10(-4) mutations per bp). These frequencies varied according to the number, type of heterologous passages and the genomic region considered. The shape of the quasispecies structure was dramatically modified by passages in mice, in which the mutation frequencies increased by 12-31 x 10(-4) mutations per bp, depending on the region considered. Non-synonymous mutations were preponderant particularly in the glycoprotein gene, stressing the importance of positive selection in the maintenance and fixation of substitutions. Two mechanisms of genomic evolution of the rabies virus quasispecies, while adapting to environmental changes, have been identified: a limited accumulation of mutations with no replacement of the original master sequence and a less frequent but rapid selective overgrowth of favoured variants.

Genetic and fitness changes accompanying adaptation of an arbovirus to vertebrate and invertebrate cells

The alternating host cycle and persistent vector infection may constrain the evolution of arboviruses. To test this hypothesis, eastern equine encephalitis virus was passaged in BHK or mosquito cells, as well as in alternating (both) host cell passages. High and low multiplicities were used to examine the effect of defective interfering particles. Clonal BHK and persistent mosquito cell infections were also evaluated. Fitness was measured with one-step growth curves and competition assays, and mutations were evaluated by nucleotide sequencing and RNA fingerprinting. All passages and assays were done at 32 degrees C to eliminate temperature as a selection factor. Viruses passaged in either cell type alone exhibited fitness declines in the bypassed cells, while high-multiplicity and clonal passages caused fitness declines in both types of cells. Bypassed cell fitness losses were mosquito and vertebrate specific and were not restricted to individual cell lines. Fitness increases occurred in the cell line used for single-host-adaptation passages and in both cells for alternately passaged viruses. Surprisingly, single-host-cell passage increased fitness in that cell type no more than alternating passages. However, single-host-cell adaptation resulted in more mutations than alternating cell passages. Mosquito cell adaptation invariably resulted in replacement of the stop codon in nsP3 with arginine or cysteine. In one case, BHK cell adaptation resulted in a 238-nucleotide deletion in the 3' untranslated region. Many nonsynonymous substitutions were shared among more than one BHK or mosquito cell passage series, suggesting positive Darwinian selection. Our results suggest that alternating host transmission cycles constrain the evolutionary rates of arboviruses but not their fitness for either host alone.

Serial passage of human immunodeficiency virus type 1 generates misalignment deletions in non-essential accessory genes

Human immunodeficiency virus type 1 (HIV-1) derived from an infectious molecular clone pNL432 was extensively passaged in tissue culture by repeated rounds of acute infection. We previously showed the natural occurrence of a nonsense mutation in the vpr gene during continued passage of this virus. In this report, we show that two forms of large deletions (561 and 518 base pairs containing short direct repeats at the deletion junctions) occur after passage 50 in the region that spans the vif and vpr open reading frames. One model to explain the occurrence of these deletion regions is that such mutations result from misalignment of the growing point at a limited number of nucleotide positions. Infection of CD4+ T-cells with a recombinant HIV-1 construct containing the same vif to vpr deletion showed virtually no cytopathogenic phenotype. Thus, misalignment deletions at non-essential accessory genes of HIV-1 might be induced during replication, which result in the generation of virus with a low cytopathogenic potential.

Differences in mutagenesis during minus strand, plus strand and strand transfer (recombination) synthesis of the HIV-1 gene in vitro

We have developed an HIV nef-Escherichia coli lacZ fusion system in vitro that allows the detection of low frequency mutations, including frameshifts, deletions and insertions. A portion of the nef gene that encompasses a hypervariable region was fused in-frame with a downstream lacZalpha peptide coding region. The resulting lacZalpha peptide fusion protein remained functional. Any frameshift mutations in the nef insert would put the downstream lacZ alpha peptide gene out of frame, eliminating alpha complementation. With this system we compared the error rates of frameshift mutations that arise during DNA-directed and RNA-directed DNA synthesis. Results showed that DNA-directed and RNA-directed DNA synthesis did not contribute equally to the generation of mutations. DNA-directed DNA synthesis generated frameshift mutations at a frequency approximately 10-fold higher than those arising from RNA-directed DNA synthesis. RNA-directed DNA synthesis in the presence of acceptor templates showed an increase in mutation rate and differences in the mutation spectrum. The enhancement of mutation rate was caused by the appearance of mutations at three new locations that correlated with likely recombination sites. Results indicate that recombination is another source of mutations during viral replication.

Generation and role of defective proviruses in cytopathic feline leukemia virus (FeLV-FAIDS) infections

Cytopathic feline leukemia virus (FeLV) infections of feline T-cell line (FeT-cell) cultures led to the accumulation and maintenance of threefold more proviruses with deletions within the polymerase gene (pol) than minimally cytopathic FeLV infections. Over 60% of the viral DNA from cytopathic infections bore deletions in pol. Characterization of DNA sequences adjoining the deletions revealed that the junctions were most often flanked by RNA splice donor and acceptor consensus motifs. A thymidine-to-cytidine mutation introduced at the +2 position of one RNA splice donor-like motif inhibited formation of the two most prevalent viral DNA species with deletions, confirming the origin of many proviruses with deletions from reverse transcription of aberrantly spliced viral RNA species. An example of deletion by misalignment was also characterized. Viral inocula obtained from cells recovered after cytopathic infections were attenuated in their ability to cause cytopathic effects (CPE) and were able to confer superinfection resistance to naïve FeT-cells, despite maintaining envelope gene (env) sequences with full cytopathic potential. This suggested that viral genomes with deletions, rather than being required for cytopathicity, play a role in protecting cells from CPE. Indeed, expression of a molecularly cloned provirus bearing one of the characterized deletions attenuated CPE in FeT-cells caused by superinfecting cytopathic virus.

Virus-induced immunosuppression: kinetic analysis of the selection of a mutation associated with viral persistence

Infection of neonatal mice with lymphocytic choriomeningitis virus (LCMV) strain Armstrong (ARM) results in a lifelong persistent infection. Viral variants (cytotoxic T lymphocyte [CTL] negative, persistence positive [CTL- P+]) can be isolated from the lymphoid tissues of such mice. Adult mice inoculated with these CTL- P+ viruses fail to generate sufficient cytotoxic T lymphocytes to clear the acute infection and become persistently infected. By contrast, inoculation of a similar dose of the parental ARM virus (CTL+ P-) into adult mice leads to the generation of a vigorous virus-specific CTL response that clears the infection. Sequence analysis revealed a phenylalanine (Phe)-to-Leucine (Leu) change at amino acid 260 of the viral glycoprotein (GP) as a marker for variant viruses with the CTL- P+ phenotype. An RNA PCR assay that detects the variant GP sequence and thus allows kinetic studies of the selection of the Leu at position 260 was developed. We found that although CTL- P+ viruses are known to be lymphotropic, mature T and B cells were not required for the generation and selection of the Leu at GP amino acid 260. Kinetically, in mice infected at birth with LCMV ARM, as early as 3 weeks postinfection the Phe-to-Leu change was found in virus in the serum. By 5 weeks, viral nucleic acid obtained from peritoneal macrophages, spleen, lymph nodes, and liver showed the Phe-to-Leu change. At 2 months postinfection, the Leu change was detected in virus from the thymus, heart, lung, and kidney. By contrast, virus replicating in the central nervous system showed only minimal levels of the Leu change by 4 months and as long as 1 year postinfection. In vitro studies showed that the parental LCMV ARM CTL+ P- virus replicates more efficiently and outcompetes CTL- P+ virus in a cultured neuronal cell line, indicating that differential growth properties in neurons are likely the basis for the selection of the parental virus over the CTL- P+ variant in the brain.

New viruses and new disease: mutation, evolution and ecology

Given the extraordinarily high mutation rate of viruses, particularly those with RNA genomes, it is not surprising that new viruses are continually evolving. However, the symptomatology of old viral diseases has remained stable for centuries. The combination of genetic and ecological factors that constrain as well as facilitate the emergence of new viruses is analyzed.

Quantitation of relative fitness and great adaptability of clonal populations of RNA viruses

We describe a sensitive, internally controlled method for comparing the genetic adaptability and relative fitness of virus populations in constant or changing host environments. Certain monoclonal antibody-resistant mutants of vesicular stomatitis virus can compete equally during serial passages in mixtures with the parental wild-type clone from which they were derived. These genetically marked "surrogate wild-type" neutral mutants, when mixed with wild-type virus, allow reliable measurement of changes in virus fitness and of virus adaptation to different host environments. Quantitative fitness vector plots demonstrate graphically that even clones of an RNA virus are composed of complex variant populations (quasispecies). Variants of greater fitness (competitive replication ability) were selected within very few passages of virus clones in new host cells or animals. Even clones which were well adapted to BHK21 cells gained further fitness during repeated passages in BHK21 cells.

RNA virus quasispecies populations can suppress vastly superior mutant progeny

A variant clone of vesicular stomatitis virus recovered from a high-passage, evolving virus population replicated rapidly and produced remarkably high yields of virus, but these variants never dominated during further passages. We show that this clone is highly competitive, but it can overwhelm its progenitor population only when seeded above threshold level during dilute passages.

Polymerase errors accumulating during natural evolution of the glycoprotein gene of vesicular stomatitis virus Indiana serotype isolates

We report the entire glycoprotein (G) gene nucleotide sequences of 26 vesicular stomatitis virus Indiana serotype (VSV IND) type 1 isolates from North and Central America. These sequences are also compared with partial G gene sequences of VSV IND type 2 (Cocal) and type 3 (Alagoas) viruses and the complete G gene sequences of the more distantly related VSV New Jersey (NJ) and Chandipura viruses. Phylogenetic analysis of the G gene sequences by maximum parsimony revealed four major lineages or subtypes within the classical VSV IND (type 1) viruses, each with a distinct geographic distribution. A high degree of VSV genetic diversity was found in Central America, with several virus subtypes of both VSV IND and NJ serotypes existing in this mainly enzootic disease region. Nineteen percent sequence variation but no deletions or insertions were evident within the 5' noncoding and the coding regions of the VSV IND type 1 G genes. In addition to numerous base substitutions, the 3' noncoding regions of these viruses also contained numerous base insertions and deletions. This resulted in striking variation in G gene sizes, with gene lengths ranging from 1,652 to 1,868 nucleotides. As the VSV IND type 1 subtypes have diverged from the common ancestor with the NJ subtypes, their G mRNAs have accumulated more 3' noncoding sequence inserts, ranging up to 303 nucleotides in length. These primarily consist of an imprecise reiteration of the sequence UUUUUAA, apparently generated by a unique polymerase stuttering error. Analysis of the deduced amino acid sequence differences among VSV IND type 1 viruses revealed numerous substitutions within defined antigenic epitopes, suggesting that immune selection may play a role in the evolution of these viruses.

Stability of the nucleotide sequence of the phosphoprotein gene of measles virus during lytic infections

Three clones with cDNA inserts encoding large portions of the measles virus phosphoprotein mRNA were characterized and compared with a previously published sequence of the Edmonston strain of measles virus. The two cloned viruses were separated by more than 100 passages. Only one out of 1477 nucleotides differed in the two sequences reflecting a very low mutation rate of the phosphoprotein gene during dilute lytic passages. The discovery that a third reading frame in the phosphoprotein gene may code for a novel peptide chain in addition to the P and C peptides may explain some of the high stability of the gene. The new reading frame was accessed by a translational shift caused by insertion of one extra G at a particular site in one of three otherwise identical cDNA sequences. A discrepancy was also found between the presumably high error rate of viral RNA polymerases and the stability of nucleotides in which mutations would not lead to amino acid substitutions. A few errors in the previously published sequence were discovered and the corrections are presented.

Very high frequency of reversion to guanidine resistance in clonal pools of guanidine-dependent type 1 poliovirus

We have carefully examined the frequency of guanidine-resistant revertants in six different clonal pools of guanidine-dependent mutants of type 1 poliovirus. The mutation frequency was (6.5 +/- 6.3) x 10(-4) (with all amino acid substitutions occurring at position 227). The minimal corrected base substitution frequency per single nucleotide site in the codon for amino acid 227 was (2.1 +/- 1.9) x 10(-4).

Virus mutation frequencies can be greatly underestimated by monoclonal antibody neutralization of virions

Monoclonal antibody-resistant mutants have been widely used to estimate virus mutation frequencies. We demonstrate that standard virion neutralization inevitably underestimates monoclonal antibody-resistant mutant genome frequencies of vesicular stomatitis virus, due to phenotypic masking-mixing when wild-type (wt) virions are present in thousandsfold greater numbers. We show that incorporation of antibody into the plaque overlay medium (after virus penetration at 37 degrees C) can provide accurate estimates of genome frequencies of neutral monoclonal antibody-resistant mutant viruses in wt clones. By using this method, we have observed two adjacent G----A base transition frequencies in the I3 epitope to be of the order of 10(-4) in a wt glycine codon. This appears to be slightly lower than the frequencies observed at other sites for total (viable and nonviable) virus genomes when using a direct sequence approach.

Genetic heterogeneity of the RNA genome population of the plant virus U5-TMV

The genetic heterogeneity in a population of the U5 strain of tobacco mosaic virus (U5-TMV) was studied. The T1 fingerprint characterizing a cloned population did not vary after a new cloning step in the local lesion host Nicotiana tabacum Xanthi-nc, nor during four series of 20 passages in the systemic host N. tabacum Samsum. No heterogeneity was observed among 10 clones derived from the cloned populations, while 1 of 18 clones derived from a 20-fold passaged population differed from the rest in 1 of 55 oligonucleotides. A higher heterogeneity was found in an uncloned field isolate in which 2 of 10 clones differed in type in 1 and 2 oligonucleotides, respectively. These data agree with those reported for bacterial and animal RNA viruses and are compatible with the quasi-species model for RNA populations. On the other hand, the intrapopulational heterogeneities found for U5-TMV are considerably smaller than those reported for other RNA viruses, our data showing a high genetic stability for U5-TMV.

Partial RNA sequencing of eight supposed derivatives of type 3 poliovirus/USA/Saukett/50 reveals remarkable differences between three apparent substrains

Eight supposed derivatives of type 3 poliovirus/USA/Saukett/50 could be divided in three subgroups differing from each other as much as from the independent P3/Sabin strain as judged by partial genomic sequences covering about 25% of the portion of RNA coding for the structural proteins. This suggests that strains designated as Saukett in different laboratories are derived from three separate but related American isolates of type 3 poliovirus. Deduced amino acid sequence of the "Saukett" strains revealed amino acid substitutions at all known major antigenic sites compared with P3/Sabin or P3/Finland/23127/84 strain, but also between individual Saukett strains. These substitutions may be responsible for the known antigenic differences between the studied strains.

High nucleotide substitution error frequencies in clonal pools of vesicular stomatitis virus

Nucleotide substitution error frequencies were determined for several specific guanine base positions in the genomes of cloned vesicular stomatitis virus populations. Predetermined sites were examined in coding regions for the N, M, and L proteins and at a site in the genome 5'-end regulatory region. Misincorporation frequencies were estimated to be on the order of 10(-3) to 10(-4) at all positions analyzed. Isolates taken from virus populations after disruption of equilibrium conditions displayed replicase fidelity similar to that of cloned wild-type vesicular stomatitis virus. These mutation frequencies apply to all virus genomes present, including viruses rendered nonviable by lethal mutations. At one selected site in the N gene, two of three G----N base substitutions generated lethal nonsense mutations, yet their frequency was also very high. Biological implications for rapid virus evolution are discussed.

Extreme heterogeneity in populations of vesicular stomatitis virus

Vesicular stomatitis virus (VSV) sequence evolution and population heterogeneity were examined by T1 oligonucleotide mapping. Individual clones isolated from clonal pools of wild-type Indiana serotype VSV displayed identical T1 maps. This was observed even after one passage at high concentrations of the potent viral mutagen 5-fluorouracil. Under low-multiplicity passage conditions, the consensus T1 fingerprint of this virus remained unchanged after 523 passages. Interestingly, however, individual clones from this population (passage 523) differed significantly from each other and from consensus sequence. When virus population equilibria were disrupted by high-multiplicity passage (in which defective interfering particle interference is maximized) or passage in the presence of mutagenic levels of 5-fluorouracil, rapid consensus sequence evolution occurred and extreme population heterogeneity was observed (with some members of these population differing from others at hundreds of genome positions). A limited sampling of clones at one stage during high-multiplicity passages suggested the presence of at least several distinct master sequences, the related subpopulations of which exhibit at least transient competitive fitness within the total virus population (M. Eigen and C.K. Biebricher, p. 211-245, in E. Domingo, J.J. Holland, P. Ahlquist, ed., RNA Genetics, vol. 3, 1988). These studies further demonstrate the important role of selective pressure in determining the genetic composition of RNA virus populations. This is true under equilibrium conditions in which little consensus sequence evolution is observed owing to stabilizing selection as well as under conditions in which selective pressure is driving rapid RNA virus genome evolution.

Glycoprotein evolution of vesicular stomatitis virus New Jersey

A T1 ribonuclease fingerprinting study of a large number of virus isolates had previously demonstrated that considerable genetic variability existed among natural isolates of the vesicular stomatitis virus (VSV) New Jersey (NJ) serotype [S.T. Nichol (1988) J. Virol. 62, 572-579]. Based on these results, 34 virus isolates were chosen as representing the extent of genetic diversity within the VSV NJ serotype. We report the entire glycoprotein (G) gene nucleotide sequence and the deduced amino acid sequence for each of these viruses. Up to 19.8% G gene sequence differences could be seen among NJ serotype isolates. Analysis of the distribution of nucleotide substitutions relative to nucleotide codon position revealed that third position changes were distributed randomly throughout the gene. Third base changes constituted 84% of the observed nucleotide substitutions and affected 89% of the third base positions located in the G gene. Only three short oligonucleotide stretches of complete sequence conservation were observed. The remaining nucleotide changes located in the first and second positions were not distributed randomly, indicating that most of the amino acids coded by the G gene cannot be altered without reducing the fitness of the VSV NJ serotype viruses. Despite these constraints, up to 8.5% amino acid differences were observed between virus isolates. These differences were located throughout the G protein including regions adjacent to defined major antibody neutralization epitopes. Apparent clusters of amino acid substitutions were present in the hydrophobic signal sequence, transmembrane domain, and within the cytoplasmic domain of the G protein. A maximum parsimony analysis of the G gene nucleotide sequences allowed construction of a phylogram indicating the evolutionary relationship of these viruses. The VSV NJ serotype appears to contain at least three distinct lineages or subtypes. All recent virus isolates from the United States and Mexico are within subtype I and appear to have evolved from an ancestor more closely related to the Hazelhurst historic strain than other older strains. The implications of these findings for the evolution, epizootiology, and classification of these viruses are discussed.

RNA virus evolution and the control of viral disease

RNA viruses and other RNA genetic elements must be viewed as organized distributions of sequences termed quasi-species. This means that the viral genome is statistically defined but individually indeterminate. Stable distributions may be maintained for extremely long time periods under conditions of population equilibrium. Perturbation of equilibrium results in rapid distribution shifts. This genomic organization has many implications for viral pathogenesis and disease control. This review has emphasized the problem of selection of viral mutants resistant to antiviral drugs and the current difficulties encountered in the design of novel synthetic vaccines. Possible strategies for antiviral therapy and vaccine development have been discussed.

Determination of the mutation rate of a retrovirus

The mutation rate of Rous sarcoma virus (RSV) was measured. Progeny descended from a single virion were collected after one replication cycle, and seven regions of the genome were analyzed for mutations by denaturing-gradient gel electrophoresis. In all, 65,250 nucleotides were screened, yielding nine mutations, and the RSV mutation rate was calculated as 1.4 x 10(-4) mutations per nucleotide per replication cycle. These results indicate that RSV is an extremely mutable virus. We speculate that the mutation rate of a virus may correlate inversely with the effectiveness of vaccination against a given virus and suggest that prevention of retrovirus-mediated disease via vaccination may prove difficult.

RNA genome stability of Toscana virus during serial transovarial transmission in the sandfly Phlebotomus perniciosus

We have carried out a T1 ribonuclease fingerprinting analysis of the RNA genomes of Toscana virus isolates from successive generations of an experimentally virus-infected laboratory colony of Phlebotomus perniciosus sandflies. This analysis detected no virus RNA genome changes during transovarial transmission of the virus over 12 sandfly generations (a period of almost 2 years). These results demonstrate that although RNA viruses can exhibit high rates of mutational change under a variety of conditions, Toscana virus RNA genomes can be maintained in a stable manner during repeated transovarial virus transmission in the natural insect host. The implications of these results for insect RNA virus evolution are discussed.

Evolution of sex in RNA viruses

The distribution of deleterious mutations in a population of organisms is determined by the opposing effects of two forces, mutation pressure and selection. If mutation rates are high, the resulting mutation-selection balance can generate a substantial mutational load in the population. Sex can be advantageous to organisms experiencing high mutation rates because it can either buffer the mutation-selection balance from genetic drift, thus preventing any increases in the mutational load (Muller, 1964: Mut. Res. 1, 2), or decrease the mutational load by increasing the efficiency of selection (Crow, 1970: Biomathematics 1, 128). Muller's hypothesis assumes that deleterious mutations act independently, whereas Crow's hypothesis assumes that deleterious mutations interact synergistically, i.e., the acquisition of a deleterious mutation is proportionately more harmful to a genome with many mutations than it is to a genome with a few mutations. RNA viruses provide a test for these two hypotheses because they have extremely high mutation rates and appear to have evolved specific adaptations to reproduce sexually. Population genetic models for RNA viruses show that Muller's and Crow's hypotheses are also possible explanations for why sex is advantageous to these viruses. A re-analysis of published data on RNA viruses that are cultured by undiluted passage suggests that deleterious mutations in such viruses interact synergistically and that sex evolved there as a mechanism to reduce the mutational load.

Genetic diversity of enzootic isolates of vesicular stomatitis virus New Jersey

The RNA genomes of 43 vesicular stomatitis virus (VSV) isolates of the New Jersey (NJ) serotype were T1-ribonuclease fingerprinted to compare the extent of genetic diversity of virus from regions of epizootic and enzootic disease activity. Forty of these viruses were obtained from Central America during 1982 to 1985. The other three were older isolates, including a 1970 isolate from Culex nigripalpus mosquitos in Guatemala, a 1960 bovine isolate from Panama, and a 1976 isolate from mosquitos (Mansonia indubitans) in Ecuador. The data indicate that extensive genetic diversity exists among virus isolates from this predominantly enzootic disease zone. Six distinct T1 fingerprint groups were identified for the Central American VSV NJ isolates from 1982 to 1985. The 1960 VSV NJ isolate from Panama and the 1976 isolate from Ecuador formed two additional distinct fingerprint groups. This finding is in sharp contrast to the relatively close genetic relationship existing among VSV NJ isolates obtained from predominantly epizootic disease areas of the United States and Mexico during the same period (S. T. Nichol, J. Virol. 61:1029-1036, 1987). In this previous study, RNA genome T1 fingerprint differences were observed among isolates from different epizootics; however, the isolates were all clearly members of one large T1 fingerprint group. The eight T1 fingerprint groups described here for Central American and Ecuadorian viruses are distinct from those characterized earlier for virus isolates from the United States and Mexico and for the common laboratory virus strains Ogden and Hazelhurst. Despite being isolated 14 years earlier, the 1970 insect isolate from Guatemala is clearly a member of one of the 1982 to 1985 Central American virus fingerprint groups. This indicates that although virus genetic diversity in the region is extensive, under certain natural conditions particular virus genotypes can be relatively stably maintained for an extended period. The implications of these findings for the evolution of VSV NJ and epizootiology of the disease are discussed.

Molecular epizootiology and evolution of vesicular stomatitis virus New Jersey

Vesicular stomatitis virus (VSV) has been shown previously to be capable of undergoing rapid mutational change during sequential experimental infections in various tissue culture cell systems (J. Holland, K. Spindler, F. Horodyski, E. Grabau, S. Nichol, and S. Vandepol, Science 215:1577-1585, 1982). The present study was undertaken to determine the degree of genetic diversity and evolution of the virus under natural infection conditions and to gain insight into the epizootiology of the disease. Between 1982 and 1985, numerous outbreaks of VSV of the New Jersey serotype were reported throughout regions of the United States and Mexico. A T1 RNase fingerprint analysis was performed on the RNA genomes of 43 virus isolates from areas of epizootic and enzootic virus activity. This indicates that virus populations were genetically relatively homogeneous within successive U.S. virus epizootics. The data included virus isolates from different epizootic stages, geographical locations, host animals, and host lesion sites. In contrast, only distant genome RNA T1 fingerprint similarities were observed among viruses of the different U.S. epizootics. However, Mexican viruses isolated before or concurrent with U.S. epizootics had very similar RNA genome fingerprints, suggesting that Mexico may have been the possible origin of virus initiating recent U.S. VSV New Jersey outbreaks. Comparison of T1 fingerprints of viruses with enzootic disease areas revealed a greater extent of virus genetic diversity in these areas relative to that observed in epizootic areas. The evolutionary significance of these findings and their relationship to experimental data on VSV evolution are discussed.

Continuing coevolution of virus and defective interfering particles and of viral genome sequences during undiluted passages: virus mutants exhibiting nearly complete resistance to formerly dominant defective interfering particles

We quantitatively analyzed the interference interactions between defective interfering (DI) particles and mutants of cloned vesicular stomatitis virus passaged undiluted hundreds of times in BHK-21 cells. DI particles which predominated at different times in these serial passages always interfered most strongly (and very efficiently) with virus isolated a number of passages before the isolation of the DI particles. Virus isolated at the same passage level as the predominant DI particles usually exhibited severalfold resistance to these DI particles. Virus mutants (Sdi- mutants) isolated during subsequent passages always showed increasing resistance to these DI particles, followed by decreasing resistance as new DI particles arose to predominate and exert their own selective pressures on the virus mutant population. It appears that such coevolution of virus and DI particle populations proceeds indefinitely through multiple cycles of selection of virus mutants resistant to a certain DI particle (or DI particle class), followed by mutants resistant to a newly predominant DI particle, etc. At the peak of resistance, virus mutants were isolated which were essentially completely resistant to a particular DI particle; i.e., they were several hundred thousand-fold resistant, and they formed plaques of normal size and numbers in the presence of extremely high multiplicities of the DI particle. However, they were sensitive to interference by other DI particles. Recurring population interactions of this kind can promote rapid virus evolution. Complete sequencing of the N (nucleocapsid) and NS (polymerase associated) genes of numerous Sdi- mutants collected at passage intervals showed very few changes in the NS protein, but the N gene gradually accumulated a series of stable nucleotide and amino acid substitutions, some of which correlated with extensive changes in the Sdi- phenotype. Likewise, the 5' termini (and their complementary plus-strand 3' termini) continued to accumulate extensive base substitutions which were strikingly confined to the first 47 nucleotides. We also observed addition and deletion mutations in noncoding regions of the viral genome at a level suggesting that they probably occur at a high frequency throughout the genome, but usually with lethal or debilitating consequences when they occur in coding regions.

Evolution of human influenza A viruses over 50 years: rapid, uniform rate of change in NS gene

Variation in influenza A viruses was examined by comparison of nucleotide sequences of the NS gene (890 bases) of 15 human viruses isolated over 53 years (1933 to 1985). Changes in the genes accumulate with time, and an evolutionary tree based on the maximum parsimony method can be constructed. The evolutionary rate is approximately 2 X 10(-3) substitution per site per year in the NS genes, which is about 10(6) times the evolutionary rate of germline genes in mammals. This uniform and rapid rate of evolution in the NS gene is a good molecular clock and is compatible with the hypothesis that positive selection is operating on the hemagglutinin (or perhaps some other viral genes) to preserve random mutations in the NS gene.

Direct method for quantitation of extreme polymerase error frequencies at selected single base sites in viral RNA

Methods are described which allow direct quantitation and sequence analysis of base substitution levels at predetermined single nucleotide positions in cloned pools of an RNA virus genome or in its RNA transcripts in vitro. Base substitution frequencies for vesicular stomatitis virus (VSV) at one highly conserved site examined were reproducible and extremely high, averaging between 10(-4) and 4 X 10(-4) substitutions per base incorporated at this single site. If polymerase error frequencies averaged as high at all other sites in the 11-kilobase VSV genome, then every member of a cloned VSV population would differ from most other genomes in that clone at a number of different nucleotide positions. The preservation of a consensus sequence in such variable RNA virus genomes then could only result from strong biological selection (in a single host or multihost environment) for the most fit and competitive representatives of extremely heterogeneous virus populations.

Evolutionary variants of Rous sarcoma virus: large deletion mutants do not result from homologous recombination

Large deletion (LD) mutants of Prague strain Rous sarcoma virus, subgroup B (PrB), derived by serial undiluted passage through chicken (C/E) cells, were isolated and characterized. Individual LD viruses were initially isolated by cloning in soft agar of infected, chemically transformed quail (QT6) cells. Two regions of the PrB genome were deleted in the formation of the LD virus. This resulted in the junction of gag sequences in p12 to env sequences in gp37, and in the loss of the src gene. DNA restriction analysis of biologically active lambda Charon 27-LD recombinant clones indicated that individual LD viruses contained similar but not identical deletion endpoints. Two LD isolates, LD25 and LD85, were further subcloned into pBR322, and the deletion junctions were examined by DNA sequencing. Although the gag-env deletion endpoints were identical in the two subclones, heterogeneity was observed across the src deletion in that both mutants analyzed had the same 5' endpoint but slightly different 3' endpoints. In all cases, only a single homologous base (always an A residue) was found at the deletion endpoint. S1 nuclease analysis of the RNA from a number of QT6-LD clones gave similar results, indicating that the LD population was composed of viruses with similar but not identical deletion endpoints. Such viruses may have been generated from errors during reverse transcription of the virion RNA with subsequent selection assuring their dominance in the population.

Genome rearrangements of bovine rotavirus after serial passage at high multiplicity of infection

After serial passage at high multiplicity of infection of standard bovine rotavirus in MA104 cells, different genome rearrangements occurred in which segment 5 was lost from the RNA profile and distinct additional bands of double-stranded (ds) RNA were found in positions on gels between segments 1 and 6. It was shown that some of the additional RNA bands contained segment 5-specific sequences. The additional RNA bands were transcribed in vitro to apparent full length. Analysis of the proteins synthesized in cells infected with viruses possessing rearranged genomes showed that in all cases the product of RNA segment 5, VP5, was missing; however, in one case an abnormal protein was observed which corresponded in size to the coding capacity of the mRNA transcribed from the additional genomic RNA band. Viruses with rearranged genomes could be plaque purified, and they grew in the absence of standard virus to titers comparable to those obtained from standard virus. In mixed infections of standard virus and virus possessing genome rearrangements, standard virus overgrew during passage at low multiplicity of infection whereas virus possessing genome rearrangements overgrew during passage at high multiplicity of infection.

Vesicular stomatitis virus defective interfering particles can contain extensive genomic sequence rearrangements and base substitutions

We sequenced the 5' and 3' RNA termini of 16 defective interfering (DI) particles of vesicular stomatitis virus (VSV) isolated at intervals from persistent infections and from a series of undiluted lytic passages. All DI RNAs exhibited complementary termini, but sequences internal to these termini were extensively rearranged in a variety of ways. Despite extensive rearrangement, these internal sequences (in addition to the termini) apparently are important for DI particle interference properties. Some of these DI particles are derived from multiple intrastrand and interstrand recombination events, and the generation of each can be explained by current replicase error models. During viral evolution in persistent and acute infections, DI particles with specific termini base substitutions are selected. One DI particle exhibits a remarkable clustering of specific A----G (and complementary U----C) substitutions, apparently as a result of repetitive misincorporations by an error-prone viral polymerase complex.

Characterization of genetic changes occurring in attenuated poliovirus 2 during persistent infection in mouse central nervous systems

Genomic changes occurring in the attenuated W-2 strain of poliovirus 2 during persistent infection of the central nervous system of immunosuppressed mice were analyzed. The RNase T1 oligonucleotide fingerprints of 34 different viruses, isolated from the brains and spinal cords of paralyzed and nonparalyzed mice during a 105-day period, were used to quantitate and compare the mutations occurring in each isolate. Although mice were inoculated with plaque-purified virus, genetically distinct viruses were recovered from the central nervous system. The number of oligonucleotide changes occurring in isolates from paralyzed mice generally was greater than that observed in isolates from nonparalyzed mice. However, differences in the extent of mutation in isolates from the two groups of mice did not appear to be related to the level of virus replication. In paralyzed mice, the number of oligonucleotide changes on average was greater in viruses isolated during the first 60 days of the infection than in the last 45 days. The number of oligonucleotide changes was essentially constant throughout the infection, however, in viruses isolated from the brains of nonparalyzed mice. In addition, several specific oligonucleotide changes were found only in viruses isolated from paralyzed animals.

Structure and origin of a novel class of defective interfering particle of vesicular stomatitis virus

The genome structure and terminal sequences of a 'copyback' defective interfering (DI) particle ST1, and a novel complexly rearranged 'snapback' DI particle ST2 of vesicular stomatitis virus have been determined. The ST1 DI genome RNA possesses 54 base long inverted complementary termini, the 5' end of which is homologous to the standard virus genome 5' end. Following this region of inverted complementarity the DI RNA 5' end continues to be homologous to standard virus RNA 5' sequences, whereas the 3' end diverges into sequences within the virus L gene internal sequences. ST2 DI genome RNA does not contain colinear covalently linked plus and minus sense RNA copies of the standard infectious virus RNA 5' terminus as predicted from the prototype snapback DI structure, but instead appears to be a hairpin copy of the ST1 DI RNA genome. This is the first evidence suggesting that DI particles may be generated from RNA templates other than the standard virus RNA. Generation models and the implications of these findings for RNA virus evolution are discussed.

Multiple genetic variants arise in the course of replication of foot-and-mouth disease virus in cell culture

The genetic heterogeneity generated upon passage of foot-and-mouth disease virus (FMDV) in cell culture has been evaluated by T1-oligonucleotide fingerprinting of genomic RNA. Plaque-purified FMDV O-S7 and C-S8 were propagated by serial low multiplicity infections of BHK-21 (c-13) or IBRS-2 (c-26) cells. In independent parallel passage of the same virus, different oligonucleotide variations were fixed in the RNAs. T1-oligonucleotide fingerprinting of RNA from 34 individual viral clones derived from two passaged populations shows an extensive heterogeneity, with some mutations present in only one of the cloned genomes analyzed. Some FMDV variants are phenotypically distinct in that they yield increased progeny in infections of cell monolayers. From the number of variant sequences it can be estimated that each infectious RNA in the population differs in two to eight mutations from the average parental sequence. Thus, passaged FMDV populations consist of a pool of variants, an observation previously made with phage Q beta (E. Domingo, D. Sabo, T. Taniguchi, and C. Weissmann, Cell 13, 735-744, 1978). The FMDV genome must be described as a fluctuating distribution of sequences due to its high mutability. This may be the basis of the extensive genetic and antigenic diversity of this virus in nature.

Analysis of genetic variation in Theiler's virus during persistent infection in the mouse central nervous system

The genetic changes occurring in the BeAn strain of Theiler's murine encephalomyelitis virus (TMEV) during persistent infection in the mouse central nervous system (CNS) were studied. RNase T1-oligonucleotide fingerprinting of the RNAs of 28 BeAn viruses isolated at various times postinfection (p.i.) demonstrated that mutation occurred throughout the infection. Although plaque-purified BeAn virus was used to inoculate mice intracerebrally, genetically different viruses were recovered from the CNS. One to three oligonucleotide changes were found up to Day 152 p.i., but all three viruses isolated at Day 180 had four to nine oligonucleotide changes. No pattern of oligonucleotide changes occurring in different virus isolates was found, yet three viruses isolated from different animals at Day 180 had the same four new oligonucleotides. Overall, the number of oligonucleotide changes represented a 0.1 to 1.2% change in the virus genome. In addition, the analytical two-dimensional gel technique of P.Z. O'Farrell, H.M. Goodman, and P.H. O'Farrell (Cell 12, 1133-1142, 1977) suggested that mutation occurred in all virus isolates. In nine isolates, one to three proteins were found to have charge changes, and in general, as many nonstructural proteins had charge changes as structural proteins. P20, a nonstructural protein probably equivalent to the protease described for encephalomyocarditis virus, was found to have shifted cathodally in six different viruses. Several virus isolates had doublet patterns, suggesting the possibility that within the CNS, subpopulations existed which had proteins of slightly different charge or that virus-specified proteins had been modified after translation. Finally, antigenic variation of neutralizing site(s) on BeAn virus isolates as a way for virus to evade immune surveillance and thereby maintain the persistent state was studied. The ability of mouse serum to neutralize persisting virus isolates was not significantly different from the ability to neutralize the infecting virus. Therefore, antigenic variation does not appear to be a factor in TMEV persistence.