Document not found! Please try again

Positive Selection of Hepatitis Delta Antigen in ... - Journal of Virology

0 downloads 0 Views 613KB Size Report
Medicine, Mackay Memorial Hospital, Taipei,5 and Department of Life .... fied by using Micro-SSP HLA DNA typing trays (One Lambda, Inc., Canoga Park,. CA). ..... 2000. Hepatitis delta virus: the molecular basis of laboratory diagnosis. Crit.
JOURNAL OF VIROLOGY, May 2007, p. 4438–4444 0022-538X/07/$08.00⫹0 doi:10.1128/JVI.02847-06 Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Vol. 81, No. 9

Positive Selection of Hepatitis Delta Antigen in Chronic Hepatitis D Patients䌤 Shen-Yung Wang,1,5 Jaw-Ching Wu,1,3* Tzen-Yuh Chiang,6 Yi-Hsiang Huang,2,3 Chien-Wei Su,2,4 and I-Jane Sheen2 Department of Medical Research and Education1 and Division of Gastroenterology, Department of Medicine,2 Taipei Veterans General Hospital, Taipei, Institute of Clinical Medicine3 and Faculty of Internal Medicine,4 National Yang-Ming University School of Medicine, Taipei, Division of Gastroenterology, Department of Medicine, Mackay Memorial Hospital, Taipei,5 and Department of Life Sciences, National Cheng-Kung University, Tainan,6 Taiwan Received 22 December 2006/Accepted 5 February 2007

Liver disease may become ameliorated in some patients with chronic hepatitis D virus (HDV) infection. We present here a study based on longitudinal sampling to investigate the viral dynamics in chronic HDV infection. We examined the HDV variants from different time points, especially those before and after the elevation of serum aminotransferase levels. The datasets from each patient were tested for positive selection by using maximum-likelihood methods with heterogeneous selective pressures along the nucleotide sequence. An average of 4.9%, ranging from 3.1 to 6.8%, of the entire delta antigen sites was regulated by a diversifying selection. Most of the positively selected sites were associated with immunogenic domains. Likelihood ratio tests revealed a significant fitness of positive selection over neutrality of the hepatitis delta antigen gene in all patients. We further adapted a neural network method to predict potential cytotoxic T ligand epitopes. Among the HLA-A*0201 cytotoxic T ligand epitopes, three consistent epitopes across all three genotypes were identified: amino acids (aa) 43 to 51, 50 to 58, and 114 to 122. Three patients (60%) had sites evolving under positive selection in the epitope from aa 43 to 51, and four patients (80%) had sites evolving under positive selection in the epitope from aa 114 to 122. The discovery of immunogenic epitopes, especially cytotoxic-T-lymphocyte ligands, associated with chronic HDV infection may be crucial for further development of novel treatments or designs in vaccine for HDV superinfection. 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

RNA virus exists in a quasispecies nature, which represents microheterogeneity among viral populations. The mutations occur associated with the lack of proofreading activity of the RNA polymerase during replication. The genomic changes may render the viral genome able to replicate more efficiently or to escape the host immune responses. Positive selection usually preserves such favorable mutations. In contrast, negative selection removes deleterious substitutions from the viral populations, thereby ensuring survival of the viral species. Through generations, viral populations evolve under various selective forces at different regions and sites that display different functional constraints. The hepatitis D virus (HDV), with its genome simplicity, provides an ideal model for elucidating this viral evolution. 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).

* Corresponding author. Mailing address: Department of Medical Research and Education, Taipei Veterans General Hospital, 201, Shih-Pai Road, Sec. 2, Taipei 11217, Taiwan. Phone: 886-2-28712121, ext. 3218. Fax: 886-2-28745074. E-mail: [email protected]. 䌤 Published ahead of print on 14 February 2007. 4438

VOL. 81, 2007

SELECTION OF HDAg IN PATIENTS WITH HDV INFECTION

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. MATERIALS AND METHODS Patients. Serum samples were collected at several time points from each of five patients with chronic hepatitis D. These patients went into remission status during longitudinal follow-up. Chronic hepatitis was diagnosed by an elevated ALT level lasting for more than 6 months. Biochemical remission was defined as normal ALT levels in at least three separate samples lasting for at least 12 months without histological or clinical evidence of cirrhosis or hepatocellular carcinoma (HCC) (22, 27). Of the patients studied, two were infected with HDV genotype I, another two were infected with genotype IV (previously named IIb) (18), while another was infected with genotype II (previously named IIa). The patients were all positive for serum hepatitis B surface antigen and antibody to HDV antigen (anti-HDV) for more than 6 months and negative for immunoglobulin M (IgM) antibody to hepatitis B core antigen (anti-HBc), antibody to hepatitis C virus (anti-HCV), and IgM antibody to hepatitis A virus (anti-HAV). The HDV variants before and after elevation of the serum ALT levels were cloned from the sera of these patients. 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

4439

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).

RESULTS Test of positive selection. The results of likelihood ratio tests and the maximum-likelihood estimates of different models’ implantation are listed in Table 1. Model M2 suggested that 3.1 to 6.8% (mean 4.9%) of the delta antigen sites evolved under a positive selection with a ␻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. Anisimova and Yang (2) reported that model M2, by adding an additional class of sites to neutral model M1, may identify purifying rather than positive selected sites (2). In the present study, the optimal peaks of model M2 all favored the hypothesis of positive selection. The likelihood ratio tests also showed strong evidences of positive selection in all five patients. Positively selected sites predicted by the Bayes theorem under model M3 with a possibility of ⬎90% are presented in Table 2.

4440

WANG ET AL.

J. VIROL. TABLE 1. Maximum-likelihood estimates and likelihood ratio tests

Patient

A

B

C

D

E

l0

Null model

M0 ␻ ⫽ 0.661

⫺1,634.227

M1 ␻0 ⫽ 0, ␻1 ⫽ 1 P0 ⫽ 0.849, P1 ⫽ 0.151

⫺1,561.448

M0 ␻ ⫽ 0.498

⫺2,676.195

M1 ␻0 ⫽ 0, ␻1 ⫽ 1 P0 ⫽ 0.712, P1 ⫽ 0.288

⫺2,565.263

M0 ␻ ⫽ 1.572

⫺1,348.454

M1 ␻0 ⫽ 0, ␻1 ⫽ 1 P0 ⫽ 0.739, P1 ⫽ 0.261

⫺1,323.364

M0 ␻ ⫽ 0.675

⫺1,156.107

M1 ␻0 ⫽ 0, ␻1 ⫽ 1 P0 ⫽ 0.629, P1 ⫽ 0.372

⫺1,151.462

M0 ␻ ⫽ 0.575

⫺1,945.590

M1 ␻0 ⫽ 0, ␻1 ⫽ 1 P0 ⫽ 0.665, P1 ⫽ 0.335

⫺1,907.835

l1

PLRT

M3 ␻0 ⫽ 0.097, ␻1 ⫽ 9.846, ␻2 ⫽ 31.399 P0 ⫽ 0.942, P1 ⫽ 0.036, P2 ⫽ 0.022 M2 ␻0 ⫽ 0, ␻1 ⫽ 1, ␻2 ⫽ 23.958 P0 ⫽ 0.797, P1 ⫽ 0.156, P2 ⫽ 0.048

⫺1,489.977

⬍0.001

⫺1,493.223

⬍0.001

M3 ␻0 ⫽ 0.051, ␻1 ⫽ 0.550, ␻2 ⫽ 9.590 P0 ⫽ 0.696, P1 ⫽ 0.242, P2 ⫽ 0.061 M2 ␻0 ⫽ 0, ␻1 ⫽ 1, ␻2 ⫽ 12.163 P0 ⫽ 0.668, P1 ⫽ 0.279, P2 ⫽ 0.053

⫺2,441.145

⬍0.001

⫺2,445.689

⬍0.001

M3 ␻0 ⫽ 0.00001, ␻1 ⫽ 2.418, ␻2 ⫽ 29.955 P0 ⫽ 0.771, P1 ⫽ 0.169, P2 ⫽ 0.060 M2 ␻0 ⫽ 0, ␻1 ⫽ 1, ␻2 ⫽ 24.871 P0 ⫽ 0.647, P1 ⫽ 0.285, P2 ⫽ 0.068

⫺1,270.337

⬍0.001

⫺1,270.970

⬍0.001

M3 ␻0 ⫽ 0.228, ␻1 ⫽ 0.228, ␻2 ⫽ 7.269 P0 ⫽ 0.0001, P1 ⫽ 0.943, P2 ⫽ 0.057 M2 ␻0 ⫽ 0, ␻1 ⫽ 1, ␻2 ⫽ 7.582 P0 ⫽ 0.667, P1 ⫽ 0.288, P2 ⫽ 0.046

⫺1,147.487

0.002

⫺1,147.751

0.024

M3 ␻0 ⫽ 0.068, ␻1 ⫽ 1.703, ␻2 ⫽ 12.036 P0 ⫽ 0.755, P1 ⫽ 0.218, P2 ⫽ 0.027 M2 ␻0 ⫽ 0, ␻1 ⫽ 1, ␻2 ⫽ 10.367 P0 ⫽ 0.627, P1 ⫽ 0.342, P2 ⫽ 0.031

⫺1,879.115

⬍0.001

⫺1,880.275

⬍0.001

Alternative model

Most positively selected sites (ranging from 50 to 83% in each patient) were located within the immunogenic domains. Moreover, patient A had only four sites (9, 116, 121, and 180), patient B had three sites (9, 121, and 180), patient C had seven sites (9, 17, 28, 149, 158, 180, and 191), patient D had four sites (6, 16, 171, and 187), and patient E had five sites (6, 16, 171, 187, and 190) that were homologous to those reported by Anisimova and Yang (2). By analyzing HDV variants before and after intense immune selections, additional sites under positive selection were disclosed. Indeed, these sites demonstrated a stronger association with the

TABLE 2. Positively selected sites predicted in the delta antigen gene Patient

immunogenic domains (ranging from 73.6 to 100%). In the present study, we also observed a higher percentage of conserved sites than previously reported (2). This suggests that most sites essential for viral structure or replication are regulated by a purifying selection. 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 HLAA*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 HLAA*0201 patients. By analyzing the published peptides’ binding affinity along with those predicted by the neural network algo-

Positively selected sitesa

A ..........................9, 34*, 47, 48, 49, 112, 113, 116, 121, 180, 188 B...........................8, 9, 43, 48, 78*, 81, 121, 170, 172, 173*, 174, 180 C...........................7, 9, 16, 17, 26, 27, 28, 29, 37, 43, 44, 48, 78, 112, 122, 124, 131, 141, 149, 154, 155, 158, 169, 180, 184, 191 D ..........................6*, 16*, 171, 172, 181, 187 E...........................6, 14, 16, 26*, 36*, 54*, 62, 77*, 86, 94, 103*, 112, 119, 131*, 141*, 171, 172, 173, 183*, 187, 190, 203 a That is, the sites under selection inferred under model M3. Sites predicted to be under positive selection with a probability of ⬎0.9 were included. Those with a probability of ⬎0.95 are indicated with an asterisk; those with a probability of ⬎0.99 are in boldface.

TABLE 3. Characteristics of patients and HDV quasispecies Patient

Genotype

A B

I I

C D E

II IV IV

HLA typing

Selected blood samplesa

A2, A11 1435, 5132 A2 1573, 2643, 2667, 2688, 2748, 2781, 3123, 3678 A2, A24 1385, 4077 A2, A11 1025, 3038 A11 D62, D78, 338, 962, 2199, 2621

Total no. of clones

22 67 14 20 64

a The numbers indicate serial numbers of blood samples in our serum bank that have been collected and labeled consecutively.

VOL. 81, 2007

SELECTION OF HDAg IN PATIENTS WITH HDV INFECTION

TABLE 4. Potential CTL epitopes on HDAg of four patients with HLA-A*0201 Sitesa (aa)

Sequence

Binding level (10⫺9 M)

26–34* 29–37 43–51* 50–58 114–122 154–162 194–202 198–206

KLEELERDL ELERDLRKV KLEDDNPWL WLGNIKGIL QLSAGGKSL PLEGGSRGA FPWDILFPA ILFPADPPF

5,199 7,419 175 3,807 2,373 9,100 5,815 3,437

26–34* 43–51* 50–58 114–122 154–162 166–174 194–202 198–206

KLEDLEREA KLEDENPWL WLGNIKGIL QLAAGGKHL PLEGGSRGA GFVPNMLSV FPWDILFPA ILFPADPPF

8,432 205 3,807 2,815 9,100 2,817 5,815 3,437

II, C

44–52* 51–59 114–122

RLEDENPWL WLGNILGII QLSSGGKNL

360 1,405 3,791

IV, D

29–37* 43–51* 50–58 113–122 153–161

ELERDLRKV RLEEDNPWL WLGNIRGII QLSSGGKNL PLDGGPRGA

7,419 226 3,789 3,791 8,765

Genotype, patient

I, A

I, B

a Predicted epitopes are included if the binding level was ⬍10⫺5 M. *, CTL epitopes in an HLA-A*0201 transgenic mice model (10).

rithm, a cutoff level of 10⫺5 M was determined. Of the potential epitopes, only those predicted binding affinity stronger than 10⫺5 M were included. Table 4 shows the potential cytotoxic T ligand (CTL) epitopes of the four HLA-A2 patients. The potential CTL locations predicted by neural network methods were highly conserved within each genotype. Two HLAA*0201 HDV-I-infected patients had an 87.5% (7 of 8) consistency; two HLA-A*1101 HDV-IV-infected patients had a 100% identity. Of the HLA-A*0201 CTL epitopes, all genotypes had three consistent ones: aa 43 to 51, 50 to 58, and 114 to 122 (annotation of genotype I). Genotypes IV and I shared another common epitope of aa 153 to 161. The predicted affinity of aa 43 to 51 topped among the predicted ones in all genotypes. Of the HLA-A*1101 CTL epitopes, the genotype IV and I shared a common epitope of aa 52 to 60. The differences on predicted epitopes between genotypes may be due to the heterogeneity in sequences of different genotypes. The prediction within each genotype was highly consistent, which suggested a high reproducibility of the algorithm. 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

4441

TABLE 5. Patients with selected sites located within immunogenic epitopes of HDAg Patient(s) with selected sites (P ⬎ 0.9)

Cell type

N- to C-terminal aa

CD4⫹ T cell

26 to 41 50 to 65 66 to 81 106 to 121

A, C, E E B, C, E A, B, C, E

B cell

2 to 7 63 to 74 86 to 91 94 to 100 159 to 172 174 to 195 197 to 207

C, D, E E E E B, C, D, E A, B, C, D, E E

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. Figure 1 shows the consequence of the selection of novel HDV dominants after acute exacerbation in patient A. Although this patient had active cirrhosis before the last episode of acute exacerbation, the disease course went into biochemical remission. The novel dominant HDV variant of sample 5132 had an alteration of amino acid residues at the CD4 epitopes, CD8 epitopes, cryptic RNA-binding domain (aa 2 to 27), coiled-coil structure (aa 31 to 52), and helix-loop-helix (aa 108 to 135) of the RNA-binding domain (aa 97 to 146). The serum ALT level gradually became normalized and remained within normal limits for 70 months at the last follow-up. Figure 2 demonstrates the emergence of a novel HDV dominant variant after acute exacerbation of chronic hepatitis D. The novel dominant variant of sample 3038 had amino acid residues changes at B-cell epitopes and cryptic RNA-binding domain (aa 2 to 27). DISCUSSION We examined the population of HDV during chronic infection. Despite the majorities of delta antigen sites being functionally constrained and thereby under purifying selection, the analysis revealed that positive selection predominates the modeling of HDV population during intense immune responses from host. An average of 4.9% HDAg sites were shaped by a positive selection in our study. Crucially, most of these positively selected sites were located within immunogenic domains. Various immunologic mechanisms contributed to the modeling of HDV population. In our previous study, DNA-based immunization could lead the endogenous small

4442

WANG ET AL.

J. VIROL.

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.

HDAg to generate a significantly higher anti-HDV response than the large HDAg (11). Epitope mapping further showed that the peptide from aa 174 to 195 could bind more strongly with antibodies induced by the DNA vaccine expressing the small HDAg than those induced by the DNA vaccine expressing the large HDAg. In another study, we also observed that HDV variants with residue changes within the epitope from aa 174 to 195 often emerged after severe hepatitis attacks in chronic HDV-infected patients (29). Interestingly, all patients in the present study had a positive selection occurred within the B-cell epitope. The humoral immunity associated with this epitope may have an important role in HDV chronicity. 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 HDAgspecific 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

VOL. 81, 2007

SELECTION OF HDAg IN PATIENTS WITH HDV INFECTION

4443

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).

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 HDVinfected 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 (ELERD LRKV) 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 ⫻ 10⫺5 M, which was higher than the 0.7419 ⫻ 10⫺5 M of ligand 29-37 or our designated 10⫺5 M. Thus, the peptide in the Huang et al. study was suggested to have a weaker binding affinity. The

4444

WANG ET AL.

J. VIROL.

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.

11. 12. 13. 14. 15.

16.

17.

18.

19. 20.

21.

ACKNOWLEDGMENTS This study was supported by grants NSC94-2314-B-010-041, NSC952314-B-010-011, and NSC95-2314-B-010-021-MY3 from the National Science Council of the Executive Yuan, Taiwan. REFERENCES 1. Accapezzato, D., R. Nisini, M. Paroli, G. Bruno, F. Bonino, M. Houghton, and V. Barnaba. 1998. Generation of an MHC class II-restricted T-cell epitope by extracellular processing of hepatitis delta antigen. J. Immunol. 160:5262–5266. 2. Anisimova, M., and Z. Yang. 2004. Molecular evolution of the hepatitis delta virus antigen gene: recombination or positive selection? J. Mol. Evol. 59: 815–826. 3. Chang, J., S. O. Gudima, and J. M. Taylor. 2005. Evolution of hepatitis delta virus RNA genome following long-term replication in cell culture. J. Virol. 79:13310–13316. 4. Chao, Y. C., C. M. Lee, H. S. Tang, S. Govindarajan, and M. M. Lai. 1991. Molecular cloning and characterization of an isolate of hepatitis delta virus from Taiwan. Hepatology 13:345–352. 5. Chung, C. T., and R. H. Miller. 1988. A rapid and convenient method for the preparation and storage of competent bacterial cells. Nucleic Acids Res. 16:3580. 6. Cole, S. M., E. J. Gowans, T. B. Macnaughton, P. D. Hall, and C. J. Burrell. 1991. Direct evidence for cytotoxicity associated with expression of hepatitis delta virus antigen. Hepatology 13:845–851. 7. Govindarajan, S., K. P. Chin, A. G. Redeker, and R. L. Peters. 1984. Fulminant B viral hepatitis: role of delta agent. Gastroenterology 86:1417–1420. 8. Gruner, N. H., T. J. Gerlach, M. C. Jung, H. M. Diepolder, C. A. Schirren, W. W. Schraut, R. Hoffmann, R. Zachoval, T. Santantonio, M. Cucchiarini, A. Cerny, and G. R. Pape. 2000. Association of hepatitis C virus-specific CD8⫹ T cells with viral clearance in acute hepatitis C. J. Infect. Dis. 181: 1528–1536. 9. Hadler, S. C., M. De Monzon, A. Ponzetto, E. Anzola, D. Rivero, A. Mondolfi, A. Bracho, D. P. Francis, M. A. Gerber, S. Thung, et al. 1984. Delta virus infection and severe hepatitis: an epidemic in the Yucpa Indians of Venezuela. Ann. Intern. Med. 100:339–344. 10. Huang, Y. H., M. H. Tao, C. P. Hu, W. J. Syu, and J. C. Wu. 2004. Identi-

22.

23. 24.

25. 26.

27.

28. 29. 30. 31. 32. 33.

fication of novel HLA-A*0201-restricted CD8⫹ T-cell epitopes on hepatitis delta virus. J. Gen. Virol. 85:3089–3098. Huang, Y. H., J. C. Wu, S. C. Hsu, and W. J. Syu. 2003. Varied immunity generated in mice by DNA vaccines with large and small hepatitis delta antigens. J. Virol. 77:12980–12985. Krushkal, J., and W. H. Li. 1995. Substitution rates in hepatitis delta virus. J. Mol. Evol. 41:721–726. Kwok, S., and R. Higuchi. 1989. Avoiding false positives with PCR. Nature 339:237–238. Modahl, L. E., and M. M. Lai. 2000. Hepatitis delta virus: the molecular basis of laboratory diagnosis. Crit. Rev. Clin. Lab. Sci. 37:45–92. Nielsen, M., C. Lundegaard, P. Worning, S. L. Lauemoller, K. Lamberth, S. Buus, S. Brunak, and O. Lund. 2003. Reliable prediction of T-cell epitopes using neural networks with novel sequence representations. Protein Sci. 12:1007–1017. Nisini, R., M. Paroli, D. Accapezzato, F. Bonino, F. Rosina, T. Santantonio, F. Sallusto, A. Amoroso, M. Houghton, and V. Barnaba. 1997. Human CD4⫹ T-cell response to hepatitis delta virus: identification of multiple epitopes and characterization of T-helper cytokine profiles. J. Virol. 71:2241–2251. Ottobrelli, A., A. Marzano, A. Smedile, S. Recchia, M. Salizzoni, C. Cornu, M. E. Lamy, J. B. Otte, B. De Hemptinne, A. Geubel, et al. 1991. Patterns of hepatitis delta virus reinfection and disease in liver transplantation. Gastroenterology 101:1649–1655. Radjef, N., E. Gordien, V. Ivaniushina, E. Gault, P. Anais, T. Drugan, J. C. Trinchet, D. Roulot, M. Tamby, M. C. Milinkovitch, and P. Deny. 2004. Molecular phylogenetic analyses indicate a wide and ancient radiation of African hepatitis delta virus, suggesting a deltavirus genus of at least seven major clades. J. Virol. 78:2537–2544. Rizzetto, M., M. G. Canese, J. L. Gerin, W. T. London, D. L. Sly, and R. H. Purcell. 1980. Transmission of the hepatitis B virus-associated delta antigen to chimpanzees. J. Infect. Dis. 141:590–602. Rizzetto, M., G. Verme, S. Recchia, F. Bonino, P. Farci, S. Arico, R. Calzia, A. Picciotto, M. Colombo, and H. Popper. 1983. Chronic hepatitis in carriers of hepatitis B surface antigen, with intrahepatic expression of the delta antigen. An active and progressive disease unresponsive to immunosuppressive treatment. Ann. Intern. Med. 98:437–441. Roingeard, P., F. Dubois, P. Marcellin, J. Bernuau, S. Bonduelle, J. P. Benhamou, and A. Goudeau. 1992. Persistent delta antigenaemia in chronic delta hepatitis and its relation with human immunodeficiency virus infection. J. Med. Virol. 38:191–194. Su, C. W., Y. H. Huang, T. I. Huo, H. H. Shih, I. J. Sheen, S. W. Chen, P. C. Lee, S. D. Lee, and J. C. Wu. 2006. Genotypes and viremia of hepatitis B and D viruses are associated with outcomes of chronic hepatitis D patients. Gastroenterology 130:1625–1635. Taylor, J. M. 2003. Replication of human hepatitis delta virus: recent developments. Trends Microbiol. 11:185–190. Thimme, R., S. Wieland, C. Steiger, J. Ghrayeb, K. A. Reimann, R. H. Purcell, and F. V. Chisari. 2003. CD8⫹ T cells mediate viral clearance and disease pathogenesis during acute hepatitis B virus infection. J. Virol. 77: 68–76. Wang, J. G., R. W. Jansen, E. A. Brown, and S. M. Lemon. 1990. Immunogenic domains of hepatitis delta virus antigen: peptide mapping of epitopes recognized by human and woodchuck antibodies. J. Virol. 64:1108–1116. Wu, J. C., C. L. Chen, M. C. Hou, T. Z. Chen, S. D. Lee, and K. J. Lo. 1994. Multiple viral infection as the most common cause of fulminant and subfulminant viral hepatitis in an area endemic for hepatitis B: application and limitations of the polymerase chain reaction. Hepatology 19:836–840. Wu, J. C., T. Z. Chen, Y. S. Huang, F. S. Yen, L. T. Ting, W. Y. Sheng, S. H. Tsay, and S. D. Lee. 1995. Natural history of hepatitis D viral superinfection: significance of viremia detected by polymerase chain reaction. Gastroenterology 108:796–802. Wu, J. C., T. Y. Chiang, and I. J. Sheen. 1998. Characterization and phylogenetic analysis of a novel hepatitis D virus strain discovered by restriction fragment length polymorphism analysis. J. Gen. Virol. 79(Pt. 5):1105–1113. Wu, J. C., T. Y. Chiang, W. K. Shiue, S. Y. Wang, I. J. Sheen, Y. H. Huang, and W. J. Syu. 1999. Recombination of hepatitis D virus RNA sequences and its implications. Mol. Biol. Evol. 16:1622–1632. Wu, J. C., K. B. Choo, C. M. Chen, T. Z. Chen, T. I. Huo, and S. D. Lee. 1995. Genotyping of hepatitis D virus by restriction-fragment length polymorphism and relation to outcome of hepatitis D. Lancet 346:939–941. Yang, Z. 1997. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13:555–556. Yang, Z., and J. P. Bielawski. 2000. Statistical methods for detecting molecular adaptation. Trends Ecol. Evol. 15:496–503. Yang, Z., R. Nielsen, N. Goldman, and A. M. Pedersen. 2000. Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 155:431–449.