Virus Genes (2008) 36:299–305 DOI 10.1007/s11262-008-0205-5
Complete genomic sequence and phylogenetic relatedness of hepatitis B virus isolates in Cambodia Tran Thien Tuan Huy Æ Amadou Alpha Sall Æ Jean Marc Reynes Æ Kenji Abe
Received: 15 September 2007 / Accepted: 21 January 2008 / Published online: 9 February 2008 Ó Springer Science+Business Media, LLC 2008
Abstract Although it is known that Cambodia is one of the high endemic area of hepatitis B virus (HBV) infection, molecular characterization of HBV circulating in this country has not been reported. In this study, pre-S gene of HBV from 12 Cambodian patients was sequenced. Phylogenetic analysis based on the pre-S gene sequence revealed that 8 out of 12 isolates (66.7%) belonged to HBV/C1 and remaining four (33.3%) were HBV/B4. Furthermore, complete genomic sequences were also determined for three Cambodian HBV isolates. They all comprised of 3,215 bp long and two of them belonged to subgenotype B4, which had recombination event with genotype C in the precore/core gene confirmed by SimPlot and BootScanning analyses. Our results showed that both HBV strains belonged to subgenotypes B4 and C1, which are circulating in this country. This is the first report on molecular characterization of the HBV prevalent in Cambodia. Keywords Hepatitis B virus (HBV) HBV genome HBV genotypes HBV Asian strain Cambodia
T. T. T. Huy K. Abe (&) Department of Pathology, National Institute of Infectious Diseases, Toyama 1-23-1, Shinjuku-ku, Tokyo 162-8640, Japan e-mail:
[email protected] A. A. Sall J. M. Reynes Virology Unit, Pasteur Institute of Cambodia, Phnom Penh, Cambodia Present Address: A. A. Sall Institut Pasteur de Dakar, Dakar, Senegal Present Address: J. M. Reynes Institut Pasteur de Madagascar, Antananarivo, Madagascar
Introduction Hepatitis B virus (HBV) is an important etiologic agent of acute or chronic hepatitis and hepatocellular carcinoma, with over 350 million chronic infected patients around the world [1]. From the genetic variation of HBV strains isolated from various geographic regions, HBV is now classified into at least eight different genotypes, from A to H [2–5]. This HBV genotype classification has a distinct geographical distribution and has been shown to affect the clinical manifestation and course of HBV infection [6, 7]. In the Asian continent, HBV genotypes B (HBV/B) and C (HBV/C) are the prevalent types and are called the major Asian strain. HBV/C seems to be more pathogenic, in comparison with HBV/B [8]. Thus, characterization of the HBV genome could be important for understanding not only the molecular epidemiology, but also the relevant clinical fields. Cambodia, a Southeast Asian country, has been known to have endemic HBV infection, with an HBsAg carrier rate of approximately 2–14% of children and 8% of adults [9]. However, there have been no reports on HBV sequences, or the characteristics of the HBV genotype in this area. In the present study, full genome and pre-S gene sequences of HBV were investigated to determine the viral characterization of HBV in this country.
Materials and methods A total of 12 serum samples positive for HBsAg from Cambodian HBV chronic carriers (all males, aged 19– 32 years) were collected in the year 2003, from the Pasteur Institute of Cambodia, Cambodia. These sera were stored at -80°C until use. Serum HBsAg, anti-HCV antibody and
123
300
anti-HIV antibody were screened by enzyme linked immunosorbent assay (ELISA), using commercially available kits (axsym HBsAg for HBsAg, axsym HCV for anti-HCV antibody and axsym HIV Ag/Ab for anti-HIV antibody from Abbott Laboratories, North Chicago, IL). These sera were negative for anti-HCV and anti HIV-1 antibody. Informed consent for participation in this study was obtained from each individual. Viral DNA was extracted from 100 ll of serum using the DNA/RNA extraction Kit (SepaGene RV-R, Sanko Junyaku Co., Ltd., Tokyo, Japan). The resulting pellet was eluted in 50 ll of RNase-free water and kept in -20°C until use.
Amplification of HBV DNA Partial gene of HBV covered the pre-S region (from beginning of pre-S1 to the end of pre-S2 (522 bases)) in 12 isolates was amplified by hemi-nested PCR. Primer pairs PS2–1 (50 -TCCTGCTGGTGGCTCCAGTTC-30 , 55–75) and S6R (50 -TGCRTCAGCAAACACTTGGCA-30 , 1,194– 1,174, with R: A/G) were used for the first round; and primer pairs PS2-1, S7R (50 -GGCCTTRTAAGTTGGCGARAA-30 , 1,116–1,096, with R: A/G) for the second round PCR. After a denaturing step at 95°C for 2 min, PCR reactions were performed in 40 cycles (95°C: 20 s, 55°C: 20 s, and 72°C: 1 min) followed by an extension step at 72°C for 7 min. Amplification for full-length HBV genome was performed as described previously with modification [10]. In brief, five overlapping fragments were obtained by nested or hemi-nested PCR using the following primer combinations for the second round of PCR: HBV-1(50 -CCGG AAGAATTCTTTTTCACCTCTGCCTAATCA-30 , 1,821– 1,841; underlining of the primer sequence indicates the EcoRI site)/BG1R (50 -ATAGGGGCATTTGGTGGTCT-30 , 2,316–2,297) for fragment A; HBc1 (50 -AGTGTGGAT TCGCACTCCT-30 , 2,269–2,287)/PS8R (50 -ARGCCCTGA GCCTGAGGGCTC-30 , 3,098–3,078, with R: A/G) for fragment B; P1 (50 -TCACCATATTCTTGGGAACAAGA30 , 2,817–2,839)/S1-2 (50 -CGAACCACTGAACAAAT GGC-30 , 704–685) for fragment C; S1-1 (50 -TCGTGTTACAGGCGGGGTTT-30 , 192–211)/HBx2 (50 -ACGTGCAGA GGTGAAGCGAAG-30 , 1,604–1,584) for fragment D and HBx1 (50 -GTCCCCTTCTTCATCTGCCATCTGCCGT-30 , 1,487–1,507)/HBV-2 (50 -CCGGAAGAATTCAAAAAGT TGCATGGTGCTGG-30 , 1,825–1,806; underlining of the primer sequence indicates the EcoRI site) for fragment E. The first round of PCR was done with the primer combinations of HBV-1/HBV-2 which was described previously [11]. The first round of PCR was carried out for 40 cycles (95°C for 15 s, 55°C for 45 s, and 72°C for 3.5 min). The
123
Virus Genes (2008) 36:299–305
second round was carried out for 35 cycles (94°C for 20 min, 55°C for 20 min, and 72°C for 1 min) followed by extension at 72°C for 7 min. PCR products were separated by 1% agarose gel electrophoresis and recovered using the QIAquick gel extraction kit (Qiagen Inc., Chatsworth, CA, USA). DNA sequencing Purified PCR products were subjected to direct sequencing using the ABI PRISMTM Big Dye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems, Foster City, CA, USA). The inner primer pairs were used as sequencing primers for pre-S gene. Sequences of amplified DNA were determined using automated DNA sequencer (ABI model 3130, Applied Biosystems). Phylogenetic analysis Nucleotide sequences were multiple-aligned, analyzed using Genetyx for Windows ver. 6.0 software (GENETYX, Tokyo, Japan), and corrected manually by visual inspection. Genetic distances were calculated by using the Kimura’s twoparameter method. Phylogenetic trees were constructed by the neighbor-joining method [12]. To confirm the reliability of the pairwise comparison and phylogenetic tree analysis, bootstrap resampling was carried out in 1,000 data sets. Evolutionary and phylogenetic tree analyses were conducted using MEGA Version 2.1 [13]. Recombination analysis Recombination of viral genomes was investigated using SimPlot, BootScanning program (SimPlot software, Version 3.0, distributed kindly by Stuart C. Ray, Division of Infectious Diseases, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA). Nucleotide sequence accession numbers The nucleotide sequence data reported in this article have been submitted to the DDBJ/GenBank/EMBL databases under Accession No. AB214636–AB214647 (for pre-S gene), and AB115551, AB117758, AB117759 (for fulllength genome).
Results Pre-S gene (pre-S1/pre-S2, 522 nucleotides) isolated from 12 Cambodian patients were directly sequenced and were
Virus Genes (2008) 36:299–305
301
highly conserve amino acid sequence to HBV/B4 or HBV/ C1. Genotype C strains possessed Gln51, Ala54, Val60, Ser62, which were known to be conserved in HBV/C1. Deletion of ten amino acids (30 nucleotides) was detected in one strain (CB16) with no frame-shift of the sequence. This strain also had the highest substitution rate, up to nine amino acids. Genotype B strains were found rather conserved especially in the pre-S2 region, despite of the higher rate of substitution detected among genotype C strains.
compared to those of 27 genotype A–H strains retrieved from GenBank. Phylogenetic tree with neighbor-joining method showed that Cambodian strains were grouped into genotypes B and C branches with bootstrap value 100% and 92%, respectively (Fig. 1). Furthermore, 8 out of 12 isolates (66.7%) classified to subtype C1 of genotype C (HBV/C1) and remaining four isolates (33.3%) belonged to subtype B4 of genotype B (HBV/B4). In the pre-S gene, genotypes B and C strains of Cambodian HBV showed Fig. 1 Phylogram generated by neighbor-joining analysis of genetic distances in the pre-S1/pre-S2 gene of HBV. The percentages of bootstrap replicates supporting these branches are shown in number
Pre-S gene 100
G E D
78
F H
F-AB166850_Bolivia H-AY090454_Nicaragua
G-AB056513_USA E-AB205190_Ghana D-AB126581_Russia B1-AB010291_Japan B1-D23678_Japan B1-D00329_Japan B1-AB205121_Japan 99
B1
B2-AF121249_Vietnam B2-AF121250_Vietnam
B2
CB15_AB115551_Cambodia
100
B
B4-AB212626_Vietnam B4-AB205122_Vietnam B4-AB212625_Vietnam B4-AB031266_Vietnam
100
B4
CB14_AB117759_Cambodia CB11_AB214637_Cambodia CB17_AB214643_Cambodia B3-AB033554_Indonesia B3-D00331_Indonesia B5-AB241116_Philippines 100 B5-AB241117_Philippines A1-AB076678_Malawi A1_M57663_Philippines 100 A2-AJ012207_Germany A2-AB126580_Russia C4-AB0487048_Aus. Aborigine 100 C4-AB048705_Aus. Aborigine C3-X75656_Polynesia C3-X75665_Polynesia C5-AB241109_Philippines 100 C5-EU410080_Philippines C5-EU410081_Philippines C2-AF182802_China C2-AB205123_China C2-AB049610_Japan C2-AB205124_Japan C2-AB049609_Japan C2-AY247032_Korea C1-AB205125_Vietnam C1-AB112065_Vietnam C1-AB112063_Vietnam C1-AB111946_Vietnam 100 C1-AB112472_Thailand C1-AB112348_Myanmar C1-AB112066_Myanmar C1-AB112471_Thailand
B3
100
98
100 97
A
100
89 93
92
C
98
70
0.01
B5
C4 C3 C5
C2
CB16_AB214642_Cambodia CB23_AB214645_Cambodia CB12_AB214638_Cambodia CB13_AB214639_Cambodia CB20_AB214644_Cambodia CB24_AB214646_Cambodia CB25_AB214647_Cambodia CB10_AB117758_Cambodia
C1
C1-AB112408_Myanmar
123
302
Virus Genes (2008) 36:299–305
3,215 nucleotides, which were the characteristic of both genotypes B and C and designated HBV-CB10, HBVCB14, and HBV-CB15, respectively. Phylogenetic analysis of these complete genomes, compared with isolates retrieved from the database, confirmed these Cambodian strains belonged to HBV/B4 (CB14 and CB15) and HBV/
However, the virus attachment protein at amino acid positions 21–47 within the pre-S2 gene was conserved in all Cambodian strains. To further investigate the characteristic of HBV in Cambodia, complete genomic sequences from three Cambodian isolates were determined. They all comprised of Fig. 2 Phylogram generated by neighbor-joining analysis of genetic distances in the fulllength sequence of HBV. The percentages of bootstrap replicates supporting these branches are shown in number
Full genome F H
100
F-AB166850_Bolivia H-AY090454_Nicaragua
G
G-AB056513_USA
D 98 E
D-AB126581_Russia E-AB205190_Ghana A1-AB076678_Malawi
100
A
A1-M57663_Philippine A2-AJ012207_Germany
100 100
A2-AB126580_Russia 100
B1-D00329_Japan
100
B1-AB205121_Japan B1-AB010291_Japan
81 100
B2-AF121250_Vietnam
100
99
100
100
B2
B2-AB205120_Japan
CB14_AB117759_Cambodia
100
B
B1
B1-D23678_Japan B2-AF121249_Vietnam
B4-AB212625_Vietnam
100
CB15_AB115551_Cambodia
B4
B4-AB205122_Vietnam B4-AB212626_Vietnam B3-AB033554_Indonesia
100
B3-D00331_Indonesia
100
99
100 100
B5-AB241116_Philippines B5-AB241117_Philippines C4-AB0487048_Aus. Aborigine C4-AB048705_Aus. Aborigine
C3-X75656_Polynesia
100
100
C3-X75665_Polynesia
C
100 68
96
B5 C4 C3 C5
C2-AB049610_Japan
99 67
C5-AB241109_Philippines C5-EU410080_Philippines
B3
C2-AB205123_China C2-AB049609_Japan
97 98
C2
C2-AB205124_Japan 89
C1-AB112472_Thailand C1-AB205125_Vietnam C1-AB112063_Vietnam
100
C1-AB111946_Vietnam 100
C1-AB112066_Myanmar C1-AB112348_Myanmar
79 0.01
123
C1-AB112471_Thailand
CB10_AB117758_Cambodia
C1
Virus Genes (2008) 36:299–305
303
Table 1 Mean percentage and range of nucleotide divergence over the complete genome among three Cambodian isolates of HBV with different genotypes Genotype
Cambodia HBV/B4
Cambodia HBV/B4
0.026
Cambodia HBV/C1
A
B
C
D
E
F
Cambodia HBV/C1
0.088
A
0.096
0.090
B
0.046
0.097
0.094
C
0.089
0.038
0.088
0.097
D
0.110
0.107
0.101
0.109
0.105
E F
0.110 0.148
0.107 0.146
0.094 0.145
0.108 0.148
0.104 0.143
0.076 0.149
0.139
G
0.127
0.128
0.113
0.132
0.126
0.120
0.112
0.155
H
0.144
0.145
0.144
0.149
0.144
0.151
0.146
0.086
G
H
0.155
Bold italic indicates values obtained by comparison between Cambodian strains isolated in this study and bold indicate values between Cambodian isolates and distinct isolates retrieved from the database
C1 (CB10) with both of bootstrap value 100% (Fig. 2). Pairwise distance analysis showed that HBV/B4 strains from Cambodia had 97.4% similarity to each other and 91.2% to Cambodian HBV/C1 strain (Table 1). Cambodian HBV/B4 strains showed 85.2–95.4% homology with other reported HBV genotypes from database. Similarly, Cambodian HBV/C1 strain showed 85.4–96.2% with other reported HBV genotypes. In addition, SimPlot and BootScanning analyses were conducted to determine the possible recombination sites in these Cambodian isolates. Three Cambodian HBV isolates with full-length sequence Fig. 3 The location of recombinant event in HBVCB14 isolate was determined using the SimPlot program and bootscanning analysis. The HBV-CB14 was compared with genotypes B and C on the full genome. The dotted screen shows breakpoint of recombination (estimated at nt1824–2318) with genotype C
were compared with distinct isolates belonging to eight genotypes. By this analysis, the recombination event was observed in two (CB14 and CB15) HBV/B4 strains, but not in CB10 of HBV/C1. These two isolates had recombination event with genotype C in the precore/core gene (Fig. 3). The position of the recombination breakpoints were estimated at nt1824–2318. In the S gene, based on the deduced amino acid at positions 122, 127, and 160, the subtypes of the Cambodian strains were identified. Genotype C strain (CB10) was identified as subtype adrq+, and genotype B strains
BootScan - Query: HBV-CB14 (HBV/B4) preC/C P Pre-S2/S
P X
Pre-S1 1824
2318
100 95 90 85 80 75
% of Permuted Trees
70 65 60 55 50
45 40 35 30 25 20 15 10 5 0 0
100 200 300 400 500 600 700 800 900 1,000 1,100 1,200 1,300 1,400 1,500 1,600 1,700 1,800 1,900 2,000 2,100 2,200 2,300 2,400 2,500 2,600 2,700 2,800 2,900 3,000 3,100 3,200
Nucleotide position
Genotype B_D00329
Genotype C_AB049609
123
304
(CB14 and CB15) belonged to subtype adw2 and ayw1, respectively. For the ‘a’ determinant domain at positions 122–147 in the HBsAg region, isolate CB14 had amino acid substitutions at Q129R, G130N, and T131N. The other two strains showed conserved amino acid sequence in both of the first and the second loops of this domain. In the precore/core ORF, no deletion or insertion was found. Within the core promoter region¸ one isolate (CB10) had A to T substitution at position 1762 (A1762T) and possessed precore stop codon mutation from G to A at position 1,896. Other two remained the wild type of both the basal core promoter and precore sequence. In the polymerase ORF, three Cambodian strains showed no deletion or insertion. The YMDD motif, from residues 549–552 of the P gene, which placed in the reverse transcriptase (RT) region, was conserved for all Cambodian strains isolated in this study.
Discussion Our data presented here showed that genotypes B and C strains of HBV were detected in this country. To our knowledge, this is the first report on molecular characterization of the HBV prevalent in Cambodia. It is known that genotypes B and C of HBV are Asian strains since they are prevalent mainly throughout Asian countries. The coexistence of genotypes B and C in this country were concordant with the distribution of HBV genotypes in the neighboring countries in Southeast Asia, such as Thailand and Vietnam [14]. Asian strains of HBV sometimes have genetic variations and can be classified into subgroups, such as subgenotypes B1 to B5 and subgenotypes C1 to C5 [15], respectively. Subgenotype C1 is distributed mainly in Southeast Asia, whereas subgenotypes B1 and C2 are prevalent throughout Far East Asia such as Japan, Korea, and China [10, 14, 16–18]. HBV isolated in this study belonged to subgenotypes B4 and C1, respectively, by phylogenetic analysis. Interestingly, subgenotype B4 shows recombination with genotype C and this event always occurs at the same position in the precore/core region from genotype C [19]. The characteristics of three Cambodian HBV strains were further investigated by full genome phylogenetic analysis. Two HBV/B4 strains from Cambodia exhibited high homology (97.4%) and also showed recombination event with genotype C within the precore/core gene. The HBV/C1 strain also had conserved amino acids Lys85, Gln143, Ser249, which is a characteristic of the subgroup of genotype C in Southeast Asia [10]. The HBV subtypes based on the deduced amino acid in the ‘‘a’’ determinant region also confirmed their relatedness to the distribution of genotypes B and C. Subtypes adw2 and ayw1 were
123
Virus Genes (2008) 36:299–305
found in genotype B and subtype adrq+ in genotype C in the Cambodian strains isolated in this study. This subtypegenotype relationship is well documented in Southeast Asia countries like Vietnam and Thailand [unpublished data]. Within the ‘‘a’’ determinant region located within the S gene, an amino acid substitution was detected in one isolate. The ‘‘a’’ determinant is the major immune target of the polyclonal antibody to HBsAg. In patients infected with the HBV harboring a mutant in this region, the mutant HBsAg is not detectable by commonly used HBsAg assays. More importantly, the amino acid mutation at several points within the ‘‘a’’ determinant might lead to vaccine escape [20, 21]. Unfortunately, we could not obtain sufficient information on this patient (CB14) whether he had vaccination of HBV. In conclusion, we reported the genomic characterization of HBV isolated in Cambodia. The results revealed that HBV variants are distributed in this region and the information obtained may enable further characterization of the HBV Asian strain. Further investigations are needed to clarify the relatedness of the HBV variants and clinical manifestations including progression to hepatocellular carcinoma, vaccine escape mutants and the influence on lamivudine therapy in this country where actual situation on hepatitis virus-related liver diseases is not clear. Acknowledgments We thank Dr. Stephan Schaefer, Department of Virology, Institute for Medical Microbiology, Virology and Hygiene, Rostock University, Rostock, Germany, for kindly providing the sequence information of HBV/C5, and Dr. Tetsutaro Sata, National Institute of Infectious Diseases, Japan, for his continuous encouragement during this study. This study was supported in part by grantsin-aid for Science Research of the Ministry of Education, Culture, Sports, Science and Technology of Japan and the Ministry of Health, Labor and Welfare of Japan.
References 1. H. Okamoto, F. Tsuda, H. Sakugawa, R.I. Sastrosoewignjo, M. Imai, Y. Miyakawa, M. Mayumi, J. Gen. Virol. 69, 2575–2583 (1988) 2. H. Norder, A.M. Courouce, L.O. Magnius, Virology 198, 489– 503 (1994) 3. W.M. Lee, N. Engl. J. Med. 337, 1733–1745 (1997) 4. P. Arauz-Ruiz, H. Norder, B.H. Robertson, L.O. Magnius, J. Gen. Virol. 83, 2059–2073 (2002) 5. L. Stuyver, S. De Gendt, C. Van Geyt, F. Zoulim, M. Fried, R.F. Schinazi, R. Rossau, J. Gen. Virol. 81, 67–74 (2000) 6. J.H. Kao, P.J. Chen, M.Y. Lai, D.S. Chen, Gastroenterology 118, 554–559 (2000) 7. K. Kidd-Ljunggren, Y. Miyakawa, A.H. Kidd, J. Gen. Virol. 83, 1267–1280 (2002) 8. T.T. Huy, H. Ushijima, K.M. Win, P. Luengrojanakul, P.K. Shrestha, Z.H. Zhong, A.V. Smirnov, T.C. Taltavull, T. Sata, K. Abe, J. Clin. Microbiol. 41, 5449–5455 (2003) 9. E.G. Thuring, H.I. Joller-Jemelka, H. Sareth, U. Sokhan, C. Reth, P. Grob, Southeast Asian J. Trop. Med. Public Health 24, 239– 249 (1993)
Virus Genes (2008) 36:299–305 10. T.T. Huy, H. Ushijima, V.X. Quang, K.M. Win, P. Luengrojanakul, K. Kikuchi, T. Sata, K. Abe, J. Gen. Virol. 85, 283–292 (2004) 11. S. Gunther, B.C. Li, S. Miska, D.H. Kruger, H. Meisel, H. Will, J. Virol. 69, 5437–5444 (1995) 12. N. Saitou, M. Nei, Mol. Biol. Evol. 4, 406–425 (1987) 13. S. Kumar, K. Tamura, I.B. Jakobsen, M. Nei, Bioinformatics 17, 1244–1245 (2001) 14. T.T. Huy, K. Abe, Pediatr. Int. 46, 223–230 (2004) 15. S. Schaefer, World J. Gastroenterol. 13, 14–21 (2007) 16. A. Banerjee, S. Datta, P.K. Chandra, S. Roychowdhury, C.K. Panda, R. Chakravarty, World J. Gastroenterol. 12, 5964– 5971 (2006) 17. W.-C. Liu, P.H. Phiet, T.-Y. Chiang, K.-T. Sun, K.-H. Hung, K.-C. Young, I.-C. Wu, P.-N. Cheng, T.-T. Chang, Virus Res. 129, 212–223 (2007)
305 18. Z. Wang, Y. Tanaka, Y. Huang, F. Kurbanov, J. Chen, G. Zeng, B. Zhou, M. Mizokami, J. Houl, J. Clin. Microbiol. 45, 1491– 1496 (2007) 19. F. Sugauchi, E. Orito, T. Ichida, H. Kato, H. Sakugawa, S. Kakumu, T. Ishida, A. Chutaputti, C.L. Lai, R. Ueda, Y. Miyakawa, M. Mizokami, J. Virol. 76, 5985–5992 (2002) 20. W.F. Carman, A.R. Zanetti, P. Karayiannis, J. Waters, G. Manzillo, E. Tanzi, A.J. Zuckerman, H.C. Thomas, Lancet 336, 325–329 (1990) 21. J.A. Waters, M. Kennedy, P. Voet, P. Hauser, J. Petre, W. Carman, H.C. Thomas, J. Clin. Invest. 90, 2543–2547 (1992)
123