- Research article
- Open Access
- Open Peer Review
Subgenotype reclassification of genotype B hepatitis B virus
BMC Gastroenterologyvolume 12, Article number: 116 (2012)
Nine subgenotypes from genotype B have been identified for hepatitis B virus (HBV). However, these subgenotypes were less conclusive as they were often designated based on a few representative strains. In addition, subgenotype B6 was designated twice for viruses of different origin.
All complete genome sequences of genotype B HBV were phylogenetically analyzed. Sequence divergences between different potential subgenotypes were also assessed.
Both phylogenetic and sequence divergence analyses supported the designation of subgenotypes B1, B2, B4, and B6 (from Arctic). However, sequence divergences between previously designated B3, B5, B7, B8, B9 and another B6 (from China) were mostly less than 4%. In addition, subgenotype B3 did not form a monophyly.
Current evidence failed to classify original B5, B7, B8, B9, and B6 (from China) as subgenotypes. Instead, they could be considered as a quasi-subgenotype B3 of Southeast Asian and Chinese origin. In addition, previously designated B6 (from Arctic) should be renamed as B5 for continuous numbering. This novel classification is well supported by both the phylogeny and sequence divergence of > 4%.
By comparing 18 HBV genomes, Okamoto et al. proposed that sequence divergence of > 8% over the entire genome should be used for HBV genotyping in 1988 . Based on this, four genotypes were identified as A, B, C, and D. So far, at least eight genotypes, from A to H, have been identified and widely accepted. Genotypes B and C are prevalent in Asia and co-infection and/or super-infection lead to frequent recombination between these two genotypes [2, 3].
Based on the rule that different subgenotypes should diverge by at least 4% over the genome , genotypes A , B , C , D , and F  have been reported to have evolved into various subgenotypes . Genotype B was initially divided into two subgenotypes: Bj (j for Japan) and Ba (a for Asia) . Viruses of subgenotype Bj are not recombinants, while those of subgenotype Ba are B/C recombinants, with their preC-C genes coming from genotype C . In 2004, Bj was renamed as B1 and Ba was renamed as B2 . In this report, subgenotypes B3 and B4 were also described. Subgenotype B3 was composed of four strains from Indonesia, while subgenotype B4 mostly comprised strains from Vietnam and France . Subgenotype B5 was initially reported in 2006 from the Philippines . A few months later, Sakamoto et al. found a few viruses from the Philippines that differed from subgenotypes B1 to B4, and also designated it as a novel subgenotype, B5 . Subgenotype B6 was identified in 2007 from arctic indigenous populations . In particular, just as with subgenotype B1, viruses of B6 were not recombinants. This observation led Sakamoto et al. to classify genotype B into two types, non-recombinant (B1 and B6) and recombinant (B2 to B5) . Viruses of B7 were isolated from the Nusa Tenggara islands in Eastern Indonesia from 2008 . Subsequently, subgenotype B8 was also identified in Indonesia by analyzing a large cohort of patient samples . More recently, viruses of subgenotype B9 were isolated from the same region in Indonesia where subgenotype B8 was isolated . This suggested that the distribution of HBV subgenotypes might relate to the ethnic origin of the infected patients .
However, previous studies just used a number of selected representative strains to perform phylogenetic analysis and estimate sequence divergence. Hence, the designation of novel HBV subgenotypes was less conclusive. In this study, we reanalyzed all of the full-length genome sequences of genotype B using phylogenetic analysis. This analysis led to the proposal of a novel and consistent classification for HBV of genotype B.
We have previously analyzed a total of 3471 full length genome sequences of HBV using a phylogenetic approach . Our results showed that 860 sequences belonged to genotype B. These sequences were selected to generate a new dataset for further analysis. In addition, a sequence of genotype C (GenBank:EU939604) was also included in the new dataset and used as an outgroup. Information of these sequences, such as subgenotype and recombination, was also extracted from GenBank annotations. For sequences with a reference available in Pubmed, we performed an extensive literature review to obtain their subgenotype and recombination information. This information was used to define the subgenotypes.
Phylogenetic analysis of the genotype B strains was carried out using RAxML  under the GTRCAT approximation  and random starting trees. Three thousand rapid bootstrap replicates were performed and all other parameters were set to default. Trees were visualized and analyzed using Dendroscope . The tree is available as Additional file 1.
The mean nucleotide divergence (mean ± SD) between different subgenotypes was calculated using Mega 5  under the Kimura 2-parameter model . Five hundred bootstrap replicates were applied to obtain consistent and reliable sequence divergence values.
Phylogenetic analysis of all HBV genotype B full-length genome sequences revealed five distinct major clusters. They were named as cluster 1 to 5 from top to bottom (Figure 1). Most of the clusters were well supported with high bootstrap values. Cluster 1 was composed of 27 viruses belonging to subgenotype B6 , with a bootstrap value of 100%. 37 viruses of cluster 2 were isolated from Japan and were designated as subgenotype B1 [2, 11]. The bootstrap value for this cluster was 92%. 46 viruses of subgenotype B4  constituted cluster 4, with a bootstrap value of 98%.
Phylogenetic analysis grouped 655 subgenotype B2 viruses into cluster 5 [2, 11], representing approximately 76% of all genotype B strains. All viruses of this cluster were B/C recombinants, with the great majority of the strains isolated from Asia (Additional file 1). However, the bootstrap value for this cluster was only 46%.
Apart from the sequence divergence between B4 and B2 (3.9 ± 0.3%), sequence divergences between any other two subgenotypes (B6, B1, B4 and B2) were greater than 4% (Table 1). Furthermore, the within subgenotype divergences of these four subgenotypes were below 4% (Table 1).
Cluster 3 was also well supported with a bootstrap value of 72%. Unlike the described clusters composed of a single subgenotype, cluster 3 included several previously reported subgenotypes, such as B5, B3, B9, B8 and B7 (Figure 1). In addition, a few Chinese strains named as B6  were included into cluster 3. Due to the existence of another subgenotype B6 , the Chinese strains were tentatively renamed as B6’. Apart from the previously designated subgenotypes, two sequences isolated from China (GenBank:AY800391 and GenBank:AY800392) have not been designated a subgenotype and were tentatively named Bx. In particular, viruses of subgenotype B3 did not form a monophyly. Instead, subgenotype B3 viruses were scattered in three parts, B3-1, B3-2 and B3-3 (Figure 1).
To resolve the aforementioned problems, we calculated sequence divergences between B5, B3-1, B3-2, B3-3, B9, B8, B7, B6’ and Bx (Table 2). Apart from Bx, sequence divergence between any two of these potential subgenotypes was usually below 4% (Table 2). Although sequence divergences between Bx and other potential subgenotypes were mostly greater than 4%, the phylogeny did not support Bx to be a novel subgenotype in that it did not form a monophyly compared to other potential subgenotypes in this cluster (Figure 1). Therefore, sequence divergences did not support the designation of subgenotypes B5, B9, B8, B7, B6’, and Bx.
Previous studies have shown that different HBV genotypes and subgenotypes may cause differences in disease progression, response to anti-viral treatment regimens and in clinical outcomes [4, 10, 23]. Therefore, the accurate classification of HBV is important. To resolve the problems in HBV subgenotyping and to propose a consistent classification for genotype B HBV, we analyzed 860 complete genome sequences of genotype B using phylogenetic analysis.
Phylogenetic analysis showed that genotype B HBV has evolved into five major clusters. Four of them corresponded to B6 (cluster 1), B1 (cluster 2), B4 (cluster 4) and B2 (cluster 5). Apart from subgenotype B2, the remaining three subgenotypes were well supported by both high bootstrap values and sequence divergences of > 4%. Therefore, subgenotypes B6, B1 and B4 were properly designated. Although the bootstrap value for subgenotype B2 was only 46% and the sequence divergence between subgenotypes B4 and B2 was less than 4%, since subgenotype B2 has been widely accepted, it should be maintained in order to avoid more confusion.
However, our results failed to support the designation of subgenotypes B5, B3, B9, B8, B7 and B6’. First, subgenotype B3 was not a monophyly, but scattered in the tree. This was against the rule that an HBV genotype and subgenotype should be a monophyly . Second, sequence divergences between the above subgenotypes were mostly less than 4%. This also did not support them to be separate subgenotypes.
Instead, considering that viruses in this cluster are all isolated from Southeast Asia and China and this cluster is well supported by a high bootstrap value of 72%, we proposed that a quasi-subgenotype B3 should be used for cluster 3 according to previous reports, in which the definition of quasi-subgenotype has been used to resolve the inconsistency in the subgenotyping of genotype A HBV [24, 25]. Sequence divergences between the quasi-subgenotype B3 and B6, B1, B4 and B2 were 6.1 ± 0.4%, 6.8 ± 0.3%, 5.1 ± 0.3% and 5.1 ± 0.3%, respectively (Table 3). In addition, the within subgenotype divergence of quasi-subgenotype B3 was 3.1 ± 0.1%, less than 4% to differentiate a subgenotype.
Taken together, we corrected the incongruence in the classification of genotype B HBV and proposed a consistent classification for genotype B based on a phylogenetic analysis of all genotype B HBV complete genome sequences (Figure 2). In this classification, originally designated B3, B5, B7, B8, B9 and B6’ (sequences from China) comprised a quasi-subgenotype B3. For continuous numbering, the subgenotype B6 has been renamed as new B5. This classification is well supported by both phylogenetic analysis and sequence divergences. To avoid misclassification of HBV subgenotypes, we insist that the rules [24, 26] proposed to define novel genotypes and subgenotypes should be strictly obeyed. In addition, if possible, novel genotype/subgenotype should be made following a complete comparison of all relevant sequences rather than with a few representative sequences.
Our results revealed that subgenotypes B1, B2, B4, and B6 (from Arctic) have been properly designated and should be maintained in the new classification. However, sequence divergences between previously designated B3, B5, B7, B8, B9 and another B6 (from China) were mostly less than 4%. In addition, subgenotype B3 did not form a monophyly. Therefore, current evidence failed to classify original B5, B7, B8, B9, and B6 (from China) as subgenotypes. Instead, they could be considered as a quasi-subgenotype B3 of Southeast Asian and Chinese origin. Moreover, previously designated B6 (from Arctic) should be renamed as B5 for continuous numbering. This novel classification is well supported by both the phylogeny and sequence divergence of > 4%.
Hepatitis B virus.
Okamoto H, Tsuda F, Sakugawa H, et al: Typing hepatitis B virus by homology in nucleotide sequence: comparison of surface antigen subtypes. J Gen Virol. 1988, 69: 2575-2583. 10.1099/0022-1317-69-10-2575.
Sugauchi F, Orito E, Ichida T, et al: Hepatitis B virus of genotype B with or without recombination with genotype C over the precore region plus the core gene. J Virol. 2002, 76: 5985-5992. 10.1128/JVI.76.12.5985-5992.2002.
Sakamoto T, Tanaka Y, Simonetti J, et al: Classification of hepatitis B virus genotype B into 2 major types based on characterization of a novel subgenotype in Arctic indigenous populations. J Infect Dis. 2007, 196: 1487-1492. 10.1086/523111.
Kramvis A, Kew MC: Relationship of genotypes of hepatitis B virus to mutations, disease progression and response to antiviral therapy. J Viral Hepat. 2005, 12: 456-464. 10.1111/j.1365-2893.2005.00624.x.
Hubschen JM, Mbah PO, Forbi JC, et al: Detection of a new subgenotype of hepatitis B virus genotype A in Cameroon but not in neighbouring Nigeria. Clin Microbiol Infect. 2011, 17: 88-94. 10.1111/j.1469-0691.2010.03205.x.
Thedja MD, Muljono DH, Nurainy N, Sukowati CH, Verhoef J, Marzuki S: Ethnogeographical structure of hepatitis B virus genotype distribution in Indonesia and discovery of a new subgenotype, B9. Arch Virol. 2011, 156: 855-868. 10.1007/s00705-011-0926-y.
Mulyanto A, Depamede SN, Wahyono A, et al: Analysis of the full-length genomes of novel hepatitis B virus subgenotypes C11 and C12 in Papua, Indonesia. J Med Virol. 2011, 83: 54-64. 10.1002/jmv.21931.
Abdou Chekaraou M, Brichler S, Mansour W, et al: A novel hepatitis B virus (HBV) subgenotype D (D8) strain, resulting from recombination between genotypes D and E, is circulating in Niger along with HBV/E strains. J Gen Virol. 2010, 91: 1609-1620. 10.1099/vir.0.018127-0.
Huy TTT, Ushijima H, Sata T, Abe K: Genomic characterization of HBV genotype F in Bolivia: genotype F subgenotypes correlate with geographic distribution and T(1858) variant. Arch Virol. 2006, 151: 589-597. 10.1007/s00705-005-0671-1.
Cao GW: Clinical relevance and public health significance of hepatitis B virus genomic variations. World J Gastroenter. 2009, 15: 5761-5769. 10.3748/wjg.15.5761.
Norder H, Courouce AM, Coursaget P, et al: Genetic diversity of hepatitis B virus strains derived worldwide: genotypes, subgenotypes, and HBsAg subtypes. Intervirology. 2004, 47: 289-309. 10.1159/000080872.
Nagasaki F, Niitsuma H, Cervantes JG, et al: Analysis of the entire nucleotide sequence of hepatitis B virus genotype B in the Philippines reveals a new subgenotype of genotype B. J Gen Virol. 2006, 87: 1175-1180. 10.1099/vir.0.81525-0.
Sakamoto T, Tanaka Y, Orito E, et al: Novel subtypes (subgenotypes) of hepatitis B virus genotypes B and C among chronic liver disease patients in the Philippines. J Gen Virol. 2006, 87: 1873-1882. 10.1099/vir.0.81714-0.
Nurainy N, Muljono DH, Sudoyo H, Marzuki S: Genetic study of hepatitis B virus in Indonesia reveals a new subgenotype of genotype B in east Nusa Tenggara. Arch Virol. 2008, 153: 1057-1065. 10.1007/s00705-008-0092-z.
Mulyanto A, Depamede SN, Surayah K, et al: A nationwide molecular epidemiological study on hepatitis B virus in Indonesia: identification of two novel subgenotypes, B8 and C7. Arch Virol. 2009, 154: 1047-1059. 10.1007/s00705-009-0406-9.
Shi W, Carr MJ, Dunford L, Zhu CD, Hall WW, Higgins DG: Identification of novel inter-genotypic recombinants of human hepatitis b viruses by large-scale phylogenetic analysis. Virology. 2012, 427: 51-59. 10.1016/j.virol.2012.01.030.
Stamatakis A, Ludwig T, Meier H: RAxML-III: a fast program for maximum likelihood-based inference of large phylogenetic trees. Bioinformatics. 2005, 21: 456-463. 10.1093/bioinformatics/bti191.
Stamatakis A: Phylogenetic models of rate heterogeneity: a high performance computing perspective. 2006, Rhodos, Greece: Proceedings of 20th IEEE/ACM International Parallel and Distributed Processing Symposium (IPDPS2006), High Performance Computational Biology Workshop
Huson DH, Richter DC, Rausch C, Dezulian T, Franz M, Rupp R: Dendroscope: An interactive viewer for large phylogenetic trees. BMC Bioinformatics. 2007, 8: 460-10.1186/1471-2105-8-460.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.
Shen T, Gao JM, Zou YL, Dong H, Yan XM: Novel hepatitis B virus subgenotype in the southern Yunnan Province of China. Intervirology. 2009, 52: 340-346. 10.1159/000252912.
Schaefer S: Hepatitis B, virus: significance of genotypes. J Viral Hepat. 2005, 12: 111-124. 10.1111/j.1365-2893.2005.00584.x.
Pourkarim MR, Amini-Bavil-Olyaee S, Lemey P, Maes P, Van Ranst M: Are hepatitis B virus "subgenotypes" defined accurately?. J Clin Virol. 2010, 47: 356-360. 10.1016/j.jcv.2010.01.015.
Pourkarim MR, Amini-Bavil-Olyaee S, Lemey P, Maes P, Van Ranst M: HBV subgenotype misclassification expands quasi-subgenotype A3. Clin Microbiol Infect. 2011, 17: 947-949. 10.1111/j.1469-0691.2010.03374.x.
Schaefer S, Magnius L, Norder H: Under construction: classification of hepatitis B virus genotypes and subgenotypes. Intervirology. 2009, 52: 323-325. 10.1159/000242353.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-230X/12/116/prepub
This study was partly supported by Science Foundation Ireland (PI grant 07/IN.1/B1783). We thank Dr. XY Lang and Dr. XN Wang in the Supercomputing Center of Chinese Academy of Sciences (CAS) for their help in installing and optimizing RAxML on the SCIGRID.
The authors declare that they have no competing interests.
Conceived and designed the experiments: WS, ZZ. Performed the experiments: WS, CZ, WZ. Analyzed the data: WS. Wrote the paper: WS. Revised the paper: MJC, DGH and ZZ. All the authors read and approved the final manuscript.