J Gen Virol Email Content Delivery
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 89 (2008), 1829-1839; DOI 10.1099/vir.0.83660-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tallo, T.
Right arrow Articles by Norder, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tallo, T.
Right arrow Articles by Norder, H.
Agricola
Right arrow Articles by Tallo, T.
Right arrow Articles by Norder, H.

D2: major subgenotype of hepatitis B virus in Russia and the Baltic region

Tatjana Tallo1,2,3, Valentina Tefanova1, Ljudmilla Priimägi1, Jelena Schmidt4, Olga Katargina1, Michail Michailov5, Sergey Mukomolov6, Lars Magnius2,3 and Heléne Norder2,3

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
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Complete or almost complete hepatitis B virus (HBV) genomes were sequenced for 13 genotype A and 42 genotype D strains from the former USSR. The strains were classifiable within subgenotypes A2, D1, D2 and D3. Comparison of the deduced gene products for the four ORFs of 89 genotype D strains revealed 27 subgenotype-specific residues, and a region spanning residues 58–128 in the spacer region of the P gene could be used to distinguish between D1 and D4. This enabled the allocation to subgenotype of strains with partially sequenced genomes. D2 was dominating, while D3 was found in low frequency in the whole region. D1 was most prevalent in the Middle Asian Republics. Mean inter-subgenotype divergences between D1 and D2, D1 and D3 and D2 and D3 were 2.7, 3.4 and 3.4 %, respectively. The intra-subgenotype divergence was 0.4, 1.1, 1.0 and 1.8 % for A2, D1, D2 and D3, respectively. All D1 and D3 strains encoded subtype ayw2, whereas most D2 strains encoded ayw3. Two D2 strains encoded ayw4. Strains with identical S genes were closely related at the level of complete genomes and formed geographically specific clades with low intraclade divergences, possibly indicating past iatrogenic spread. It is not clear whether the finding of four subgenotypes in the area corresponds to separate introductions of the virus or to previous population migrations into the area. An earlier introduction of D3 compared with D2 was supported by its higher intra-subgenotype divergence, while the lower divergence within D1 is probably due to a more recent emergence.

The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper are EU594382-EU594436.


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Hepatitis B virus (HBV) is a partially double-stranded DNA virus of the family Hepadnaviridae with a genome of 3.2 kb containing four partially overlapping open reading frames designated S, C, P and X. Based on serological heterogeneity of the small S gene product, the hepatitis B surface antigen (HBsAg) has been classified into nine different subtypes (Le Bouvier, 1971Down; Bancroft et al., 1972Down; Magnius et al., 1975Down; Couroucé et al., 1976Down). This heterogeneity is based on single amino acid substitutions of the HBsAg with the most important residues being 122, 127 and 160 (Okamoto et al., 1987Down; Norder et al., 1992Down). On the basis of genetic variability of complete genomes, HBV was more recently classified into genotypes designated A–H (Okamoto et al., 1988Down; Naumann et al., 1993Down; Norder et al., 1994Down; Stuyver et al., 2000Down; Arauz-Ruiz et al., 2002Down). Within the genotypes, a number of subgenotypes designated with arabic numerals are now described (Kramvis et al., 2002Down, 2008Down; Norder et al., 2004Down; Sugauchi et al., 2004Down; Kurbanov et al., 2005Down; Banerjee et al., 2006Down; Huy et al., 2006Down; Sakamoto et al., 2006Down).

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., 2002Down; Norder et al., 2004Down; Sugauchi et al., 2004Down). Recently A3 was described in native populations of West and Central Africa (Kurbanov et al., 2005Down; Makuwa et al., 2006Down). Definite geographical predilections have also been shown for genotypes B and C (Norder et al., 2004Down; Sakamoto et al., 2006Down; Liu et al., 2007Down). Genotype D is the most widespread and divides into five subgenotypes (Norder et al., 2004Down; Banerjee et al., 2006Down), 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., 2004Down). 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
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Sera selection.
One hundred and twenty six HBV strains, 26 genotype A and 100 genotype D strains obtained from patients and blood donors, were selected for sequence analyses of complete HBV genomes from the 205 strains previously published (Tallo et al., 2004Down). Fifteen of the selected genotype A strains had identical S gene sequences and 59 selected genotype D strains belonged to seven clades, each formed by strains with identical small S gene sequences (Table 1Down; Tallo et al., 2004Down). Criteria for selection were inclusion of strains with identical and divergent small S genes encoding subtypes ayw2, ayw3, ayw4 and adw2 originating from European, Siberian, Middle Asia and Far-Eastern regions of the former USSR (Table 1Down).


View this table:
[in this window]
[in a new window]

 
Table 1. Number of complete HBV genomes obtained after selection of strains from eight clades characterized by identical small S gene sequences and other representative strains from Tallo et al. (2004)Down

 
DNA extraction, polymerase chain reaction and sequencing.
HBV DNA was extracted by the standard phenol-chloroform method. Briefly, 50 µl serum was incubated for 2 h at 42 °C with 250 mg proteinase K ml–1 (Sigma), 1 % sodium dodecyl sulphate, 2.5 mM disodium EDTA (Sigma), 25 mM sodium acetate and 0.25 µg tRNA ml–1. The DNA was subsequently extracted with phenol and phenol : chloroform, 1 : 1 mixture, precipitated in absolute ethanol with 0.15 M potassium acetate and dissolved in 40 µl MilliQ H2O.

A modification of the technique and primers described by Gunther et al. (1995)Down 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., 2005Down) 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., 1996Down; Kramvis et al., 2005Down), 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 2Down). 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.


View this table:
[in this window]
[in a new window]

 
Table 2. Wild-type HBV complete genomes used for analysis of signature motifs of HBV subgenotypes

 
Amino acid analysis.
The 80 genotype D strains without pre-core stop or core promoter mutations were divided into subgenotypes according to the phylogenetic analysis of the complete genomes. The deduced amino acid sequences of each of the four ORFs were aligned and compared within each subgenotype. When more than 90 % of the strains within a subgenotype expressed the same amino acid residue at a specific position, this residue was considered a consensus residue for this subgenotype at that position.


   RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
The genomes could be amplified for 19 of the 26 genotype A and for 64 of the 100 genotype D strains selected, although complete genomes were only obtained for 46 of these (Table 1Up). The complete or almost complete P genes could be sequenced for another nine strains, which could be used for subgenotyping of these strains. Two of these had identical S genes encoding ayw4 within clade III (Tallo et al., 2004Down; Table 1Up).

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. 1aDown). 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 %).


Figure 1
View larger version (44K):
[in this window]
[in a new window]

 
Fig. 1. Phylogenetic trees obtained by the neighbour-joining method (NJ), based on 123 complete HBV genomes representing (a) the tree formed by 34 genotype A strains and (b) the tree formed by 89 genotype D strains. Strains from this study are shown in bold. Sequences of 77 HBV isolates from GenBank are indicated by their accession numbers. The designation and origin of each strain is indicated at the nodes of the branches. Bootstrap values higher than 75 % obtained from 1000 replicates are given at the internal nodes. Strains with identical small S genes are indicated with * for genotype A in Fig. 1(a)Up and with the clade designation according to Tallo et al. (2004)Down for genotype D strains in Fig. 1(b)Up.

 
The 33 complete genotype D genomes segregated into three clades each with strains previously shown to belong to subgenotypes D1, D2 and D3 (Fig. 1bUp). The majority, 27 (82 %), of these strains could be classified as D2; all encoded ayw3. One strain each from Uzbekistan and Kazakhstan belonged to D1, and four belonged to D3. Three of these D3 strains derived from the European part of the region and encoded ayw2 (Table 3Down, Fig. 1bUp). The mean intra-group divergence for the strains studied within subgenotypes D1, D2 and D3 were 1.1, 1.0 and 1.8 %, respectively.


View this table:
[in this window]
[in a new window]

 
Table 3. Number of strains with complete genomes with indicated subgenotype of HBV by area of origin [within parentheses databased on small S gene sequences according to Tallo et al. (2004)Down]

One genotype C strain and one genotype D strain both from Khabarovsk and four genotype D strains with undetermined subtypes, originating from Estonia, the Volga region, Kazakhstan and Chita, are not included in the table. The complete P sequence was obtained for seven strains. For two strains encoding ayw4 in clade III, 2871 nucleotides could be used for subgenotyping.

 
All genotype D strains with identical S gene sequences could be distinguished, although they were still closely related, when complete genomes were compared (Fig. 1bUp). The strains corresponding to S gene clades I and V were both primitive as compared with II and IV, that seem to have evolved from the former. Genotype A strains with identical S genes were also found to be closely related when their complete genomes were compared. However, one strain each from Germany and the USA were also related to these strains. Two strains, 4600-97 and 4625-97, from blood donors in Yakutsk, remained identical even when complete genomes were compared (Fig. 1aUp).

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-dDown). The X gene for subgenotype D1 could not be distinguished from that of D2 (Fig. 2dDown), although the other subgenotypes, D3–D5, also showed distinct clustering in this region (Fig. 2dDown).


Figure 2
View larger version (37K):
[in this window]
[in a new window]

 
Fig. 2. NJ trees based on the phylogenetic analysis of four open reading frames of 89 HBV isolates belonging to genotype D (a) for polymerase, positions 2307–3182, 1–1623, (b) for preS1/S2/S gene, positions 2848–3182, 1–835, (c) for preC/C region of genome, positions 1814–2452, and (d) for the X gene, positions 1374–1838. Enumeration is given according to strain Ehime12 with GenBank accession no. AB110075. Designations of the strains are indicated at the nodes of the branches. Bootstrap values higher than 75 % obtained from 1000 replicates are given at the internal nodes.

 
PreS/S, P, PreC/C and X proteins
The deduced primary structures for the gene products of all four ORFs were compared for subgenotypes D1–D5, and the amino acid residues expressed by more than 90 % of the D1–D4 strains and those specific for the D5 strain were identified (Table 4Down). These signature amino acid residues were also compared with the residues at the corresponding positions for the other genotypes (Table 4Down).


View this table:
[in this window]
[in a new window]

 
Table 4. Conserved amino acid residues in the polymerase, preS/S, preC/C and X genes in 80 wild-type genotype D strains, 27 within D1, 36 within D2, 13 within D3, 3 within D4 and one D5 strain

Substitutions present in more than 90 % of the strains within each subgenotype are shown. Subgenotype-specific substitutions are boxed, genotype-specific amino acid residues are shown in bold italic. Subgenotype-specific substitutions shared by strains belonging to other genotypes are shown in bold and underlined. The specific residues for subgenotype D5 are not marked, since there is only one complete genome available of wild-type D5.

tp, terminal protein; sp, spacer; rt, reverse transcriptase; rh, ribonuclease (RNase) H.

 
Most signature motif residues, 25, distinguishing one subgenotype from the other were found in the P gene product, when the D5-specific residues were also included. Ten of these were unique, i.e. not shared with any other genotype or D subgenotype (Table 4Up). There were another nine genotype D specific residues that were not shared by all D subgenotypes. In addition, there were 15 subgenotype-specific residues that were also present in strains belonging to other genotypes (Table 4Up). In all, D2 strains had two specific residues, D3 strains five, D4 strains ten and the D5 strain had seven unique residues in the P gene (Table 4Up). By combining two or three residues in this region the D1 subgenotype could thus be assessed. Two residues, Asp32 and Glu178, in the terminal protein region (tp) were found in the deduce amino acid sequence in the D2 strains from this study. Interestingly, subgenotype D1–D3 strains had Phe104 in the spacer region (sp) and Tyr256 in the reverse transcriptase region (rt), whereas D4 and D5 expressed Ser104 and Trp256 found in the strains of the other genotypes. The deduced amino acid sequence of the seven genotype D strains, for which only the P gene was sequenced, revealed that all seven strains had the three unique and the two specific D2 residues. For two ayw4 strains with incompletely sequenced P gene, due to gaps from residue 58 in the sp region to residue 53 in the rt region, one unique subgenotype D2 residue was identified, Arg at position rt125.

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 4Up). One signature residue each was identified for D1, D2 and D3, two for D4 and six for D5 (Table 4Up). 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 4Up). 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 4Up). 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
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
The amplification strategy using the primers P1/P2 for obtaining full genomes will mainly amplify genomes from HBsAg-positive sera, which are optimal for phylogenetic analyses. This strategy will mostly not allow the amplification of genomes from sera with low concentration of HBV DNA. In agreement with this just a few full genomes had promoter or pre-core mutations. However, it was possible to obtain complete genomes for 46 of the 126 selected strains. There were regions encompassing the core promoter and pre-core regions that could not be sequenced for a number of strains, possibly due to highly structured regions within this region. Interestingly, for strains encoding ayw4 regions upstream, the preS region was also difficult to sequence.

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., 2002Down; Norder et al., 2004Down). 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., 2006Down). 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., 2006Down). 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. 2005Down; Bozdayi et al., 2005Down). 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., 2004Down; De Maddalena et al., 2007Down). 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., 2007Down; Lazarevic et al., 2007Down).

When comparing the dendrogram based on the small S gene (Tallo et al., 2004Down) with that based on complete genomes, it was possible to classify all strains from the major cluster encoding ayw3 or ayw4 within D2 (Table 4Up). 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., 2004Down).

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., 2004Down). The strains with identical S genes often originated from regions distantly located from each other (Tallo et al., 2004Down). 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., 2006Down) and 1.6 % for subgenotype C5 strains in the Philippines (Sakamoto et al., 2006Down). 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., 2007Down). 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
 
The study was supported by grants from the New Visby Program and from the Swedish Institute by grant no. 01543/2006, from the Swedish Research Council grant VR521-2006-2753 (K2007-58X-20363-01-3), and partially by grant no. 5961 from the Estonian Science Foundation.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Amini-Bavil-Olyaee, S., Sarrami-Forooshani, R., Adeli, A., Sabahi, F., Abachi, M., Azizi, M. & Mahboudi, F. (2005). Complete genomic sequence and phylogenetic relatedness of hepatitis B virus isolates from Iran. J Med Virol 76, 318–326.[CrossRef][Medline]

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.[Abstract/Free Full Text]

Bancroft, W. H., Mundon, F. K. & Russell, P. K. (1972). Detection of additional antigenic determinants of hepatitis B antigen. J Immunol 109, 842–848.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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.



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tallo, T.
Right arrow Articles by Norder, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tallo, T.
Right arrow Articles by Norder, H.
Agricola
Right arrow Articles by Tallo, T.
Right arrow Articles by Norder, H.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS