Previous Article | Next Article ![]()
Journal of Virology, September 2004, p. 9782-9789, Vol. 78, No. 18
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.18.9782-9789.2004
Deborah Taylor,1 Duncan R. Morhardt,1 Kathleen Mihalik,1 Montserrat Puig,1 Charles M. Rice,2 Stephen M. Feinstone,1 and Marian E. Major1*
Laboratory of Hepatitis Viruses, Division of Viral Products, Center for Biologics Evaluation and Research, Food and Drug Administration, Bethesda, Maryland,1 Center for the Study of Hepatitis C, Rockefeller University, New York, New York2
Received 30 January 2004/ Accepted 18 May 2004
|
|
|---|
1 x 103 substitutions/site/year. Analysis of PCR products using single-stranded conformational polymorphism indicated a low level of heterogeneity in the viral genome. The results of these studies confirm that the persistence of infection is not solely due to changes in HVR1 or heterogeneity and that the majority of variants observed in natural infections could not arise simply through mutation during the time period most humans and chimpanzees are observed. These data also indicate that immune pressure and selection continue throughout the chronic phase. |
|
|---|
The mechanisms leading to viral persistence, which is associated with the more severe forms of liver disease, are as yet undefined. Any single HCV isolate exists as a quasispecies with sequence variability throughout the RNA genome (3, 27). This variation could lead to evasion of the host immune response through the selection of neutralizing antibody or cytotoxic T-lymphocyte escape mutants and thereby the establishment of persistent infection. Evidence for both types of escape mutants have been reported in HCV infections (4, 7, 13, 31). Reports have indicated that hypervariable region 1 (HVR1), located in the N terminus of the E2 protein, evolves more rapidly in vivo than the rest of the viral genome (15, 18) and that it plays a major role in the maintenance of persistent infections (reviewed in reference 17). Previous studies have hypothesized that a higher complexity of virus species provide an indicator of progression to chronicity, particularly in HVR1 (8, 24). However, RNA transcribed from an infectious cDNA clone lacking HVR1 caused a persistent infection in a chimpanzee, indicating that this region is not essential for infection or persistence (10).
The quasispecies nature of natural isolates makes it impossible to distinguish true de novo mutations. In this study, we have examined the molecular evolution of HCV over a 4-year period in two chimpanzees infected with a virus consisting of a single sequence as the starting population. The predominant circulating virus at different times after infection was analyzed by direct sequencing of PCR amplicons. We have been able to monitor the true accumulation of mutations in the HCV genome in chimpanzees over time relative to host responses and viral kinetics. We propose that the majority of variant sequences observed during infections with quasispecies isolates, particularly multiple substitutions observed in HVR1, arise primarily through selection of sequence variants already present in the pool of quasispecies rather than by mutation during the time period most humans and chimpanzees have been observed. Our studies indicate that even over several years only single-amino-acid mutations become fixed at distinct sites and that many more years of infection would be required before significant sequence changes would be achieved.
|
|
|---|
RNA extractions and real-time RT-PCR. Total RNA was prepared from 100 µl of a serum sample or a section of a liver biopsy sample with an area of approximately 1 mm3 using TRIzol (Life Technologies, Gaithersburg, Md.) as previously described (16). RNA pellets were resuspended in 10-µl portions of RNasin-dithiothreitol-water (0.2 U of RNasin per µl of water, 10 mM dithiothreitol) (Promega, Madison, Wis.) and stored at 80°C until use. Negative controls, in the form of serum samples from Ch1535 and Ch1536 prior to inoculation or serum samples from uninoculated chimpanzees, were included in all extractions. RNA levels in serum samples were quantified by real-time reverse transcription-PCR (RT-PCR) (22).
RT-PCR. RT of purified RNA was performed using the first-strand cDNA synthesis kit (Pharmacia, Uppsala, Sweden) according to the manufacturer's instructions, and 10 pg of random hexamers in a 15- or 33-µl reaction mixture volume. Nested PCRs, consisting of 40 cycles each, were performed on RT products using the Expand High Fidelity PCR system (Boehringer Mannheim, Indianapolis, Ind.) as previously described (16). The HCV genome was amplified as overlapping 400- to 1,100-bp fragments. Primer sequences were based on the cDNA consensus clone (designated p90) used for the synthesis of the original RNA transcripts inoculated into Ch1535 and Ch1536.
Sequence analysis of PCR products. Following PCR, specific fragments were isolated from ethidium bromide-stained gels using Genelute agarose spin columns (Sigma, St. Louis, Mo.) according to the manufacturer's instructions. Sequencing reactions were performed by SeqWright (Houston, Tex.) using HCV-specific primers, and the data were analyzed using the DNAstar (Madison, Wis.) sequence analysis package.
Enzyme-linked immunosorbent assay for E1E2- and HVR1-specific antibodies. Sera were tested for HCV-specific antibodies using purified E1E2 protein or biotinylated peptides covering amino acids 384 (aa 384) to 410 of the HCV polyprotein as previously described (16). Mean optical density (OD) values were expressed as P/N ratios calculated by dividing the OD at 405 nm for test sera after infection by that obtained for preimmune serum. The cutoff value was taken as P/N = 2.
SSCP analysis of PCR products. Single-stranded conformational polymorphism (SSCP) analysis was performed using 20 ng (5 µl) of 164- to 176-bp PCR products denatured with 5 µl of 95% formamide. Products were generated by amplification of cDNA obtained by random-primed RT of viral RNA extracted from serum. DNA was resolved on 8% polyacrylamide-1x Tris-borate-EDTA gels (Novex; Invitrogen). Gels were prerun for 30 min at 100 V and 18°C, then loaded, and run for an additional 2 h at 100 V. Bands were visualized by silver staining.
|
|
|---|
![]() View larger version (29K): [in a new window] |
FIG. 1. Clinical and virologic responses up to 216 weeks p.i. in Ch1535 (A) and Ch1536 (B) after inoculation with RNA transcripts representing the HCV infectious cDNA.
|
HCV genome analysis. The predominant circulating virus from each animal was analyzed for mutations at 26 (Ch1535) or 22 (Ch1536), 60, 130, and 216 weeks postinfection (p.i.). Analysis of the viral genome sequence isolated from Ch1536 at week 22 indicated no substitutions; therefore, a sample from week 26 was analyzed from this animal. Overlapping fragments were amplified by PCR and directly sequenced. The observed nucleotide and amino acid changes are shown for each animal in Tables 1 and 2.
|
View this table: [in a new window] |
TABLE 1. Nucleotide and amino substitutions observed in the predominant circulating virus from serum samples from Ch1535a
|
|
View this table: [in a new window] |
TABLE 2. Nucleotide and amino substitutions observed in the predominant circulating virus from serum samples from Ch1536a
|
Positive versus negative or random selection. The substitution rates per codon site per year after 60 weeks were calculated as 1.57 x 103 (Ch1535) and 1.48 x 103 (Ch1536) substitutions/site at the nucleotide level and 2.59 x 103 (Ch1535) and 2.88 x 103 (Ch1536) substitutions/site at the amino acid level. A method to distinguish positive selection from random accumulation of mutations is to compare the nonsynonymous (protein sequence-altering changes) substitution rate (dN) with the synonymous (silent) substitution rate (dS). These comparisons are shown in Fig. 2A for each animal at each sampling time point. It is significant that for both animals the mutations observed early were 100% nonsynonymous, suggesting that the greatest level of selective pressure occurs during this early phase. When the substitution rates were calculated at the later time points, dN decreased steadily (Fig. 2A) from 3.32 x 103 (Ch1535) and 3.14 x 103 (Ch1536) substitutions/site/year to 1.36 x 103 (Ch1535) and 1.6 x 103 (Ch1536) substitutions/site/year. This indicates that fewer amino acid substitutions were acquired and maintained by the virus as the persistent infection progressed, possibly because they no longer confer a replicative advantage due to diminished immune pressure during the chronic phase.
![]() View larger version (29K): [in a new window] |
FIG. 2. (A) Comparison of nonsynonymous and synonymous fixation of mutation rates for Ch1535 and Ch1536 at each sampling time point during the course of the chronic infection. (B) Nonsynonymous substitution rate/synonymous substitution rate ratios for Ch1535 and Ch1536.
|
Envelope and HVR1 substitutions. Approximately 50% of the amino acid substitutions were observed in the envelope region; this proportion remained consistent throughout the study period and was higher than for other regions of the genome. The appearance of substitutions seemed to coincide with increased anti-E1E2 antibody titers, suggesting a response to immune pressure. Figure 3 shows antibody responses to the p90 HVR1 peptide and to recombinant E1E2 antigen. The antibody profiles against both antigens are very similar for each animal, with an increase in antibody to HVR1 correlating with an increase to the whole E1E2 protein. Ch1536 developed higher levels of antibody to the envelope region than Ch1535. This increase coincides with a higher amino acid substitution rate in E2 for Ch1536, 6.61 x 103 and 4.65 x 103 substitutions/site at weeks 130 and 216, respectively, versus 4.41 x 103 and 2.66 x 103 substitutions/site for the same dates in Ch1535. The concomitant higher antibody level and the dN/dS ratios of 6 and 3.5 at these time points for E2 substitutions in Ch1536 suggest that selective immune pressure took place.
![]() View larger version (26K): [in a new window] |
FIG. 3. Comparison of anti-E1E2 and anti-HVR1 antibody levels in Ch1535 (A) and Ch1536 (B) after inoculation with RNA transcripts representing the HCV infectious cDNA. Responses are represented as P/N ratios with a cutoff value of 2.
|
Genome heterogeneity. It is possible that although the predominant circulating virus incorporates a low number of stable mutations over time, the overall population of viruses present in serum shows a greater variation.
Figure 4 shows the results of SSCP analysis of specific regions from the viral genomes of Ch1535 and Ch1536 compared to the same regions amplified from the p90 cDNA clone and from a 6-week p.i. plasma sample from Ch1397, an animal that received H77 plasma. HVR1 and E2 regions amplified from serum samples from Ch1397 show a more complex band pattern than these same regions from the serum samples of Ch1535 and Ch1536, suggesting a greater diversity of these viral sequences in this animal. The NS3 regions show less complexity for all animals, suggesting that these regions are not areas of extensive heterogeneity even in natural isolates. The different band patterns for Ch1535 and Ch1536 are reflective of different viral sequences in these animals (Tables 1 and 2). Analysis of cloned PCR fragments covering the E2 (including HVR1), NS3, and NS5B regions of virus from Ch1535 and Ch1536 and the H77 inoculum (ranging from 335 to 719 bp in length) also confirmed this reduced heterogeneity in the animals that received the infectious clone (data not shown). Between 80 and 100% of the clones contained the mutations observed from PCR sequencing; the clones that did not carry the mutated nucleotide contained the wild-type sequence. These data confirm a low level of sequence heterogeneity after more than 4 years of HCV persistence in the sera of animals that received clonal infectious RNA.
![]() View larger version (114K): [in a new window] |
FIG. 4. SSCP analysis of HCV regions. The HCV regions were amplified from the p90 cDNA clone (lanes 1), serum from Ch1535 at week 216 (lanes 2), serum from Ch1536 at week 216 (lanes 3), and Ch1397 that received H77 plasma (lanes 4). Arrowheads indicate major bands observed in each sample.
|
|
|
|---|
Although the virus isolated from the two animals in this study have only one common mutation, several of the amino acid changes are located in previously identified HCV epitopes recognized by human or chimpanzee class I-restricted CD8+ T cells (reviewed in reference 22). There have been several reports in the literature of viral immune escape at the neutralizing antibody and T-cell levels in both chimpanzees and humans (7, 11, 26, 30), giving weight to the argument that these events play a major role in the establishment and maintenance of persistent infections. The observation that higher antienvelope antibody titers in Ch1536 correlated with a greater number of amino acid substitutions in E2 and a higher dN/dS ratio for this gene product suggests that escape from neutralizing antibody occurs during the chronic phase. The low antienvelope antibody levels during the first 20 to 30 weeks of infection suggest that this mechanism occurs more often in the chronic phase than in the acute phase. Antibody to the surface proteins has been shown to have some neutralizing effect in several studies (5, 9, 10, 21), and this hypothesis of escape can now be addressed using the retrovirus pseudoparticle systems recently developed (1, 12). Interestingly, two amino acid changes within the E2 region (aa 403 and 444) occurring late in infection (at week 130) lie within recently identified monoclonal antibody epitopes shown to have potent neutralizing ability using retrovirus pseudotype particles bearing HCV glycoproteins (12). The substitution at residue 444 was the only shared mutation between both animals.
One pattern that is emerging from this and other studies is that HVR1 does not accumulate mutations rapidly in the chimpanzee even with high titers of HVR1 antibody (23). Previous reports that show rapid change of HVR1 over time have used infections that are initiated from quasispecies populations. We consider these to represent selected changes in the dominant sequence from the preexisting population as opposed to de novo mutations. The changes that were incorporated into HVR1 in the virus isolated from Ch1535 and Ch1536 did not appear to represent antibody escape mutants and, given the mutation rates in these two animals, substantially longer periods of time would be required to generate HVR1 regions that are no longer recognized by preexisting antibody. It has been previously suggested that the rate of HCV evolution in chimpanzees and humans is different (23) with weaker positive selection in chimpanzees. This may account for our observations, although other studies have shown higher rates of HCV substitution in chimpanzees (19, 29). These discrepancies using natural isolates in chimpanzees may be due to subtype differences or different levels of complexity. Infections in humans usually result from higher doses of HCV than the <200 infectious doses normally used in chimpanzee studies, which would therefore result in a greater variant population for subsequent selection.
It has been reported that there is a restricted pattern of amino acid replacement within HVR1 (20, 28), with some positions having little or no variation among all genotypes. The substitutions observed in Ch1535 and Ch1536 follow these previously identified patterns of constraint (28). In this study, we also showed very little sequence heterogeneity of the virus after more than 4 years of infection even in HVR1, which would suggest further that there is no causal relationship between HVR1 variation and HCV persistence.
The comparisons between synonymous and nonsynonymous mutations at different stages of disease have shown that the early mutations arising immediately after control of virus replication during the acute phase are 100% nonsynonymous. This result would suggest that the greatest level of immune pressure occurs during this early phase. Analysis of the viral kinetics in these two animals indicates that significant control of viral replication occurs at weeks 10 to 15; this period coincides with ALT elevations and seroconversion, indicators of induction of the adaptive immune response. As the period of persistence increases, the calculated dN/dS ratio becomes <1 for Ch1535 while remaining higher for Ch1536. This again suggests that there is less immune pressure occurring later during the chronic phase of Ch1535. The difference in dN/dS for Ch1535 and Ch1536 could be due to a slightly higher level of replication in Ch1535 (titers are approximately 0.5 log unit higher), which would allow for the introduction of a greater number of random mutations. Alternatively, the replication rate in Ch1536 may be lower due to a relatively stronger immune response, which would in turn result in more nonsynonymous mutations due to selective pressure. This increase in synonymous mutations in Ch1535 was observed at week 130 p.i. after an increase in viral replication of almost 2 log units (Fig. 1A). This event would argue for a reduced or less effective immune response in this animal that allowed higher levels of viral replication and in turn led to a higher incorporation of random mutations.
Previous studies that have examined molecular evolution of HCV have begun with natural isolates consisting of a quasispecies population. Shifts in sequence due to mutations cannot be differentiated from shifts in virus populations. It is possible in any system that amino acid-changing substitutions can occur by chance and confer no selective advantage or disadvantage to the virus. Using infectious clones, the functional significance of mutations can be more accurately studied. Such analyses are currently being addressed using recently developed in vitro systems such as retrovirus pseudotypes (12) or subgenomic HCV replicons (2).
A steady-state level of virus titers in any persistent infection suggests that virus is being produced and cleared at the same rate, suggesting that some immune response is occurring in these chimpanzees during the chronic phase to suppress the virus, albeit less than that occurring during the acute phase. This may be at the T-cell level and/or due to neutralizing antibody production. The substitution data presented in this study support this hypothesis, and exploitation of these immune responses may be the key to combating persistent infections of HCV.
These studies were supported in part by internal FDA funds, a grant from the National Vaccine Program Office, and a grant from the National Cancer Institute (CA85883).
Present address: Ceregene, Inc., San Diego, CA 92121. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»