Skip to main content
Download PDF

Complete genome analysis of hepatitis B virus in Qinghai-Tibet plateau: the geographical distribution, genetic diversity, and co-existence of HBsAg and anti-HBs antibodies

Virology Journal volume 17, Article number: 75 (2020) Cite this article

Abstract

Background

The genetic variation and origin of Hepatitis B Virus (HBV) in Qinghai-Tibet Plateau were poorly studied. The coexistence of HBsAg and anti-HBs has been described as a puzzle and has never been reported in the indigenous population or in recombinant HBV sequences. This study aimed to report geographical distribution, genetic variability and seroepidemiology of HBV in southwest China.

Methods

During 2014–2017, 1263 HBsAg positive serum were identified and 183 complete genome sequences were obtained. Serum samples were collected from community-based populations by a multistage random sampling method. Polymerase chain reaction (PCR) was used to amplify the HBV complete genome sequences. Then recombination, genetic variability, and serological analysis were performed.

Results

(1) Of the 1263 HBsAg positive serum samples, there were significant differences between the distribution of seromarkers in Tibet and Qinghai. (2) Of 183 complete genome sequences, there were 130 HBV/CD1 (71.0%), 49 HBV/CD2 (26.8%) and four HBV/C2 isolates (2.2%). Serotype ayw2 (96.1%) was the main serological subtype. (3) Several nucleotide mutations were dramatically different in CD1 and CD2 sequences. Clinical prognosis-related genetic variations such as nucleotide mutation T1762/A1764 (27.93%), A2189C (12.85%), G1613A (8.94%), T1753C (8.38%), T53C (4.47%) T3098C (1.68%) and PreS deletion (2.23%) were detected in CD recombinants. (4) From the inner land of China to the northeast boundary of India, different geographical distributions between CD1 and CD2 were identified. (5) Twenty-seven (2.14%) HBsAg/HBsAb coexistence serum samples were identified. S protein amino acid mutation and PreS deletion were with significant differences between HBsAg/HBsAb coexistence group and control group.

Conclusions

HBV/CD may have a mixed China and South Asia origin. Based on genetic variations, the clinical prognosis of CD recombinant seems more temperate than genotype C strains in China. The HBsAg/HBsAb coexistence is a result of both PreS deletion and aa variation in S protein. Several unique mutations were frequently detected in HBV/CD isolates, which could potentially influence the clinical prognosis.

Background

Hepatitis B Virus (HBV) is considered to be a major global health problem with more than 250 million chronic HBV (CHB) carriers and more than one million HBV-associated human deaths per year [1]. In the nationwide investigation in 2006, Hepatitis B surface antigen (HBsAg) was identified in 7.2% of the whole population in China [2]. In Qinghai-Tibet Plateau, due to high altitude environment, religion issue, delayed vaccine inoculation or other unknown reasons, the HBV prevalence is over 10% according to our recent study (unpublished), which is much higher than in other areas of China. These carriers of HBV are at increased risk of developing liver cirrhosis (LC) and hepatocellular carcinoma (HCC) [3].

Qinghai-Tibet Plateau covers more than 2.5 million square kilometers, which is the second- largest plateau in the world and connecting south and northeast Asian [4]. In this area, a special recombinant of HBV genotype C and genotype D (HBV/CD) was reported [4]. However, the nature and origin of the recombinant are still poorly studied. At the same time, the coexistence of HBsAg and anti-HBs (HBsAb) in CHB patients has been occasionally reported [5, 6]. Due to the low incidence of HBsAg/ HBsAb coexistence in HBV carriers, the sample sizes in many previous studies [5, 7] were small (less than 20) and lack of adequate control subjects [8]. The mechanism underlying the emergence of anti-HBs in CHB patients remains unclear and related information has not been reported in the indigenous population or in recombinant HBV sequences.

To gain deeper insights into HBV genomic diversity of the special recombinants, 1263 HBsAg positive serum samples were obtained from eight areas of Qinghai-Tibet Plateau and 183 complete genome sequences of HBV isolates were under further analysis. Twenty-seven HBsAg/HBsAb coexistence serum samples were identified and studied with genome variation.

Methods

Procedures to detect point mutations and recombination of HBV by PCR and direct sequencing analysis

Sera samples from Qinghai-Tibet plateau → HBV infection serological markers detection→ DNA purification → HBV complete genome amplification → PCR-DNA purification → Automatic sequencing → Sequences assembling → Comparing with reference sequences.

Subject sample collection

During 2014–2017, the subjects were collected from community-based populations from eight regions of Qinghai-Tibet Plateau based on the population density and covered most of the indigenous habitations, including Hainan, Lhasa, Shannan, Nyingchi, Ali, Nakqu, Chamdo and Rikaze. In this study, the multistage random sampling method was used to ensure the sample representativeness of the whole area. Firstly, two or three counties were selected at random from eight areas, respectively. Secondly, two villages were selected from every county. Thirdly, populations of 18–59 years old were selected from every village. Basic information was recorded on the questionnaire prepared beforehand, including name, gender, birth date, address, and medical information, and 5 mL of venous blood was taken from each participant. HBV infection markers including hepatitis B surface antigen (HbsAg), anti-hepatitis B surface antibody (HBsAb), anti-hepatitis B core antibody (HBcAb), hepatitis B e antigen (HBeAg) and anti-hepatitis B e antibody (HBeAb) were detected by chemiluminescent assays (AXSYM; Abbott Laboratories, North Chicago, IL, USA).

HBV DNA extraction, whole-genome amplification, and sequencing

HBV DNA was extracted from 200 μL serum using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany). Full-length HBV DNA (about 3.2 kb) was amplified by Nested PCR finally performed in seven fragments. Based on the previously reported methods [9] used in our preliminary study [10], the primers and thermal profile were optimized in this study to adapt to HBV/CD isolates and identify more complete genomes. The basal core promoter (BCP) region was generated by two rounds. First round was conducted using the primer combination of BcpF1 and BcpR1 in a 25 μL reaction volume containing 5 μL extracted DNA and 12.5 μL premix Taq polymerase. The whole- genome (without BCP region) consisting of six fragments was also generated by two rounds of PCR. First round was conducted using the primer combination of HBV1799FLong and HBV1801RLong in a 50 μL reaction volume containing 15 μL extracted DNA and 25 μL premix Taq polymerase. All the primers and thermal profiles were listed in Table 1.

Table 1 List of primers used to amplify the different regions of the HBV genome and their respective thermal profile
Full size table

After purification of the PCR products with a QIA Gel Extraction Kit (Qiagen, Valencia, CA), the sequences were determined using the Sanger dideoxy terminator sequencing method with DNA sequencer ABI 3700 (PE Applied Biosystem). Then the sequences were assembled using SeqMan II software (DNAStar Inc.), and the correct nucleotide position of the complete HBV sequence was revised through alignment with the reference sequence in the GenBank database.

The whole-genome nucleotide sequences reported in this article have been deposited in the National Center for Biotechnology Information GenBank database with accession numbers MN683570- MN683729, MN657315- MN657318, KX660674- KX660690.

Recombination and serological analysis

Genotypes and recombination in genome sequences of HBV were investigated using SimPlot v3.5.1 software and JPHMM (jumping profile Hidden Markov Model) method [11]. MEGA 7.0 software was used to calculate genetic distance and pairwise distance comparisons. The HBsAg serotypes were deduced from the sequences of the S gene region by the identification of amino acids at positions 122 (Lys-Arg for d-y determinants), 160 (Lys-Arg for w-r determinants), 127 (Pro-Thr-Leu/Ile for w2-w3-w4), and in the case of Arg122 Pro127 Lys160, also at positions 159 (Ala-not Ala for ayw1-ayw2 and ayw4) and 140 (not Ser-Ser for ayw2-ayw4) [12].

Mutation analysis

Generally, the mutation definition and analysis were performed as previously described [13]. The whole genome, including S gene, PreC/C gene, reverse transcriptase (RT), X gene and basal core promoter (BCP) were under nucleotide mutation analysis. Base on the structure of HBV/CD recombinants, genotype C and genotype D reference sequences were used as parental sequences for mutation detection [14, 15]. Nucleotide sequences of datasets in this study and reference HBV sequences were analyzed using the MEGA 7.0 software and Mutation Reporter Tool [16].

The nucleotide PreS/HBsAg sequences obtained were translated into amino acid sequences, aligned and compared with reference sequences. Amino acid variability was defined as the frequencies of residue substitutions at each position [17]. For analysis, HBsAg was divided into subregions corresponding to structural and/or functional domains: the N- terminal region (aa 1 to 99), the major hydrophilic region (MHR, aa 100 to 169) and the C- terminal region (aa 170 to 226). The “a” determinant in MHR (aa 124 to 147), the first loop of “a” determinant (aa 124 to 137) and second loop of “a” determinant (aa 139 to 147) were also analyzed respectively.

Statistical analysis

Statistical analyses were performed using SPSS 22.0 (IBM Corp., Armonk, NY, USA). Chi-squared test and two-tailed Student’s t-test were used for analysis, as appropriate. A P-Value of ≤0.05 was considered as statistically significant.

Ethical approval

The study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki and was approved by the Ethics Committee of the Chinese Center for Disease Control and Prevention. The purpose of the study and the right to information were explained to the participants by research staff. Written informed consent was obtained from each participant before the interview and venous blood collection.

Results

Patient characteristics

During 2014–2017, 1263 HBsAg positive serum samples were identified from the community population of Qinghai-Tibet plateau. There were no differences in age (29.29 ± 17.93 vs 30.56 ± 18.98, P = 0.211) and gender (449/404 vs 201/210, P = 0.206) distribution between study population in Tibet and Qinghai. Of all the 1263 positive serum samples, twenty-seven serum samples (2.14%) were HBsAg/HBsAb coexistence. There were significant differences between the distribution of seromarkers in Tibet and Qinghai. Details of serological information are shown in Supplement Table 1.

One hundred and eighty-three samples were selected from the 1263 HBsAg-positive serum samples and then full-length sequences were amplified for analysis of the HBV genome. A total of 87 samples were HBeAg-negative, and 96 samples were HBeAg-positive. Most isolates had a genome size of 3215 nt, except eight sequences with evidence of deletion. Basic information of HBV/CD complete genome sequences was listed in Table 2.

Table 2 Basic information of HBV CD1 and CD2 isolates in Qinghai-Tibet Plateau
Full size table

Recombination and serological results

There were two types of HBV/CD recombinants detected in this study. HBV/CD1 was with around nt10–800 from genotype D integrated into genotype C to form the recombinant viral strain; HBV/CD2 was with a region around nt10–1500 from genotype D integrated into genotype C to form the recombinant viral strain. Graphical display of JPHMM and Simplot recombination analysis results were shown in Fig. 1. Of all the 183 HBV complete sequences, there were 130(71.0%) HBV/CD1 isolates, 49 (26.8%) HBV/CD2 isolates and four (2.2%) HBV/C2 isolates. As the results of evolutionary divergence over sequences showed, the genotype D fragment of HBV/CD is closest to HBV subgenotype D4, while the genotype C fragment of HBV/CD is closest to the subgenotype C2 (Supplement Table 2 and 3).

Fig. 1
figure1

Genome recombination analyses in two types of HBV/CD isolates (XZ4-CD1, XZ109-CD2) using two approaches. a Schematic diagram of jpHMM results showing the circular HBV genome with alternate genotype C and genotype D regions. b Graphical representations of SimPlot and Bootscaning analysis demonstrate the recombination of genotypes C and D

Full size image

The main serological subtype of HBV/CD is ayw2 (172/179, 96.1%). Of the 130 HBV/CD1 isolates, 126 (96.9%) were predicted to belong to serotype ayw2, two (1.5%) to ayr, one (0.8%) to ayw1 and one (0.8%) to ayw4. Of the 49 HBV/CD2 isolates, 46 (93.9%) belonged to ayw2 and three (6.1%) to adw2. All four HBV/C2 isolates belonged to adr serotype.

Geographical distribution of HBV subgenotypes in the plateau

HBV/CD1 recombinant sequences were predominantly identified in all the geographic regions analyzed in this study. However, most of the HBV/CD2 isolates (93.9%, 46/49) were identified in two regions (Shannan & Rikaze), which was significantly more than the other six regions in the plateau (P = 0). Four genotype C isolates were all identified in the eastern part of the plateau. The distribution of HBV subgenotypes in the plateau is shown in Fig. 2 and Supplement Table 4.

Fig. 2
figure2

Distribution of HBV genotypes in eight study regions of Qinghai-Tibet Plateau. The proportion of the HBV recombinant CD1, CD2 and genotype C in each region is shown. The relative distance to the northeast of Indian is shown by the color of region in the map

Full size image

Nucleotide mutations

Four genotype C isolates in this study were not enough for statistical analysis, so only the HBV/CD recombinants (130 CD1 and 49 CD2) were under mutation analysis. In the level of nucleotide, the double mutations A1762T/G1764A in the BCP region of HBV (nt1742–1849) were frequently observed in HBV/CD sequences (50/179, 27.9%). In addition, T1753C, another mutation in BCP, was found at a frequency of 4.6% in CD1 and 18.4% in CD2 (P = 0.003). Compare to genotype C and genotype D reference sequences, there were different mutation patterns respectively detected in CD1 and CD2 isolates. Nucleotide mutations, including A1984G, T150AC, A2735C, C2523G, A37G, C2712AT, and C1653T were dramatically different in CD1 and CD2, though these mutations were in the same recombination regions (genotype C or genotype D fragment). Many nucleotide mutations cause amino acid change which could potentially change the clinical prognosis of the recombinants. Other mutations, including stop codon mutation G1896A, were identically detected in CD1 and CD2 isolates. Details of main nucleotide mutations and corresponding amino acid changes identified in the complete genome sequences of HBV/CD1 and HBV/CD2 were shown in Table 3 and Fig. 3.

Table 3 Main nucleotide and amino acid mutations of HBV CD1 and CD2 isolates in Qinghai-Tibet Plateau
Full size table
Fig. 3
figure3

Distribution of wild type and nucleotide mutations (amino acid substitutions) in HBV/CD1 and HBV/CD2 genome. Each bar represents the percentage of isolates with mutated nucleotide (amino acid residues) in CD1 and CD2 recombinants

Full size image

Compare to reference sequences of genotype D and genotype C, several nucleotide (amino acid) positions changed in nearly all the HBV/CD1 and HBV/CD2 sequences, such as A942T(aaL613QH for HBV/CD1 and aaH613K for HBV/CD2), T1485A and T3210A(aaS272TN) in P gene, T1485C(aaS38P) in X gene.

Amino acid substitution in PreS/S region

One hundred and seventy-nine HBV CD recombinants with complete genome sequences were under analyses of amino acid substitution in PreS/S region. Amino acid substitution of 27 HBsAg+/HBsAb+ strains (Group I) were compared with 152 HBsAg+/HbsAb- strains (Group II). The distribution of different recombination type (HBV/CD1 and HBV/CD2, P = 0.677) and HBeAg status (P = 0.213) had no significant difference between Group I and Group II. DNA level of the 179 serum samples were all above 5 log10 and had no significant difference between Group I and Group II (P = 0.355). Significant aa substitution diversity was observed within S gene of HBV between Group I and Group II (1.03 vs. 0.39, for substitution per 100 aa, the same below, P < 0.001). Moreover, the aa variabilities in MHR (P < 0.001), “a” determinant (P < 0.001), the first loop (P < 0.001) and the second loop (P < 0.001) of “a” determinant were all more variable in Group I than in Group II. The frequency of PreS deletion was 2.23% (4/179) and also with significant differences between the two groups. Details were listed in Table 4 and Fig. 4.

Table 4 Incidence of amino acid mutation/ deletion between HBsAg+/HBsAb+ group (Group I) and HBsAg+/HBsAb- group (Group II)
Full size table
Fig. 4
figure4

Frequencies of residue substitutions within the S protein. Isolates from HBs Ag/anti-HBs patients (Group I, black bars, n = 27) and solely HBsAg-positive patients (Group II, gray bars, n = 152) were analyzed in intervals of 10 amino acids each. Each bar represents the percentage of patients with mutated amino acid residues in each group at each interval of 10 amino acids per group. Positions where the proportion of sequences harboring mutations was significant between two groups are marked with an asterisk (*, P < 0.05)

Full size image

Discussion

HBV genotypes are related to the severity of liver disease and response to clinical therapy [18]. Compare to other genotypes, HBV genotype C and genotype D carry a higher lifetime risk of liver cirrhosis and hepatocellular carcinoma development [19]. It is believed that recombination can exert an influence on clinical important properties more dramatically than the steady accumulation of natural mutations, which suggests the potential pathogen significance of the HBV/CD recombinants [20]. As far as we know, there was no detailed molecular epidemiology or genetic variability study carried out based on a large number of HBV/CD recombinant complete genome sequences.

In this study, HBV/CD recombinant was the main genotype (179/183, 97.81%) isolated in the plateau. This result is different from recent but smaller sample size reports in Tibet that genotype C and genotype D were the most dominant HBV genotypes [21, 22]. The HBV genotypes showed a distinct geographical distribution all over the world [19]. In both HBV CD1 and CD2 recombinant genomes, the ‘C fragment’ and ‘D fragment’ were genetically close to subgenotype C2 and D4, respectively. HBV/D4 isolates were mainly reported in South America and the Pacific islands [14]. The only report of HBV/D4 in Asia was in northeast India [23]. Specifically, subtype C2 was the main genotype in central and north China [24], so Qinghai-Tibet Plateau was found at the presumed geographical junction to the distribution of these two HBV subgenotypes (Fig. 2). In this study, recombinants CD1 and CD2 had significant different geographic distribution in the plateau. It could also partly explain the different distribution of seromarkers in Tibet and Qinghai. HBV/CD2 isolates were mainly identified in two regions (Shannan & Rikaze) which are close to the boundary of northeast India (Fig. 2). Genotype D fragment was in around nt10–800 of HBV/CD1 and around nt10–1500 of HBV/CD2 (Fig. 1). Then both CD1 and CD2 have the same initial recombination site at nt10. Based on the facts mentioned above, it can be inferred that these two types of recombinations might have a mixed China (subgenotype C2) and Indian (subgenotype D4) origin.

Though lots of HBV genotype D or genotype C isolates were reported in India [23, 25,26,27], no HBV/CD recombinant strain was reported in south Asia until now. Actually, nearly all the HBV/CD recombinants in this study and in previous studies scattered in indigenous of high altitude areas, such as Qinghai-Tibet Plateau, Yunnan-Kweichow Plateau [28], Loess Plateau [29] and Mongolia Plateau (Accession Number: AB270534, AB270534). A possible reason for this phenomenon might be due to the genetic background or migration history of highlanders.

As the diagnosis information in these rural areas is still not detailed enough to analyze the clinical feature of these HBV recombinants, the complete genome features in this study offer us a way to estimate the clinical prognosis of these HBV strains in the indigenous population. Serotype awy2 was the predominant serotype of HBV/CD in this study, as same as the results of genotype D [30]. It indicates that the main serological character had not dramatically changed after the recombination. Previous studies reported that A1762T/G1764A double mutation in BCP (nt.1742–1849) was the strongest viral factor associated with the development of liver disease [20, 31]. Other mutations, such as T53C, G1613A, C1653T, T1753C, A2189C, T3098C and PreS deletions were also reported associated with clinical progress [32,33,34,35,36,37,38]. In this study, the A1762T/G1764A double mutations were observed in 27.93% in HBV CD recombinant sequences. These frequencies were less than previous reports of genotype C in CHB carriers [39], which indicates a lower risk of HCC in the population. The same result was also found in frequencies of mutation A2189C (12.85%), G1613A (8.94%), T1753C (8.38%), C1653T (5.59%), T53C (4.47%), T3098C (1.68%) and PreS deletion (2.23%, 4/179), which were also lower than mutation frequencies of genotype C in CHB patients in Asia [32, 34,35,36,37,38,39,40]. This indicates the clinical progress of CD recombinant seems genetically more temperate than genotype C which had caused the most liver disease and related death in China.

There were several mutations frequently identified in PreC/C and PreS/S regions, which potentially influence the clinical prognosis. However, many of these prevalent nucleotide mutations or amino acid substitutions were unique in the CD recombinants, which were not previously reported in genotype C or other genotypes of HBV genome. The function and influence of these mutations need further analysis.

The prevalence of anti-HBs coexistence was 2.14% in the population of HBsAg-positive patients. The coexistence frequency is lower than previous reports [5,6,7, 41], suggesting possible geographical variability or special recombinant characteristics in these statistical differences. HBeAg status and DNA level of the serum were the focus of the argument about comparability of the HBsAg/HBsAb coexistence studies [8]. In this study, no differences were identified in the distribution of DNA level, HBeAg positivity or recombination types between Group I and Group II. And the DNA level of serum samples in Group I was all above 5 log10 copies/mL. The aa variability of S region in Group I was significantly higher than that in Group II (Table 3, P < 0.001), which was consistent with previous reports [41]. In further analysis of each segment in the S region, the aa variability in Group I was significantly increased in N-terminal, MHR and C-terminal, compared with Group II (P < 0.001, P < 0.001, P = 0.013, Fig. 3).

The “a” determinant in MHR is the hot zone of HBsAg/HBsAb coexistence studies, which is located in the hydrophilic region between aa124 and aa147 and act as the most important antigen- binding site in the S region of each HBV serotypes [17, 41]. In this study, compared with Group II, the aa variability in the “a” determinant of Group I increased significantly (1.54% vs 0.06%, P < 0.001). This region is composed of two stem- loop structures, which act as two binding sites of monoclonal antibodies and are of great significance for the effectiveness of hepatitis B vaccine and the clinical detection of antigens. In this study, the aa variation of the first and second loop of the “a” determinant were both statistically correlated with HBsAg/HBsAb coexistence, which consistent with previous studies [5, 17]. There were different opinions about the effect of the second loop [42] or even the entire S protein variablity [7] on the HBsAg/HBsAb coexistence. However, former studies were about genotype B and genotype C [42], while HBV/CD recombinant was the main genotype in this study and the “a” determinant was located in genotype D fragment (Fig. 1). It may be due to the amino acid variation background in different genotypes. Moreover, due to the low incidence of HBsAg/ HBsAb coexistence in all HBV genotypes, the sample sizes in previous studies were less than 20 [5, 7], and the insufficient number of samples may affect the stability of statistical results.

The association between aa variation of S protein and HBsAg/HBsAb coexistence has been reported in many studies [5, 6, 17, 41, 42]. However, the aa mutations in any part of S protein is not a sufficient and necessary condition for the explanation of HBsAg/HBsAb coexistence enigma [8], so there should be at least one auxiliary or secondary condition in the emergence of the coexistence. Previous studies suggested that mutations or deletions in the PreS region were the reason for the coexistence of HBsAg/ HBsAb [43]. The significant difference was also found in the distribution of deletions in PreS region between Group I and Group II (Table 3), suggesting the change of PreS region is also associated with coexistence. Interestingly, the frequency of PreS deletion in HBV/CD recombinants was 2.23%, which is lower than the frequency in other genotypes (4.9%) [7]. Combine with the fact that the frequency of HBsAg/HBsAb coexistence is 2.14% in this study, which is also lower than other genotypes [5,6,7, 41]. PreS deletion may act as a genetic variation background which occasionally affects the interaction between antigens and antibodies, together with the aa mutation in MHR to cause the coexistence.

This study has several limitations. First of all, multiple infections, minor populations of immune escape variants in viral quasispecies may not be identified in this study by PCR product sequencing. Secondly, the combined action of PreS deletion and aa variation in MHR need further support. Finally, the result would be more convincing if we could have a larger sample size of HBsAg/HBsAb coexistence, though 27 HBsAg/HBsAb coexistence in this study is more than most of the previous studies.

Conclusions

In summary, this study describes the geographical distribution, genetical variability and HBsAg/HBsAb coexistence phenomena of HBV isolates in Qinghai-Tibet plateau. HBV/CD recombinant has become the predominant genotype in Qinghai-Tibet Plateau. There were signs that HBV/CD had a mixed China and South Asia origin. PreS deletion and aa variation in S protein may cause the HBsAg/HBsAb coexistence together. Several unique nucleotide mutations were frequently detected in HBV/CD isolates, which could potentially influence the clinical prognosis.

Availability of data and materials

The datasets generated in the current study are based on complete genome sequences which are uploaded by the author in the National Center for Biotechnology Information GenBank database with accession numbers MN683570- MN683729, MN657315- MN657318, KX660674- KX660690.

Abbreviations

HBV:

Hepatitis B Virus

HBsAg:

Hepatitis B surface antigen

anti-HBs:

Antibody to HBsAg

CHB:

Chronic HBV

nt:

Nucleotides

HBcAb:

Anti-hepatitis B core antibody

HBeAg:

Hepatitis B e antigen

HBeAb:

Anti-hepatitis B e antibody

PCR:

Polymerase chain reaction

BCP:

Basal core promoter

MHR:

Major hydrophilic region

aa:

Amino acid

LC:

Liver cirrhosis

HCC:

Hepatocellular carcinoma

References

  1. 1.

    World Health Organization. Global hepatitis report 2017. Geneva: Licence: CC BY-NC-SA 3.0 IGO; 2017. p. 2017.

    Google Scholar 

  2. 2.

    Liang X, Bi S, Yang W, Wang L, Cui G, Cui F, Zhang Y, Liu J, Gong X, Chen Y, et al. Epidemiological serosurvey of hepatitis B in China—declining HBV prevalence due to hepatitis B vaccination. Vaccine. 2009;27:6550–7.

    PubMed  Article  Google Scholar 

  3. 3.

    de Martel C, Maucort-Boulch D, Plummer M, Franceschi S. World-wide relative contribution of hepatitis B and C viruses in hepatocellular carcinoma. Hepatology. 2015;62:1190–200.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  4. 4.

    Ye Z. Climate characteristics of ozone over Tibetan plateau. J Geophys Res. 2003;108:4654.

    Article  Google Scholar 

  5. 5.

    Colson P, Borentain P, Motte A, Henry M, Moal V, Botta-Fridlund D, Tamalet C, Gerolami R. Clinical and virological significance of the co-existence of HBsAg and anti-HBs antibodies in hepatitis B chronic carriers. Virology. 2007;367:30–40.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Liu W, Hu T, Wang X, Chen Y, Huang M, Yuan C, Guan M. Coexistence of hepatitis B surface antigen and anti-HBs in Chinese chronic hepatitis B virus patients relating to genotype C and mutations in the S and P gene reverse transcriptase region. Arch Virol. 2012;157:627–34.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Zhang JM, Xu Y, Wang XY, Yin YK, Wu XH, Weng XH, Lu M. Coexistence of hepatitis B surface antigen (HBsAg) and heterologous subtype-specific antibodies to HBsAg among patients with chronic hepatitis B virus infection. Clin Infect Dis. 2007;44:1161–9.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Gerlich WH. The enigma of concurrent hepatitis B surface antigen (HBsAg) and antibodies to HBsAg. Clin Infect Dis. 2007;44:1170–2.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Sumathi R, Xiangjun Z, Hong T, Campo DS, Guoliang X, Ganova-Raeva LM, Jan D, Khudyakov YE. Evaluation of intra-host variants of the entire hepatitis B virus genome. PLoS One. 2011;6:e25232.

  10. 10.

    Liu T, Wang F, Zhang S, Wang F, Meng Q, Zhang G, Cui F, Dunzhu D, Yin W, Bi S, Shen L. Whole-gene analysis of two groups of hepatitis B virus C/D inter-genotype recombinant strains isolated in Tibet, China. PLoS One. 2017;12:e0179846.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  11. 11.

    Anne-Kathrin S, Ingo B, Mariama AC, Emmanuel G, Burkhard M, Fabien Z, Paul D. Mario S: jpHMM: recombination analysis in viruses with circular genomes such as the hepatitis B virus. Nucleic Acids Res. 2012;40:193–8.

    Article  CAS  Google Scholar 

  12. 12.

    Purdy MA, Talekar G, Swenson P, Araujo A, Fields H. A new algorithm for deduction of hepatitis B surface antigen subtype determinants from the amino acid sequence. Intervirology. 2007;50:45–51.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Liu BM, Li T, Xu J, Li XG, Dong JP, Yan P, Yang JX, Yan L, Gao ZY, Li WP, et al. Characterization of potential antiviral resistance mutations in hepatitis B virus reverse transcriptase sequences in treatment-naive Chinese patients. Antivir Res. 2010;85:512–9.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Spitz N, Mello FCA, Moreira AS, Gusatti CS, Martins RMB, Gomes SA, Bello G, Araujo NM. Reconstruction of the spatial and temporal dynamics of hepatitis B virus genotype D in the Americas. PLoS One. 2019;14:e0220342.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Zhu HL, Wang CT, Xia JB, Li X, Zhang ZH. Establishment of reference sequences of hepatitis B virus genotype C subgenotypes. Genet Mol Res. 2015;14:16521–34.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Bell TG, Kramvis A. Mutation reporter tool : an online tool to interrogate loci of interest, with its utility demonstrated using hepatitis B virus. Virol J. 2013;10:1–8.

    Article  Google Scholar 

  17. 17.

    Lada O, Benhamou Y, Poynard T, Thibault V. Coexistence of hepatitis B surface antigen (HBs Ag) and anti-HBs antibodies in chronic hepatitis B virus carriers: influence of “a” determinant variants. J Virol. 2006;80:2968–75.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Fujiwara K, Matsunami K, Iio E, Nojiri S, Joh T. Novel non-canonical genetic rearrangements termed “complex structural variations” in HBV genome. Virus Res. 2017;238:84–93.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Lin CL, Kao JH. Natural history of acute and chronic hepatitis B: the role of HBV genotypes and mutants. Best Pract Res Clin Gastroenterol. 2017;31:249–55.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Ghosh S, Banerjee P, Deny P, Mondal RK, Nandi M, Roychoudhury A, Das K, Banerjee S, Santra A, Zoulim F, et al. New HBV subgenotype D9, a novel D/C recombinant, identified in patients with chronic HBeAg-negative infection in eastern India. J Viral Hepat. 2013;20:209–18.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Zhang Q, Liao Y, Chen J, Cai B, Su Z, Ying B, Lu X, Tao C, Wang L. Epidemiology study of HBV genotypes and antiviral drug resistance in multi-ethnic regions from Western China. Sci Rep. 2015;5:17413.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Du L, Liu M, Liu M, Zhou T, Cheng X, Liu C, Tang H. The prevalence and replication capacity of a Tibetan dominant HBV strain, C/D Recombinant. Biomed Res Int. 2017;2017:8415907.

    PubMed  PubMed Central  Google Scholar 

  23. 23.

    Banerjee P, Mondal RK, Nandi M, Ghosh S, Khatun M, Chakraborty N, Bhattacharya S, RoyChoudhury A, Banerjee S, Santra A, et al. A rare HBV subgenotype D4 with unique genomic signatures identified in North-Eastern India--an emerging clinical challenge? PLoS One. 2014;9:e109425.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  24. 24.

    Wang Z, Hou J, Zeng G, Wen S, Tanaka Y, Cheng J, Kurbanov F, Wang L, Jiang J, Naoumov NV, et al. Distribution and characteristics of hepatitis B virus genotype C subgenotypes in China. J Viral Hepat. 2007;14:426–34.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Ismail AM, Samuel P, Eapen CE, Kannangai R, Abraham P. Antiviral resistance mutations and genotype-associated amino acid substitutions in treatment-Naïve hepatitis B virus-infected individuals from the Indian subcontinent. Intervirology. 2012;55:36–44.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Narayanan K, Thayumanavan L, Mariakuttikan J. Predominance of HBV genotype D in southern part of India. BMC Infect Dis. 2014;14:P22.

    PubMed Central  Article  Google Scholar 

  27. 27.

    Saikia A, Bose M, Barman NN, Deka M, Thangkhiew RS, Bose S. Molecular epidemiology of HBV infection in chronic hepatitis B virus infected patients in Northeast India. J Med Virol. 2015;87:1539–48.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Wang B, Feng Y, Li Z, Duan H, Zhao T, Zhang A, Liu L, Baloch Z, Xia X. Distribution and diversity of hepatitis B virus genotypes in Yunnan, China. J Med Virol. 2014;86:1675–82.

    PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Zhou B, Xiao L, Wang Z, Chang ET, Chen J, Hou J. Geographical and ethnic distribution of the HBV C/D recombinant on the Qinghai-Tibet plateau. PLoS One. 2011;6:e18708.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Bahramali G, Sadeghizadeh M, Amini-Bavil-Olyaee S, Alavian SM, Behzad-Behbahani A, Adeli A, Aghasadeghi MR, Amini S, Mahboudi F. Clinical, virologic and phylogenetic features of hepatitis B infection in Iranian patients. World J Gastroenterol. 2008;14:5448–53.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Fang ZL, Sabin CA, Dong BQ, Ge LY, Wei SC, Chen QY, Fang KX, Yang JY, Wang XY, Harrison TJ. HBV A1762T, G1764A mutations are a valuable biomarker for identifying a subset of male HBsAg carriers at extremely high risk of hepatocellular carcinoma: a prospective study. Am J Gastroenterol. 2008;103:2254–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Jiang SZ, Gao ZY, Tong LI, Ya-Juan LI, Chen XM, Wang L, Feng-Min LU, Zhuang H. T3098C and T53C mutations of HBV genotype C is associated with HBV infection Progress1. Biomed Environ Sci. 2009;22:511–7.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Li MS, Lau TC, Chan SK, Wong CH, Ng PK, Sung JJ, Chan HL, Tsui SK. The G1613A mutation in the HBV genome affects HBeAg expression and viral replication through altered core promoter activity. PLoS One. 2011;6:e21856.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Kim JK, Chang HY, Lee JM, Baatarkhuu O, Yoon YJ, Park JY, Kim DY, Han KH, Chon CY, Ahn SH. Specific mutations in the enhancer II/core promoter/precore regions of hepatitis B virus subgenotype C2 in Korean patients with hepatocellular carcinoma. J Med Virol. 2009;81:1002–8.

    PubMed  Article  CAS  Google Scholar 

  35. 35.

    Shinkai N, Tanaka Y, Ito K, Mukaide M, Hasegawa I, Asahina Y, Izumi N, Yatsuhashi H, Orito E, Joh T, Mizokami M. Influence of hepatitis B virus X and core promoter mutations on hepatocellular carcinoma among patients infected with subgenotype C2. J Clin Microbiol. 2007;45:3191–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Zhu Y, Jin Y, Guo X, Bai X, Chen T, Wang J, Qian G, Groopman JD, Gu J, Li J, Tu H. Comparison study on the complete sequence of hepatitis B virus identifies new mutations in core gene associated with hepatocellular carcinoma. Cancer Epidemiol Biomark Prev. 2010;19:2623–30.

    CAS  Article  Google Scholar 

  37. 37.

    Jiao F, Long L, Ding S, Xie X, Jia L, Lu F. Complete genome sequencing and clinical analysis of intrahepatic hepatitis B virus cccDNA from HCC. Microb Pathog. 2017;109:49–55.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Li F, Li X, Yan T, Liu Y, Cheng Y, Xu Z, Shao Q, Liao H, Huang P, Li J, et al. The preS deletion of hepatitis B virus (HBV) is associated with liver fibrosis progression in patients with chronic HBV infection. Hepatol Int. 2018;12:107–17.

    PubMed  Article  Google Scholar 

  39. 39.

    Li X, Liu Y, Xu Z, Wan Z, Bai S, Mao P, Mao Y, Xin S, Xu D. A complete genomic analysis of hepatitis B virus isolated from 516 Chinese patients with different clinical manifestations. J Med Virol. 2013;85:1698–704.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Zhang K, Imazeki F, Fukai K, Arai M, Kanda T, Mikata R, Yokosuka O. Analysis of the complete hepatitis B virus genome in patients with genotype C chronic hepatitis in relation to HBeAg and anti-HBe. J Med Virol. 2007;79:683–93.

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Chen Y, Qian F, Yuan Q, Li X, Wu W, Guo X, Li L. Mutations in hepatitis B virus DNA from patients with coexisting HBsAg and anti-HBs. J Clin Virol. 2011;52:198–203.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Wang L, Liu H, Ning X, Gao F. Sequence analysis of the S gene region in HBV DNA from patients positive for both HBsAg and HBsAb tests. Hepatol Res. 2010;40:1212–8.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Huang X, Qin Y, Zhang P, Tang G, Shi Q, Xu J, Qi F, Shen Q. PreS deletion mutations of hepatitis B virus in chronically infected patients with simultaneous seropositivity for hepatitis-B surface antigen and anti-HBS antibodies. J Med Virol. 2010;82:23–31.

    CAS  PubMed  Article  Google Scholar 

Download references

Acknowledgements

The authors thank Jingyuan Cao, Qiudong Su, Wenjiao Yin and Wenting Zhou for providing us with technical support and data analysis guidance for this research.

Funding

This study was supported by funds from the Chinese National Science and Technology Major Project (Grant No. 2017ZX10105015–001-002; Recipient: Liping Shen). The funders had no role in study design, data collection and analysis and preparation of the manuscript.

Author information

Author notes

  1. He Liu and Liping Shen contribute the same in this article.

Affiliations

  1. NHC Key Laboratory of Medical Virology and Viral Diseases, National Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Changbai Rd 155#, Changping District, Beijing, China

    He Liu, Liping Shen, Shuang Zhang, Feng Wang, Feng Qiu & Shengli Bi

  2. Department of Infectious Disease, Tianjin Center for Disease Control and Prevention, Tianjin, People’s Republic of China

    He Liu

  3. National Immunization Program, Chinese Center for Disease Control and Prevention, Beijing, People’s Republic of China

    Guomin Zhang, Zundong Yin & Fuzhen Wang

  4. Chinese Center for Disease Control and Prevention, Beijing, People’s Republic of China

    Xiaofeng Liang

Authors
  1. He Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Liping Shen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Shuang Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Feng Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Guomin Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Zundong Yin
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Feng Qiu
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xiaofeng Liang
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Fuzhen Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Shengli Bi
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LH and SL designed and carried out the study and wrote the manuscript. ZS, WF, ZG, YZ, LT and QF carried out the serological and molecular biology assays. LX, WF and BS designed the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Shengli Bi.

Ethics declarations

Ethics approval and consent to participate

The study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki and was approved by the Ethics Committee of the Chinese Center for Disease Control and Prevention. The purpose of the study and the right to information were explained to the participants by research staff. Written informed consent was obtained from each participant before the interview and venous blood collection.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Supplementary Table 1.

This table shows the serum markers of 1263 HBV isolates in this study.

Additional file 2: Supplementary Table 2.

This table shows estimates of Evolutionary Divergence (%) over Sequence Pairs between CD recombinants and subgenotypes D1-D11(nt10–800).

Additional file 3: Supplementary Table 3.

This table shows estimates of Evolutionary Divergence (%) over Sequence Pairs between CD recombinants and subgenotypes C1-C7(nt1500–9).

Additional file 4: Supplementary Table 4.

This table shows distribution of HBV subgenotypes in different regions of Qinghai-Tibet Plateau.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liu, H., Shen, L., Zhang, S. et al. Complete genome analysis of hepatitis B virus in Qinghai-Tibet plateau: the geographical distribution, genetic diversity, and co-existence of HBsAg and anti-HBs antibodies. Virol J 17, 75 (2020). https://doi.org/10.1186/s12985-020-01350-w

Download citation

Keywords

  • Hepatitis B virus
  • Mutation/mutation rate
  • Recombination
  • Hepatitis B surface antigen
  • Antibody to HBsAg

Virology Journal

ISSN: 1743-422X