Previous Article | Next Article ![]()
Journal of Virology, July 2006, p. 6469-6477, Vol. 80, No. 13
0022-538X/06/$08.00+0 doi:10.1128/JVI.00245-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Microbiology and Immunology, Georgetown University, Washington, D.C. 20057,1 Division of Molecular Virology and Immunology, Georgetown University Medical Center, Rockville, Maryland 20850,2 Department of Clinical Sciences, New York State College of Veterinary Medicine, Cornell University, Ithaca, New York3
Received 2 February 2006/ Accepted 12 April 2006
|
|
|---|
|
|
|---|
There are two modes of HDV infection: superinfection of an individual with chronic HBV infection and coinfection with both viruses of an individual who has not been exposed previously to HBV. While the latter type of exposure typically leads to recovery from both viral infections, the former results in chronic HDV infection with a frequency of 70% to 80% (13). The determinants of the outcome of HDV superinfection are not known. Previous studies of viral genetic changes during infection have analyzed HDV RNA obtained either from the sera of chronically infected patients during the course of HDV infection (7, 14, 19) or from the liver of a woodchuck infected with virus that had been passaged several times (21). These studies did not address whether genetic changes in the virus could be important for the establishment of chronic infection.
In this study, we analyze genetic changes occurring in the HDV genome following the infection of woodchucks with an inoculum derived from a molecular clone of HDV genotype I. We found that a limited number of sequence changes occur during chronic infection. Some of these modifications alter the consensus sequence of hepatitis delta antigen (HDAg) and could indicate a mechanism for avoidance of host immune responses.
|
|
|---|
Transfection of a woodchuck with HDV cDNA. All animal experiments for the current study were conducted under protocols reviewed and approved by the Cornell University Institutional Animal Care and Use Committee. Woodchuck 4928, a chronic WHV carrier, was transfected with 50 µg of the HDV construct pGDLx1.2 in a total volume of 0.5 ml of phosphate-buffered saline by direct injection at several separate sites in the surgically exposed liver. Sera were obtained weekly and analyzed by reverse transcription (RT) and semiquantitative PCR for HDV RNA, as described previously (22, 31). Highly positive sera from weeks 7 to 15 were combined to create the woodchuck hepatitis delta virus type 1 (WHDV-1) pool. This serum pool contained 1 x 109 genome equivalents of HDV RNA per ml, as determined by blot hybridization.
Infection of WHV woodchucks with woodchuck-derived HDV. WHV carrier woodchucks were produced experimentally by administering the standardized cWHV8P1 inoculum (8) to neonatal woodchucks born to WHV-negative females. The WHV carrier status of woodchucks was verified by detection of persistent WHV viremia and antigenemia (8). Chronic WHV carrier woodchucks were infected at 18 months of age with woodchuck-derived HDV by intravenous administration of a 0.5-ml volume of a 1:5 dilution of the WHDV-1 pool. Serum samples were obtained weekly and analyzed for HDV RNA by reverse transcription and semiquantitative PCR (22, 31).
Sequence analysis of HDV RNA. Serum samples were subjected to reverse transcription, followed by PCR with one of two sets of primers (Table 1). Primer pair 5414 and 6657 amplifies a 775-bp product from positions 886 to 1660 (numbering according to GenBank accession no. L22066) that encompasses the entire HDAg coding region; primer pair 3415A and 6465 amplifies a 746-bp product that includes the amino-terminal part of the HDAg coding region and the noncoding portion of the genome that is base paired with it in the unbranched rod structure (Fig. 1). PCR amplification was performed using Pfu polymerase (Stratagene, La Jolla, CA) for 35 cycles. PCR products were either sequenced directly (MWG Biotech, High Point, NC) or cloned into the plasmid pZero (Invitrogen, Carlsbad, CA). When the products were cloned, two separate RT-PCRs were performed to minimize the potential of analyzing sibling clones. Changes at the amber/W site were not included in the sequence comparisons because this site is modified as part of the normal HDV replication cycle (25).
|
View this table: [in a new window] |
TABLE 1. Oligonucleotide primers used for RT-PCR amplification
|
![]() View larger version (6K): [in a new window] |
FIG. 1. Schematic of the HDV RNA genome showing PCR primers and products analyzed. The flattened circular line represents the HDV genome, the open rectangle the HDAg gene, the filled triangle the location of the genomic ribozyme, and the open triangle the location of the antigenomic ribozyme. Arrows indicate the locations and orientations of primers, and the dashed lines show the PCR products.
|
|
|
|---|
Analysis of HDV genetic stability in the transfected woodchuck. To analyze the genetic diversity of the WHDV-1 pool, 11 clones of the HDAg coding region were obtained by RT-PCR, using primers 6657 and 5415 (Fig. 1), and then sequenced. Clones were obtained from two independent RT-PCR amplifications to minimize the chances of obtaining sibling clones. Compared with the sequence of the transfected cDNA, a total of 18 sequence changes occurred in the 695-nucleotide (nt) region sequenced, or 0.24% (Table 2). Remarkably, 11 of the 18 changes occurred at position 1008 in the genome, which was modified from A to G in all 11 clones. The change at this position alters the predicted HDAg sequence at amino acid 198 from isoleucine to threonine. This amino acid change has been described previously after multiple passages of HDV in woodchucks (9). The remaining seven nucleic acid sequence changes were observed at different positions in different clones. In similar analyses of cDNA clones obtained by RT-PCR, we observed that the error rate due to Pfu polymerase was less than 0.1% (26). Thus, while the modification at position 1008 is clearly not likely to be due to misincorporation during RT-PCR amplification, it is impossible to determine to what extent the other seven changes occurred during the course of infection or during reverse transcription and PCR amplification.
|
View this table: [in a new window] |
TABLE 2. Number and types of sequence changes observed in cDNA clones
|
![]() View larger version (13K): [in a new window] |
FIG. 2. Time course of HDV RNA detected in the sera of animals inoculated with WHDV-1. Sera were obtained weekly and analyzed for HDV RNA by reverse transcription and semiquantitative PCR (22). Quantitation was based on a comparison of PCR product yields with standard dilutions of an HDV source of known concentration.
|
In order to examine how rapidly sequence divergence occurred, six HDV cDNA clones were obtained at 8 weeks postinoculation from two of the animals (animals 4543 and 4569) that eventually became chronically infected. Analysis of the sequences of these 12 clones indicated a much lower degree of nucleotide substitution than at 73 weeks postinoculation (Table 2). Subtracting the frequency of sequence changes observed at 8 weeks from that at 73 weeks gives an overall rate of 4 x 103 substitutions per site per year in the HDAg coding region, similar to results previously reported for infected patients (7, 14, 19).
Types of sequence changes. In the clones of HDV RNA from week 73, adenosine-to-guanosine changes occurred with the highest frequency (1.4% of adenosines modified), followed by uridine to cytosine (1% of uridines modified) (Table 2). Both of these modifications are consistent with editing by host RNA adenosine deaminase (ADAR) activity. ADAR deaminates adenosines in RNA to inosines, which are subsequently transcribed as guanosines; adenosines that are edited on the genome will appear as A-to-G changes, while those on the antigenome will appear as U-to-C changes on the genome. A-to-G and U-to-C transitions accounted for 67% (129/191) of the sequence substitutions observed in the clones from week 73. These were also the predominant sequence changes observed in the clones obtained during the peak of viremia (Table 2). The high relative frequency of these changes may indicate that ADAR is responsible for much of the genetic divergence that arises during the course of acute and chronic HDV infection.
Previous analyses of HDV sequence divergence among different isolates have indicated that the noncoding region of the genome is, in general, more variable than the HDAg coding region (22, 28). Moreover, studies of the evolution rate of HDV during the course of chronic infection have shown that sequence changes accumulate more rapidly in the noncoding region (7, 19). To compare the extent of sequence changes in the coding and noncoding regions, the region 1267 to 1681/1 to 308 was amplified with primers 3415A and 6465 (Fig. 1), using serum obtained at 73 weeks postinoculation from animal 4569. Analysis of six cDNA clones of this region indicated a rate of sequence divergence of 1 x 102 changes per nucleotide per year (Table 2). This rate was considerably higher than that observed in the coding region over the same period of time. Notably, sequence modifications in the noncoding region included insertions and deletions (mostly single nucleotide), and these accounted for the difference in the rates of sequence divergence between the coding and noncoding regions analyzed. No insertions or deletions were found in any of the clones from the HDAg coding region (Table 2).
Changes in the HDAg consensus sequence. Comparison of the HDAg sequence of the inoculum with those of the clones obtained at 73 weeks postinoculation from the five chronically infected animals indicated that the distribution of the modifications was not random (Fig. 3). In all five animals, at least one change occurred in the consensus HDAg coding sequence (defined here as the sequence in six or more of the eight clones) (Fig. 3 and Table 3). Most of these consensus changes were in the N-terminal region of the protein. In animal 4543, the consensus was altered at two positions; in the other four animals, there was a single modification in the consensus sequence (Table 3). The consensus changes observed in animals 4543 and 4547 were unique to those two animals. At least six of eight clones from woodchucks 4553, 4559, and 4569 exhibited a change in the consensus at nucleotide position 1480 that changed amino acid 41 from isoleucine to valine.
![]() View larger version (11K): [in a new window] |
FIG. 3. Histogram of HDAg amino acid sequence changes in HDV cDNA clones obtained from chronically infected animals at 73 weeks postinfection. Sequence changes from the inoculum are indicated by vertical bars; bar heights indicate the number of clones changed at different positions. The dashed horizontal line indicates one clone.
|
|
View this table: [in a new window] |
TABLE 3. Number of specific sequence changes observed in eight cDNA clones obtained at 73 weeks
|
V); one clone was also modified at amino acid 51 (L
M). Finally, the two clones from animal 4559 that were not modified at amino acid 41 exhibited I
T changes at position 37. As can be seen in Fig. 4, these modifications were unique to particular animals. Thus, for each of the chronically infected animals, the HDAg sequence was modified in all clones, either at a single position or at two or more positions within a 15-amino-acid segment (Fig. 4). Furthermore, while there were similarities in the sequence changes observed in some of the woodchucks, the overall changes in each animal were unique.
![]() View larger version (7K): [in a new window] |
FIG. 4. Comparison of HDAg sequences in cDNA clones obtained from chronically infected animals at 73 weeks postinfection. The consensus sequence of the inoculum is shown at the top of each panel. Sequence changes are shown by single-letter amino acid codes; dots indicate no change.
|
Kinetics of sequence changes in animals that either recovered from HDV infection or became chronically infected. To examine the dynamics of sequence changes, RNA was isolated from sera obtained at weeks 8, 17, and 27 postinoculation. Products of RT-PCR amplifications were sequenced directly to determine the consensus sequences at these times. In the animals that eventually became chronically infected, no changes in the majority sequence of HDAg were observed at 8 weeks postinoculation and there was no evidence of minor populations (ca. 10 to 20%) with sequence differences at positions that changed at later times. By 27 weeks postinoculation, the majority sequence had changed in just one of the animals that eventually became chronically infected (animal 4543) (Table 4). In two others (animals 4553 and 4569), there was evidence of some sequence heterogeneity, but the abundance of the new sequence was 20% or less (Table 4).
|
View this table: [in a new window] |
TABLE 4. Changes at week 27
|
![]() View larger version (28K): [in a new window] |
FIG. 5. Electropherogram of HDV RT-PCR products obtained at weeks 8, 17, and 27 postinfection from animal 4543. Black, G; green, A; blue, C; red, U. Positions 1554 and 1556 are shown by dashed lines. Note that because of the primer used for sequencing, the sequence shown is the reverse complement of the genome.
|
The HDAg amino acid sequence changes found at week 27 in the woodchucks that recovered from infection differed somewhat from the changes that occurred by week 73 in the chronically infected animals (Fig. 6). Only six different modifications were observed that changed the HDAg consensus sequence, and three of these (E15D, I16V, and I41V) were found in two or more animals. An I16V substitution occurred in two of the animals that recovered from infection, and I41V was common among three of the five chronically infected animals. On the other hand, the E15D modification was observed in both animal 4548, which recovered from HDV infection, and animal 4543, which became chronically infected. The E15D change was also detected at the 50% level in animal 4543 at 27 weeks postinfection (Table 4). Thus, it is not clear to what extent the apparent segregation of sequence changes between animals that became chronically infected and those that recovered from infection is due to the outcome of infection, to the timing of the sequence change, or to other factors.
![]() View larger version (5K): [in a new window] |
FIG. 6. Consensus HDAg sequence changes in animals that either recovered from HDV infection or became chronically infected. Sequences for animals whose infection resolved were obtained by direct sequencing of RT-PCR products obtained from HDV RNA isolated at 27 weeks postinfection (Table 4); two consensus changes were observed for animal 4548 and none for animal 4556. Sequences for animals that became chronically infected were obtained by sequencing of eight clones of RT-PCR products obtained from HDV RNA isolated at 73 weeks postinfection (Fig. 3). Locations of sequence modifications are indicated in their approximate positions along the HDAg coding sequence.
|
|
|
|---|
Editing on the HDV RNA by ADAR1 is highly specific for the amber/W site. During a 2-week course of replication in cultured cells, non-amber/W sites were edited at very low levels (26). However, editing at non-amber/W sites is readily detectable in vitro, with some sites modified more efficiently than others (26). In the clones of the week 73 PCR products, most of the sites at which substitutions consistent with ADAR editing occurred were found to be modified in only a single clone (data not shown), but some, including those at which four of the six HDAg consensus sequence changes occurred, were found to be modified in many clones. The abundance of changes consistent with ADAR editing in the viral population during the course of infection is likely to have been determined by both the susceptibility of different sites to ADAR editing and any long-term selective growth advantage conferred by the sequence changes.
Remarkably, a limited number of modifications in the consensus sequence of HDAg occurred in seven of eight animals infected with the molecularly derived inoculum (Fig. 4 and Table 4). Moreover, a consensus change was observed in virus isolated from the woodchuck transfected with the HDV cDNA clone, even though the overall rate of sequence divergence was low in clones isolated from this animal. A previous study of the genetic stability of HDV in woodchucks (21) reported no consensus sequence changes during serial passages over a total length of 253 days. However, only the C-terminal half of the HDAg coding region was analyzed in that study and some changes could have been missed. Indeed, most of the consensus changes that we observed were in the N-terminal half of the protein (Fig. 3).
It is difficult to distinguish between the effects of positive and negative selective pressures on changes in the virus sequence. Most likely, both negative selective pressure due to host immune responses and positive selective advantages due to increased replication efficiency play a role. Indeed, there is evidence of both processes in our results. The dominant A-to-G change at position 1008 (amino acid 198) may be the result of adaptation of the virus to the woodchuck from the native human host. This change, which was observed in another passage of HDV to woodchucks (9), was observed in all 11 clones obtained during the acute phase of viremia in animal 4928 following transfection and was found in all sequences obtained from woodchucks that were subsequently infected. It thus appears that this change confers a significant replication advantage in woodchucks. The A1008G sequence modification replaces I198 with T in large HDAg (HDAg-L); I198 is among the additional amino acids present in HDAg-L that are required for particle formation with the hepadnavirus surface antigen. O'Malley and Lazinski found that mutation of I198 to T was fully tolerated for packaging of the mutant HDAg-L by HBV envelope proteins (24). Perhaps the I198-to-T change produces a more favorable interaction between HDAg-L and the WHV surface antigen. That this sequence change was not observed in a serial passage of HDV in WHV-infected woodchucks (21) might be explained by the use of different WHV strains.
The results obtained for the woodchucks infected with the molecularly derived inoculum are consistent with a strong influence of host immune responses on the selection of viral genetic changes. The woodchucks used in this study were not inbred; thus, somewhat diverse immune responses and subsequent evasive genetic changes in the viral sequence could be expected. Particular changes that affected the amino acid sequence were clearly segregated among different animals; however, two other frequent changes, at positions 1358 and 1617, that did not affect HDAg sequence appeared to be distributed randomly among clones from different animals (Table 3).
While consensus changes at some HDAg amino acid positions in woodchucks 4543 and 4569 were due to modifications in all clones, the amino acid at the position changed in the consensus sequence was not modified in all clones in woodchucks 4547, 4553, and 4559 (Fig. 4). However, for these three woodchucks, all clones not containing modifications at the position changed in the consensus sequence were modified at one or more sites within 15 amino acids of this position (Fig. 4). Such changes, whether restricted to a single amino acid or clustered within a 15-amino-acid region, could readily affect the immune response by affecting the association of specific HDAg-derived peptide fragments with major histocompatibility complex class I or class II molecules (23). The evolution of sequence changes in animal 4543 is consistent with this interpretation. Initial changes in the HDAg consensus sequence at position 1554 in animal 4543 at 17 weeks postinoculation may indicate a virus-evasive reaction to immune pressure that either did not fully avoid the host immune response or imposed a replication disadvantage. Thus, the subsequent reversion at this position combined with the change at position 1556 would be consistent with either a more successful avoidance of host immune responses (in the former case) or a better combination of immune avoidance and viral replicative capacity (in the latter case).
HDV sequence changes observed at 73 weeks postinoculation in the chronically infected animals predominated at a limited number of sites, most of which are in the N-terminal segment of the HDAg coding region (Fig. 3 and 4; Tables 3 and 4). Previous analyses of humoral antibody responses to HDV infection identified an immunodominant region between amino acids 52 to 93 (2); relatively few sequence changes were observed in this region (Fig. 3). However, high percentages of HDV-infected woodchucks also exhibited humoral antibody responses to peptides from the amino-terminal region of HDAg (amino acids 1 to 18 and 36 to 52) (2), and sequences in this region were observed to be modified frequently in our study (Fig. 3). The significance of this association is not clear because humoral antibody to HDAg is not likely to be protective; however, the correlation suggests the possibility of immunological cross talk between the development of humoral responses and other immune responses that are involved in the clearance of virus infection, i.e., cytotoxic T-lymphocyte (CTL) responses. Dissection of the immune response to HDV infection in future studies may help to directly determine the extent to which HDAg sequence changes are the result of immunologic pressure.
The determinants of the outcome of HDV superinfection are not known. Analyses of the viral genetics of hepatitis C virus have indicated that progression to chronic infection is correlated with sequence changes and increased sequence diversity in the hypervariable region of the envelope protein (11). Previous studies of genetic changes that occur during the course of HDV infection have analyzed sequence modifications that occur over time in HDV RNA isolated either from the sera of patients with established chronic infection (7, 14, 19) or from the liver of a woodchuck infected with HDV that had been passaged in woodchucks several times previously (21). However, none of these studies specifically addressed the role of viral genetic changes in the establishment of chronic infection.
In our study, the timing and degree of sequence changes appeared to be correlated with the outcome of infection. Five of five chronically infected woodchucks exhibited consensus sequence changes at 73 weeks postinfection, but only one of the five exhibited sequence substitutions at earlier times (Fig. 4 and Table 4). In contrast, direct sequencing of RT-PCR products indicated that majority sequence changes were apparent by 27 weeks in two of three woodchucks that recovered from HDV infection, and in the third, there was evidence of multiple sequence changes that had not yet dominated the population (Table. 4). Moreover, the number of positions changed was generally greater in the animals that recovered from HDV infection than in those that developed chronic infection (Fig. 4 and Table 4). Although the correlation between sequence changes and the outcome of infection is based on a small number of animals, we note that this correlation is also supported by analysis of HDV infection in woodchuck 4928, which was transfected with the cDNA clone and eventually recovered from HDV infection; like the three animals that recovered from infection initiated by intravenous inoculation, this animal also produced virus in which a consensus change was detected early, between 8 and 15 weeks posttransfection.
One explanation of the differences observed between woodchucks that recovered and those that became chronically infected is that the virus sequence changed earlier in the woodchucks that recovered because of a more vigorous host immune response. In the animals that recovered, either the changes in HDAg sequence were insufficient to avoid immune pressure or the immune system was able to recognize additional epitopes that facilitated viral clearance. In the animals that became chronically infected, the delayed appearance and limited nature of the changes in consensus sequence may indicate that the immune response was weak or delayed and may have been limited to only one or two epitopes. Thus, chronic HDV infection in woodchucks may result from a delayed and weak immune response that is limited to a small number of HDAg epitopes. Indeed, analyses of proliferative peripheral blood mononuclear cell responses to HDAg in infected humans and woodchucks have indicated weak responses that are typically limited in most individuals to just one to three epitopes (10, 12, 23).
If the explanation is correct that the observed genetic changes are determined by the immune response, then stimulation of a potent immune response, which is inferred by the sequence changes observed at 27 weeks postinfection in the animals that recovered, could be achieved by vaccination and such a response could result in viral clearance. Most likely, the alteration of the course of infection by some vaccines (1, 10, 12, 15, 16) is due to the ability to elicit appropriate cell-mediated immune responses. However, no correlation has been observed between vaccine efficacy and anti-HD T-cell proliferative responses (12). Perhaps, CTL activity is the more important contributor to clearance of HDV-infected cells. Future vaccine studies may benefit from measurement of CTL activity induced by HDV infection and/or vaccines as well as analysis of viral genetics. The ability of the virus sequence to change at a limited number of sites could be a critical factor in the ability to develop a successful vaccine against HDV infection.
We gratefully acknowledge the expert assistance of Betty Baldwin, Lou Ann Graham, Richard Moore, and Chris Bellezza (Cornell University) and of Frances Wells and James Nupp (Georgetown University).
|
|
|---|
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»