Viral RNA mutations are region specific and increased by ribavirin in a full-length hepatitis C virus replication system

High rates of genetic variation ensure the survival of RNA viruses. Although this variation is thought to result from error-prone replication, RNA viruses must also maintain highly conserved genomic segments. A balance between conserved and variable viral elements is especially important in order for viruses to avoid "error catastrophe." Ribavirin has been shown to induce error catastrophe in other RNA viruses. We therefore used a novel hepatitis C virus (HCV) replication system to determine relative mutation frequencies in variable and conserved regions of the HCV genome, and we further evaluated these frequencies in response to ribavirin. We sequenced the 5' untranslated region (5' UTR) and the core, E2 HVR-1, NS5A, and NS5B regions of replicating HCV RNA isolated from cells transfected with a T7 polymerase-driven full-length HCV cDNA plasmid containing a cis-acting hepatitis delta virus ribozyme to control 3' cleavage. We found quasispecies in the E2 HVR-1 and NS5B regions of untreated replicating viral RNAs but not in conserved 5' UTR, core, or NS5A regions, demonstrating that important cis elements regulate mutation rates within specific viral segments. Neither T7-driven replication nor sequencing artifacts produced these nucleotide substitutions in control experiments. Ribavirin broadly increased error generation, especially in otherwise invariant regions, indicating that it acts as an HCV RNA mutagen in vivo. Similar results were obtained in hepatocyte-derived cell lines. These results demonstrate the potential utility of our system for the study of intrinsic factors regulating genetic variation in HCV. Our results further suggest that ribavirin acts clinically by promoting nonviable HCV RNA mutation rates. Finally, the latter result suggests that our replication model may be useful for identifying agents capable of driving replicating virus into error catastrophe.     

Specific interaction of hepatitis C virus protease/helicase NS3 with the 3'-terminal sequences of viral positive- and negative-strand RNA

The hepatitis C virus (HCV)-encoded protease/helicase NS3 is likely to be involved in viral RNA replication. We have expressed and purified recombinant NS3 (protease and helicase domains) and Delta pNS3 (helicase domain only) and examined their abilities to interact with the 3'-terminal sequence of both positive and negative strands of HCV RNA. These regions of RNA were chosen because initiation of RNA synthesis is likely to occur at or near the 3' untranslated region (UTR). The results presented here demonstrate that NS3 (and Delta pNS3) interacts efficiently and specifically with the 3'-terminal sequences of both positive- and negative-strand RNA but not with the corresponding complementary 5'-terminal RNA sequences. The interaction of NS3 with the 3'-terminal negative strand [called 3'(-) UTR(127)] was specific in that only homologous (and not heterologous) RNA competed efficiently in the binding reaction. A predicted stem-loop structure present at the 3' terminus (nucleotides 5 to 20 from the 3' end) of the negative-strand RNA appears to be important for NS3 binding to the negative-strand UTR. Deletion of the stem-loop structure almost totally impaired NS3 (and Delta pNS3) binding. Additional mutagenesis showed that three G-C pairs within the stem were critical for helicase-RNA interaction. The data presented here also suggested that both a double-stranded structure and the 3'-proximal guanosine residues in the stem were important determinants of protein binding. In contrast to the relatively stringent requirement for 3'(-) UTR binding, specific interaction of NS3 (or Delta pNS3) with the 3'-terminal sequences of the positive-strand RNA [3'(+) UTR] appears to require the entire 3'(+) UTR of HCV. Deletion of either the 98-nucleotide 3'-terminal conserved region or the 5' half sequence containing the variable region and the poly(U) and/or poly(UC) stretch significantly impaired RNA-protein interaction. The implication of NS3 binding to the 3'-terminal sequences of viral positive- and negative-strand RNA in viral replication is discussed.     

HCV RNA-dependent RNA polymerase replicates in vitro the 3' terminal region of the minus-strand viral RNA more efficiently than the 3' terminal region of the plus RNA.

The NS5B protein, or RNA-dependent RNA polymerase of the hepatitis virus type C, catalyzes the replication of the viral genomic RNA. Little is known about the recognition domains of the viral genome by the NS5B. To better understand the initiation of RNA synthesis on HCV genomic RNA, we used in vitro transcribed RNAs as templates for in vitro RNA synthesis catalyzed by the HCV NS5B. These RNA templates contained different regions of the 3' end of either the plus or the minus RNA strands. Large differences were obtained depending on the template. A few products shorter than the template were synthesized by using the 3' UTR of the (+) strand RNA. In contrast the 341 nucleotides at the 3' end of the HCV minus-strand RNA were efficiently copied by the purified HCV NS5B in vitro. At least three elements were found to be involved in the high efficiency of the RNA synthesis directed by the HCV NS5B with templates derived from the 3' end of the minus-strand RNA: (a) the presence of a C residue as the 3' terminal nucleotide; (b) one or two G residues at positions +2 and +3; (c) other sequences and/or structures inside the following 42-nucleotide stretch. These results indicate that the 3' end of the minus-strand RNA of HCV possesses some sequences and structure elements well recognized by the purified NS5B.     

Production of infectious chimeric hepatitis C virus genotype 2b harboring minimal regions of JFH-1

To establish a cell culture system for chimeric hepatitis C virus (HCV) genotype 2b, we prepared a chimeric construct harboring the 5' untranslated region (UTR) to the E2 region of the MA strain (genotype 2b) and the region of p7 to the 3' UTR of the JFH-1 strain (genotype 2a). This chimeric RNA (MA/JFH-1.1) replicated and produced infectious virus in Huh7.5.1 cells. Replacement of the 5' UTR of this chimera with that from JFH-1 (MA/JFH-1.2) enhanced virus production, but infectivity remained low. In a long-term follow-up study, we identified a cell culture-adaptive mutation in the core region (R167G) and found that it enhanced virus assembly. We previously reported that the NS3 helicase (N3H) and the region of NS5B to 3' X (N5BX) of JFH-1 enabled replication of the J6CF strain (genotype 2a), which could not replicate in cells. To reduce JFH-1 content in MA/JFH-1.2, we produced a chimeric viral genome for MA harboring the N3H and N5BX regions of JFH-1, combined with a JFH-1 5' UTR replacement and the R167G mutation (MA/N3H+N5BX-JFH1/R167G). This chimeric RNA replicated efficiently, but virus production was low. After the introduction of four additional cell culture-adaptive mutations, MA/N3H+N5BX-JFH1/5am produced infectious virus efficiently. Using this chimeric virus harboring minimal regions of JFH-1, we analyzed interferon sensitivity and found that this chimeric virus was more sensitive to interferon than JFH-1 and another chimeric virus containing more regions from JFH-1 (MA/JFH-1.2/R167G). In conclusion, we established an HCV genotype 2b cell culture system using a chimeric genome harboring minimal regions of JFH-1. This cell culture system may be useful for characterizing genotype 2b viruses and developing antiviral strategies.     

Specific interaction between the classical swine fever virus NS5B protein and the viral genome.

The NS5B protein of the classical swine fever virus (CSFV) is the RNA-dependent RNA polymerase of the virus and is able to catalyze the viral genome replication. The 3' untranslated region is most likely involved in regulation of the Pestivirus genome replication. However, little is known about the interaction between the CSFV NS5B protein and the viral genome. We used different RNA templates derived from the plus-strand viral genome, or the minus-strand viral genome and the CSFV NS5B protein obtained from the Escherichia coli expression system to address this problem. We first showed that the viral NS5B protein formed a complex with the plus-strand genome through the genomic 3' UTR and that the NS5B protein was also able to bind the minus-strand 3' UTR. Moreover, it was found that viral NS5B protein bound the minus-strand 3' UTR more efficiently than the plus-strand 3' UTR. Further, we observed that the plus-strand 3' UTR with deletion of CCCGG or 21 continuous nucleotides at its 3' terminal had no binding activity and also lost the activity for initiation of minus-strand RNA synthesis, which similarly occurred in the minus-strand 3' UTR with CATATGCTC or the 21 nucleotide fragment deleted from the 3' terminal. Therefore, it is indicated that the 3' CCCGG sequence of the plus-strand 3' UTR, and the 3' CATATGCTC fragment of the minus-strand are essential to in vitro synthesis of the minus-strand RNA and the plus-strand RNA, respectively. The same conclusion is also appropriate for the 3' 21 nucleotide terminal site of both the 3' UTRs.     

In vitro selection of RNA aptamers against the HCV NS3 helicase domain

Nonstructural protein 3 (NS3) of hepatitis C virus (HCV) has two distinct domains, protease and helicase, that are essential for HCV proliferation. Therefore, NS3 is considered a target for anti-HCV treatment. To study RNA aptamers of the NS3 helicase domain, we carried out in vitro selection against the HCV NS3 helicase domain. RNA aptamers obtained after eight generations possessed 5' extended single-stranded regions and the conserved sequence (5'-GGA(U/C)GGAGCC-3') at stem-loop regions. Aptamer 5 showed strong inhibition of helicase activity in vitro. Deletion and mutagenesis analysis clarified that the conserved stem-loop is important and that the whole structure is needed for helicase inhibition. We compared the inhibition of helicase activity between aptamer 5 and 3'+-UTR of HCV.     

The helicase activity associated with hepatitis C virus nonstructural protein 3 (NS3).

To assess the RNA helicase activity of hepatitis C virus (HCV) nonstructural protein 3 (NS3), a polypeptide encompassing amino acids 1175 to 1657, which cover only the putative helicase domain, was expressed in Escherichia coli by a pET expression vector. The protein was purified to near homogeneity and assayed for RNA helicase activity in vitro with double-stranded RNA substrates prepared from a multiple cloning sequence and an HCV 5' nontranslated region (5'-NTR) or 3'-NTR. The enzyme acted successfully on substrates containing both 5' and 3' single-stranded regions (standard) or on substrates containing only the 3' single-stranded regions (3'/3') but failed to act on substrates containing only the 5' single-stranded regions (5'/5') or on substrates lacking the single-stranded regions (blunt). These results thus suggest 3' to 5' directionality for HCV RNA helicase activity. However, a 5'/5' substrate derived from the HCV 5'-NTR was also partially unwound by the enzyme, possibly because of unique properties inherent in the 5' single-stranded regions. Gel mobility shift analyses demonstrated that the HCV NS3 helicase could bind to either 5'- or 3'-tailed substrates but not to substrates lacking a single-stranded region, indicating that the polarity of the RNA strand to which the helicase bound was a more important enzymatic activity determinant. In addition to double-stranded RNA substrates, HCV NS3 helicase activity could displace both RNA and DNA oligonucleotides on a DNA template, suggesting that HCV NS3 too was disposed to DNA helicase activity. This study also demonstrated that RNA helicase activity was dramatically inhibited by the single-stranded polynucleotides. Taken altogether, our results indicate that the HCV NS3 helicase is unique among the RNA helicases characterized so far.     

Cell culture-adaptive mutations promote viral protein-protein interactions and morphogenesis of infectious hepatitis C virus

Recent genetic studies suggested that viral nonstructural (NS) proteins play important roles in morphogenesis of flaviviruses, particularly hepatitis C virus (HCV). Adaptive and compensatory mutations occurring in different NS proteins were demonstrated to promote HCV production in cell culture. However, the underlying molecular mechanism of NS proteins in HCV morphogenesis is poorly understood. We have isolated a cell culture-adapted HCV of genotype 2a (JFH1) which grew to an infectious titer 3 orders of magnitude higher than that of wild-type virus. Sequence analysis identified a total of 16 amino acid mutations in core (C), E1, NS2, NS3, NS5A, and NS5B, with the majority of mutations clustered in NS5A. Reverse genetic analysis of these mutations individually or in different combinations demonstrated that amino acid mutations in NS2 and NS5A markedly enhanced HCV production. Additionally, mutations in C, E1, NS3, and NS5B synergistically promoted HCV production in the background of NS2 and NS5A mutations. Adaptive mutations in NS5A domains I, II, and III independently enhanced HCV production, suggesting that all three domains of NS5A are important for HCV morphogenesis. More importantly, adaptive mutations greatly enhanced physical interactions among HCV structural and NS proteins, as determined by studies with coimmunoprecipitation and mammalian two-hybrid assays. Collectively, these findings demonstrate that adaptive mutations can enhance specific protein-protein interactions among viral structural and NS proteins and therefore promote the assembly of infectious HCV particles.     

Identification of determinants involved in initiation of hepatitis C virus RNA synthesis by using intergenotypic replicase chimeras

The 5' nontranslated region (NTR) and the X tail in the 3' NTR are the least variable parts of the hepatitis C virus (HCV) genome and play an important role in the initiation of RNA synthesis. By using subgenomic replicons of the HCV isolates Con1 (genotype 1) and JFH1 (genotype 2), we characterized the genotype specificities of the replication signals contained in the NTRs. The replacement of the JFH1 5' NTR and X tail with the corresponding Con1 sequence resulted in a significant decrease in replication efficiency. Exchange of the X tail specifically reduced negative-strand synthesis, whereas substitution of the 5' NTR impaired the generation of progeny positive strands. In search for the proteins involved in the recognition of genotype-specific initiation signals, we analyzed recombinant nonstructural protein 5B (NS5B) RNA polymerases of both isolates and found some genotype-specific template preference for the 3' end of positive-strand RNA in vitro. To further address genotype specificity, we constructed a series of intergenotypic replicon chimeras. When combining NS3 to NS5A of Con1 with NS5B of JFH1, we observed more-efficient replication with the genotype 2a X tail, indicating that NS5B recognizes genotype-specific signals in this region. In contrast, a combination of the NS3 helicase with NS5A and NS5B was required to confer genotype specificity to the 5' NTR. These results present the first genetic evidence for an interaction between helicase, NS5A, and NS5B required for the initiation of RNA synthesis and provide a system for the specific analysis of HCV positive- and negative-strand syntheses.     

Replication of a hepatitis C virus replicon clone in mouse cells

Background

The duration of treatment for HCV infection is partly indicated by the genotype of the virus. For studies of disease transmission, vaccine design, and surveillance for novel variants, subtype-level classification is also needed. This study used the Shimodaira-Hasegawa test and related statistical techniques to compare phylogenetic trees obtained from coding and non-coding regions of a whole-genome alignment for the reliability of subtyping in different regions.

Results

Different regions of the HCV genome yield inconsistent phylogenies, which can lead to erroneous conclusions about classification of a given infection. In particular, the highly conserved 5' untranslated region (UTR) yields phylogenetic trees with topologies that differ from the HCV polyprotein and complete genome phylogenies. Phylogenetic trees from the NS5B gene reliably cluster related subtypes, and yield topologies consistent with those of the whole genome and polyprotein.

Conclusion

These results extend those from previous studies and indicate that, unlike the NS5B gene, the 5' UTR contains insufficient variation to resolve HCV classifications to the level of viral subtype, and fails to distinguish genotypes reliably. Use of the 5' UTR for clinical tests to characterize HCV infection should be replaced by a subtype-informative test.     

Identification of a naturally occurring recombinant genotype 2/6 hepatitis C virus

Hepatitis C viruses (HCVs) display a high level of sequence diversity and are currently classified into six genotypes and an increasing number of subtypes. Most likely, this heterogeneity is caused by genetic drift; evidence for recombination is scarce. To study the molecular heterogeneity of HCV in Vietnam, we analyzed 58 HCV RNA-positive sera from Vietnamese blood donors by sequence analysis of the CORE and NS5B regions. Phylogenetic analyses revealed the presence of genotype 1 (38%), genotype 2 (10.3%), and genotype 6 viruses (51.7%). All samples showed concordant results except for two (D3 and D54). Sample D54 was a mixed infection of genotype 2i and 6h viruses. Whole-genome analysis and bootscan analysis of sample D3, on the other hand, revealed a recombinant virus with genotype 2i and genotype 6p sequences at the 5' and 3' ends, respectively. The crossover point was located between nucleotide positions 3405 to 3464 (numbering according to prototype strain HCV-H, M67463) at the NS2/NS3 junction. The identification of this naturally occurring recombinant virus strengthens the concept that recombination may play a role in HCV epidemiology and evolution. Furthermore, the location of the recombination breakpoint may be relevant for constructing infectious chimeric viruses.     

Genetic variation in the 5' end and NS5B regions of classical swine fever virus genome among Japanese isolates.

Sixteen clinical strains of classical swine fever virus (CSFV) isolated in Japan were subjected to analyses of nucleotide sequence variations in the 5' end and NS5B regions of the genome. These isolates were divided into three genovars, CSFV-1, CSFV-2 and CSFV-3, based on palindromic nucleotide substitutions at the three variable loci in the 5' untranslated region (UTR). Phylogenetic trees constructed from nucleotide sequences in the 5'-UTR and NS5B gene indicated that the CSFV strains were divided into three clusters, I, II and III. CSFV strains included in clusters I, II and III were identical to those in the CSFV-1, CSFV-2 and CSFV-3 genovars, respectively.     

Determination of Sustained Virological Response in Hepatitis C Virus Genotypes by the Number of Mutations in the E2 and NS5A-ISDR Regions: A Meta-Analysis

Since last few decades hepatitis C virus (HCV) has been a major cause of death due to the involvement of acute and chronic type of liver diseases throughout the world. Genotype variability and mutations occurring at different regions of HCV genome provides a critical parameter for the study of sustained virological response (SVR) against mono and combinational therapies. Most of these mutations occurring in E2 and NS5A-ISDR regions in HCV genotypes play a significant role in SVR against Interferon-monotherapy and combination therapy. Therefore, the aim of the present study is to evaluate the SVR in various genotypes and the role of mutations in specific regions. In line with this, the NS5A and E2 proteins of HCV genotype 1 were found to suppress the double-stranded (ds) RNA-dependent protein kinase (PKR), which in turn is entailed in the cellular antiviral response stimulated by interferon (IFN). The response to IFN therapy varies between genotypes, with response rates among patients infected with types 2 and 3 nearly two-three-fold greater than in patients infected with type 1. Surprisingly, a considerable percentage of HCV genotype 3a infected patients do not react to treatment at all. In Japan, a link was observed between the numbers of mutations in an “interferon sensitivity determining region” (ISDR) and the result of interferon treatment in genotype 1b infected patients. Therefore, we published data on E2 (PePHD) and NS5A-ISDR regions including our data on SVR of different HCV genotypes, the relationship between the number and patterns of mutations in the E2 (PePHD) and NS5A-ISDR regions and responsiveness to IFN therapy.     

An optimal cis-replication stem-loop IV in the 5' untranslated region of the mouse coronavirus genome extends 16 nucleotides into open reading frame 1

The 288-nucleotide (nt) 3' untranslated region (UTR) in the genome of the bovine coronavirus (BCoV) and 339-nt 3' UTR in the severe acute respiratory syndrome (SARS) coronavirus (SCoV) can each replace the 301-nt 3' UTR in the mouse hepatitis coronavirus (MHV) for virus replication, thus demonstrating common 3' cis-replication signals. Here, we show that replacing the 209-nt MHV 5' UTR with the ～63%-sequence-identical 210-nt BCoV 5' UTR by reverse genetics does not yield viable virus, suggesting 5' end signals are more stringent or possibly are not strictly 5' UTR confined. To identify potential smaller, 5'-common signals, each of three stem-loop (SL) signaling domains and one inter-stem-loop domain from the BCoV 5' UTR was tested by replacing its counterpart in the MHV genome. The SLI/II domain (nucleotides 1 to 84) and SLIII domain (nucleotides 85 to 141) each immediately enabled near-wild-type (wt) MHV-like progeny, thus behaving similarly to comparable 5'-proximal regions of the SCoV 5' UTR as shown by others. The inter-stem-loop domain (nt 142 to 173 between SLs III and IV) enabled small plaques only after genetic adaptation. The SLIV domain (nt 174 to 210) required a 16-nt extension into BCoV open reading frame 1 (ORF1) for apparent stabilization of a longer BCoV SLIV (nt 174 to 226) and optimal virus replication. Surprisingly, pleiomorphic SLIV structures, including a terminal loop deletion, were found among debilitated progeny from intra-SLIV chimeras. The results show the inter-stem-loop domain to be a potential novel species-specific cis-replication element and that cis-acting SLIV in the viral genome extends into ORF1 in a manner that stabilizes its lower stem and is thus not 5' UTR confined.     

Specific Interaction between the Hepatitis C Virus NS5B RNA Polymerase and the 3′ End of the Viral RNA

Hepatitis C virus (HCV) NS5B protein is the viral RNA-dependent RNA polymerase capable of directing RNA synthesis. In this study, an electrophoretic mobility shift assay demonstrated the interaction between a partially purified recombinant NS5B protein and a 3' viral genomic RNA with or without the conserved 98-nucleotide tail. The NS5B-RNA complexes were specifically competed away by the unlabeled homologous RNA but not by the viral 5' noncoding region and very poorly by the 3' conserved 98-nucleotide tail. A 3' coding region with conserved stem-loop structures rather than the 3' noncoding region of the HCV genome is critical for the specific binding of NS5B. Nevertheless, no direct interaction between the 3' coding region and the HCV NS5A protein was detected. Furthermore, two independent RNA-binding domains (RBDs) of NS5B were identified, RBD1, from amino acid residues 83 to 194, and RBD2, from residues 196 to 298. Interestingly, the conserved motifs of RNA-dependent RNA polymerase for putative RNA binding (220-DxxxxD-225) and template/primer position (282-S/TGxxxTxxxNS/T-292) are present in the RBD2. Nevertheless, the RNA-binding activity of RBD2 was abolished when it was linked to the carboxy-terminal half of the NS5B. These results provide some clues to understanding the initiation of HCV replication.     

cis-acting RNA signals in the NS5B C-terminal coding sequence of the hepatitis C virus genome

The cis-replicating RNA elements in the 5' and 3' nontranslated regions (NTRs) of the hepatitis C virus (HCV) genome have been thoroughly studied before. However, no cis-replicating elements have been identified in the coding sequences of the HCV polyprotein until very recently. The existence of highly conserved and stable stem-loop structures in the RNA polymerase NS5B coding sequence, however, has been previously predicted (A. Tuplin, J. Wood, D. J. Evans, A. H. Patel, and P. Simmonds, RNA 8:824-841, 2002). We have selected for our studies a 249-nt-long RNA segment in the C-terminal NS5B coding region (NS5BCR), which is predicted to form four stable stem-loop structures (SL-IV to SL-VII). By deletion and mutational analyses of the RNA structures, we have determined that two of the stem-loops (SL-V and SL-VI) are essential for replication of the HCV subgenomic replicon in Huh-7 cells. Mutations in the loop and the top of the stem of these RNA elements abolished replicon RNA synthesis but had no effect on translation. In vitro gel shift and filter-binding assays revealed that purified NS5B specifically binds to SL-V. The NS5B-RNA complexes were specifically competed away by unlabeled homologous RNA, to a small extent by 3' NTR RNA, and only poorly by 5' NTR RNA. The other two stem-loops (SL-IV and SL-VII) of the NS5BCR domain were found to be important but not essential for colony formation by the subgenomic replicon. The precise function(s) of these cis-acting RNA elements is not known.     

The 3' untranslated regions of influenza genomic sequences are 5'PPP-independent ligands for RIG-I

Retinoic acid inducible gene-I (RIG-I) is a key regulator of antiviral immunity. RIG-I is generally thought to be activated by ssRNA species containing a 5'-triphosphate (PPP) group or by unphosphorylated dsRNA up to ~300 bp in length. However, it is not yet clear how changes in the length, nucleotide sequence, secondary structure, and 5' end modification affect the abilities of these ligands to bind and activate RIG-I. To further investigate these parameters in the context of naturally occurring ligands, we examined RNA sequences derived from the 5' and 3' untranslated regions (UTR) of the influenza virus NS1 gene segment. As expected, RIG-I-dependent interferon-β (IFN-β) induction by sequences from the 5' UTR of the influenza cRNA or its complement (26 nt in length) required the presence of a 5'PPP group. In contrast, activation of RIG-I by the 3' UTR cRNA sequence or its complement (172 nt) exhibited only a partial 5'PPP-dependence, as capping the 5' end or treatment with CIP showed a modest reduction in RIG-I activation. Furthermore, induction of IFN-β by a smaller, U/A-rich region within the 3' UTR was completely 5'PPP-independent. Our findings demonstrated that RNA sequence, length, and secondary structure all contributed to whether or not the 5'PPP moiety is needed for interferon induction by RIG-I.