|
|
||||||||
1 Department of Virology, National Institute for Health Development, Tallinn, Estonia
2 Department of Virology, Immunology and Vaccinology, Swedish Institute for Infectious Disease Control, Stockholm, Sweden
3 Department of Microbiology, Tumor and Cell Biology, Karolinska Institute, Stockholm, Sweden
4 Department for Infectious Disease, Ida-Viru Central Hospital, Kohtla-Järve, Estonia
5 MP Chumakov Institute of Poliomyelitis and Viral Encephalitis, Moscow, Russia
6 St Petersburg Pasteur Institute, St Petersburg, Russia
Correspondence
Heléne Norder
helene.norder{at}smi.ki.se
| ABSTRACT |
|---|
|
|
|---|
The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper are EU594382-EU594436.
| INTRODUCTION |
|---|
|
|
|---|
In general, HBV genotypes as well as subgenotypes show a distinct geographical distribution, although for the latter information for different countries is still scarce. Thus, subgenotype A1 is predominant in South Asia and sub-Saharan Africa, and subgenotype A2 prevalent mainly in European and North-American countries (Kramvis et al., 2002
; Norder et al., 2004
; Sugauchi et al., 2004
). Recently A3 was described in native populations of West and Central Africa (Kurbanov et al., 2005
; Makuwa et al., 2006
). Definite geographical predilections have also been shown for genotypes B and C (Norder et al., 2004
; Sakamoto et al., 2006
; Liu et al., 2007
). Genotype D is the most widespread and divides into five subgenotypes (Norder et al., 2004
; Banerjee et al., 2006
), although the geographical distribution of these subgenotypes are so far less defined than those for genotypes A through C and F.
In a previous study, 205 HBV strains from Estonia and 14 other regions of the former Soviet Union (USSR) were genotyped by sequencing the small S genes. Genotype D was found to be dominating. However, an unexpectedly high number of strains were found to have identical S gene sequences, even though they were collected in geographically distant regions over a 13 years period (Tallo et al., 2004
). The aim of this study was to characterize the complete genomes of representative strains from the different regions to classify them into subgenotypes and to further investigate the subclades of strains with identical S genes.
| METHODS |
|---|
|
|
|---|
|
A modification of the technique and primers described by Gunther et al. (1995)
was used for amplification of complete HBV genomes. Briefly, 5 µl DNA was amplified with Expand Long Template PCR System (Roche Applied Science) in a 50 µl mixture containing 0.3 µM modified primers P1* (1821–1843), 5'-TTTTTCACCTCGCCTAATCA-3', and P2* (1825–1801), 5'-AAAAAGTTGCATGRTGMTGG-3', (Kramvis et al., 2005
) and 350 µM dNTPs. Amplification was performed for 40 cycles; the thermal profile of the reaction was denaturation at 94 °C for 40 s, annealing at 60 °C for 1 min and elongation at 68 °C for 3 min with an increment of 2 min after every 10 cycles. To obtain complete genomes, the core promoter region was amplified directly from extracted DNA with 0.5 µM primers Hep1230S and Hep56AS (position 2262–2244), 5'-AGTGCGAATCCACACTCCG-3', 1.75 mM MgCl2, 200 µM dNTPs and 1.25 U AmpliTaq DNA polymerase (Applied Biosystems) for 40 cycles, denaturation at 92 °C for 20 s, annealing at 58 °C for 25 s and elongation at 72 °C for 90 s. PCR fragments were purified using the GFXTM PCR DNA and Gel Band Purification kit (GE Healthcare). Purified products P1*/P2* (110 ng) and 1230/56 (40 ng) were used as templates in the sequencing reaction using the dideoxynucleotide chain-termination method with ABI PRISM TM BigDye Terminator Cycle Sequencing Reaction kit version 3.1 (Applied Biosystems) with a set of 12 sequencing primers: Hep35S, Hep38AS, Hep39S, Hep50AS, Hep56AS, Hep58AS, HepFAS, Hep1876S, Hep61AS, and Hep230S previously described (Norder et al., 1996
; Kramvis et al., 2005
), Hep55AS (position 3081–3064), 5'-TGAGCCTGAGGGCTCCAC-3', and Hep53AS (position 1430–1447), 5'-TCCCGTCGGCGCTGAATC-3'. Positions are given according to strain pHBV3200 (GenBank accession no. X02763
[GenBank]
). Primers Hep53AS, Hep1230S and Hep56AS were used for sequencing the fragment amplified with primers Hep1230S and Hep56AS. The ABI PRISM 3100 Genetic Analyzer (Applied Biosystems) was used for electrophoresis and data collection.
Sequence analysis.
The sequences obtained were edited using the SeqMan program in the Lasergene package (DNAStar) and the BioEdit program (http://www.mbio.ncsu.edu/BioEdit/BioEdit.html). The sequences obtained were aligned with 21 complete genotype A genomes belonging to subgenotypes A1–A3, and 55 wild-type genotype D sequences all within subgenotypes D1–D5. These sequences were retrieved from GenBank (Table 2
). All genotype D sequences in GenBank with core promoter and/or pre-core stop mutations were excluded. Phylogenetic analysis was carried out with the PHYLIP package version 3.67 (http://evolution.genetics.washington.edu/phylip/getme.html). Evolutionary distances were estimated with the DNADIST program using the F84 algorithm with a transition : transversion ratio of 1.05. Phylogenetic trees were constructed using the neighbour-joining (NJ) method in the PHYLIP package. Bootstrap analysis for 1000 replicas was performed with the SEQBOOT and CONSENSE programs.
|
| RESULTS |
|---|
|
|
|---|
Subgenotypes, genetic distance and analysis of complete and partial genome sequences
In the phylogenetic trees a number of clades supported by significant bootstrap values were identified within both genotypes A and D. The 13 genotype A strains encoding adw2 were in one clade on the branch formed by subgenotype A2 strains (Fig. 1a
). The mean intra-subgenotype divergence for these A2 strains was 0.4 %, while the intra-subgenotype divergence for all strains within A2 was 1.7 % (0–7 %).
|
|
Four strains, 4605-97, 4651-97, 320-95 and 321-96, shared the A1762T mutation within the basal core promoter. Another two strains, 4630-97 and 4496-97, had the G1896A stop mutation in the preC/C region, and one strain, 1901-99, had both these mutations. All these strains belonged to subgenotype D2 apart from 320-95, which belonged to subgenotype D3.
Subgenotypes and phylogenetic analysis of four ORFs
The sizes of preS1/S2/S, P, preC/C and X ORFs were as expected: 1203, 2538, 645 and 465 nt for genotype A strains and 1169, 2499, 639 and 465 nt for genotype D strains. All five subgenotypes of D were clearly distinguished when the P, preS1/S2/S and preC/C ORFs were analysed (Fig. 2a-d
). The X gene for subgenotype D1 could not be distinguished from that of D2 (Fig. 2d
), although the other subgenotypes, D3–D5, also showed distinct clustering in this region (Fig. 2d
).
|
|
When the preS/S gene products were compared, signature residues were identified for all the subgenotypes. There were five residues unique for the D subgenotypes, one each for D2, D3 and D4, and two for D5. In addition, there were five residues specific for genotype D that were not shared by all D subgenotypes, and at five positions there were subgenotype-specific residues shared by strains belonging to other genotypes (Table 4
). One signature residue each was identified for D1, D2 and D3, two for D4 and six for D5 (Table 4
). The D2 strains in this study expressed Val118 and Val128 in the S gene. In the preS region D1–D3 strains expressed Ala39 and Asn115 not found in other strains, while D4 strains had Asn39, and D5 strain expressed Arg39 and Asp115, the latter shared by the other genotypes. The D4 strains had a unique Thr115 not found in other strains. The deduced amino acid sequence for the seven strains with only the P-gene sequenced revealed a subgenotype D2-specific amino acid sequence, with Ala39, Ile85, Ser96, Ala158, Ile161, Arg167 and Arg173 in the preS region. In the S gene all nine strains had the subgenotype D2 unique substitutions at residue 188, either Val or Ala.
In the preC/C product there was one specific residue each for D1 and D5. Two residues were shared by only D4 and D5, thereby distinguishing them from each other and from the other subgenotypes. D2 could not be distinguished from D3 in this region. All these signature residues were shared with strains belonging to other genotypes (Table 4
). D1–D3 strains expressed Ile116 in the core region, whereas D4 and D5 strains had Leu116, shared by all other genotypes. The two strains, expressing ayw4, shared a combination of three substitutions, Glu40, Val74 and Ile116, characteristic for subgenotypes D2 and D3.
In the X product, strains belonging to D1 could not be distinguished from those belonging to D2, while D3 had one signature residue that was shared with genotype E strains. D4 and D5 had two signature residues each shared by strains belonging to other genotypes. In addition, D5 had two unique substitutions not found in any other strain (Table 4
). D1 and D2 strains shared a Ser46, while the other strains belonging to D3–D5 and to all other genotypes had Pro46. The deduced amino acid sequence of this region of the two ayw4 strains showed a combination of subgenotype D1- and D2-specific residues.
| DISCUSSION |
|---|
|
|
|---|
All genotype A strains were classified within subgenotype A2 by phylogenetic analysis of complete HBV genomes, as expected for strains from European populations (Kramvis et al., 2002
; Norder et al., 2004
). Strains belonging to genotype D were subdivided into subgenotypes D1–D3, a subdivision that was also obtained by phylogenetic analysis of the polymerase, preS1/S2/S and preC/C genes, but not for the X gene. Most genotype D strains belonged to D2 and encoded ayw3. These strains were derived from the European region, including Estonia and the Siberian and Far-Eastern parts of the former USSR. D2 has also been found in a small area in Western Japan, although genotypes B and C are prevalent in Japan (Michitaka et al., 2006
). It was suggested that the Japanese–Russian War from 1904 to 1905 was the most likely event that led to the transmission of D2 into Japan (Michitaka et al., 2006
). However, comparison of D2 strains from this study and those from Japan could not exclude this possibility, although the corresponding subclades differed by 2.3 % and the mean intra-subclade divergences were larger for the Japanese strains, 1.3 vs 1.0 %. The strains from the Middle Asia part of the former USSR belonged to D1. This subgenotype has also been reported as the most prevalent one in countries in the Middle East such as Turkey and Iran (Amini-Bavil-Olyaee et al. 2005
; Bozdayi et al., 2005
). Several strains from the European part of the investigated region belonged to subgenotype D3, all encoding ayw2 and none of them encoding ayw3, which is characteristic for the addict D3 clade (Norder et al., 2004
; De Maddalena et al., 2007
). Recently D3 strains encoding ayw2 were reported as the most prevalent subgenotype in Serbia and Italy, although D3 strains belonging to the addict ayw3 clade were also identified (De Maddalena et al., 2007
; Lazarevic et al., 2007
).
When comparing the dendrogram based on the small S gene (Tallo et al., 2004
) with that based on complete genomes, it was possible to classify all strains from the major cluster encoding ayw3 or ayw4 within D2 (Table 4
). The strains from two distinct clades encoding ayw2 specificity could be classified as subgenotypes D1 and D3. Interestingly, some strains encoding ayw4 specificity most probably belonged to D3 according to the S gene. Due to the insufficient sequence data it was not possible to classify five aberrant strains from the previous study (Tallo et al., 2004
).
The geographical distribution of genotype D subgenotypes showed that D2 seems to be the most prevalent subgenotype in Northern Eurasia, while D1 is highly prevalent in Middle East with extensions into the Republics of Middle Asia. D3 encoding ayw2 circulated mostly in low frequency in the studied region, but seems more widely distributed in other parts of Europe such as Italy and Serbia where it is the dominating subgenotype.
In a previous study, approximately 50 % of the HBV strains from the former USSR had identical S gene sequences (Tallo et al., 2004
). The strains with identical S genes often originated from regions distantly located from each other (Tallo et al., 2004
). In this study, sequencing of complete HBV genomes allowed the division of all strains with identical S genes, apart from one A2 strain, into distinct strains that separated mainly according to geographical origin. Thus, within subgroup D2 a number of strains from Estonia and Khabarovsk formed distinct subclades. This was also observed for subgenotype A2 strains from Yakutsk.
Analysis of molecular signature motifs demonstrated the specific combination of amino acid substitutions distinguishing the D1–D5 subgenotypes in the P and S ORFs. All subgenotypes could also be distinguished from each other in the two other ORFs, apart from D2 and D3, in the preC/C gene and D1 and D2 in the X gene. Interestingly, the subdivision into subgenotypes was also observed at the nucleotide level within three ORFs: P, preS/S and preC/C. The finding of signature substitutions for the subgenotypes will mostly enable the classification of HBV strains into subgenotypes by limited sequencing, in particular within the preS region. Interestingly, D1–D3 strains shared five unique residues, two in the P, two in the S, and one in the X product, while the D4 and D5 strains in these positions shared residues with strains belonging to the other genotypes. This indicates that D4 and D5 strains evolved before D1–D3. Interestingly, D2 had more unique substitutions than D1; that would argue for D2 rather than D1 being the last split within genotype D. Based on the deduced amino acid sequences, the seven incompletely sequenced strains could be assigned to subgenotype D2.
There were three regions with subgenotype-specific residues suitable for subgenotyping genotype D strains. All three were within the preS/S gene, which overlaps with the spacer and reverse transcriptase regions of the P gene. The shortest amino acid sequence containing specific residues distinguishing the subgenotypes was located between residues 58 and 128 in the spacer region of the P gene. This region overlaps with the preS residues 39 of preS1 to residue 42 of preS2 (preS161), which also contained residues discriminating between the subgenotypes. The subgenotypes could also be distinguished by specific substitutions in both reading frames between preS161 and residue 118 of the small S gene, which overlaps with the sp128 and rt121 of the P gene. In a third region the subgenotypes could only be discriminated by the amino acid sequence of the P gene reading frame between rt122 and rt247, which overlaps with s117 to s242, encompassing the a-determinant of the small S gene.
All genotype D strains with identical S gene sequences remained closely related when their complete genomes were compared, and the strains in clades II and IV seemed to have evolved from S gene clades I and V. Thus, there was no evidence for convergent evolution or recombination as an explanation for the occurrence of strains with identical S genes.
The mean intra-group divergence for the studied strain within subgenotypes D1, D2 and D3 was low, ranging from 1.0 to 1.8 %. The divergence 0–1.3 % was reported for subgenotype C2 strains in India (Banerjee et al., 2006
) and 1.6 % for subgenotype C5 strains in the Philippines (Sakamoto et al., 2006
). In most cases, the low genetic variability of HBV strains was explained by relatively recent introduction into the population and/or close geographical origin of strains. In contrast, a high genetic heterogeneity was reported for genotype D strains in chronically infected patients in Italy, and analysis of the P gene showed the intra-subgenotypic divergences were 3.1, 2.1, 3.0 and 1.6 % for subgenotypes D1, D2, D3 and D4, respectively, although for subgenotypes A1 and A2, the intra-subgenotypic divergences were in general lower (De Maddalena et al., 2007
). The full genome amplification used for sequencing might explain the low genetic diversity of HBV in our study, since this strategy will amplify HBV genomes from sera with high level HBV DNA which consequently are mostly HBsAg positive. The virus in such sera has a low genetic variability because of lack of immune pressure.
It is not clear to which extent the finding of the four HBV subgenotypes A2, D1, D2 and D3 in the territory of the former USSR corresponds to that number of separate introductions of the virus or to which extent it reflects previous population migrations into the area. A2 might be the original subgenotype in the area, while there have been later separate introductions of subgenotypes D3, D2 and D1, in that order. The earlier introduction of D3 is supported by its highest intra-subgenotype divergence, 1.8 %, while the relatively lower level of divergence within subtype D1 is probably due to the more restricted distribution of this subgenotype in the Asian part of the former USSR. However, the finding of almost identical strains within subgenotypes A2 and D2 in patients and blood donors from distantly located regions might suggest more recent common epidemiological links, possibly of iatrogenic nature. This notion was further strengthened by the clades formed by strains with identical S genes; all belonged to a larger clade within D2, where almost all strains derived from the former USSR.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Arauz-Ruiz, P., Norder, H., Robertson, B. H. & Magnius, L. O. (2002). Genotype H: a new Amerindian genotype of hepatitis B virus revealed in Central America. J Gen Virol 83, 2059–2073.
Bancroft, W. H., Mundon, F. K. & Russell, P. K. (1972). Detection of additional antigenic determinants of hepatitis B antigen. J Immunol 109, 842–848.
Banerjee, A., Kurbanov, F., Datta, S., Chandra, P. K., Tanaka, Y., Mizokami, M. & Chakravarty, R. (2006). Phylogenetic relatedness and genetic diversity of hepatitis B virus isolates in Eastern India. J Med Virol 78, 1164–1174.[CrossRef][Medline]
Bozdayi, G., Turkyilmaz, A. R., Idilman, R., Karatayli, E., Rota, S., Yurdaydin, C. & Bozdayi, A. M. (2005). Complete genome sequence and phylogenetic analysis of hepatitis B virus isolated from Turkish patients with chronic HBV infection. J Med Virol 76, 476–481.[CrossRef][Medline]
Couroucé, A. M., Holland, P. V., Muller, J. Y. & Soulier, J. P. (1976). HBsAg subtypes. In Bibliotheca Hematologica, vol. 42, p.133. Basel: Karger.
De Maddalena, C., Giambelli, C., Tanzi, E., Colzani, D., Schiavini, M., Milazzo, L., Bernini, F., Ebranati, E., Cargnel, A. & other authors (2007). High level of genetic heterogeneity in S and P genes of genotype D hepatitis B virus. Virology 365, 113–124.[CrossRef][Medline]
Gunther, S., Li, B. C., Miska, S., Kruger, D. H., Meisel, H. & Will, H. (1995). A novel method for efficient amplification of whole hepatitis B virus genomes permits rapid functional analysis and reveals deletion mutants in immunosuppressed patients. J Virol 69, 5437–5444.[Abstract]
Huy, T. T., Ushijima, H., Sata, T. & Abe, K. (2006). Genomic characterization of HBV genotype F in Bolivia: genotype F subgenotypes correlate with geographic distribution and T (1858) variant. Arch Virol 151, 589–597.[CrossRef][Medline]
Kramvis, A., Weitzmann, L., Owiredu, W. K. & Kew, M. C. (2002). Analysis of the complete genome of subgroup A' hepatitis B virus isolates from South Africa. J Gen Virol 83, 835–839.
Kramvis, A., Restorp, K., Norder, H., Botha, J. F., Magnius, L. O. & Kew, M. C. (2005). Full genome analysis of hepatitis B virus genotype E strains from South-Western Africa and Madagascar reveals low genetic variability. J Med Virol 77, 47–52.[CrossRef][Medline]
Kramvis, A., Arakawa, K., Yu, M. C., Nogueira, R., Stram, D. O. & Kew, M. C. (2008). Relationship of serological subtype, basic core promoter and precore mutations to genotypes/subgenotypes of hepatitis B virus. J Med Virol 80, 27–46.[CrossRef][Medline]
Kurbanov, F., Tanaka, Y., Fujiwara, K., Sugauchi, F., Mbanya, D., Zekeng, L., Ndembi, N., Ngansop, C., Kaptue, L. & other authors (2005). A new subtype (subgenotype) Ac (A3) of hepatitis B virus and recombination between genotypes A and E in Cameroon. J Gen Virol 86, 2047–2056.
Lazarevic, I., Cupic, M., Delic, D., Svirtlih, N. S., Simonovic, J. & Jovanovic, T. (2007). Distribution of HBV genotypes, subgenotypes and HBsAg subtypes among chronically infected patients in Serbia. Arch Virol 152, 2017–2025.[CrossRef][Medline]
Le Bouvier, G. L. (1971). The heterogeneity of Australia antigen. J Infect Dis 123, 671–675.[Medline]
Liu, W. C., Phiet, P. H., Chiang, T. Y., Sun, K. T., Hung, K. H., Young, K. C., Wu, I. C., Cheng, P. N. & Chang, T. T. (2007). Five subgenotypes of hepatitis B virus genotype B with distinct geographic and virological characteristics. Virus Res 129, 212–223.[CrossRef][Medline]
Magnius, L. O., Kaplan, L., Vyas, G. N. & Perkins, H. A. (1975). A new virus-specified determinant of hepatitis B surface antigen. Acta Pathol Microbiol Scand [B] 83, 295–297.[Medline]
Makuwa, M., Souquiere, S., Telfer, P., Apetrei, C., Vray, M., Bedjabaga, I., Mouinga-Ondeme, A., Onanga, R., Marx, P. A. & other authors (2006). Identification of hepatitis B virus subgenotype A3 in rural Gabon. J Med Virol 78, 1175–1184.[CrossRef][Medline]
Michitaka, K., Tanaka, Y., Horiike, N., Duong, T. N., Chen, Y., Matsuura, K., Hiasa, Y., Mizokami, M. & Onji, M. (2006). Tracing the history of hepatitis B virus genotype D in western Japan. J Med Virol 78, 44–52.[CrossRef][Medline]
Naumann, H., Schaefer, S., Yoshida, C. F., Gaspar, A. M., Repp, R. & Gerlich, W. H. (1993). Identification of a new hepatitis B virus (HBV) genotype from Brazil that expresses HBV surface antigen subtype adw4. J Gen Virol 74, 1627–1632.
Norder, H., Couroucé, A. M. & Magnius, L. O. (1992). Molecular basis of hepatitis B virus serotype variations within the four major subtypes. J Gen Virol 73, 3141–3145.
Norder, H., Couroucé, A. M. & Magnius, L. O. (1994). Complete genomes, phylogenetic relatedness, and structural proteins of six strains of the hepatitis B virus, four of which represent two new genotypes. Virology 198, 489–503.[CrossRef][Medline]
Norder, H., Ebert, J. W., Fields, H. A., Mushahwar, I. K. & Magnius, L. O. (1996). Complete sequencing of a gibbon hepatitis B virus genome reveals a unique genotype distantly related to the chimpanzee hepatitis B virus. Virology 218, 214–223.[CrossRef][Medline]
Norder, H., Courouce, A. M., Coursaget, P., Echevarria, J. M., Lee, S. D., Mushahwar, I. K., Robertson, B. H., Locarnini, S. & Magnius, L. O. (2004). Genetic diversity of hepatitis B virus strains derived worldwide: genotypes, subgenotypes, and HBsAg subtypes. Intervirology 47, 289–309.[CrossRef][Medline]
Okamoto, H., Imai, M., Tsuda, F., Tanaka, T., Miyakawa, Y. & Mayumi, M. (1987). Point mutation in the S gene of hepatitis B virus for a d/y or w/r subtypic change in two blood donors carrying a surface antigen of compound subtype adyr or adwr. J Virol 61, 3030–3034.
Okamoto, H., Tsuda, F., Sakugawa, H., Sastrosoewinjo, R., Imai, M., Miyakawa, Y. & Mayumi, M. (1988). Typing hepatitis B virus by homology in nucleotide sequence: comparison of surface antigen subtypes. J Gen Virol 69, 2575–2583.
Sakamoto, T., Tanaka, Y., Orito, E., Co, J., Clavio, J., Sugauchi, F., Ito, K., Ozasa, A. & Quino, A. (2006). Novel subtypes (subgenotypes) of hepatitis B virus genotypes B and C among chronic liver disease patients in the Philippines. J Gen Virol 87, 1873–1882.
Stuyver, L., De Gendt, S., Van Geyt, C., Zoulim, F., Freid, M., Schinazi, R. F. & Rossau, R. (2000). A new genotype of hepatitis B virus: complete genome and phylogenetic relatedness. J Gen Virol 81, 67–74.
Sugauchi, F., Kumada, H., Acharya, S. A., Shrestha, S. M., Gamutan, M. T., Khan, M., Gish, R. G., Tanaka, Y. & Kato, T. (2004). Epidemiological and sequence differences between two subtypes (Ae and Aa) of hepatitis B virus genotype A. J Gen Virol 85, 811–820.
Tallo, T., Norder, H., Tefanova, V., Krispin, T., Priimagi, L., Mukomolov, S., Mikhailov, M. & Magnius, L. O. (2004). Hepatitis B virus genotype D strains from Estonia share sequence similarity with strains from Siberia and may specify ayw4. J Med Virol 74, 221–227.[CrossRef][Medline]
Received 3 December 2007;
accepted 1 April 2008.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |