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

Stochastic model of vesicular stomatitis virus replication reveals mutational effects on virion production

We present the first complete stochastic model of vesicular stomatitis virus (VSV) intracellular replication. Previous models developed to capture VSV's intracellular replication have either been ODE-based or have not represented the complete replicative cycle, limiting our ability to understand the impact of the stochastic nature of early cellular infections on virion production between cells and how these dynamics change in response to mutations. Our model accurately predicts changes in mean virion production in gene-shuffled VSV variants and can capture the distribution of the number of viruses produced. This model has allowed us to enhance our understanding of intercellular variability in virion production, which appears to be influenced by the duration of the early phase of infection, and variation between variants, arising from balancing the time the genome spends in the active state, the speed of incorporating new genomes into virions, and the production of viral components. Being a stochastic model, we can also assess other effects of mutations beyond just the mean number of virions produced, including the probability of aborted infections and the standard deviation of the number of virions produced. Our model provides a biologically interpretable framework for studying the stochastic nature of VSV replication, shedding light on the mechanisms underlying variation in virion production. In the future, this model could enable the design of more complex viral phenotypes when attenuating VSV, moving beyond solely considering the mean number of virions produced.

Cooperative Virus-Virus Interactions: An Evolutionary Perspective

Despite extensive evidence of virus-virus interactions, not much is known about their biological significance. Importantly, virus-virus interactions could have evolved as a form of cooperation or simply be a by-product of other processes. Here, we review and discuss different types of virus-virus interactions from the point of view of social evolution, which provides a well-established framework for interpreting the fitness costs and benefits of such traits. We also classify interactions according to their mechanisms of action and speculate on their evolutionary implications. As in any other biological system, the evolutionary stability of viral cooperation critically requires cheaters to be excluded from cooperative interactions. We discuss how cheater viruses exploit cooperative traits and how viral populations are able to counteract this maladaptive process.

Newcastle disease virus genotype VII gene expression in experimentally infected birds

Newcastle disease virus genotype VII (NDV-GVII) is a highly contagious pathogen responsible for pandemics that have caused devastating economic losses in the poultry industry. Several features in the transcription of NDV mRNA, including differentially expressed genes across the viral genome, are shared with that for other single, non-segmented, negative-strand viruses. Previous studies measuring viral gene expression using northern blotting indicated that the NDV transcription produced non-equimolar levels of viral mRNAs. However, deep high-throughput sequencing of virus-infected tissues can provide a better insight into the patterns of viral transcription. In this report, the transcription pattern of virulent NDV-GVII was analysed using RNA-seq and qRT-PCR. This study revealed the transcriptional profiling of these highly pathogenic NDV-GVII genes: NP:P:M:F:HN:L, in which there was a slight attenuation at the NP:P and HN:L gene boundaries. Our result also provides a fully comprehensive qPCR protocol for measuring viral transcript abundance that may be more convenient for laboratories where accessing RNA-seq is not feasible.

A mutant vesicular stomatitis virus with reduced cytotoxicity and enhanced anterograde trans-synaptic efficiency

Understanding the connecting structure of brain network is the basis to reveal the principle of the brain function and elucidate the mechanism of brain diseases. Trans-synaptic tracing with neurotropic viruses has become one of the most effective technologies to dissect the neural circuits. Although the retrograde trans-synaptic tracing for analyzing the input neural networks with recombinant rabies and pseudorabies virus has been broadly applied in neuroscience, viral tools for analyzing the output neural networks are still lacking. The recombinant vesicular stomatitis virus (VSV) has been used for the mapping of synaptic outputs. However, several drawbacks, including high neurotoxicity and rapid lethality in experimental animals, hinder its application in long-term studies of the structure and function of neural networks. To overcome these limitations, we generated a recombinant VSV with replication-related N gene mutation, VSV-N_R7A, and examined its cytotoxicity and efficiency of trans-synaptic spreading. We found that by comparison with the wild-type tracer of VSV, the N_R7A mutation endowed the virus lower rate of propagation and cytotoxicity in vitro, as well as significantly reduced neural inflammatory responses in vivo and much longer animal survival when it was injected into the nucleus of the mice brain. Besides, the spreading of the attenuated VSV was delayed when injected into the VTA. Importantly, with the reduced toxicity and extended animal survival, the number of brain regions that was trans-synaptically labeled by the mutant VSV was more than that of the wild-type VSV. These results indicated that the VSV-N_R7A, could be a promising anterograde tracer that enables researchers to explore more downstream connections of a given brain region, and observe the anatomical structure and the function of the downstream circuits over a longer time window. Our work could provide an improved tool for structural and functional studies of neurocircuit.

Phenotypic Consequence of Rearranging the N Gene of RABV HEP-Flury

Nucleoprotein (N) is a key element in rabies virus (RABV) replication. To further investigate the effect of N on RABV, we manipulated an infectious cDNA clone of the RABV HEP-Flury to rearrange the N gene from its wild-type position of 1 (N-P-M-G-L) to 2 (P-N-M-G-L), 3 (P-M-N-G-L), or 4 (P-M-G-N-L), using an approach that left the viral nucleotide sequence unaltered. Subsequently, viable viruses were recovered from each of the rearranged cDNA and examined for their gene expression levels, growth kinetics in cell culture, pathogenicity in suckling mice and protection in mice. The results showed that gene rearrangement decreased N mRNA transcription and vRNA replication. As a result, all viruses with rearranged genomes showed worse replication than that of rHEP-Flury in NA cells at a MOI of 0.01, but equivalent or slightly better replication levels at a MOI of 3. Consequently, the lethality in suckling mice infected with N4 was clearly attenuated compared with rHEP-Flury. However, the protection to mice was not enhanced. This study not only gives us insight into the understanding of the phenotype of RABV N gene rearrangement, but also helps with rabies vaccine candidate construction.

QSAR and molecular docking studies on 4-quinoline carboxylic acid derivatives as inhibition of vesicular stomatitis virus replication

The current study describes the development of in silico models based on quantitative structure-activity relationship (QSAR) analysis has been performed on 4-quinoline carboxylic acid derivatives as inhibition capacity of vesicular stomatitis virus replication in Madin Darby canine kidney epithelial cells. A highly descriptive and predictive QSAR model was obtained through the calculation of alignment-independent descriptors using MOE 2009.10 software. For a training set of 20 compounds, the partial least squares analyses result in a model which displays a squared correlation coefficient (r2) of 0.913. Validation of this model was performed using leave-one-out (q2) of 0.842. This model gives (r2pre) of 0.889 for a test set of five compounds. Docking studies were performed for 25 compounds to investigate the mode of interaction between 4-quinoline carboxylic acid derivatives and the active site of the human dihydroorotate dehydrogenase.

Monoclonal antibody resistant mutant of Peste des petits ruminants vaccine virus

The available vaccines for control of Peste des petits ruminants do not favour differentiation of infected and vaccinated animals (DIVA). Hence, the present study was aimed to isolate and characterize monoclonal antibody resistant mutant of an Indian strain of vaccine virus "PPRV-Sungri/96" under selection pressure of virus neutralizing monoclonal antibody '4B11' specific to haemagglutinin (H) protein. We successfully isolated five monoclonal antibody resistant (mAr) mutants (PPRV-RM5, PPRV-RM6, PPRV-RM7, PPRV- E6 and PPRV- E7). The mAr mutants did not react with the anti-H mAb 4B11 whereas reacted with control anti-nucleoprotein mAb 4G6, similar to the parent vaccine virus "PPRV-Sungri/96" in indirect ELISA, cell ELISA and indirect immunofluorescence test. Cytometry analysis of mAr mutants revealed loss of binding to mAb 4B11 while maintaining binding to mAb 4G6, more or less similar to "PPRV-Sungri/96". The sequence analysis of the H-protein gene of the mAr mutants resulted in identification of two nucleotide changes leading to amino acid substitutions at position 263 and 502 (L263P and R502P) of the H protein indicating that the epitope of mAb 4B11 could be conformational in nature. Though, mAr mutant grew to a similar titre as parent vaccine virus (PPRV-Sungri/96), the in vivo work in goats to study the mAr mutant as possible negative marker vaccine candidate could not be successfully proved with mAb 4B11 based competitive ELISA. However, one of the nucleotide change (T-C) at position 788, unique to mAr mutant virus resulted in abolition of a restriction enzyme recognition site (BglII). This could be used to differentiate mAr mutant vaccine virus from other available vaccine and field strains using restriction fragment length polymorphism. However, the mAr mutant PPRV-E6 cannot be used as a candidate strain for DIVA vaccine as immune response against it cannot be differentiated based on serology.

Multiple Barriers to the Evolution of Alternative Gene Orders in a Positive-Strand RNA Virus

The order in which genes are organized within a genome is generally not conserved between distantly related species. However, within virus orders and families, strong conservation of gene order is observed. The factors that constrain or promote gene-order diversity are largely unknown, although the regulation of gene expression is one important constraint for viruses. Here we investigate why gene order is conserved for a positive-strand RNA virus encoding a single polyprotein in the context of its authentic multicellular host. Initially, we identified the most plausible trajectory by which alternative gene orders could evolve. Subsequently, we studied the accessibility of key steps along this evolutionary trajectory by constructing two virus intermediates: (1) duplication of a gene followed by (2) loss of the ancestral gene. We identified five barriers to the evolution of alternative gene orders. First, the number of viable positions for reordering is limited. Second, the within-host fitness of viruses with gene duplications is low compared to the wild-type virus. Third, after duplication, the ancestral gene copy is always maintained and never the duplicated one. Fourth, viruses with an alternative gene order have even lower fitness than viruses with gene duplications. Fifth, after more than half a year of evolution in isolation, viruses with an alternative gene order are still vastly inferior to the wild-type virus. Our results show that all steps along plausible evolutionary trajectories to alternative gene orders are highly unlikely. Hence, the inaccessibility of these trajectories probably contributes to the conservation of gene order in present-day viruses.

Perfusion Pressure Is a Critical Determinant of the Intratumoral Extravasation of Oncolytic Viruses

Antitumor efficacy of oncolytic virotherapy is determined by the density and distribution of infectious centers within the tumor, which may be heavily influenced by the permeability and blood flow in tumor microvessels. Here, we investigated whether systemic perfusion pressure, a key driver of tumor blood flow, could influence the intratumoral extravasation of systemically administered oncolytic vesicular stomatitis virus (VSV) in myeloma tumor-bearing mice. Exercise was used to increase mean arterial pressure, and general anesthesia to decrease it. A recombinant VSV expressing the sodium iodide symporter (NIS), which concentrates radiotracers at sites of infection, was administered intravenously to exercising or anesthetized mice, and nuclear NIS reporter gene imaging was used to noninvasively track the density and spatial distribution of intratumoral infectious centers. Anesthesia resulted in decreased intratumoral infection density, while exercise increased the density and uniformity of infectious centers. Perfusion state also had a significant impact on the antitumor efficacy of the VSV therapy. In conclusion, quantitative dynamic radiohistologic imaging was used to noninvasively interrogate delivery of oncolytic virotherapy, highlighting the critical importance of perfusion pressure as a driver of intratumoral delivery and efficacy of oncolytic viruses.

Genome rearrangement affects RNA virus adaptability on prostate cancer cells

Gene order is often highly conserved within taxonomic groups, such that organisms with rearranged genomes tend to be less fit than wild type gene orders, and suggesting natural selection favors genome architectures that maximize fitness. But it is unclear whether rearranged genomes hinder adaptability: capacity to evolutionarily improve in a new environment. Negative-sense non-segmented RNA viruses (order Mononegavirales) have specific genome architecture: 3' UTR - core protein genes - envelope protein genes - RNA-dependent RNA-polymerase gene - 5' UTR. To test how genome architecture affects RNA virus evolution, we examined vesicular stomatitis virus (VSV) variants with the nucleocapsid (N) gene moved sequentially downstream in the genome. Because RNA polymerase stuttering in VSV replication causes greater mRNA production in upstream genes, N gene translocation toward the 5' end leads to stepwise decreases in N transcription, viral replication and progeny production, and also impacts the activation of type 1 interferon mediated antiviral responses. We evolved VSV gene-order variants in two prostate cancer cell lines: LNCap cells deficient in innate immune response to viral infection, and PC-3 cells that mount an IFN stimulated anti-viral response to infection. We observed that gene order affects phenotypic adaptability (reproductive growth; viral suppression of immune function), especially on PC-3 cells that strongly select against virus infection. Overall, populations derived from the least-fit ancestor (most-altered N position architecture) adapted fastest, consistent with theory predicting populations with low initial fitness should improve faster in evolutionary time. Also, we observed correlated responses to selection, where viruses improved across both hosts, rather than suffer fitness trade-offs on unselected hosts. Whole genomics revealed multiple mutations in evolved variants, some of which were conserved across selective environments for a given gene order.

Evolutionary reversion of live viral vaccines: Can genetic engineering subdue it?

Attenuated, live viral vaccines have been extraordinarily successful in protecting against many diseases. The main drawbacks in their development and use have been reliance on an unpredictable method of attenuation and the potential for evolutionary reversion to high virulence. Methods of genetic engineering now provide many safer alternatives to live vaccines, so if live vaccines are to compete with these alternatives in the future, they must either have superior immunogenicity or they must be able to overcome these former disadvantages. Several live vaccine designs that were historically inaccessible are now feasible because of advances in genome synthesis. Some of those methods are addressed here, with an emphasis on whether they enable predictable levels of attenuation and whether they are stable against evolutionary reversion. These new designs overcome many of the former drawbacks and position live vaccines to be competitive with alternatives. Not only do new methods appear to retard evolutionary reversion enough to prevent vaccine-derived epidemics, but it may even be possible to permanently attenuate live vaccines that are transmissible but cannot evolve to higher virulence under prolonged adaptation.

Increased replicative fitness of a dengue virus 2 clade in native mosquitoes: potential contribution to a clade replacement event in Nicaragua

The four dengue virus (DENV) serotypes (DENV serotype 1 [DENV-1] to DENV-4) are transmitted by Aedes aegypti and A. albopictus mosquitoes, causing up to 390 million DENV infections worldwide each year. We previously reported a clade replacement of the DENV-2 Asian-American genotype NI-1 clade by the NI-2B clade in Managua, Nicaragua. Here, we describe our studies of the replicative ability of NI-1 and NI-2B viruses in an A. aegypti cell line (Aag2) and A. aegypti mosquitoes reared from eggs collected in Managua. In coinfection experiments, several different pairs of NI-1 and NI-2B clinical isolates were used to infect Aag2 cells or blood-fed A. aegypti mosquitoes. Results consistently showed a significant replicative advantage of NI-2B over NI-1 viruses early after infection in vitro, and in mosquitoes, NI-2B viruses attained a higher replicative index than NI-1 isolates 3 to 7 days postinfection (dpi). At 7 dpi, NI-2B viruses displayed a significantly higher replicative index in legs and salivary glands; however, this advantage was lost by 14 and 21 dpi. We also found that the percentage of mosquitoes in which NI-2B viruses were dominant was significantly higher than that in which NI-1 viruses were dominant on day 7 but not at later time points. Taken together, these data demonstrate that clade NI-2B holds a replicative advantage over clade NI-1 early in infection that wanes at later time points. This early fitness advantage of NI-2B viruses over NI-1 viruses in the native vector, A. aegypti, suggests a shorter extrinsic incubation period for NI-2B viruses, which could have contributed to the clade replacement event in Managua.

IMPORTANCE

Dengue virus (DENV), one of the most medically important arthropod-borne viruses, is transmitted to humans by Aedes aegypti and A. albopictus mosquitoes in tropical and subtropical regions worldwide. Dengue epidemics continue to increase in frequency, geographic range, and severity and are a major public health concern. This is due to globalization, unplanned urbanization, and climate change, as well as host genetics and immune responses and viral genetic changes. DENV consists of four serotypes, in turn composed of genotypes and genetically distinct clades. What drives the frequent replacement of a previously circulating DENV clade by another is unclear. Here, we investigate the replicative fitness of two clades of DENV serotype 2 in Aedes aegypti cells and mosquitoes collected from the region where the viruses circulated and conclude that increased replicative fitness could have contributed to a DENV clade replacement event in Nicaragua. These findings provide insight into vector-driven evolution of DENV epidemics.

Consequences of in vitro host shift for St. Louis encephalitis virus

Understanding the potential for host range shifts and expansions of RNA viruses is critical to predicting the evolutionary and epidemiological paths of these pathogens. As arthropod-borne viruses (arboviruses) experience frequent spillover from their amplification cycles and are generalists by nature, they are likely to experience a relatively high frequency of success in a range of host environments. Despite this, the potential for host expansion, the genetic correlates of adaptation to novel environments and the costs of such adaptations in originally competent hosts are still not characterized fully for arboviruses. In the studies presented here, we utilized experimental evolution of St. Louis encephalitis virus (SLEV; family Flaviviridae, genus Flavivirus) in vitro in the Dermacentor andersoni line of tick cells to model adaptation to a novel invertebrate host. Our results demonstrated that levels of adaptation and costs in alternate hosts are highly variable among lineages, but also that significant fitness increases in tick cells are achievable with only modest change in consensus genetic sequence. In addition, although accumulation of diversity may at times buffer against phenotypic costs within the SLEV swarm, an increased proportion of variants with an impaired capacity to infect and spread on vertebrate cell culture accumulated with tick cell passage. Isolation and characterization of a subset of these variants implicates the NS3 gene as an important host range determinant for SLEV.

Evolutionarily stable attenuation by genome rearrangement in a virus

Live, attenuated viruses provide many of the most effective vaccines. For the better part of a century, the standard method of attenuation has been viral growth in novel environments, whereby the virus adapts to the new environment but incurs a reduced ability to grow in the original host. The downsides of this approach were that it produced haphazard results, and even when it achieved sufficient attenuation for vaccine production, the attenuated virus was prone to evolve back to high virulence. Using bacteriophage T7, we apply a synthetic biology approach for creating attenuated genomes and specifically study their evolutionary stability. Three different genome rearrangements are used, and although some initial fitness recovery occurs, all exhibit greatly impaired abilities to recover wild-type fitness over a hundred or more generations. Different degrees of stable attenuation appear to be attainable by different rearrangements. Efforts to predict fitness recovery using the extensive background of T7 genetics and biochemistry were only sometimes successful. The use of genome rearrangement thus offers a practical mechanism of evolutionary stable viral attenuation, with some progress toward prediction.

Understanding and altering cell tropism of vesicular stomatitis virus

Vesicular stomatitis virus (VSV) is a prototypic nonsegmented negative-strand RNA virus. VSV's broad cell tropism makes it a popular model virus for many basic research applications. In addition, a lack of preexisting human immunity against VSV, inherent oncotropism and other features make VSV a widely used platform for vaccine and oncolytic vectors. However, VSV's neurotropism that can result in viral encephalitis in experimental animals needs to be addressed for the use of the virus as a safe vector. Therefore, it is very important to understand the determinants of VSV tropism and develop strategies to alter it. VSV glycoprotein (G) and matrix (M) protein play major roles in its cell tropism. VSV G protein is responsible for VSV broad cell tropism and is often used for pseudotyping other viruses. VSV M affects cell tropism via evasion of antiviral responses, and M mutants can be used to limit cell tropism to cell types defective in interferon signaling. In addition, other VSV proteins and host proteins may function as determinants of VSV cell tropism. Various approaches have been successfully used to alter VSV tropism to benefit basic research and clinically relevant applications.

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.

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

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

Level of gene expression is a major determinant of protein evolution in the viral order Mononegavirales

Although the rate at which proteins change is a key parameter in molecular evolution, its determinants are poorly understood in viruses. A variety of factors, including gene length, codon usage bias, protein abundance, protein function, and gene expression level, have been shown to affect the rate of protein evolution in a diverse array of organisms. However, the role of these factors in viral evolution has yet to be addressed. The polar 3'-5' stepwise attenuation of transcription in the Mononegavirales, a group of single-strand negative-sense RNA viruses, provides a unique system to explore the determinants of protein evolution in viruses. We analyzed the relative importance of a variety of factors in shaping patterns of sequence variation in full-length genomes from 13 Mononegavirales species. Our analysis suggests that the level of gene expression, and by extension the relative genomic position of each gene, is a key determinant of the protein evolution in these viruses. This appears to be the consequence of selection for translational robustness, but not for translational accuracy, in highly expressed genes. The small genome size and number of proteins encoded by these viruses allowed us to identify other protein-specific factors that may also play a role in virus evolution, such as host-virus interactions and functional constraints. Finally, we explored the evolutionary pressures acting on noncoding regions in Mononegavirales genomes and observed that, despite being less constrained than coding regions, their evolutionary rates are also associated with genomic position.

Why do RNA viruses recombine?

Recombination occurs in many RNA viruses and can be of major evolutionary significance. However, rates of recombination vary dramatically among RNA viruses, which can range from clonal to highly recombinogenic. Here, we review the factors that might explain this variation in recombination frequency and show that there is little evidence that recombination is favoured by natural selection to create advantageous genotypes or purge deleterious mutations, as predicted if recombination functions as a form of sexual reproduction. Rather, recombination rates seemingly reflect larger-scale patterns of viral genome organization, such that recombination may be a mechanistic by-product of the evolutionary pressures acting on other aspects of virus biology.

Analysis of respiratory syncytial virus preclinical and clinical variants resistant to neutralization by monoclonal antibodies palivizumab and/or motavizumab

BACKGROUND

Palivizumab is a US Food and Drug Administration-approved monoclonal antibody for the prevention of respiratory syncytial virus (RSV) lower respiratory disease in high-risk infants. Motavizumab, derived from palivizumab with enhanced antiviral activity, has recently been tested in humans. Although palivizumab escape mutants have been generated in the laboratory, the development of resistant RSV in patients receiving palivizumab has not been reported previously.

METHODS

We generated palivizumab and motavizumab escape mutants in vitro and examined the development of resistant mutants in RSV-breakthrough patients receiving immunoprophylaxis. The effect of these mutations on neutralization by palivizumab and motavizumab and in vitro fitness was studied.

RESULTS

Antibody-resistant RSV variants selected in vitro had mutations at position 272 of the fusion protein, from lysine to asparagine, methionine, threonine, glutamine, or glutamate. Variants containing mutations at positions 272 and 275 were detected in breakthrough patients. All these variants were resistant to palivizumab, but only the glutamate variant at position 272 demonstrated resistance to motavizumab. Mixtures of wild-type and variant RSV soon lost the resistant phenotype in the absence of selection.

CONCLUSIONS

Resistant RSV variants were detected in a small subset (∼ 5%) of RSV breakthrough cases. The fitness of these variants was impaired, compared to wild-type RSV.

Analysis of parainfluenza virus-5 hemagglutinin-neuraminidase protein mutants that are blocked in internalization and degradation

The PIV-5 hemagglutinin-neuraminidase (HN) protein is a multifunctional protein with sialic acid binding, neuraminidase and fusion promotion activity. HN is internalized by clathrin-mediated endocytosis and degraded. HN lacks internalization signals in its cytoplasmic tail but a single glutamic acid present at residue 37 at the putative transmembrane/ectodomain boundary is critical. We rescued rPIV-5 with mutations E37D or E37K, which have been shown to impair or abolish HN internalization, respectively. These viruses exhibited growth properties similar to wild-type (wt) virus but are impaired for fitness in tissue culture. Biochemical analysis of HN activities showed differences between HN E37D and HN E37K in fusion promotion and incorporation of HN and F into virions. Furthermore, oligomeric analyses indicate that HN E37 mutants perturb the tetrameric organization of HN, probably by destabilizing the dimer-of-dimers interface.

Safety studies on intrahepatic or intratumoral injection of oncolytic vesicular stomatitis virus expressing interferon-beta in rodents and nonhuman primates

Toxicology studies were performed in rats and rhesus macaques to establish a safe starting dose for intratumoral injection of an oncolytic vesicular stomatitis virus expressing human interferon-beta (VSV-hIFNbeta) in patients with hepatocellular carcinoma (HCC). No adverse events were observed after administration of 7.59 x 10(9) TCID(50) (50% tissue culture infective dose) of VSV-hIFNbeta into the left lateral hepatic lobe of Harlan Sprague Dawley rats. Plasma alanine aminotransferase and alkaline phosphatase levels increased and platelet counts decreased in the virus-treated animals on days 1 and 2 but returned to pretreatment levels by day 4. VSV-hIFNbeta was also injected into normal livers or an intrahepatic McA-RH7777 HCC xenograft established in Buffalo rats. Buffalo rats were more sensitive to neurotoxic effects of VSV; the no observable adverse event level (NOAEL) of VSV-hIFNbeta in Buffalo rats was 10(7) TCID(50). Higher doses were associated with fatal neurotoxicity and infectious virus was recovered from tumor and brain. Compared with VSV-hIFNbeta, toxicity of VSV-rIFNbeta (recombinant VSV expressing rat IFN-beta) was greatly diminished in Buffalo rats (NOAEL, >10(10) TCID(50)). Two groups of two adult male rhesus macaques received 10(9) or 10(10) TCID(50) of VSV-hIFNbeta injected directly into the left hepatic lobe under computed tomographic guidance. No neurological signs were observed at any time point. No abnormalities (hematology, clinical chemistry, body weights, behavior) were seen and all macaques developed neutralizing anti-VSV antibodies. Plasma interleukin-6, tumor necrosis factor-alpha, and hIFN-beta remained below detection levels by ELISA. On the basis of these studies, we will be proposing a cautious approach to dose escalation in a phase I clinical trial among patients with HCC.

Infidelity of SARS-CoV Nsp14-exonuclease mutant virus replication is revealed by complete genome sequencing

Most RNA viruses lack the mechanisms to recognize and correct mutations that arise during genome replication, resulting in quasispecies diversity that is required for pathogenesis and adaptation. However, it is not known how viruses encoding large viral RNA genomes such as the Coronaviridae (26 to 32 kb) balance the requirements for genome stability and quasispecies diversity. Further, the limits of replication infidelity during replication of large RNA genomes and how decreased fidelity impacts virus fitness over time are not known. Our previous work demonstrated that genetic inactivation of the coronavirus exoribonuclease (ExoN) in nonstructural protein 14 (nsp14) of murine hepatitis virus results in a 15-fold decrease in replication fidelity. However, it is not known whether nsp14-ExoN is required for replication fidelity of all coronaviruses, nor the impact of decreased fidelity on genome diversity and fitness during replication and passage. We report here the engineering and recovery of nsp14-ExoN mutant viruses of severe acute respiratory syndrome coronavirus (SARS-CoV) that have stable growth defects and demonstrate a 21-fold increase in mutation frequency during replication in culture. Analysis of complete genome sequences from SARS-ExoN mutant viral clones revealed unique mutation sets in every genome examined from the same round of replication and a total of 100 unique mutations across the genome. Using novel bioinformatic tools and deep sequencing across the full-length genome following 10 population passages in vitro, we demonstrate retention of ExoN mutations and continued increased diversity and mutational load compared to wild-type SARS-CoV. The results define a novel genetic and bioinformatics model for introduction and identification of multi-allelic mutations in replication competent viruses that will be powerful tools for testing the effects of decreased fidelity and increased quasispecies diversity on viral replication, pathogenesis, and evolution.

A temperature sensitive VSV identifies L protein residues that affect transcription but not replication

To investigate the polymerase components selectively involved in transcription versus replication of vesicular stomatitis virus (VSV), we sequenced the polymerase gene of a conditionally RNA defective, temperature sensitive VSV: ts(G)114, which has a phenotype upon shift from permissive to non-permissive temperature of shut-down of mRNA transcription and unaffected genome replication. Sequence analysis of the ts(G)114 L gene identified three altered amino acid residues in the L protein. These three changes were specifically engineered individually and in combinations into a functional cDNA clone encoding the VSV genome and tested for association with the temperature sensitive and RNA defective phenotypes in the background of recovered engineered viruses. The data presented in this study show a specific amino acid substitution in domain II of the VSV L protein that significantly affects total RNA synthesis, but when in combination with two additional amino acid substitutions identified in the ts(G)114 L protein, leads to a specific reduction in mRNA transcription, but not replication.

In vitro and in vivo fitness of respiratory syncytial virus monoclonal antibody escape mutants

Respiratory syncytial virus (RSV) is the only infectious disease for which a monoclonal antibody (MAb) is used in humans. Palivizumab (PZ) is a humanized murine MAb to the F protein of RSV. PZ-resistant viruses appear after in vitro and in vivo growth of RSV in the presence of PZ. Fitness for replication could be a determinant of the likelihood of dissemination of resistant viruses. We assessed the fitness of two PZ-resistant viruses (F212 and MP4). F212 grew less well in cell culture than the parent A2 virus and was predicted to be less fit than A2. Equal amounts of F212 and A2 were mixed and passaged in cell culture. F212 disappeared from the viral population, indicating it was less fit than the A2 virus. The MP4 virus grew as well as A2 in culture and in cotton rats. A2/MP4 virus input ratios of 1:1, 10:1, 100:1, and 1,000:1 were compared in competitive replication. For all input ratios except 1,000:1, the MP4 virus became dominant, supplanting the A2 virus. The MP4 virus also dominated the A2 virus during growth in cotton rats. Thus, the mutant MP4 virus was more fit than A2 virus in both in vitro and in vivo competitive replication. Whether this fitness difference was due to the identified nucleotide substitutions in the F gene or to mutations elsewhere in the genome is unknown. Understanding the mechanisms by which mutant virus fitness increased or decreased could prove useful for consideration in attenuated vaccine design efforts.

Neurovirulence properties of recombinant vesicular stomatitis virus vectors in non-human primates

Although vesicular stomatitis virus (VSV) neurovirulence and pathogenicity in rodents have been well studied, little is known about VSV pathogenicity in non-human primates. To address this question, we measured VSV viremia, shedding, and neurovirulence in macaques. Following intranasal inoculation, macaques shed minimal recombinant VSV (rVSV) in nasal washes for 1 day post-inoculation; viremia was not detected. Following intranasal inoculation of macaques, wild type (wt) VSV, rVSV, and two rVSV-HIV vectors showed no evidence of spread to CNS tissues. However, macaques inoculated intrathalamically with wt VSV developed severe neurological disease. One of four macaques receiving rVSV developed clinical and histological signs similar to the wt group, while the remaining three macaques in this group and all of the macaques in the rVSV-HIV vector groups showed no clinical signs of disease and reduced severity of histopathology compared to the wt group. The implications of these findings for rVSV vaccine development are discussed.

Virus fitness: concept, quantification, and application to HIV population dynamics

Viral fitness has been broadly studied during the past three decades, mainly to test evolutionary models and population theories difficult to analyze and interpret with more complex organisms. More recent studies, however, are focused in the role of fitness on viral transmission, pathogenesis, and drug resistance. Here, we used human immunodeficiency virus (HIV) as one of the most relevant models to evaluate the importance of viral quasispecies and fitness in HIV evolution, population dynamics, disease progression, and potential clinical implications.

Rational attenuation of a morbillivirus by modulating the activity of the RNA-dependent RNA polymerase

Negative-strand RNA viruses encode a single RNA-dependent RNA polymerase (RdRp) which transcribes and replicates the genome. The open reading frame encoding the RdRp from a virulent wild-type strain of rinderpest virus (RPV) was inserted into an expression plasmid. Sequences encoding enhanced green fluorescent protein (EGFP) were inserted into a variable hinge of the RdRp. The resulting polymerase was autofluorescent, and its activity in the replication/transcription of a synthetic minigenome was reduced. We investigated the potential of using this approach to rationally attenuate a virus by inserting the DNA sequences encoding the modified RdRp into a full-length anti-genome plasmid from which a virulent virus (rRPV(KO)) can be rescued. A recombinant virus, rRPV(KO)L-RRegfpR, which grew at an indistinguishable rate and to an identical titer as rRPV(KO) in vitro, was rescued. Fluorescently tagged polymerase was visible in large cytoplasmic inclusions and beneath the cell membrane. Subcutaneous injection of 10(4) TCID(50) of the rRPV(KO) parental recombinant virus into cattle leads to severe disease symptoms (leukopenia/diarrhea and pyrexia) and death by 9 days postinfection. Animals infected with rRPV(KO)L-RRegfpR exhibited transient leukopenia and mild pyrexia, and the only noticeable clinical signs were moderate reddening of one eye and a slight ocular-nasal discharge. Viruses that expressed the modified polymerase were isolated from peripheral blood lymphocytes and eye swabs. This demonstrates that a virulent morbillivirus can be attenuated in a single step solely by modulating RdRp activity and that there is not necessarily a correlation between virus growth in vitro and in vivo.

Gene order constrains adaptation in bacteriophage T7

The order of genes in the genome is commonly thought to have functional significance for gene regulation and fitness but has not heretofore been tested experimentally. We adapted a bacteriophage T7 variant harboring an ectopically positioned RNA polymerase gene to determine whether it could regain the fitness of the wild type. Two replicate lines maintained the starting gene order and showed only modest recovery of fitness, despite the accumulation of over a dozen mutations. In both lines, a mutation in the early terminator signal is responsible for the majority of the fitness recovery. In a third line, the phage evolved a new gene order, restoring the wild-type position of the RNA polymerase gene but also displacing several other genes to ectopic locations. Due to the recombination, the fitness of this replicate was the highest obtained but it falls short of the wild type adapted to the same growth conditions. The large benefits afforded by the terminator mutation and the recombination are explicable in terms of T7 biology, whereas several mutations with lesser benefits are not easily accounted for. These results support the premise that gene order is important to fitness and that wild-type fitness is not rapidly re-evolved in reorganized genomes.

A role for virus promoters in determining the pathogenesis of Rinderpest virus in cattle

Rinderpest virus (RPV) is a morbillivirus that causes cattle plague, a disease of large ruminants. The viral genome is flanked at the 3′ and 5′ genome termini by the genome promoter (GP) and antigenome promoter (AGP), respectively. These promoters play essential roles in directing replication and transcription as well as RNA encapsidation and packaging. It has previously been shown that individual changes to the GP of RPV greatly affect promoter activity in a minigenome assay and it was therefore proposed that individual nucleotide changes in the GP and AGP might also have significant effects on the ability of the virus to replicate and cause disease in cattle. The Plowright vaccine strain of RPV has been derived by tissue-culture passage from the virulent Kabete ‘O’ isolate (KO) and is highly attenuated for all ruminant species in which it has been used. Here, it was shown that swapping the GP and the first 76 nt of the AGP between virulent and avirulent strains affected disease progression. In particular, it was shown that flanking the virulent strain with the vaccine GP and AGP sequences, while not appreciably affecting virus growth in vitro, led to attenuation in vivo. The reverse was not true, since the KO promoters did not alter the vaccine's attenuated nature. The GP/AGP therefore play a role in attenuation, but are not the only determinants of attenuation in this vaccine.

The Plowright vaccine strain of Rinderpest virus has attenuating mutations in most genes

The currently used vaccine strain of Rinderpest virus was derived by serial passage of the highly virulent Kabete ‘O’ strain (KO). A full-length cDNA copy of the KO strain was made from which a virus identical in pathogenicity to the wild-type virus was rescued. A series of chimeric viruses was prepared in which the coding sequences for the N, P, F, H or L proteins were replaced with the corresponding sequences from the vaccine strain. The KO-based virus with the vaccine strain H gene and that with the carboxy-terminal half of the L gene replaced with the corresponding sequence from the vaccine strain retained all or almost all of the virulence of the original KO virus. Animals infected with the KO-based virus containing the vaccine strain N, P or F gene, or the amino-terminal half of the L gene, developed high and prolonged pyrexia and leukopenia, but with reduced or absent lesions and other clinical signs; although partially attenuated, none was nearly as attenuated as the vaccine strain itself. These data indicate that the high attenuation and stability of the current vaccine are due to the accumulation of a number of separate mutations, none of which is itself so sufficiently debilitating that there is strong selective pressure in favour of the revertant.