Ribavirin-Resistant Mutants of Human Enterovirus 71 Express a High Replication Fidelity Phenotype during Growth in Cell Culture

It has been shown in animal models that ribavirin-resistant poliovirus with a G64S mutation in its 3D polymerase has high replication fidelity coupled with attenuated virulence. Here, we describe the effects of mutagenesis in the human enterovirus 71 (HEV71) 3D polymerase on ribavirin resistance and replication fidelity. Seven substitutions were introduced at amino acid position 3D-G64 of a HEV71 full-length infectious cDNA clone (26M). Viable clone-derived virus populations were rescued from the G64N, G64R, and G64T mutant cDNA clones. The clone-derived G64R and G64T mutant virus populations were resistant to growth inhibition in the presence of 1,600 μM ribavirin, whereas the growth of parental 26M and the G64N mutant viruses were inhibited in the presence of 800 μM ribavirin. Nucleotide sequencing of the 2C and 3D coding regions revealed that the rate of random mutagenesis after 13 passages in the presence of 400 μM ribavirin was nearly 10 times higher in the 26M genome than in the mutant G64R virus genome. Furthermore, random mutations acquired in the 2C coding regions of 26M and G64N conferred resistance to growth inhibition in the presence of 0.5 mM guanidine, whereas the G64R and G64T mutant virus populations remained susceptible to growth inhibition by 0.5 mM guanidine. Interestingly, a S264L mutation identified in the 3D coding region of 26M after ribavirin selection was also associated with both ribavirin-resistant and high replication fidelity phenotypes. These findings are consistent with the hypothesis that the 3D-G64R, 3D-G64T, and 3D-S264L mutations confer resistance upon HEV71 to the antiviral mutagen ribavirin, coupled with a high replication fidelity phenotype during growth in cell culture.     

High fidelity of murine hepatitis virus replication is decreased in nsp14 exoribonuclease mutants

Replication fidelity of RNA virus genomes is constrained by the opposing necessities of generating sufficient diversity for adaptation and maintaining genetic stability, but it is unclear how the largest viral RNA genomes have evolved and are maintained under these constraints. A coronavirus (CoV) nonstructural protein, nsp14, contains conserved active-site motifs of cellular exonucleases, including DNA proofreading enzymes, and the severe acute respiratory syndrome CoV (SARS-CoV) nsp14 has 3'-to-5' exoribonuclease (ExoN) activity in vitro. Here, we show that nsp14 ExoN remarkably increases replication fidelity of the CoV murine hepatitis virus (MHV). Replacement of conserved MHV ExoN active-site residues with alanines resulted in viable mutant viruses with growth and RNA synthesis defects that during passage accumulated 15-fold more mutations than wild-type virus without changes in growth fitness. The estimated mutation rate for ExoN mutants was similar to that reported for other RNA viruses, whereas that of wild-type MHV was less than the established rates for RNA viruses in general, suggesting that CoVs with intact ExoN replicate with unusually high fidelity. Our results indicate that nsp14 ExoN plays a critical role in prevention or repair of nucleotide incorporation errors during genome replication. The established mutants are unique tools to test the hypothesis that high replication fidelity is required for the evolution and stability of large RNA genomes.     

Evolution of drug-resistant viral populations during interruption of antiretroviral therapy

Analysis of a large number of HIV-1 genomes at multiple time points after antiretroviral treatment (ART) interruption allows determination of the evolution of drug-resistant viruses and viral fitness in vivo in the absence of drug selection pressure. Using a parallel allele-specific sequencing (PASS) assay, potential primary drug-resistant mutations in five individual patients were studied by analyzing over 18,000 viral genomes. A three-phase evolution of drug-resistant viruses was observed after termination of ART. In the first phase, viruses carrying various combinations of multiple-drug-resistant (MDR) mutations predominated with each mutation persisting in relatively stable proportions while the overall number of resistant viruses gradually increased. In the second phase, viruses with linked MDR mutations rapidly became undetectable and single-drug-resistant (SDR) viruses emerged as minority populations while wild-type viruses quickly predominated. In the third phase, low-frequency SDR viruses remained detectable as long as 59 weeks after treatment interruption. Mathematical modeling showed that the loss in relative fitness increased with the number of mutations in each viral genome and that viruses with MDR mutations had lower fitness than viruses with SDR mutations. No single viral genome had seven or more drug resistance mutations, suggesting that such severely mutated viruses were too unfit to be detected or that the resistance gain offered by the seventh mutation did not outweigh its contribution to the overall fitness loss of the virus. These data provide a more comprehensive understanding of evolution and fitness of drug-resistant viruses in vivo and may lead to improved treatment strategies for ART-experienced patients.     

Quantitative SNP-detection method for estimating HIV-1 replicative fitness: application to protease inhibitor-resistant viruses

We have improved the methods for the standard competitive growth assay of human immunodeficiency virus type 1 (HIV-1). The cloning step for the mixed viral population and subsequent genotype analysis for arbitrary numbers of clones were excluded from procedures. Instead, a single nucleotide polymorphism (SNP)-detection step was devised for the determination of viral populations. The quantitative SNP-detection method can rapidly estimate the proportion of wild-type and mutant populations with high reproducibility. Consequently, this method allows manipulation of many samples within a short period. Using this new competitive growth assay, replicative fitness of drug-resistant HIV-1 containing an M46I amino acid mutation in the protease was assessed in the presence or absence of indinavir. Without indinavir, replicative fitness of wild-type HIV-1 surpassed that of M46I-mutated HIV-1, and the fraction of mutated virus was reduced to about 10% at passage #9. In contrast, the fraction of M46I-mutated virus increased to >90% at passage #5 in the presence of 26.4 nM indinavir. Almost identical results were obtained for L90M-mutated HIV-1 with or without saquinavir. HIV-1 can survive under indinavir pressure by acquiring M46I mutation, as with acquisition of the L90M mutation under saquinavir pressure. However, these mutations damage viral replicative fitness under natural conditions without any drugs. Subtle differences between wild-type and mutant viruses are thus easily detected using the improved method.     

A Multi-Step Process of Viral Adaptation to a Mutagenic Nucleoside Analogue by Modulation of Transition Types Leads to Extinction-Escape

Resistance of viruses to mutagenic agents is an important problem for the development of lethal mutagenesis as an antiviral strategy. Previous studies with RNA viruses have documented that resistance to the mutagenic nucleoside analogue ribavirin (1-β-D-ribofuranosyl-1-H-1,2,4-triazole-3-carboxamide) is mediated by amino acid substitutions in the viral polymerase that either increase the general template copying fidelity of the enzyme or decrease the incorporation of ribavirin into RNA. Here we describe experiments that show that replication of the important picornavirus pathogen foot-and-mouth disease virus (FMDV) in the presence of increasing concentrations of ribavirin results in the sequential incorporation of three amino acid substitutions (M296I, P44S and P169S) in the viral polymerase (3D). The main biological effect of these substitutions is to attenuate the consequences of the mutagenic activity of ribavirin —by avoiding the biased repertoire of transition mutations produced by this purine analogue—and to maintain the replicative fitness of the virus which is able to escape extinction by ribavirin. This is achieved through alteration of the pairing behavior of ribavirin-triphosphate (RTP), as evidenced by in vitro polymerization assays with purified mutant 3Ds. Comparison of the three-dimensional structure of wild type and mutant polymerases suggests that the amino acid substitutions alter the position of the template RNA in the entry channel of the enzyme, thereby affecting nucleotide recognition. The results provide evidence of a new mechanism of resistance to a mutagenic nucleoside analogue which allows the virus to maintain a balance among mutation types introduced into progeny genomes during replication under strong mutagenic pressure.     

Role of the human immunodeficiency virus type 1 envelope gene in viral fitness

A human host offers a variety of microenvironments to the infecting human immunodeficiency virus type 1 (HIV-1), resulting in various selective pressures, most of them directed against the envelope (env) gene. Therefore, it seems evident that the replicative capacity of the virus is largely related to viral entry. In this study we have used growth competition experiments and TaqMan real-time PCR detection to measure the fitness of subtype B HIV-1 primary isolates and autologous env-recombinant viruses in order to analyze the contribution of wild-type env sequences to overall HIV-1 fitness. A significant correlation was observed between fitness values obtained for wild-type HIV-1 isolates and those for the corresponding env-recombinant viruses (r = 0.93; P = 0.002). Our results suggest that the env gene, which is linked to a myriad of viral characteristics (e.g., entry into the host cell, transmission, coreceptor usage, and tropism), plays a major role in fitness of wild-type HIV-1. In addition, this new recombinant assay may be useful for measuring the contribution of HIV-1 env to fitness in viruses resistant to novel antiretroviral entry inhibitors.     

A mutation present in the amino terminus of Sabin 3 poliovirus VP1 protein is attenuating.

The attenuated phenotype of Sabin 3 poliovirus compared with its neurovirulent progenitor strain has been largely accounted for by mutations in the genome at positions 472 and 2034 (G. D. Westrop, K. A. Wareham, D. M. A. Evans, G. Dunn, P. D. Minor, D. I. Magrath, F. Taffs, S. Marsden, M. A. Skinner, G. C. Schild, and J. W. Almond, J. Virol. 63:1338-1344, 1989). By sequencing vaccine virus RNA, we recently identified another Sabin 3-specific mutation at position 2493 (U----C), which predicts an Ile----Thr change at the sixth residue of VP1 (C. Weeks-Levy, J. M. Tatem, S. J. DiMichele, W. Waterfield, A. F. Georgiu, and S. J. Mento, Virology 185:934-937, 1991). Viruses generated by using cDNAs which represent the vaccine sequence (LED3) and a derivative (VR318) possessing a single base change to the wild-type nucleotide (U) at 2493 were used to determine the impact of the 2493 mutation on virus phenotype. The VP1 proteins of LED3 and VR318 viruses were distinguishable by denaturing electrophoretic analysis. LED3 produced smaller plaques in Vero cells than VR318 virus did. Neurovirulence testing of these cDNA-derived viruses in monkeys demonstrated that the 2493 mutation in LED3 virus is attenuating.     

Adaptation of reovirus to growth in the presence of protease inhibitor E64 segregates with a mutation in the carboxy terminus of viral outer-capsid protein sigma3

Reovirus virions are internalized into cells by receptor-mediated endocytosis. Within the endocytic compartment, the viral outer capsid undergoes acid-dependent proteolysis leading to degradation of sigma3 protein and proteolytic cleavage of micro1/micro1C protein. E64 is a specific inhibitor of cysteine-containing proteases that blocks disassembly of reovirus virions. To identify domains in reovirus proteins that influence susceptibility to E64-mediated inhibition of disassembly, we selected variant viruses by serial passage of strain type 3 Dearing (T3D) in murine L929 cells treated with E64. E64-adapted variant viruses (D-EA viruses) produced 7- to 17-fold-greater yields than T3D did after infection of cells treated with 100 microM E64. Viral genes that segregate with growth of D-EA viruses in the presence of E64 were identified by using reassortant viruses isolated from independent crosses of E64-sensitive strain type 1 Lang and two prototype D-EA viruses. Growth of reassortant viruses in the presence of E64 segregated with the S4 gene, which encodes outer-capsid protein sigma3. Sequence analysis of S4 genes of three D-EA viruses isolated from independent passage series revealed a common tyrosine-to-histidine mutation at amino acid 354 in the deduced amino acid sequence of sigma3. Proteolysis of D-EA virions by endocytic protease cathepsin L occurred with faster kinetics than proteolysis of wild-type T3D virions. Treatment of D-EA virions, but not T3D virions, with cathepsin D resulted in proteolysis of sigma3, a property that also was found to segregate with the D-EA S4 gene. These results indicate that a region in sigma3 protein containing amino acid 354 influences susceptibility of sigma3 to proteolysis during reovirus disassembly.     

Sequence-Specific Fidelity Alterations Associated with West Nile Virus Attenuation in Mosquitoes

High rates of error-prone replication result in the rapid accumulation of genetic diversity of RNA viruses. Recent studies suggest that mutation rates are selected for optimal viral fitness and that modest variations in replicase fidelity may be associated with viral attenuation. Arthropod-borne viruses (arboviruses) are unique in their requirement for host cycling and may necessitate substantial genetic and phenotypic plasticity. In order to more thoroughly investigate the correlates, mechanisms and consequences of arbovirus fidelity, we selected fidelity variants of West Nile virus (WNV; Flaviviridae, Flavivirus) utilizing selection in the presence of a mutagen. We identified two mutations in the WNV RNA-dependent RNA polymerase associated with increased fidelity, V793I and G806R, and a single mutation in the WNV methyltransferase, T248I, associated with decreased fidelity. Both deep-sequencing and in vitro biochemical assays confirmed strain-specific differences in both fidelity and mutational bias. WNV fidelity variants demonstrated host-specific alterations to replicative fitness in vitro, with modest attenuation in mosquito but not vertebrate cell culture. Experimental infections of colonized and field populations of Cx. quinquefaciatus demonstrated that WNV fidelity alterations are associated with a significantly impaired capacity to establish viable infections in mosquitoes. Taken together, these studies (i) demonstrate the importance of allosteric interactions in regulating mutation rates, (ii) establish that mutational spectra can be both sequence and strain-dependent, and (iii) display the profound phenotypic consequences associated with altered replication complex function of flaviviruses.     

Determinants of RNA-dependent RNA polymerase (in)fidelity revealed by kinetic analysis of the polymerase encoded by a foot-and-mouth disease virus mutant with reduced sensitivity to ribavirin

A mutant poliovirus (PV) encoding a change in its polymerase (3Dpol) at a site remote from the catalytic center (G64S) confers reduced sensitivity to ribavirin and forms a restricted quasispecies, because G64S 3Dpol is a high-fidelity enzyme. A foot-and-mouth disease virus (FMDV) mutant that encodes a change in the polymerase catalytic site (M296I) exhibits reduced sensitivity to ribavirin without restricting the viral quasispecies. In order to resolve this apparent paradox, we have established a minimal kinetic mechanism for nucleotide addition by wild-type (WT) FMDV 3Dpol that permits a direct comparison to PV 3Dpol as well as to FMDV 3Dpol derivatives. Rate constants for correct nucleotide addition were on par with those of PV 3Dpol, but apparent binding constants for correct nucleotides were higher than those observed for PV 3Dpol. The A-to-G transition frequency was calculated to be 1/20,000, which is quite similar to that calculated for PV 3Dpol. The analysis of FMDV M296I 3Dpol revealed a decrease in the calculated ribavirin incorporation frequency (1/8,000) relative to that (1/4,000) observed for the WT enzyme. Unexpectedly, the A-to-G transition frequency was higher (1/8,000) than that observed for the WT enzyme. Therefore, FMDV selected a polymerase that increases the frequency of the misincorporation of natural nucleotides while specifically decreasing the frequency of the incorporation of ribavirin nucleotide. These studies provide a mechanistic framework for understanding FMDV 3Dpol structure-function relationships, provide the first direct analysis of the fidelity of FMDV 3Dpol in vitro, identify the β9-α11 loop as a (in)fidelity determinant, and demonstrate that not all ribavirin-resistant mutants will encode high-fidelity polymerases.     

Attenuating mutations of the matrix gene of influenza A/WSN/33 virus

The matrix protein (M1) of influenza virus plays an essential role in viral replication. Our previous studies have shown that basic amino acids 101RKLKR105 of M1 are involved in RNP binding and nuclear localization. For the present work, the functions of 101RKLKR105 were studied by introducing mutations into the M gene of influenza virus A/WSN/33 by reverse genetic methods. Individual substitution, R101S or R105S, had a minimal effect on viral replication. In contrast, the double mutation R101S-R105S was synergistic and resulted in temperature sensitivity reflected by reduced viral replication at a restrictive temperature. To investigate the in vivo effect on infection, BALB/c mice were infected with either A/WSN/33 wild-type (Wt) or mutant viruses and assessed for signs of illness, viral replication in the lungs, and survival rates. The results from mouse studies indicated that the R101S-R105S double mutant virus was strongly attenuated, while single mutant viruses R101S and R105S were minimally attenuated compared to A/WSN33 Wt under the same conditions. In challenge studies, mice immunized by infection with R101S-R105S were fully protected from lethal challenge with A/WSN/33. The replication and attenuating properties of R101S-R105S suggest its potential in development of live influenza virus vaccines.     

Rapid viral escape at an immunodominant simian-human immunodeficiency virus cytotoxic T-lymphocyte epitope exacts a dramatic fitness cost

Escape from specific T-cell responses contributes to the progression of human immunodeficiency virus type 1 (HIV-1) infection. T-cell escape viral variants are retained following HIV-1 transmission between major histocompatibility complex (MHC)-matched individuals. However, reversion to wild type can occur following transmission to MHC-mismatched hosts in the absence of cytotoxic T-lymphocyte (CTL) pressure, due to the reduced fitness of the escape mutant virus. We estimated both the strength of immune selection and the fitness cost of escape variants by studying the rates of T-cell escape and reversion in pigtail macaques. Near-complete replacement of wild-type with T-cell escape viral variants at an immunodominant simian immunodeficiency virus Gag epitope KP9 occurred rapidly (over 7 days) following infection of pigtail macaques with SHIVSF162P3. Another challenge virus, SHIVmn229, previously serially passaged through pigtail macaques, contained a KP9 escape mutation in 40/44 clones sequenced from the challenge stock. When six KP9-responding animals were infected with this virus, the escape mutation was maintained. By contrast, in animals not responding to KP9, rapid reversion of the K165R mutation occurred over 2 weeks after infection. The rapidity of reversion to the wild-type sequence suggests a significant fitness cost of the T-cell escape mutant. Quantifying both the selection pressure exerted by CTL and the fitness costs of escape mutation has important implications for the development of CTL-based vaccine strategies.     

A single mutation in the carboxy terminus of reovirus outer-capsid protein sigma 3 confers enhanced kinetics of sigma 3 proteolysis,resistance to inhibitors of viral disassembly,and alterations in sigma 3 structure

Mammalian reoviruses undergo acid-dependent proteolytic disassembly within endosomes, resulting in formation of infectious subvirion particles (ISVPs). ISVPs are obligate intermediates in reovirus disassembly that mediate viral penetration into the cytoplasm. The initial biochemical event in the reovirus disassembly pathway is the proteolysis of viral outer-capsid protein sigma 3. Mutant reoviruses selected during persistent infection of murine L929 cells (PI viruses) demonstrate enhanced kinetics of viral disassembly and resistance to inhibitors of endocytic acidification and proteolysis. To identify sequences in sigma 3 that modulate acid-dependent and protease-dependent steps in reovirus disassembly, the sigma 3 proteins of wild-type strain type 3 Dearing; PI viruses L/C, PI 2A1, and PI 3-1; and four novel mutant sigma 3 proteins were expressed in insect cells and used to recoat ISVPs. Treatment of recoated ISVPs (rISVPs) with either of the endocytic proteases cathepsin L or cathepsin D demonstrated that an isolated tyrosine-to-histidine mutation at amino acid 354 (Y354H) enhanced sigma 3 proteolysis during viral disassembly. Yields of rISVPs containing Y354H in sigma3 were substantially greater than those of rISVPs lacking this mutation after growth in cells treated with either acidification inhibitor ammonium chloride or cysteine protease inhibitor E64. Image reconstructions of electron micrographs of virus particles containing wild-type or mutant sigma 3 proteins revealed structural alterations in sigma 3 that correlate with the Y354H mutation. These results indicate that a single mutation in sigma 3 protein alters its susceptibility to proteolysis and provide a structural framework to understand mechanisms of sigma 3 cleavage during reovirus disassembly.     

Effect of D222G mutation in the hemagglutinin protein on receptor binding, pathogenesis and transmissibility of the 2009 pandemic H1N1 influenza virus

Influenza viruses isolated during the 2009 H1N1 pandemic generally lack known molecular determinants of virulence associated with previous pandemic and highly pathogenic avian influenza viruses. The frequency of the amino acid substitution D222G in the hemagglutinin (HA) of 2009 H1N1 viruses isolated from severe but not mild human cases represents the first molecular marker associated with enhanced disease. To assess the relative contribution of this substitution in virus pathogenesis, transmission, and tropism, we introduced D222G by reverse genetics in the wild-type HA of the 2009 H1N1 virus, A/California/04/09 (CA/04). A dose-dependent glycan array analysis with the D222G virus showed a modest reduction in the binding avidity to human-like (α2-6 sialylated glycan) receptors and an increase in the binding to avian-like (α2-3 sialylated glycan) receptors in comparison with wild-type virus. In the ferret pathogenesis model, the D222G mutant virus was found to be similar to wild-type CA/04 virus with respect to lethargy, weight loss and replication efficiency in the upper and lower respiratory tract. Moreover, based on viral detection, the respiratory droplet transmission properties of these two viruses were found to be similar. The D222G virus failed to productively infect mice inoculated by the ocular route, but exhibited greater viral replication and weight loss than wild-type CA/04 virus in mice inoculated by the intranasal route. In a more relevant human cell model, D222G virus replicated with delayed kinetics compared with wild-type virus but to higher titer in human bronchial epithelial cells. These findings suggest that although the D222G mutation does not influence virus transmission, it may be considered a molecular marker for enhanced replication in certain cell types.     

Viability of poliovirus/rhinovirus VPg chimeric viruses and identification of an amino acid residue in the VPg gene critical for viral RNA replication

Picornaviral RNA replication utilizes a small virus-encoded protein, termed 3B or VPg, as a primer to initiate RNA synthesis. This priming step requires uridylylation of the VPg peptide by the viral polymerase protein 3D(pol), in conjunction with other viral or host cofactors. In this study, we compared the viral specificity in 3D(pol)-catalyzed uridylylation reactions between poliovirus (PV) and human rhinovirus 16 (HRV16). It was found that HRV16 3D(pol) was able to uridylylate PV VPg as efficiently as its own VPg, but PV 3D(pol) could not uridylylate HRV16 VPg. Two chimeric viruses, PV containing HRV16 VPg (PV/R16-VPg) and HRV16 containing PV VPg (R16/PV-VPg), were constructed and tested for replication capability in H1-HeLa cells. Interestingly, only PV/R16-VPg chimeric RNA produced infectious virus particles upon transfection. No viral RNA replication or cytopathic effect was observed in cells transfected with R16/PV-VPg chimeric RNA, despite the ability of HRV16 3D(pol) to uridylylate PV VPg in vitro. Sequencing analysis of virion RNA isolated from the virus particles generated by PV/R16-VPg chimeric RNA identified a single residue mutation in the VPg peptide (Glu(6) to Val). Reverse genetics confirmed that this mutation was highly compensatory in enhancing replication of the chimeric viral RNA. PV/R16-VPg RNA carrying this mutation replicated with similar kinetics and magnitude to wild-type PV RNA. This cell culture-induced mutation in HRV16 VPg moderately increased its uridylylation by PV 3D(pol) in vitro, suggesting that it might be involved in other function(s) in addition to the direct uridylylation reaction. This study demonstrated the use of chimeric viruses to characterize viral specificity and compatibility in vivo between PV and HRV16 and to identify critical amino acid residue(s) for viral RNA replication.     

Efficient transmission of pandemic H1N1 influenza viruses with high-level oseltamivir resistance

The limited availability of approved influenza virus antivirals highlights the importance of studying the fitness and transmissibility of drug-resistant viruses. S247N is a novel, naturally occurring N1 neuraminidase mutation that reduces oseltamivir sensitivity and greatly potentiates oseltamivir resistance in the context of the H275Y mutation. Here we show that highly oseltamivir-resistant viruses containing both the S247N and H275Y mutations transmit efficiently in the guinea pig transmission model.