Extreme heterogeneity in populations of vesicular stomatitis virus

Intracellular Diversity of WNV within Circulating Avian Peripheral Blood Mononuclear Cells Reveals Host-Dependent Patterns of Polyinfection

Arthropod-borne virus (arbovirus) populations exist as mutant swarms that are maintained between arthropods and vertebrates. West Nile virus (WNV) population dynamics are host-dependent. In American crows, purifying selection is weak and population diversity is high compared to American robins, which have 100- to 1000-fold lower viremia. WNV passed in robins leads to fitness gains, whereas that passed in crows does not. Therefore, we tested the hypothesis that high crow viremia allows for higher genetic diversity within individual avian peripheral blood mononuclear cells (PBMCs), reasoning that this could have produced the previously observed host-specific differences in genetic diversity and fitness. Specifically, we infected cells and birds with a molecularly barcoded WNV and sequenced viral RNA from single cells to quantify the number of WNV barcodes in each. Our results demonstrate that the richness of WNV populations within crows far exceeds that in robins. Similarly, rare WNV variants were maintained by crows more frequently than by robins. Our results suggest that increased viremia in crows relative to robins leads to the maintenance of defective genomes and less prevalent variants, presumably through complementation. Our findings further suggest that weaker purifying selection in highly susceptible crows is attributable to this higher viremia, polyinfections and complementation.

Intracellular diversity of WNV within circulating avian peripheral blood mononuclear cells reveals host-dependent patterns of polyinfection

Error-prone replication of RNA viruses generates the genetic diversity required for adaptation within rapidly changing environments. Thus, arthropod-borne virus (arbovirus) populations exist in nature as mutant swarms that are maintained between arthropods and vertebrates. Previous studies have demonstrated that West Nile virus (WNV) population dynamics are host dependent: In American crows, which experience extremely high viremia, purifying selection is weak and population diversity is high compared to American robins, which have 100 to 1000-fold lower viremia. WNV passed in robins experiences fitness gains, whereas that passed in crows does not. Therefore, we tested the hypothesis that high crow viremia allows higher genetic diversity within individual avian peripheral-blood mononuclear cells (PBMCs), reasoning that this could have produced the previously observed host-specific differences in genetic diversity and fitness. Specifically, we infected cells and birds with a novel, barcoded version of WNV and sequenced viral RNA from single cells to quantify the number of WNV barcodes that each contained. Our results demonstrate that the richness of WNV populations within crows far exceeds that in robins. Similarly, rare WNV variants were maintained by crows more frequently than by robins. Our results suggest that increased viremia in crows relative to robins leads to maintenance of defective genomes and less prevalent variants, presumably through complementation. Our findings further suggest that weaker purifying selection in highly susceptible crows is attributable to this higher viremia, polyinfections and complementation. These studies further document the role of particular, ecologically relevant hosts in shaping virus population structure.

Author Summary

WNV mutational diversity in vertebrates is species-dependent. In crows, low frequency variants are common, and viral populations are more diverse. In robins, fewer mutations become permanent fixtures of the overall viral population. We infected crows, robins and a chicken cell line with a genetically marked (barcoded) WNV. Higher levels of virus led to multiple unique WNV genomes infecting individual cells, even when a genotype was present at low levels in the input viral stock. Our findings suggest that higher levels of circulating virus in natural hosts allow less fit viruses to survive in RNA virus populations through complementation by more fit viruses. This is significant as it allows less represented and less fit viruses to be maintained at low levels until they potentially emerge when virus environments change. Overall our data reveal new insights on the relationships between host susceptibility to high viremia and virus evolution.

Oral Rabies Vaccine Strain SPBN GASGAS: Genetic Stability after Serial In Vitro and In Vivo Passaging

Oral vaccination of wildlife has shown to be a very effective management tool in rabies control. Evaluation of the genetic stability of vaccine viruses before distributing vaccine baits in the environment is essential because all available oral rabies vaccines, including the genetically engineered rabies virus vaccine strain SPBN GASGAS (Rabitec), are based on replication-competent viruses. To evaluate the genetic stability of this vaccine strain, five serial passages of the Master Seed Virus (MSV) in the production cell line BHK21 Cl13 were performed. Furthermore, to test possible reversion to virulence, a back-passage study in suckling mouse brain (SMB) was performed. Subsequently, the pooled 5th SMB passage was inoculated intracerebrally (i.c.) in adult and suckling mice. The full genome sequences of the isolated 5th passage, in vivo and in vitro, were compared at both the consensus and the quasispecies level with the MSV. Additionally, the full genome sequence of the 6th SMB passage from the individual animals was determined and compared. Full-length integration of the double glycoprotein and modified base substitutions at amino acid position 194 and 333 of the glycoprotein could be verified in all 5th and 6th passage samples. Overall, 11 single nucleotide polymorphisms (SNPs) were detected in the 5th pooled SMB passage, 4 with frequency between 10 and 20%, and 7 with between 2.5 and 10%. SNPs that resulted in amino acid exchange were found in genes: N (one SNP), G (four SNPs), and L (three SNPs). However, none of these SNPs were associated with reversion to virulence since all adult mice inoculated i.c. with this material survived. In the individual samples of the 6th SMB passage 24 additional SNPs (>2.5%) were found, of which only 1 SNP (L-gene, position 6969) had a prevalence of >50% in 3 of 17 samples. The obtained results confirmed the stable expression of genetic modifications and the genetic stability of the consensus strain after serial in vivo and in vitro passaging.

Coxsackievirus B3-Its Potential as an Oncolytic Virus

Oncolytic virotherapy represents one of the most advanced strategies to treat otherwise untreatable types of cancer. Despite encouraging developments in recent years, the limited fraction of patients responding to therapy has demonstrated the need to search for new suitable viruses. Coxsackievirus B3 (CVB3) is a promising novel candidate with particularly valuable features. Its entry receptor, the coxsackievirus and adenovirus receptor (CAR), and heparan sulfate, which is used for cellular entry by some CVB3 variants, are highly expressed on various cancer types. Consequently, CVB3 has broad anti-tumor activity, as shown in various xenograft and syngeneic mouse tumor models. In addition to direct tumor cell killing the virus induces a strong immune response against the tumor, which contributes to a substantial increase in the efficiency of the treatment. The toxicity of oncolytic CVB3 in healthy tissues is variable and depends on the virus strain. It can be abrogated by genetic engineering the virus with target sites of microRNAs. In this review, we present an overview of the current status of the development of CVB3 as an oncolytic virus and outline which steps still need to be accomplished to develop CVB3 as a therapeutic agent for clinical use in cancer treatment.

Kinetics of Oncolytic Reovirus T3D Replication and Growth Pattern in Mesenchymal Stem Cells

Objective

Currently, application of oncolytic-virus in cancer treatment of clinical trials are growing. Oncolytic-reovirus is an attractive anti-cancer therapeutic agent for clinical testing. Many studies used mesenchymal stem cells (MSCs) as a carrier cell to enhance the delivery and quality of treatment with oncolytic-virotherapy. But, biosynthetic capacity and behavior of cells in response to viral infections are different. The infecting process of reoviruses takes from two-hours to one-week, depends on host cell and the duration of different stages of virus replication cycle. The latter includes the binding of virus particle, entry, uncoating, assembly and release of progeny-viruses. We evaluated the timing and infection cycle of reovirus type-3 strain Dearing (T3D), using one-step replication experiment by molecular and conventional methods in MSCs and L929 cell as control.

Materials and Methods

In this experimental study, L929 and adipose-derived MSCs were infected with different multiplicities of infection (MOI) of reovirus T3D. At different time points, the quantity of progeny viruses has been measured using virus titration assay and quantitative real-time polymerase chain reaction (qRT-PCR) to investigate the ability of these cells to support the reovirus replication. One-step growth cycle were examined by 50% cell culture infectious dose (CCID₅₀) and qRT-PCR.

Results

The growth curve of reovirus in cells shows that MOI: 1 might be optimal for virus production compared to higher and lower MOIs. The maximum quantity of virus production using MOI: 1 was achieved at 48-hours post-infection. The infectious virus titer became stationary at 72-hours post-infection and then gradually decreased. The virus cytopathic effect was obvious in MSCs and this cells were susceptible to reovirus infection and support the virus replication.

Conclusion

Our data highlights the timing schedule for reovirus replication, kinetics models and burst size. Further investigation is recommended to better understanding of the challenges and opportunities, for using MSCs loaded with reovirus in cancer-therapy.

A unique intra-molecular fidelity-modulating mechanism identified in a viral RNA-dependent RNA polymerase

Typically not assisted by proofreading, the RNA-dependent RNA polymerases (RdRPs) encoded by the RNA viruses may need to independently control its fidelity to fulfill virus viability and fitness. However, the precise mechanism by which the RdRP maintains its optimal fidelity level remains largely elusive. By solving 2.1-2.5 Å resolution crystal structures of the classical swine fever virus (CSFV) NS5B, an RdRP with a unique naturally fused N-terminal domain (NTD), we identified high-resolution intra-molecular interactions between the NTD and the RdRP palm domain. In order to dissect possible regulatory functions of NTD, we designed mutations at residues Y471 and E472 to perturb key interactions at the NTD-RdRP interface. When crystallized, some of these NS5B interface mutants maintained the interface, while the others adopted an 'open' conformation that no longer retained the intra-molecular interactions. Data from multiple in vitro RdRP assays indicated that the perturbation of the NTD-RdRP interactions clearly reduced the fidelity level of the RNA synthesis, while the processivity of the NS5B elongation complex was not affected. Collectively, our work demonstrates an explicit and unique mode of polymerase fidelity modulation and provides a vivid example of co-evolution in multi-domain enzymes.

Seasonal Genetic Drift of Human Influenza A Virus Quasispecies Revealed by Deep Sequencing

After a pandemic wave in 2009 following their introduction in the human population, the H1N1pdm09 viruses replaced the previously circulating, pre-pandemic H1N1 virus and, along with H3N2 viruses, are now responsible for the seasonal influenza type A epidemics. So far, the evolutionary potential of influenza viruses has been mainly documented by consensus sequencing data. However, like other RNA viruses, influenza A viruses exist as a population of diverse, albeit related, viruses, or quasispecies. Interest in this quasispecies nature has increased with the development of next generation sequencing (NGS) technologies that allow a more in-depth study of the genetic variability. NGS deep sequencing methodologies were applied to determine the whole genome genetic heterogeneity of the three categories of influenza A viruses that circulated in humans between 2007 and 2012 in France, directly from clinical respiratory specimens. Mutation frequencies and single nucleotide polymorphisms were used for comparisons to address the level of natural intrinsic heterogeneity of influenza A viruses. Clear differences in single nucleotide polymorphism profiles between seasons for a given subtype also revealed the constant genetic drift that human influenza A virus quasispecies undergo.

Kinetic Modeling of Virus Growth in Cells

When a virus infects a host cell, it hijacks the biosynthetic capacity of the cell to produce virus progeny, a process that may take less than an hour or more than a week. The overall time required for a virus to reproduce depends collectively on the rates of multiple steps in the infection process, including initial binding of the virus particle to the surface of the cell, virus internalization and release of the viral genome within the cell, decoding of the genome to make viral proteins, replication of the genome, assembly of progeny virus particles, and release of these particles into the extracellular environment. For a large number of virus types, much has been learned about the molecular mechanisms and rates of the various steps. However, in only relatively few cases during the last 50 years has an attempt been made-using mathematical modeling-to account for how the different steps contribute to the overall timing and productivity of the infection cycle in a cell. Here we review the initial case studies, which include studies of the one-step growth behavior of viruses that infect bacteria (Qβ, T7, and M13), human immunodeficiency virus, influenza A virus, poliovirus, vesicular stomatitis virus, baculovirus, hepatitis B and C viruses, and herpes simplex virus. Further, we consider how such models enable one to explore how cellular resources are utilized and how antiviral strategies might be designed to resist escape. Finally, we highlight challenges and opportunities at the frontiers of cell-level modeling of virus infections.

Resistant mutations and quasispecies complexity of hepatitis B virus during telbivudine treatment

Ultra-deep pyrosequencing (UDPS) was used to analyse the dynamics of quasispecies and resistant mutations during telbivudine (LDT) treatment of hepatitis B patients. Twenty-six HBeAg-positive chronic hepatitis B patients were treated with LDT for a period of 104 weeks and were characterized as 16 responders, six partial responders and four viral breakthrough patients based on hepatitis B virus (HBV) DNA levels. The plasma samples were subjected to UDPS of the reverse transcriptase (RT) region of HBV. Mutations rtM204I, rtL80I and rtL80V were detected in at least three of the four viral breakthrough patients, indicating the significant roles of the mutations in resistance to LDT. The degree of complexity of viral quasispecies remained in a steady state in the absence of selection pressure, but increased after the LDT treatment. The complexity in the responder group at week 12 was significantly higher than that in the group comprising partial responders and viral breakthrough patients. In vitro replication efficiency analyses showed that the RT mutations had different impacts on HBV replication, with a tendency of rtM204I>rtL80V>rtL80I. Furthermore, double mutations rtL80I/M204I and rtL80V/M204V had replication efficiency similar to that of rtL80I and rtL80V, respectively. Consistent with previous studies, mutation rtM204I was found to be highly resistant to LDT. However, in contrast with their sensitivity to lamivudine, rtL80I and rtL80V were moderately resistant to LDT. Our results indicated that rtL80I and rtL80V may not only serve as replication complementary mutations to rtM204I, but also directly contribute to the LDT resistance.

Statistical properties and error threshold of quasispecies on single-peak Gaussian-distributed fitness landscapes

The stochastic Eigen model proposed by Feng et al. (2007) (Journal of Theoretical Biology, 246, 28) showed that error threshold is no longer a phase transition point but a crossover region whose width depends on the strength of the random fluctuation in an environment. The underlying cause of this phenomenon has not yet been well examined. In this article, we adopt a single peak Gaussian distributed fitness landscape instead of a constant one to investigate and analyze the change of the error threshold and the statistical property of the quasi-species population. We find a roughly linear relation between the width of the error threshold and the fitness fluctuation strength. For a given quasi-species, the fluctuation of the relative concentration has a minimum with a normal distribution of the relative concentration at the maximum of the averaged relative concentration, it has however a largest value with a bimodal distribution of the relative concentration near the error threshold. The above results deepen our understanding of the quasispecies and error threshold and are heuristic for exploring practicable antiviral strategies.

Effect of truncation of the N-terminal region of the viral hemorrhagic septicemia virus (VHSV) P protein on viral replication

The phosphoprotein (P) of viral hemorrhagic septicemia virus (VHSV) plays an essential role in viral replication by interconnecting the L protein and the N protein-RNA complex. In this study, to investigate the role of the N-terminal region of the P protein in viral replication, we mutated the first or the first and second or the first, second, and third ATG codon into TGA stop codons. The respective mutants were named P1, P2, and P3. Recombinant VHSVs containing each mutated P gene (rVHSV-P1, -P2, and -P3) were successfully generated by supplying the intact P protein in trans. The rVHSV-P2 and -P3 were not generated from cells expressing truncated P proteins (P1, P2 or P3 protein), but the rVHSV-P1 produced infectious viruses, even in cells without any P-protein-expressing plasmids. Nucleotide sequence analysis of the P gene of rVHSV-P1 showed that a mutation had occurred that resulted in the fourth amino acid (isoleucine, ATT) being changed to methionine (ATG) without a frameshift (P0.5), suggesting that strong selection pressure might facilitate mutations that are advantageous or essential for virus replication. Infectious rVHSV-P2 and -P3 were produced in cells expressing the P0.5 protein, suggesting that the first three amino acids of the P protein of VHSV are dispensable for viral replication. Furthermore, although the P1 protein was shorter than the P0.5 protein by only two amino acid residues, no viruses were produced when the P1 protein was supplied indicating that the fourth and the fifth amino acid residues are indispensable for normal P protein functions involved in viral replication.

Error catastrophe in populations under similarity-essential recombination

Organisms are often more likely to exchange genetic information with others that are similar to themselves. One of the most widely accepted mechanisms of RNA virus recombination requires substantial sequence similarity between the parental RNAs and is termed similarity-essential recombination. This mechanism may be considered analogous to assortative mating, an important form of non-random mating that can be found in animals and plants. Here we study the dynamics of haplotype frequencies in populations evolving under similarity-essential recombination. Haplotypes are represented by a genome of B biallelic loci and the Hamming distance between individuals is used as a criterion for recombination. We derive the evolution equations for the haplotype frequencies assuming that recombination does not occur if the genetic distance is larger than a critical value G and that mutation occurs at a rate μ per locus. Additionally, uniform crossover is considered. Although no fitness is directly associated to the haplotypes, we show that frequency-dependent selection emerges dynamically and governs the haplotype distribution. A critical mutation rate μc can be identified as the error threshold transition, beyond which this selective information cannot be stored. For μ<μc the distribution consists of a dominant sequence surrounded by a cloud of closely related sequences, characterizing a quasispecies. For μ>μc the distribution becomes uniform, with all haplotypes having the same frequency. In the case of extreme assortativeness, where individuals only recombine with others identical to themselves (G=0), the error threshold results μc=1/4, independently of the genome size. For weak assortativity (G=B-1)μc=2(-(B+1)) and for the case of no assortativity (G=B) μc=0. We compute the mutation threshold for 0<G<B and show that, for large B, it depends only on the ratio G/B. We discuss the consequences of these results for recombination in viruses and for speciation.

Estimating Fitness of Viral Quasispecies from Next-Generation Sequencing Data

The quasispecies model is ubiquitous in the study of viruses. While having lead to a number of insights that have stood the test of time, the quasispecies model has mostly been discussed in a theoretical fashion with little support of data. With next-generation sequencing (NGS), this situation is changing and a wealth of data can now be produced in a time- and cost-efficient manner. NGS can, after removal of technical errors, yield an exceedingly detailed picture of the viral population structure. The widespread availability of cross-sectional data can be used to study fitness landscapes of viral populations in the quasispecies model. This chapter highlights methods that estimate the strength of selection in selective sweeps, assesses marginal fitness effects of quasispecies, and finally infers the fitness landscape of a viral quasispecies, all on the basis of NGS data.

Design requirements for interfering particles to maintain coadaptive stability with HIV-1

Defective interfering particles (DIPs) are viral deletion mutants lacking essential transacting or packaging elements and must be complemented by wild-type virus to propagate. DIPs transmit through human populations, replicating at the expense of the wild-type virus and acting as molecular parasites of viruses. Consequently, engineered DIPs have been proposed as therapies for a number of diseases, including human immunodeficiency virus (HIV). However, it is not clear if DIP-based therapies would face evolutionary blocks given the high mutation rates and high within-host diversity of lentiviruses. Divergent evolution of HIV and DIPs appears likely since natural DIPs have not been detected for lentiviruses, despite extensive sequencing of HIVs and simian immunodeficiency viruses (SIVs). Here, we tested if the apparent lack of lentiviral DIPs is due to natural selection and analyzed which molecular characteristics a DIP or DIP-based therapy would need to maintain coadaptive stability with HIV-1. Using a well-established mathematical model of HIV-1 in a host extended to include its replication in a single cell and interference from DIP, we calculated evolutionary selection coefficients. The analysis predicts that interference by codimerization between DIPs and HIV-1 genomes is evolutionarily unstable, indicating that recombination between DIPs and HIV-1 would be selected against. In contrast, DIPs that interfere via competition for capsids have the potential to be evolutionarily stable if the capsid-to-genome production ratio of HIV-1 is >1. Thus, HIV-1 variants that attempt to "starve" DIPs to escape interference would be selected against. In summary, the analysis suggests specific experimental measurements that could address the apparent lack of naturally occurring lentiviral DIPs and specifies how therapeutic approaches based on engineered DIPs could be evolutionarily robust and avoid recombination.

Second-site mutations selected in transcriptional regulatory sequences compensate for engineered mutations in the vesicular stomatitis virus nucleocapsid protein

The active template for RNA synthesis for vesicular stomatitis virus (VSV) and other negative-strand viruses is the RNA genome in association with the nucleocapsid (N) protein. The N protein molecules sequester the genomic RNA and are linked together by a network of noncovalent interactions. We previously demonstrated that mutations predicted to weaken interactions between adjacent N protein molecules altered the levels of RNA synthesis directed from subgenomic ribonucleoprotein (RNP) templates. To determine if these mutations affect virus replication, recombinant viruses containing single-amino-acid substitutions in the N protein were recovered. Four mutations altered transcription and genome replication levels, perturbed viral protein synthesis, and inhibited virus replication. Selective pressure for improved virus replication was applied by eight sequential passages. After five passages, virus replication improved and RNA synthesis recovered concomitantly with the restoration of the protein molar ratios to near-wild-type levels. Genome sequences were compared before and after passage to determine whether compensatory mutations were selected and to potentially identify interactions between N protein molecules or between the RNP template and the viral polymerase. Improved virus replication correlated with the selection of additional mutations located in cis-acting transcriptional regulatory sequences at the gene junctions of the genome rather than in coding sequences, with one exception. The engineered N gene mutations perturbed mRNA and protein expression levels, but the selection of modified transcriptional regulatory sequences with passage facilitated the restoration of wild-type protein expression by modulating transcription levels, reflecting the adaptability and versatility of gene regulation by transcriptional control.

Exploring the fitness landscape of an RNA virus by using a universal barcode microarray

Studies of viral pathogenesis have relied heavily on analyses of specific clones and their genetic determinants of virulence. It is sometimes difficult to apply this reductionist approach to the study of RNA viruses, which by virtue of their very high mutation rates, exist as a complex mixture of mutants. While quasispecies theory has provided an intellectual framework for exploring the relationship between the viral population structure and phenotype, experimental studies have been limited by the relatively poor resolution of traditional sequencing-based approaches. We have addressed this problem by developing a molecular barcoding strategy in which viral subpopulations are tagged with unique 20-nucleotide sequences. The behavior of these subpopulations can be monitored using a universal barcode microarray. We demonstrate the performance of our barcode microarray platform using poliovirus, a model RNA virus. Using this platform, we explored the fitness landscape occupied by an artificial quasispecies consisting of 48 randomly mutagenized clones. We were able to rapidly derive precise fitness measurements for a majority of these clones and identified a neutral space surrounding the wild type. The experimental paradigm presented here is readily adaptable to other viral systems and can potentially be used to track thousands of variants in a cost-effective manner.

Quasispecies-like behavior observed in catalytic RNA populations evolving in a test tube

BACKGROUND

During the RNA World, molecular populations were probably very small and highly susceptible to the force of strong random drift. In conjunction with Muller's Ratchet, this would have imposed difficulties for the preservation of the genetic information and the survival of the populations. Mechanisms that allowed these nascent populations to overcome this problem must have been advantageous.

RESULTS

Using continuous in vitro evolution experimentation with an increased mutation rate imposed by MnCl2, it was found that clonal 100-molecule populations of ribozymes clearly exhibit certain characteristics of a quasispecies. This is the first time this has been seen with a catalytic RNA. Extensive genotypic sampling from two replicate lineages was gathered and phylogenetic networks were constructed to elucidate the structure of the evolving RNA populations. A common distribution was found in which a mutant sequence was present at high frequency, surrounded by a cloud of mutant with lower frequencies. This is a typical distribution of quasispecies. Most of the mutants in these clouds were connected by short Hamming distance values, indicating their close relatedness.

CONCLUSIONS

The quasispecies nature of mutant RNA clouds facilitates the recovery of genotypes under pressure of being removed from the population by random drift. The empirical populations therefore evolved a genotypic resiliency despite a high mutation rate by adopting the characteristics of quasispecies, implying that primordial RNA pools could have used this strategy to avoid extinction.

The quasispecies nature and biological implications of the hepatitis C virus

Many RNA viruses exist as a cloud of closely related sequence variants called a quasispecies, rather than as a population of identical clones. In this article, we explain the quasispecies nature of RNA viral genomes, and briefly review the principles of quasispecies dynamics and the differences with classical population genetics. We then discuss the current methods for quasispecies analysis and conclude with the biological implications of this phenomenon, focusing on the hepatitis C virus.

Growth of an RNA virus in single cells reveals a broad fitness distribution

Genetic and environmental factors will influence the growth of an RNA virus, but their relative contributions are challenging to resolve because standard culture methods mask how virus particles interact with individual host cells. Here, single particles of vesicular stomatitis virus, a prototype RNA virus, were used to infect individual BHK cells. Infected cells produced 50 to 8000 progeny virus particles, but these differences were lost upon subsequent culture, suggesting the diversity of yields reflected cell-to-cell differences rather than viral genetic variation. Cells infected at different phases of their cell cycle produced from 1400 (early S) to 8700 (G(2)M) infectious virus particles, coinciding with the middle-to-upper range of the observed distribution. Fluctuations in virus and cell compositions and noisy gene expression may also contribute to the broad distribution of virus yields. These findings take a step toward quantifying how environmental variation can impact the fitness distribution of an RNA virus.

Field Evaluation of a Multiplex Real-Time Reverse Transcription Polymerase Chain Reaction Assay for Detection of Vesicular Stomatitis Virus

Sporadic outbreaks of vesicular stomatitis (VS) in the United States result in significant economic losses for the U.S. livestock industries because VS is a reportable disease that clinically mimics foot-and-mouth disease. Rapid and accurate differentiation of these 2 diseases is critical because their consequences and control strategies differ radically. The objective of the current study was to field validate a 1-tube multiplexed real-time reverse transcription polymerase chain reaction (real-time RT-PCR) assay for the rapid detection of Vesicular stomatitis New Jersey virus and Vesicular stomatitis Indiana virus strains occurring in Mexico and North and Central America. A comprehensive collection of 622 vesicular lesion samples obtained from cattle, horses, and swine from throughout Mexico and Central America was tested by the real-time RT-PCR assay and virus isolation. Overall, clinical sensitivity and specificity of the real-time RT-PCR were 83% and 99%, respectively. Interestingly, VS virus isolates originating from a specific region of Costa Rica were not detected by real-time RT-PCR. Sequence comparisons of these viruses with the real-time RT-PCR probe and primers showed mismatches in the probe and forward and reverse primer regions. Additional lineage-specific primers and a probe corrected the lack of detection of the missing genetic lineage. Thus, this assay reliably identified existing Mexican and Central American VS viruses and proved readily adaptable as new VS viruses were encountered. An important secondary result of this research was the collection of hundreds of new VS virus isolates that provide a foundation from which many additional studies can arise.

Characterization and genetic diversity of sugarcane streak mosaic virus causing mosaic in sugarcane

Sixty-three sugarcane leaf samples were collected from fifty-eight sugarcane varieties, evolved from eleven major sugarcane growing states in India, Australia, South Africa and USA. In RT-PCR, using gene specific primers for sugarcane streak mosaic virus (SCSMV)-CP, 58 of 63 sugarcane samples were found positive to the virus infection and rest of the five samples were negative. Partial CP gene sequences of 42 SCSMV isolates including an isolate from aphid colony (Melanaphis indosacchari) infested on sugarcane variety from this study were characterized after cloning and sequencing for selective isolates represented by at least one isolate from each location. The new sequences identified in the study were named as SCSMV-CB isolates. Fifty two sequences including the 10 database sequences (complete CP cds) deposited earlier from this institute were compared with each other as well as GenBank database sequences of Potyviridae members viz., Rymovirus, Potyvirus, Ipomovirus, Tritimovirus and eight sequences of SCSMV reported from elsewhere. Among the SCSMV-CB isolates sequenced in the study, 85.7-100% (nucleotide) and 89.9-100% (amino acid) sequence identities were observed and with the other data base sequences of SCSMV, the respective identities were 82.2-97.5 and 89.7-98.6%. Grouping of the isolates by the maximum likelihood with molecular clock model, distributed 60 SCSMV sequences including the eight database sequences deposited by other SCSMV working groups from India and USA in 16 different phylogenetic groups. Although the isolates of SCSMV were relatively close to Ipomovirus and Tritimovirus, they were sandwiched between Rymovirus and Ipomovirus. The sequence comparison and phylogenetic studies revealed that the relatedness of SCSMV with the potyviral related genera was comparatively low to consider it as a member of earlier described potyviral genera, hence the genus "Susmovirus" (sugarcane streak mosaic virus) has been proposed, with SCSMV as the sole species to be included. The 52 SCSMV-CB isolates from this institute were distributed in 14 phylogenetic groups and the grouping pattern revealed that the virus isolates could not be grouped based on geographical origin of the host varieties or longevity of the host variety.

Selection for mutational robustness in finite populations

We investigate the evolutionary dynamics of a finite population of RNA sequences replicating on a neutral network. Despite the lack of differential fitness between viable sequences, we observe typical properties of adaptive evolution, such as increase of mean fitness over time and punctuated-equilibrium transitions, after initial mutation-selection balance has been reached. We find that a product of population size and mutation rate of approximately 30 or larger is sufficient to generate selection pressure for mutational robustness, even if the population size is orders of magnitude smaller than the neutral network on which the population resides. Our results show that quasispecies effects and neutral drift can occur concurrently, and that the relative importance of each is determined by the product of population size and mutation rate.

Adaptation of canine distemper virus to canine footpad keratinocytes modifies polymerase activity and fusogenicity through amino acid substitutions in the P/V/C and H proteins

The wild-type canine distemper virus (CDV) strain A75/17 induces a non-cytocidal infection in cultures of canine footpad keratinocytes (CFKs) but produces very little progeny virus. After only three passages in CFKs, the virus produced 100-fold more progeny and induced a limited cytopathic effect. Sequence analysis of the CFK-adapted virus revealed only three amino acid differences, of which one was located in each the P/V/C, M and H proteins. In order to assess which amino acid changes were responsible for the increase of infectious virus production and altered phenotype of infection, we generated a series of recombinant viruses. Their analysis showed that the altered P/V/C proteins were responsible for the higher levels of virus progeny formation and that the amino acid change in the cytoplasmic tail of the H protein was the major determinant of cytopathogenicity.

Insights into the selective pressures restricting Pelargonium flower break virus genome variability: Evidence for host adaptation

The molecular diversity of Pelargonium flower break virus (PFBV) was assessed using a collection of isolates from different geographical origins, hosts, and collecting times. The genomic region examined was 1,828 nucleotides (nt) long and comprised the coding sequences for the movement (p7 and p12) and the coat (CP) proteins, as well as flanking segments including the entire 3' untranslated region (3' UTR). Some constraints limiting viral heterogeneity could be inferred from sequence analyses, such as the conservation of the amino acid sequences of p7 and of the shell domain of the CP, the maintenance of a leucine zipper motif in p12, and the preservation of a particular folding in the 3' UTR. A remarkable covariation, involving five specific amino acid sites, was found in the CP of isolates largely propagated in the local lesion host Chenopodium quinoa and in the progeny of a PFBV variant subjected to serial passages in this host. Concomitant with this covariation, up to 30 nucleotide substitutions in a 1,428-nt region of the viral RNA could be attributable to C. quinoa-specific adaptation, representing one of the most outstanding cases of host-driven genome variation for a plant virus. Globally, the results indicate that the selective pressures exerted by the host play a critical role in shaping PFBV populations and that these populations are likely being selected for at both protein and RNA levels.

Ribavirin and lethal mutagenesis of poliovirus: molecular mechanisms, resistance and biological implications

Positive strand RNA virus populations are a collection of similar but genetically different viruses. They exist as viral quasispecies due to the high mutation rates of the low fidelity viral RNA-dependent RNA polymerase (RdRp). It is thought that this genomic heterogeneity is advantageous to the population, allowing for adaptation to rapidly changing environments that present varying types and degrees of selective pressure. However, one consequence of this extensive diversity is the susceptibility to mutagens that further increase sequence variation. Since RNA viruses live at the edge of maximal variability, an increase in the mutation rate is likely to force the virus beyond the tolerable mutation frequency into 'error catastrophe'. One such mutagen, ribavirin, is an antiviral nucleoside analog that is mutagenic to several RNA viruses. Ribavirin is incorporated into the viral genome causing lethal mutagenesis and a subsequent decrease in the specific infectivity. Even so, passaging poliovirus in the presence of low to intermediate concentrations of the drug leads to the emergence of a viral population resistant to the effects of ribavirin. These viruses have a point mutation in the RdRp that increases the overall polymerase fidelity. Interestingly, as predicted by the quasispecies theory, ribavirin resistant viruses are less adaptable, as they are more susceptible to other non-mutagenic antiviral drugs and are highly attenuated in vivo. Here, we review the mechanism of action of ribavirin on poliovirus and other RNA viruses, the possibility for escape via increased fidelity of the viral polymerase, the consequences of this response on viral population dynamics, and the biological implications for the therapeutic use of mutagenic antiviral agents.

Incorporation fidelity of the viral RNA-dependent RNA polymerase: a kinetic, thermodynamic and structural perspective

Positive-strand RNA viruses exist as a quasi-species due to the incorporation of mutations into the viral genome during replication by the virus-encoded RNA-dependent RNA polymerase (RdRP). Therefore, the RdRP is often described as a low-fidelity enzyme. However, until recently, a complete description of the kinetic, thermodynamic and structural basis for the nucleotide incorporation fidelity of the RdRP has not been available. In this article, we review the following: (i) the steps employed by the RdRP to incorporate a correct nucleotide; (ii) the steps that are employed by the RdRP for nucleotide selection; (iii) the structure-based hypothesis for nucleotide selection; (iv) the impact of sites remote from the active site on polymerase fidelity. Given the recent observation that RNA viruses exist on the threshold of error catastrophe, the studies reviewed herein suggest novel strategies to perturb RdRP fidelity that may lead ultimately to the development of antiviral agents to treat RNA virus infection.

Co-infection weakens selection against epistatic mutations in RNA viruses

Co-infection may be beneficial in large populations of viruses because it permits sexual exchange between viruses that is useful in combating the mutational load. This advantage of sex should be especially substantial when mutations interact through negative epistasis. In contrast, co-infection may be detrimental because it allows virus complementation, where inferior genotypes profit from superior virus products available within the cell. The RNA bacteriophage phi6 features a genome divided into three segments. Co-infection by multiple phi6 genotypes produces hybrids containing reassorted mixtures of the parental segments. We imposed a mutational load on phi6 populations by mixing the wild-type virus with three single mutants, each harboring a deleterious mutation on a different one of the three virus segments. We then contrasted the speed at which these epistatic mutations were removed from virus populations in the presence and absence of co-infection. If sex is a stronger force, we predicted that the load should be purged faster in the presence of co-infection. In contrast, if complementation is more important we hypothesized that mutations would be eliminated faster in the absence of co-infection. We found that the load was purged faster in the absence of co-infection, which suggests that the disadvantages of complementation can outweigh the benefits of sex, even in the presence of negative epistasis. We discuss our results in light of virus disease management and the evolutionary advantage of haploidy in biological populations.

Recombination in primeval genomes: a step forward but still a long leap from maintaining a sizable genome

A stubborn question in early molecular evolution is how to sustain the "minimum" informational length required for the basic features of life with a putative low-copying fidelity RNA polymerase ribozyme. Proposals to circumvent the information crisis have primarily focused on networks of cooperative molecules or compartmentalization of non-cooperative unlinked templates, but success has been very limited so far. Lehman (2003) has recently suggested that recombination--a frequently ignored player in early evolution--could have been crucial to building up primeval genomes of sizable length. Here we investigate this claim by assuming (without loss of generality) that genes were already enclosed in a compartment (vesicle, protocell). The numerical results show a quite intricate interplay among mutation, recombination, and gene redundancy. Provided that the minimum number of gene copies per protocell was enough for recombination to recreate wildtype templates--but not too high to impose an unbearable burden of mutants--informational content could have increased by at least 25% by keeping the same mutational load as that for a population without recombination. However, the upper bound of informational length would still be far from the "minimal life" provisions.

A ferret model of canine distemper virus virulence and immunosuppression

Canine distemper virus (CDV) infects many carnivores, including ferrets and dogs, and is the member of the Morbillivirus genus most easily amenable to experimentation in a homologous small-animal system. To gain insights into the determinants of CDV pathogenesis, we isolated a strain highly virulent for ferrets by repeated passaging in these animals. Sequence comparison of the genome of this strain with that of its highly attenuated precursor revealed 19 mutations distributed almost evenly in the six genes. We then recovered a virus from a cDNA copy of the virulent CDV strain's consensus sequence by using a modified reverse genetics system based on B cells. We infected ferrets with this virus and showed that it fully retained virulence as measured by the timing of rash appearance, disease onset, and death. Body temperature, leukocyte number, lymphocyte proliferation activity, and cell-associated viremia also had similar kinetics. We then addressed the question of the relative importance of the envelope and other viral constituents for virulence. Viruses in which the envelope genes (matrix, fusion, and hemagglutinin) of the virulent strain were combined with the other genes of the attenuated strain caused severe rash and fever even if the disease onset was delayed. Viruses in which the nucleocapsid, polymerase, and phosphoprotein genes (coding also for the V and C proteins) of the virulent strain were combined with the envelope genes of the attenuated strain caused milder signs of disease. Thus, virulence-inducing mutations have accumulated throughout the genome.

Probability of fixation of an advantageous mutant in a viral quasispecies

The probability that an advantageous mutant rises to fixation in a viral quasispecies is investigated in the framework of multitype branching processes. Whether fixation is possible depends on the overall growth rate of the quasispecies that will form if invasion is successful rather than on the individual fitness of the invading mutant. The exact fixation probability can be calculated only if the fitnesses of all potential members of the invading quasispecies are known. Quasispecies fixation has two important characteristics: First, a sequence with negative selection coefficient has a positive fixation probability as long as it has the potential to grow into a quasispecies with an overall growth rate that exceeds that of the established quasispecies. Second, the fixation probabilities of sequences with identical fitnesses can nevertheless vary over many orders of magnitudes. Two approximations for the probability of fixation are introduced. Both approximations require only partial knowledge about the potential members of the invading quasispecies. The performance of these two approximations is compared to the exact fixation probability on a network of RNA sequences with identical secondary structure.

Evolution of mutational robustness

We review recent advances in the understanding of the mutation-selection balance of asexual replicators. For over 30 years, population geneticists thought that an expression derived by Kimura and Maruyama in 1966 fully solved this problem. However, Kimura and Maruyama's result is only correct in the absence of neutral mutations. The inclusion of neutral mutations leads to a wealth of interesting new effects, and, in particular, to a selective pressure to evolve robustness against mutations. We cover recent literature on the population dynamics of asexual replicators on networks of neutral genotypes, on the outcompetition of fast replicators by slower ones with better mutational support, and on the probability of fixation at high mutation rates. We discuss empirical evidence for the evolution of mutational robustness, and speculate on its relevance for higher organisms.

Modeling viral genome fitness evolution associated with serial bottleneck events: evidence of stationary states of fitness

Evolution of fitness values upon replication of viral populations is strongly influenced by the size of the virus population that participates in the infections. While large population passages often result in fitness gains, repeated plaque-to-plaque transfers result in average fitness losses. Here we develop a numerical model that describes fitness evolution of viral clones subjected to serial bottleneck events. The model predicts a biphasic evolution of fitness values in that a period of exponential decrease is followed by a stationary state in which fitness values display large fluctuations around an average constant value. This biphasic evolution is in agreement with experimental results of serial plaque-to-plaque transfers carried out with foot-and-mouth disease virus (FMDV) in cell culture. The existence of a stationary phase of fitness values has been further documented by serial plaque-to-plaque transfers of FMDV clones that had reached very low relative fitness values. The statistical properties of the stationary state depend on several parameters of the model, such as the probability of advantageous versus deleterious mutations, initial fitness, and the number of replication rounds. In particular, the size of the bottleneck is critical for determining the trend of fitness evolution.

Selection for fitness versus selection for robustness in RNA secondary structure folding

I investigate the competition between two quasispecies residing on two disparate neutral networks. Under the assumption that the two neutral networks have different topologies and fitness levels, it is the mutation rate that determines which quasispecies will eventually be driven to extinction. For small mutation rates, I find that the quasispecies residing on the neutral network with the lower replication rate will disappear. For higher mutation rates, however, the faster replicating sequences may be outcompeted by the slower replicating ones if the connection density on the second neutral network is sufficiently high. The analytical results are in excellent agreement with flow-reactor simulations of replicating RNA sequences.

Structure and temporal dynamics of populations within wheat streak mosaic virus isolates

Variation within the Type and Sidney 81 strains of wheat streak mosaic virus was assessed by single-strand conformation polymorphism (SSCP) analysis and confirmed by nucleotide sequencing. Limiting-dilution subisolates (LDSIs) of each strain were evaluated for polymorphism in the P1, P3, NIa, and CP cistrons. Different SSCP patterns among LDSIs of a strain were associated with single-nucleotide substitutions. Sidney 81 LDSI-S10 was used as founding inoculum to establish three lineages each in wheat, corn, and barley. The P1, HC-Pro, P3, CI, NIa, NIb, and CP cistrons of LDSI-S10 and each lineage at passages 1, 3, 6, and 9 were evaluated for polymorphism. By passage 9, each lineage differed in consensus sequence from LDSI-S10. The majority of substitutions occurred within NIa and CP, although at least one change occurred in each cistron except HC-Pro and P3. Most consensus sequence changes among lineages were independent, with substitutions accumulating over time. However, LDSI-S10 bore a variant nucleotide (G(6016)) in NIa that was restored to A(6016) in eight of nine lineages by passage 6. This near-global reversion is most easily explained by selection. Examination of nonconsensus variation revealed a pool of unique substitutions (singletons) that remained constant in frequency during passage, regardless of the host species examined. These results suggest that mutations arising by viral polymerase error are generated at a constant rate but that most newly generated mutants are sequestered in virions and do not serve as replication templates. Thus, a substantial fraction of variation generated is static and has yet to be tested for relative fitness. In contrast, nonsingleton variation increased upon passage, suggesting that some mutants do serve as replication templates and may become established in a population. Replicated mutants may or may not rise to prominence to become the consensus sequence in a lineage, with the fate of any particular mutant subject to selection and stochastic processes such as genetic drift and population growth factors.

Evidence for the Non-quasispecies Evolution of RNA Viruses

The quasispecies model of RNA virus evolution differs from those formulated in conventional population genetics in that neutral mutations do not lead to genetic drift of the population, and natural selection acts on the mutant distribution as a whole rather than on individual variants. By computer simulation, we show that this model could be inappropriate for many RNA viruses because the neutral sequence space may be too large to allow the formation of a quasispecies distribution. This view is supported by our analysis of gene sequences from vesicular stomatitis virus, which is considered a prototype RNA virus quasispecies. Our results are relevant to the evolution of RNA systems in general.

Sequences of Citrus tristeza virus separated in time and space are essentially identical

The first Citrus tristeza virus (CTV) genomes completely sequenced (19.3-kb positive-sense RNA), from four biologically distinct isolates, are unexpectedly divergent in nucleotide sequence (up to 60% divergence). Understanding of whether these large sequence differences resulted from recent evolution is important for the design of disease management strategies, particularly the use of genetically engineered mild (essentially symptomless)-strain cross protection and RNA-mediated transgenic resistance. The complete sequence of a mild isolate (T30) which has been endemic in Florida for about a century was found to be nearly identical to the genomic sequence of a mild isolate (T385) from Spain. Moreover, samples of sequences of other isolates from distinct geographic locations, maintained in different citrus hosts and also separated in time (B252 from Taiwan, B272 from Colombia, and B354 from California), were nearly identical to the T30 sequence. The sequence differences between these isolates were within or near the range of variability of the T30 population. A possible explanation for these results is that the parents of isolates T30, T385, B252, B272, and B354 have a common origin, probably Asia, and have changed little since they were dispersed throughout the world by the movement of citrus. Considering that the nucleotide divergence among the other known CTV genomes is much greater than those expected for strains of the same virus, the remarkable similarity of these five isolates indicates a high degree of evolutionary stasis in some CTV populations.

Search for the mechanism of genetic variation in the pro gene of human immunodeficiency virus

To study the mechanism of evolution of the human immunodeficiency virus (HIV) protease gene (pro), we analyzed a database of 213 pro sequences isolated from 11 HIV type 1-infected patients who had not been treated with protease inhibitors. Variation in pro is restricted to rare variable bases which are highly diverse and differ in location among individuals; an average variable base appears in about 16% of individuals. The average intrapatient distance per individual variable site, 27%, is similar for synonymous and nonsynonymous sites, although synonymous sites are twice as abundant. The latter observation excludes selection for diversity as an important, permanently acting factor in the evolution of pro and leaves purifying selection as the only kind of selection. Based on this, we developed a model of evolution, both within individuals and along the transmission chain, which explains variable sites as slightly deleterious mutants slowly reverting to the better-fit variant during individual infection. In the case of a single-source transmission, genetic bottlenecks at the moment of transmission effectively suppress selection, allowing mutants to accumulate along the transmission chain to high levels. However, even very rare coinfections from independent sources are, as we show, able to counteract the bottleneck effect. Therefore, there are two possible explanations for the high mutant frequency. First, the frequency of coinfection in the natural host population may be quite low. Alternatively, a strong variation of the best-adapted sequence between individuals could be caused by a combination of an immune response present in early infection and coselection.

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 interaction of flavivirus nonstructural proteins NS1 and NS4A as a determinant of replicase function

Nonstructural protein 1 (NS1) of yellow fever virus (YF) is a glycoprotein localized to extracytoplasmic compartments within infected cells. We have previously shown that NS1 can be supplied in trans and is required for viral RNA replication, a process thought to occur in membrane-bound cytoplasmic complexes. Here we report that the NS1 gene from a related virus, dengue virus (DEN), is unable to function in the process of YF RNA replication. This virus-specific incompatibility leads to a lack of initial minus-strand accumulation, suggesting that DEN NS1 is unable to productively interact with the YF replicase. Based on a YF deletion mutant that requires NS1 in trans, a genetic screen for suppressor mutants was used to select virus variants able to utilize DEN NS1. In three independent selections, a single mutation was mapped to the NS4A gene, which encodes a putative transmembrane replicase component. This mutation, as well as several additional mutations, was engineered into the NS1-deficient genome and confirmed a genetic interaction between NS1 and NS4A. These findings suggest a potential mechanism for integrating NS1 into the cytoplasmic process of RNA replication.

Characterization of a natural mutation in an antigenic site on the fusion protein of measles virus that is involved in neutralization

Although measles virus is an antigenically monotypic virus, nucleotide sequence analysis of the hemagglutinin and nucleoprotein genes has permitted the differentiation of a number of genotypes. In contrast, the fusion (F) protein is highly conserved; only three amino acid changes have been reported over a 40-year period. We have isolated a measles virus strain which did not react with an anti-F monoclonal antibody (MAb) which we had previously shown to be directed against a dominant antigenic site. This virus strain, Lys-1, had seven amino acid changes compared with the Edmonston strain. We have shown that a single amino acid at position 73 is responsible for its nonreactivity with the anti-F MAb. With the same MAb, antibody-resistant mutants were prepared from the vaccine strain. A single amino acid change at position 73 (R-->W) was observed. The possibility of selecting measles virus variants in vaccinated populations is discussed.

Spontaneous mutation rate of measles virus: direct estimation based on mutations conferring monoclonal antibody resistance

High mutation rates typical of RNA viruses often generate a unique viral population structure consisting of a large number of genetic microvariants. In the case of viral pathogens, this can result in rapid evolution of antiviral resistance or vaccine-escape mutants. We determined a direct estimate of the mutation rate of measles virus, the next likely target for global elimination following poliovirus. In a laboratory tissue culture system, we used the fluctuation test method of estimating mutation rate, which involves screening a large number of independent populations initiated by a small number of viruses each for the presence or absence of a particular single point mutation. The mutation we focused on, which can be screened for phenotypically, confers resistance to a monoclonal antibody (MAb 80-III-B2). The entire H gene of a subset of mutants was sequenced to verify that the resistance phenotype was associated with single point mutations. The epitope conferring MAb resistance was further characterized by Western blot analysis. Based on this approach, measles virus was estimated to have a mutation rate of 9 x 10(-5) per base per replication and a genomic mutation rate of 1.43 per replication. The mutation rates we estimated for measles virus are comparable to recent in vitro estimates for both poliovirus and vesicular stomatitis virus. In the field, however, measles virus shows marked genetic stability. We briefly discuss the evolutionary implications of these results.

The origin and evolution of human T-cell lymphotropic virus types I and II

Studies on human T-cell lymphotropic virus types I (HTLV-I) and II (HTLV-II) are briefly reviewed from the viewpoint of molecular evolution, with special reference to the evolutionary rate and evolutionary relationships among these viruses. In particular, it appears that, in contrast to the low level of variability of HTLV-I among different isolates, individual isolates form quasispecies structures. Elucidating the mechanisms connecting these two phenomena will be one of the future problems in the study of the molecular evolution of HTLV-I and HTLV-II.

Generation of an infectious cDNA of a highly cardiovirulent coxsackievirus B3(CVB3m) and comparison to other infectious CVB3 cDNAs

An infectious cDNA of a highly myocarditic coxsackievirus B3 (CVB3m; Nancy strain) was cloned. Sequence data revealed 43 extra non-viral nucleotides upstream of the initial 5' sequence. However, the authentic 5' end sequence was maintained during replication of viral RNA transfected into HeLa cells, suggesting the RNA synthesizing complex edits the picornaviral 5' terminus sequence. Nucleotide sequences of the 5' nontranslated region and the capsid protein gene sequence of CVB3m were compared with the published sequences of five other CVB3 Nancy strains and two main lineages were found. In comparative assays for cardiovirulence, three of four CVB3 tested were cardiovirulent in adolescent male CD-1 mice. Only one of the three available CVB3 strains was neutralized with several anti-CVB3m monoclonal antibodies, suggesting that mutations in the surface epitopes of the capsid polypeptides contribute to antigenic drift within the serotype, perhaps in part through immunoselective pressures. Thus, phenotypic diversity of CVB3 within the prototype Nancy strain is an example of RNA viruses adapting to changing environments (cells, mice and humans) through mutations and selective pressure.

RNA virus evolution, population dynamics, and nutritional status

Trace elements exert a strong influence on immune function. Debilitated humoral and cellular immune responses may impair virus clearance in infected organisms, and favor the generation of virus variants with altered biological properties. The population size in evolving viral quasispecies, as well as increased mutagenesis triggered by oxidative stress, may contribute to altering the outcome of quasispecies evolution in infected hosts. The genetic plasticity of RNA viruses is one of the main obstacles for the control of the diseases they cause and probably a major force in the emergence of new viral pathogens. Recent results suggest links between nutritional deficiencies and the generation of variant viruses, a possibility that is addressed in the present article.

Genetic variation in open reading frame 2 of field isolates and laboratory strains of equine arteritis virus

The open reading frame 2 (ORF2) of three laboratory strains, the live attenuated vaccine virus, and 18 field isolates of equine arteritis virus (EAV) from Europe and North America was sequenced. The ORF2 of EAV encodes the Gs protein that is abundantly expressed in infected cells but constitutes less than 2% of the virion protein mass. Variation of ORF2 among the isolates facilitated phylogenetic analysis that largely confirmed results of an earlier study based on sequence divergence of ORF5 of the same isolates of EAV, despite exposure of the proteins encoded by ORF2 (Gs) and ORF5 (GL) to potentially different selective pressures in vivo. The data indicate that the Gs protein is highly conserved between isolates, considerably more so than the GL protein, consistent with an important role of the Gs protein in virus replication.

Adaptive evolution of human immunodeficiency virus-type 1 during the natural course of infection

The rate of progression to disease varies considerably among individuals infected with human immunodeficiency virus-type 1 (HIV-1). Analyses of semiannual blood samples obtained from six infected men showed that a rapid rate of CD4 T cell loss was associated with relative evolutionary stasis of the HIV-1 quasispecies virus population. More moderate rates of CD4 T cell loss correlated with genetic evolution within three of four subjects. Consistent with selection by the immune constraints of these subjects, amino acid changes were apparent within the appropriate epitopes of human leukocyte antigen class I-restricted cytotoxic T lymphocytes. Thus, the evolutionary dynamics exhibited by the HIV-1 quasispecies virus populations under natural selection are compatible with adaptive evolution.