Previous Article | Next Article ![]()
Journal of Virology, May 2007, p. 4438-4444, Vol. 81, No. 9
0022-538X/07/$08.00+0 doi:10.1128/JVI.02847-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Medical Research and Education,1 Division of Gastroenterology, Department of Medicine, Taipei Veterans General Hospital, Taipei,2 Institute of Clinical Medicine,3 Faculty of Internal Medicine, National Yang-Ming University School of Medicine, Taipei,4 Division of Gastroenterology, Department of Medicine, Mackay Memorial Hospital, Taipei,5 Department of Life Sciences, National Cheng-Kung University, Tainan, Taiwan6
Received 22 December 2006/ Accepted 5 February 2007
|
|
|---|
|
|
|---|
HDV is a small defective virus with a single-stranded circular RNA of 1.7 kb (14). The antigenomic strand of HDV encodes the only protein, hepatitis delta antigen (HDAg), in two molecular-weight forms. The large form carries an extra 19-amino-acid (aa) extension at the C terminus which plays a key role in the packaging of HDV and suppresses viral replication in a trans-dominant-negative manner, while the small HDAg plays an essential role in transactivating the replication of the HDV RNA (14, 23).
There are two modes of the HDV infection (19). Coinfection results from acute infection with both hepatitis B virus (HBV) and HDV, whereas superinfection results from HDV infection of patients with underlying chronic hepatitis B infection. The great majority of patients with HDV superinfection progress to chronicity, whereas few patients with coinfection do (27). The disease spectrum of HDV infection varies greatly from fulminant hepatitis, rapidly progressive disease, to a subclinical course (7, 9, 20, 26, 27). Persistent replication of HDV associated with continuous hepatic inflammation and elevated alanine aminotransferase (ALT) levels is a characteristic of chronic active hepatitis D, and infection of different HDV genotypes may influence the outcomes (27, 30).
The pathogenesis of HDV infection is still unclear. Cytotoxicity was hypothesized in some studies but was not supported by many others (6, 17). An immune mechanism is believed by most investigators to be responsible for hepatic inflammation and clearance of HDV (10, 16). B-cell and T-cell epitopes were mapped based on in vitro studies (16, 25). Recently, we reported two potential cytotoxic T-cell epitopes in animal and in vitro studies (10, 11). We reported that HDV variants in quasispecies might be changed during clinical course of chronic hepatitis D (29). However, the mechanism and selected sites for the changes of these HDV variants are still obscure.
Only a few statistical studies have ever analyzed selection forces in HDV (2, 12, 29). A positive selection can be inferred if more than synonymous substitutions occur, whereas negative selection removes most nonsynonymous changes. In our previous study, positive selection of the delta antigen sequence was inferred with pairwise comparisons of synonymous and nonsynonymous substitution rates (Ka/Ks) in two patients with chronic hepatitis D (29). The substitution rate was averaged over the entire delta antigen sequence, whereas the positively selected sites were not specified. Anisimova and Yang (2) reported an analysis that allowed selective pressures to vary among codon sites on the small delta antigen (HDAg-S) gene of 33 HDV clones. About 11% sites, mostly located within the immunogenic domains of the HDAg-S gene, were predicted to be shaped by positive selection. Nevertheless, the clones used in that analysis did not display any geographical and temporal relationships. Furthermore, the C-terminal extension of the large-form HDAg has never been studied. Serving as an isoprenylation signal required for packaging of HDV virion, the sequence of the C-terminal extension of large HDAg is conserved within each genotype but highly diverged between different genotypes. The divergence caused difficulties in evaluating the selective pressures on these sites especially since different genotypes were pooled together into analysis. Moreover, the relationship between those presumably selected sites of HDAg and the clinical course of chronic HDV infection remains unanswered.
In the present study, we investigated the viral population in chronic HDV infection. We examined the genetic variation and selective modes of HDV variants from different time points especially those before and after elevation of the serum ALT levels, which are supposed to be selectively shaped by the host immunities, using the likelihood methodology of Yang et al. (33). Models with variable selective pressures among sites were implanted to reveal the positively selected sites.
|
|
|---|
Detection of viral nucleic acids. Because HDAg is the only protein encoded by the HDV and plays an essential role in viral replication and assembly (14). This study focused on changes of sequences of the HDAg-coding region. Reverse transcription PCR using the primers 120 (homologous to a sequence from nucleotides 889 to 912) and 88 (complementary to a sequence from nucleotides 1663 to 1642) according to the HDAg-coding region in the antigenomic strand of HDV RNA was performed as reported previously (4, 28). Pfu DNA polymerase (Promega, Madison, WI) with the highest fidelity was used to minimize the PCR error rate. All of the experiments were independently repeated, and the results were reproducible. Strict procedures were followed to avoid false-positive results (13). Further PCR were done as previously reported (27).
PCR cloning and sequencing.
The amplified PCR products were ligated to a plasmid pCR2 vector (Original TA cloning kit; Invitrogen Corp., Carlsbad, CA). The ligation mixture was used to transform the competent Escherichia coli strain DH5
(Gibco-BRL/Life Technologies, Gaithersburg, MD) as reported elsewhere (5). Multiple colonies from each time point of each case were randomly selected and cultured in the Luria-Bertani medium. Plasmid DNA was extracted and subjected to the dye terminator cycle sequencing reaction (dye terminator cycle sequencing core kit 402117; Perkin-Elmer Cetus Corp., Norwalk, CT). The sequencing products were precipitated with alcohol and analyzed in an ABI 373A sequencer (Perkin-Elmer Cetus Corp.). The sequences were submitted to GenBank; the corresponding GenBank accession numbers are EF187041 to EF187227.
HLA-DNA typing and serological viral markers. HLA phenotypes were classified by using Micro-SSP HLA DNA typing trays (One Lambda, Inc., Canoga Park, CA). Assays for viral hepatitis markers, including hepatitis B surface antigen, hepatitis B e antigen (HBeAg), IgM anti-HBc, anti-HDV, and IgM anti-HAV were tested in each case. All of the markers were detected by radioimmunoassay kits (Ausria II-125, HBeAg-RIA, CORAB-M, anti-Delta, and HAVAB-M; Abbott Laboratories, Chicago, IL). Anti-HCV was measured by a second generation enzyme immunoassay (Abbott Laboratories). Serum ALT levels were measured by a sequential multi-autoanalyzer (Technicon SMAC; Technicon, Tarrytown, NY).
Analysis of quasispecies. Several immunogenic domains over delta antigen have been demonstrated (16, 25). The composition of quasispecies may be altered during the disease course (29). The present study focused on the relationships between immunogenic domains and quasispecies. A dominant variant is defined as a viral strain represented for more than half of the viral population. The datasets of each patient were compiled for the model implantation of positive selection. The predicted amino acid sequences of dominant variants from each data set were processed into the prediction of cytotoxic T ligands.
Tests for positive selection.
Analysis of selective pressures was conducted with a maximum-likelihood approach, which implants different models in evaluating variable
, equivalent to nonsynonymous/synonymous substitution ratios, among codon sites (33). Model M0 assumes a single
for all sites. Model M1, for evaluation of neutral evolution, applies two classes: a portion of conserved sites with
1 = 0 and the neutral sites with
1 = 1. Model M2, denoted as a selection model, adds an additional class of sites with
2 to be estimated from the data. Model M3 uses a discrete distribution to estimate the
value in three classes and provides a sensitive test for positive selection. Likelihood ratio tests were performed to analyze the statistical significance. Sequential to model implantation, the Bayes theorem was applied to identify positively selected sites through maximum-likelihood approach (32). All methods described here were applied on the sequence sets by using CODEML programs of PAML package v3.13 (31).
Prediction of MHC class I-restricted Ligands. In addition to the known immunogenic domains over delta antigen, the domains of major histocompatibility complex (MHC) class I-restricted T-cell epitopes remain to be elucidated. The present study estimated MHC class I-restricted ligands within delta antigen with by NetMHC, which adapts a neural network method to predict binding affinity of peptides being presented by MHC class I molecules (15). Based on the assumption that host immunity interacts mostly with the majority in the viral population, the dominant variants of the first time point from each patient were selected for analysis. Among the high-affinity list of the predicted ligands, the cutoff binding level was determined by sites tested in a HLA-A*0201 transgenic mice model (10). The analyses were done by the NetMHC program v2.1 (15).
|
|
|---|
2 value ranging from 7.582 to 23.958. Model M2 assumes a class of constrained sites with
0 = 0. More than half of the sites in delta antigen gene ranging from 62.7 to 79.7% (mean, 68.1%) were conserved. Model M2 all had a significantly better log-likelihood score in five patients, indicating a greater fitness of positive selection over neutrality. |
View this table: [in a new window] |
TABLE 1. Maximum-likelihood estimates and likelihood ratio tests
|
|
View this table: [in a new window] |
TABLE 2. Positively selected sites predicted in the delta antigen gene
|
CTL epitope prediction. Of the five patients examined in the present study, four had HLA typing of HLA-A*0201, and three of them had a mixture of another allele (two HLA-A*1101 and the other HLA-A*2403). The only non-HLA-A2 patient was HLA-A*1101 (Table 3). Huang et al. (10) identified two epitopes HDV 26-34 and HDV 43-51 in four HLA-A*0201 patients. By analyzing the published peptides' binding affinity along with those predicted by the neural network algorithm, a cutoff level of 105 M was determined. Of the potential epitopes, only those predicted binding affinity stronger than 105 M were included. Table 4 shows the potential cytotoxic T ligand (CTL) epitopes of the four HLA-A2 patients.
|
View this table: [in a new window] |
TABLE 3. Characteristics of patients and HDV quasispecies
|
|
View this table: [in a new window] |
TABLE 4. Potential CTL epitopes on HDAg of four patients with HLA-A*0201
|
Immunogenic epitopes associated with selection. Table 5 shows the relationships between selected sites and immunogenic domains of the HDAg. Based on a posterior probability of >90%, positive selection was predicted in three (60%) patients in the helper T-cell epitope from aa 26 to 41, three (60%) patients in the T-cell epitope from aa 66 to 81, and four (80%) patients in the T-cell epitope from aa 106 to 201. Positive selection was predicted in three patients (60%) in a B-cell epitope spanning from aa 2 to 7 and in five patients (100%) in a B-cell epitope spanning from aa 174 to 195. Selections of sites within HDAg in the four HLA-A*0201 patients and one HLA-A*1101 patient were correlated to the potential CTL epitopes. Three (60%) patients had sites under positive selection in the cytotoxic T-cell epitope from aa 43 to 51, and four (80%) patients had sites under positive selection in an epitope from aa 114 to 122. The C-terminal extension of the large HDAg gene had only one potential epitope of aa 198 to 206 predicted within genotype I. Four of five (80%) patients did not have any evidence of positive selection in this epitope. Only patient E infected with HDV genotype IV had residue changes at aa 203. This residue was not associated with any other known immunogenic domain and may be related to the packaging of the delta antigen during replication.
|
View this table: [in a new window] |
TABLE 5. Patients with selected sites located within immunogenic epitopes of HDAg
|
![]() View larger version (18K): [in a new window] |
FIG. 1. Consequence of selection of novel dominant HDV variants after flare-up of ALT levels in patient A. (A) Clinical course. The time points of blood sampling are marked (1435 and 5132). (B) Partial amino acid sequences of delta antigen obtained from different time points. The initially dominant HDV variant is shown at top of the panel. Dots indicate identical amino acid residues. The labeled CD4+ T is a T-cell epitope reported previously (16). The labeled CD8+ CTL are predicted to be MHC-I-restricted epitopes.
|
![]() View larger version (16K): [in a new window] |
FIG. 2. Consequence of selection of novel dominant HDV variants after flare-up of ALT levels in patient D. (A) Clinical course. The time points of blood sampling are marked (1025 and 3038). (B) Partial amino acid sequences of delta antigen obtained from different time points. The initially dominant HDV variant is shown at top of the panel. Dots indicate identical amino acid residues. The labeled B are B-cell epitopes were reported previously (25).
|
|
|
|---|
Furthermore, several studies have suggested the importance of cellular immunity in HDV infection. High HDV viral loads were detected in human immunodeficiency virus-infected chronic hepatitis D patients whose circulating T cells numbers were suppressed (21), and the activity of HDV-induced liver disease correlated with CD4+ T-cell response to HDAg (16). Nisini et al. (16) also reported four epitopes from HDAg-specific CD4+ T-cell clones. Three of these four epitopes are positively selected in at least 60% of our patients. The epitope from aa 106 to 121 has an occurrence of positive selection in 80% of our patients. This MHC class II-restricted T-cell epitope is reported to generate via extracellular processing of HDAg by serum protease (1). This peptide may function both as a "blocking" peptide, inhibiting MHC binding and the presentation of "endogenous" HDAg epitopes, and as a kind of decoy inducing the exhaustion of HDAg106-121-specific T cells and consequently facilitating the persistence of virus-infected cells. The alterations in this epitope may confer a greater fitness to the viral population for escaping attacks from the host's immune system.
Importantly, CD8+ T lymphocyte-related cytotoxic immunity has a pivotal role in the elimination of intracellular pathogens. In HBV and hepatitis C virus infections, multispecific CD8+ T-cell responses to HBV and hepatitis C virus are closely associated with viral clearance (8, 24). Our previous study revealed that aa 26 to 34 and aa 43 to 51 of HDAg are novel HLA-A*0201-restricted CTL epitopes in genotype I HDV (10). Using a neural network method, our analysis suggested more potential CTL epitopes in the genotype I. Furthermore, novel epitopes were proposed in genotypes II and IV. All genotypes of HLA-A*0201 HDV- infected patients had three consistent epitopes: aa 43 to 51, aa 50 to 58, and aa 114 to 122. In our previous studies, the CD8 cells showing positive responses to HDV 43-51 were also detectable in two HLA-A*0201 chronic delta hepatitis patients (patients A and B), whose serum HDV RNA levels gradually decreased and progressed into disease remission during follow-up (10). In the present study, the majorities of the HDV population in our patients were shaped by positive selection related to aa 43 to 51 and aa 114 to 122. Therefore, cytotoxic immune response is closely associated with viral clearance, as shown by animal and human studies. These potential epitopes can be confirmed in further studies and could lead to the development of therapies.
It is worth noting that the predicted ligand 29-37 (ELERDLRKV) in patient A is located at the same loci of HDV 29-37 (DLERDLRKI) tested by Huang et al. (10). We tested these two peptides using the NetMHC algorithm. The predicted binding level of Huang et al.'s HDV 29-37 was 1.8096 x 105 M, which was higher than the 0.7419 x 105 M of ligand 29-37 or our designated 105 M. Thus, the peptide in the Huang et al. study was suggested to have a weaker binding affinity. The predicted ligand 29-37 (ELERDLRKV) may still be a potential ligand.
A recent study inspected the nucleotide changes after 1 year of in vitro replication in cell culture (3). These authors concluded that the observed nucleotide changes were essentially neutral. By supplying a consistent HDAg-S from the integrated cDNA sequence, the nucleotide changes were diverse and did not compromise replication competence. However, HDV may encounter a much more complicated environment in the human host than in cell cultures. A heterogeneous HDV population may supply the HDAg-S. Besides, the HDAg-S may be under intensive selective pressures from the host immunities. A chronic HDV infection could lead to a more complicated relationship between host and viral genomes. Our analysis revealed that, on the contrary, the nucleotide changes in chronic HDV infection were positively selected.
Positive selection is commonplace in chronic hepatitis D during intense selection forces from host immunity. Positively selected sites identified in a longitudinal study differ from those observed in a cross-sectional one. The majority of sites under diversifying selection associate with known immunogenic domains. More cytotoxic T-cell epitopes in HDAg were predicted. aa 43 to 51, 50 to 58, and 114 to 122 of HDAg are consensus epitopes across different HDV genotypes. The discovery of immune epitopes associated with selection and clearance of HDV variants may be of value for the development of novel treatments or vaccines for HDV superinfection.
Published ahead of print on 14 February 2007. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»