|Home | About | Journals | Submit | Contact Us | Français|
Hepatitis A virus (HAV), an atypical member of the Picornaviridae, grows poorly in cell culture. To define determinants of HAV growth, we introduced a blasticidin (Bsd) resistance gene into the virus genome and selected variants that grew at high concentrations of Bsd. The mutants grew fast and had increased rates of RNA replication and translation but did not produce significantly higher virus yields. Nucleotide sequence analysis and reverse genetic studies revealed that a T6069G change resulting in a F42L amino acid substitution in the viral polymerase (3Dpol) was required for growth at high Bsd concentrations whereas a silent C7027T mutation enhanced the growth rate. Here, we identified a novel determinant(s) in 3Dpol that controls the kinetics of HAV growth.
Hepatitis A virus (HAV) is an atypical member of the Picornaviridae that replicates poorly in cell culture and generally does not cause cytopathic effect (CPE). The HAV positive-strand RNA genome of about 7.5 kb is encapsidated in a 27- to 32-nm icosahedral shell (12). The HAV genome contains a long open reading frame (ORF) that codes for a polyprotein of approximately 250 kDa, which undergoes co- and posttranslational processing by the virus-encoded protease 3Cpro into structural (VP0, VP3, and VP1-2A) and nonstructural proteins (11, 13, 14, 18). VP0 undergoes structural cleavage into VP2 and VP4, and an unknown cellular protease cleaves the VP1-2A precursor (9, 23).
HAV replicates inefficiently in cell culture and in general establishes persistent infections (3, 4, 7, 8) without causing CPE. However, some strains of HAV that replicate quickly can induce cell death (5, 19, 27). Due to the growth limitations, experimentation with HAV is difficult and the biology of this virus is poorly understood. To facilitate genetic studies, we recently introduced a blasticidin (Bsd) resistance gene at the 2A-2B junction of wild-type (wt) HAV (16). Bsd, an antibiotic that blocks translation in prokaryotes and eukaryotes and thus affects HAV translation, is inactivated by the Bsd-deaminase encoded in the Bsd resistance gene (15). Cells infected with the wt HAV construct carrying the Bsd resistance gene (HAV-Bsd) grew in the presence of Bsd. We have recently used the wt HAV-Bsd construct to select human hepatoma cell lines that support the stable growth of wt HAV (16) and to establish simple and rapid neutralization and virus titration assays (17). In this study, we developed a genetic approach to study determinants of HAV replication based on the selection of HAV-Bsd variants grown under increased concentrations of Bsd. We hypothesized that by increasing the concentration of Bsd, we would select HAV variants that grew better and allowed the survival of persistently infected cells at higher concentrations of the antibiotic. We also reasoned that we would need a robust HAV-Bsd replication system to provide enough Bsd-deaminase for cell survival. Therefore, we used attenuated HAV grown in rhesus monkey fetal kidney FRhK4 cells as an experimental system because (i) the virus grows 100-fold better in this system than wt HAV in human hepatoma cells (16), and (ii) it already contains cell culture-adapting mutations (3, 4, 7, 8) that are likely to accumulate during passage of wt HAV at high concentrations of Bsd and confound our analysis.
To test our hypothesis, we introduced a Bsd resistance gene into the genome of the attenuated HM-175 strain of HAV by overlap reverse transcription-PCR (RT-PCR) using the same strategy and primers described previously (16). Briefly, we first cloned a polylinker flanked by 3Cpro protease cleavage sites into the 2A-2B junction of the infectious cDNA of HAV in pT7HAV (29) to produce pHAVvec9 and then cloned the Bsd resistance gene into the polylinker to generate pHAVvec9-Bsd (Fig. (Fig.1A).1A). All constructs were verified by automated nucleotide sequence analysis. To rescue viruses (16), linearized plasmids were transcribed in vitro with T7 RNA polymerase and transfected into FRhK4 cells. After 2 weeks of incubation at 35°C, cells were washed, lysed by three freeze-thaw cycles, and spun down. Supernatants containing the viral stocks were stored at −80°C. Viruses rescued from pT7HAV, pHAVvec9, and pHAVvec9-Bsd transfected cells were termed HAV/7, HAVvec9, and HAVvec9-Bsd, respectively. Immunofluorescence (IF) analysis staining with a neutralizing anti-HAV monoclonal antibody (MAb) revealed that FRhK4 cells infected with HAV/7, HAVvec9, and HAVvec9-Bsd had the characteristic cytoplasmic granular fluorescence of HAV-infected cells (Fig. (Fig.1B).1B). The insertion of the polylinker had no effect in viral titers, but the Bsd resistance gene reduced 100-fold the HAV titers (Fig. (Fig.1C).1C). As expected, FRhK4 cells infected with HAVvec9-Bsd but not with HAV/7, HAVvec9, or mock-infected cells survived selection with 1 μg/ml Bsd (Fig. (Fig.2A2A).
To select viral variants capable of growing at higher levels of the antibiotic, we performed serial passages of HAVvec9-Bsd in FRhK4 cells at increasing concentrations of Bsd (Fig. (Fig.2B).2B). Briefly, cells were infected with HAVvec9-Bsd at a multiplicity of infection (MOI) of 1 50% tissue culture infective dose (TCID50)/cell and incubated for 6 weeks in the presence of Bsd. Cells were lysed, and viral stocks were prepared and used to infect naive FRhK4 cells at a higher concentration of Bsd. Serial passages under increasing concentrations of Bsd were performed at 1, 5, and 20 μg/ml Bsd but could not be continued at higher concentrations of the antibiotic because infected cells stopped growing at more than 20 μg/ml Bsd. Virus stocks produced at 1, 5, and 20 μg/ml Bsd, which were termed HAVvec9-Bsd, HAVvec9-Bsd-5, and HAVvec9-Bsd-20, respectively, had similar viral titers as assessed in FRhK4 cells (approximately 2 × 105 to 3 × 105 TCID50/ml) and contained comparable genome equivalents (approximately 2 × 106 to 3 × 106 copies/ml) as determined by quantitative RT-PCR analysis (24). To analyze the variants grown at increased concentrations of Bsd, we infected FRhK4 cells with HAVvec9-Bsd, HAVvec9-Bsd-5, or HAVvec9-Bsd-20 for 1 week in the presence of 1, 5, or 20 μg/ml blasticidin and stained surviving cells with crystal violet (Fig. (Fig.2C).2C). FRhK4-infected cells formed confluent monolayers at 1 μg/ml Bsd, whereas mock-infected cells did not survive the antibiotic selection. Only a few HAVvec9-Bsd-infected cells survived selection with 5 μg/ml Bsd (observed under the microscope), and none survived at 20 μg/ml Bsd. Cells infected with HAVvec9-Bsd-5 and -20 grew subconfluent monolayers at 5 μg/ml Bsd but only formed colonies at 20 μg/ml Bsd. A one-step growth curve analysis of the variants in FRhK4 cells using 2 μg/ml Bsd, a concentration of antibiotic that allowed the growth of all the tested viruses, showed that HAVvec9-Bsd-5 and -20 grew at similar rates but faster than HAVvec9-Bsd (Fig. (Fig.2D).2D). However, these three variants reached similar viral yields at 16 days postinfection (dpi). To determine whether the difference in the growth rate was due to the asynchronous replication of the variants (2, 6, 10), we performed IF analysis staining with an anti-HAV MAb (Fig. (Fig.2E).2E). At 4 dpi, a time point at which the variants reached the maximum titer differential, all cells in the monolayers contained viral antigens. Consequently, differences in asynchronous replication were not responsible for the increased growth rate of the variants. However, as judged by the intensity of the fluorescence signal, cells infected with HAVvev9-Bsd contained less HAV antigen than cells infected with HAVvev9-Bsd-20. Quantitative RT-PCR analysis (24) also showed that the kinetics of viral RNA replication was faster in HAVvec9-Bsd-20-infected cells than in HAVvec9-Bsd-infected cells (Fig. (Fig.3A).3A). RNA replication of the HAV variants reached a maximum plateau at 3 dpi (data not shown), which contrasts with the burst at 3 dpi and further increase in RNA replication observed in slow-growing HAV strains (6). Western blot analysis of infected cells stained with rabbit anti-VP2 antibody (29), which reacts with viral capsid proteins VP2 and VP0, revealed higher levels of translation in HAVvec9-Bsd-20-infected cells than in HAVvec9-Bsd-infected cells (Fig. (Fig.3B,3B, upper panel). Bands corresponding to a constitutive FRhK4 protein that cross-reacted with the rabbit anti-VP-2 antibody confirmed that similar amounts of total protein were loaded in each line (lower panel). Taken together, our data showed that the HAV variants selected under increased concentrations of Bsd grew faster than the parental virus. Therefore, we named these viruses fast-replicating (fr) HAV mutants.
To study the HAV mutants in detail, we obtained six virus clones of each viral stock by endpoint dilution in 96-well plates containing FRhK4 cells. We then amplified the RNA of the viral clones by RT-PCR and performed automatic nucleotide sequence analysis of the complete HAV genomes (Table (Table1).1). All the HAVvec9-Bsd clones had the same sequence, which was identical to that of the input virus. Therefore, insertion of the Bsd gene, with no further selection, permitted the survival of cells at 1 μg/ml Bsd with no changes in viral RNA sequence. The survival of a few HAVvec9-Bsd-infected cells at 5 μg/ml Bsd suggested the presence of a small proportion of mutants capable of supporting cell growth at higher levels of the antibiotic. All the HAVvec9-Bsd-5 clones contained a silent C7027T change in 3Dpol. Four of the HAVvec9-Bsd-5 clones contained a T6069G mutation in addition to the C7027T change and grew well at 1, 5, and 20 μg/ml Bsd. It should be pointed out that the T6069G mutation resulted in an F-to-L amino acid change at amino acid 42 of 3Dpol (F42L). Interestingly, the two HAVvec9-Bsd-5 clones that only had the silent C7027T change did not grow at 20 μg/ml Bsd. Five HAVvec9-Bsd-20 clones that contained the C7027T and T6069G changes plus two additional silent mutations, C5545T and C6780T, grew well at all the tested concentrations of Bsd. The remaining HAVvec9-Bsd-20 clone, which did not grow in 20 μg/ml Bsd, contained the same three silent mutations but lacked the T6069G change. These data suggested that the coding T6069G mutation but not the silent C5545T and C6780T mutations was required for the fr phenotype and resistance to 20 μg/ml Bsd. However, the role of the C7027T mutation in the fr phenotype remained elusive.
To further determine the role of the T6069G in the fr phenotype, we engineered the mutation into the infectious HAVvec9-Bsd cDNA using overlap PCR and synthetic oligonucleotides containing single-nucleotide substitutions as described previously (16). The virus was rescued in FRhK4 cells, and the presence of the introduce mutation was verified by RT-PCR and nucleotide sequence analysis. The HAVvec9-Bsd mutant containing the coding T6069G change, termed HAVvec9Bsd-6069, grew faster than HAVvec9-Bsd but slower than HAVvec9-Bsd-20 in FRhK4 cells (Fig. (Fig.3C),3C), indicating that the T6069G mutation partially restored the fr phenotype. Interestingly, HAVvec9Bsd-6069 grew at 20 μg/ml Bsd (Table (Table1),1), indicating that the increase in the replication rate allowed the survival of infected cells at higher concentrations of the antibiotic.
Similarly, we introduced the C7027T change into the infectious cDNA of HAVvec9-Bsd and analyzed the phenotype of the rescued virus. The mutant containing the C7027T change, termed HAVvec9Bsd-7027, grew slowly and at the same rate than parental HAVvec9-Bsd (Fig. (Fig.3C),3C), indicating that the C7027T mutation alone did not confer the fr phenotype. Moreover, HAVvec9Bsd-7027 did not grow at 20 μg/ml Bsd (Table (Table1).1). The mutant containing both the T6069G and C7027T, termed HAVvec9-Bsd-6069 + 7027, grew as fast as HAVvec9-Bsd-20, indicating that the silent C7027T change had a synergistic effect and enhanced the replication rate conferred by the T6069G coding mutation. As expected, HAVvec9-Bsd-6069 + 7027 grew at 20 μg/ml Bsd (Table (Table1).1). Therefore, these two mutations were required to convey the full fr phenotype observed in HAVvec9-Bsd-20.
It was of interest to map the location of the F42 residue in 3Dpol affected by the mutation in nucleotide 6069. Alignment of the first 69 amino acid residues of 3Dpol of HAV with other members of the Picornaviridae (Fig. (Fig.3D)3D) revealed that two amino acid motifs, KT and PAV/A (shown in blue), are highly conserved in the family. The HAV F42 (magenta) is 3 residues from the PAV/A motif and aligned with K44 of the PV 3Dpol, for which the crystal structure has been resolved at high resolution (26). Although the HAV 3Dpol has not been crystallized, 3Dpol structures are conserved in the Picornaviridae (20), so the HAV F42 is likely to be in the proximity of the PV K44. According to the analogy of the “right hand” shape with the structure of 3Dpol (Fig. (Fig.3E),3E), the N-terminal 69 amino acid residues of PV 3Dpol form an “index finger” (residues in blue) (26) with PV K44 (magenta) lying in the middle. The fingers of 3Dpol are very flexible and exist in a highly dynamic molten globule state at physiological temperature (25). The “index finger” plays an important role in maintaining this flexibility and the function of 3Dpol (25). Residues at the tip of the “index finger” such as F30 (red) and F34 (green) in the PV 3Dpol interact with the top of the “thumb” and are important for the proper function of the polymerase (25). The G64 residue at the top of the “index finger” in PV 3Dpol (yellow) modulates the catalytic site of 3Dpol (1, 21, 26), and the G64S change reduced the rate for nucleotide incorporation of the enzyme (20, 26). Our results suggested that the HAV F42L change had the opposite effect and increased the kinetics of RNA replication (Fig. (Fig.3A)3A) and HAV growth (Fig. (Fig.3C).3C). In humans, wt HAV has a long incubation period of 7 to 50 days that could be attributed to the slow replication rate limiting the spread of the virus. In cell culture, wt HAV also grows slowly without causing CPE (16). Therefore, it is likely that HAV contains determinants capable of limiting its replication, such as the ones we identified in T6069 and C7027, to prevent cell damage and escape immune surveillance.
In this study, we showed that a T6069G mutation that resulted in the F42L amino acid change in 3Dpol increased the rate of RNA replication and HAV-specific translation. Interestingly, the faster kinetics of HAV growth did not result in an increased virus yield, indicating that other limiting factors in the life cycle of HAV, such as the use of rare codons (22) or the availability of cellular factors, are also responsible for its poor growth. The silent C7027T mutation in the 3Dpol that synergized with the T6069G change but was not sufficient to increase the rate of growth of the virus (Fig. (Fig.3C)3C) suggested that RNA structure and/or codon utilization also played a role in the fr phenotype. Recently, a cis-acting replication element (cre) was identified between nucleotides 5948 and 6057 of the HAV genome in the 3Dpol gene (28). Since nucleotide 6069 is adjacent to the cre element, we cannot rule out the possibility that the T6069G change affected the structure of the cre, increasing the rate of RNA replication. Our attempts to rescue viruses containing a silent mutation at nucleotide 6069 have been unsuccessful. Further research will be required to determine whether the fr phenotype was due to the F42L change in 3Dpol, the effect of the T6069G change in the cre, or both. The 3Dpol mutants described herein could be used to shorten the time required to produce HAV antigen and reduce the cost of vaccine production.
This work was supported by intramural research funds from the Food and Drug Administration (FDA) to G.G.K.
The findings and conclusions in this article have not been formally disseminated by the FDA and should not be construed to represent any FDA determination or policy.
We thank Susan Zullo and Jérôme Jacques for critical reviews of the manuscript.
Published ahead of print on 9 June 2010.