Skip to main content

The dynamic changes of HBV quasispecies diversity in infancy after immunoprophylaxis failure: a prospective cohort study

Abstract

Background

Previous works have observed that younger infants with chronic hepatitis B virus (HBV) infection are more responsive to antiviral treatment. However, the underlying mechanism remains unclear. In this study, the dynamic changes of HBV quasispecies in infants with immunoprophylaxis failure were investigated to provide virological explanations for clinical management on infantile antiviral therapy.

Methods

Thirteen 7-month-old infants with immunoprophylaxis failure and their mothers were enrolled from a prospective cohort, and 8 of them were followed up to 3 years old. The sequences of HBV quasispecies were analyzed by the full-length genome clone-based sequencing, and compared among mothers and their infants at different ages.

Results

The results revealed that the complexity, mutation frequency and genetic distance of HBV quasispecies decreased significantly at full-length, partial open reading frames and regulatory regions of HBV genome at nucleotide level in 7-month-old infants comparing with their mothers, whereas increased significantly to near the maternal level when infants grew up to 3 years old. Furthermore, similar changes were also found in Core, PreS2, RT and P regions of HBV genome at amino acid level, especially for potential NAs-resistant mutants in RT region and immune-escape mutants in Core and PreS2 regions.

Conclusions

This study uncovered the evolution of HBV quasispecies in infancy after mother-to-child transmission, which may provide the virological evidence for explaning that younger children are more responsive to antiviral therapy.

Introduction

Mother-to-child transmission (MTCT) remains one of the predominant routes of hepatitis B virus (HBV) dissemination worldwide. Infants with chronic HBV infection may suffer impaired quality of life and have a higher lifetime risk of developing end-stage liver disease than adults [1]. It has been shown that the risk of hepatocellular carcinoma (HCC) is higher in Asian population with prenatally acquired HBV as compared with horizontally acquired HBV [2]. Therefore, chronicity comes as a serious threat to HBV-infected infants, and appropriate clinical managements for them are in urgent need.

Fortunately, most of the studies conducted on antiviral therapy in pediatric chronic hepatitis B (CHB) patients have shown promising results with significant improvements in the rates of viral control, HBeAg seroconversion and HBsAg loss [3,4,5,6,7,8,9]. Meanwhile, the antiviral efficacy of HBV-infected infants at immune tolerant (IT) stage also showed encouraging outcomes with 61–78% of infants achieved HBV DNA loss, 22–39% achieved HBeAg seroconversion, and 17–22% achieved HBsAg loss [10,11,12,13]. Moreover, a recent study showed that the infants receiving antiviral therapies before 1 year old obtained more benefits than those receiving antiviral therapies after 1 year old, with 83% of HBsAg loss, shorter treatment duration and lower incidence of adverse events [14]. Overall, most of the studies showed that the younger patients were more responsive to antiviral treatment, suggesting the timing of the initial treatment is crucial. However, the underlying mechanisms are unclear.

Interestingly, several studies revealed that antiviral treatment responders generally had a lower quasispecies complexity than the non-responders [15,16,17,18]. Besides that, our previous study [19] revealed a significant lower HBV quasispecies diversity in mothers of infants with immunoprophylaxis failure. However, the HBV quasispecies characteristics of these infants prenatally infected are unclear, and its correlation with antiviral efficacy are unknown. Especially, the dynamic changes of HBV quasispecies during MTCT and evolution in the early stage of infection have not been well defined.

In this study, both 7-month-old infants with immunoprophylaxis failure and their mothers were included, and some infants were followed up to 3 years old. The dynamic changes of HBV quasispecies characteristics were analyzed by the full-length HBV genome clone-based sequencing assay among mothers and their paired infants. This study aims to explore the evolution of HBV quasispecies in infants with immunoprophylaxis failure, which may provide a virological evidence for establishing the optimal approach to the clinical management of infants with chronic HBV infection.

Material and methods

Subjects

As previously reported in our prospective mother-infant paired study [20], total 1177 infants were returned for post-vaccination serologic testing at 7 months, and 20 of 1177 infants were immunoprophylaxis failure after the full course of vaccination. Among 20 infants with immunoprophylaxis failure, 15 infants were infected with genotype C HBV. The other 3 infants infected with genotype B HBV and 2 infants who were not able to identify genotype due to insufficient sera were excluded to avoid the potential impact of genotype on sequence mutation analysis. Therefore, 15 mother-infant pairs infected with genotype C HBV were enrolled in this study, whereas the full-length HBV genome clones were successfully obtained from 13 infants. The blood samples of mothers were collected before delivery. Eight of the 13 infants were followed up to 3 years old. All mothers were antiviral-naïve, and infants were received 0-1-6 vaccination program, combined with one dosage of HBIG within 12 h of birth.

Amplification, cloning and sequencing

HBV genomes were extracted from 200 µL serum samples using QIAamp DNA blood mini kit (Qiagen, Hilden, Germany). The HBV genome was amplified by PCR as described by previous study [19], and followed by clone-based sequencing. Details and primers are provided in Additional file 1. PCR products of about 3200 bp were purified and cloned into the pGEM®-T Easy Vector Systems (Promega, Beijing, China) after the addition of adenylate tails.

Sequence analyses

Sequence segments were assembled to full-length HBV genome and divided into 11 coding and 7 noncoding regions. Details and nucleotide sites of these regions were provided in Additional file 1. As reported by our previous study [19], viral quasispecies characterization was evaluated by three parameters at both nucleotide and amino acid level: mutation frequency, complexity (Shannon entropy) and diversity (mean genetic distance). Quasispecies complexity was measured using normalized Shannon entropy (Sn) [Sn =  − Σi (pi × lnpi)/lnN] where N is the total number of clones, and pi is the percentage of each clone in the viral quasispecies population. Genetic distance was calculated at the nucleotide level under Tamura 3-parameter method and at the amino acid level under Jones-Taylor-Thornton (JTT) matrix-based method. The number of synonymous substitutions per synonymous site (dS) and the number of nonsynonymous substitutions per nonsynonymous site (dN) were calculated under modified Nei-Gojobori model with Jukes-Cantor correction. The phylogenetic analyses were performed by neighbor-joining tree constructed by Tamura 3-parameter model given in the MEGAX software. Mutations were identified by a same consensus sequence synthesized by all clones from mothers. B cell and CD4+ T cell immune epitopes in Core and PreS2 regions were identified based on previous publications [21,22,23].

Serological assays

Serum HBsAg and HBeAg were tested by Abbott chemiluminescent microparticle immunoassay (Abbott Diagnostic, Chicago, IL, USA), as well as HBV DNA load was measured by Abbott real-time HBV DNA assay (Abbott Molecular, IL, USA).

Statistical analyses

Categorical variables were expressed as % (m/n) and examined by χ2/Fisher’s exact test. Non-normal distributions data were expressed as median and IQR or (range) and compared by Wilcoxon Signed Rank Test. All P values were two-tailed and a P value < 0.05 was considered significant. Statistical analyses were analyzed using SPSS software V.25.0 (Chicago, IL, USA).

Results

Demographic and virological data

Totally 13 7-month-old infants and their mothers were enrolled, and 8 of them were followed up to 3 years old. The levels of serum HBV DNA, HBsAg and HBeAg did not change significantly among mothers, 7-month-old infants and 3-year-old infants (Table 1). In addition, an average of 16 ± 4.72 full-length HBV genome clones per sample were collected, and no significant differences of the clone numbers were noticed among groups.

Table 1 The demographic, clinical and laboratory data of 13 mother-infant pairs

The phylogenetic trees were constructed between 13 mothers and their paired 7-month-old infants, 8 infants at 7 months and 3 years old, as well as 8 mothers and their paired infants at 3 years old (Additional file 1: Figure S1). All clones clustered together with the reference sequence of genotype C HBV, and the sequences of all clones from the same pair were clustered together, indicating there was no contamination during the acquisition of clones.

Comparative analysis of HBV quasispecies characteristics in mothers and their paired 7-month-old infants

The quasispecies characteristic values, including complexity, mutation frequency and genetic distance, at full-length genome level, 11 coding regions and 7 noncoding regions were analyzed between mothers and 7-month-old infants. As shown in Additional file 1: Table S1 and Fig. 1A–C, the complexity, mutation frequencies and genetic distances of full-length genome and most specific regions at nucleotide level in 7-month-old infants were significantly lower than those in mothers.

Fig. 1
figure 1

Comparative analysis of HBV quasispecies characteristics during MTCT. The quasispecies complexity, mutation frequency and genetic distance were analyzed between mothers and their paired 7-month-old infants at nucleotide level (A, B and C) and amino acid level (D, E and F), respectively. The comparison of synonymous substitution rate (dS) (G) and non-synonymous substitution rate (dN) (H). *represents P < 0.05. **represents P < 0.01. “nt” represents the nucleotide level. “aa” represents the amino acid level. NTCP-BD: sodium taurocholate cotransporting polypetide binding domain. MHR: major hydrophilic region. Quasispecies complexity was measured using normalized Shannon entropy (Sn). Genetic distance was calculated at the nucleotide level under Tamura 3-parameter method and at the amino acid level under Jones–Taylor–Thornton matrix-based method. dS and dN were calculated under modified Nei–Gojobori model with Jukes–Cantor correction

For coding regions, the mutation frequencies of Core, PreS2, RT, NTCP-BD, X and PreS1 regions, as well as the complexities and genetic distances of Core, PreS2, P and RT regions at amino acid level in 7-month-old infants were also significantly lower than those in mothers (Additional file 1: Table S2 and Figs. 1D–F). Further, the synonymous substitution rates of most regions and the non-synonymous substitution rates of Core, PreS2, P and RT regions were dropped significantly after MTCT (Additional file 1: Table S3 and Fig. 1G, H).

The phylogenetic tree and mutation rate of single nucleotide site analyses were performed on each pair of infants and their mothers. Phylogenetic trees of HBV sequences in 10 infants (76.92%, 10/13) segregated from those of mothers, indicating the evolve selection was present during MTCT of HBV (Additional file 1: Figure S2). The mutation rates of some nucleotide sites were significantly different in 10 pairs of infants and their mothers (Table 2). In one of the other 3 pairs with sequences mixed in phylogenetic trees, C1826T and A1827C mutations were present in 40% of HBV sequences in infant, which is significantly higher than that of mother (Table 2).

Table 2 Detailed nucleotide site with mutation rate changed significantly during MTCT

Dynamics of HBV quasispecies characteristics from 7 months to 3 years old

The quasispecies characteristic values of full-length decreased significantly after MTCT and increased to near maternal level from 7 months to 3 years old, while the levels of serum HBV DNA, HBsAg and HBeAg did not change obviously (Fig. 2A).

Fig. 2
figure 2

Comparative analysis of serum markers and HBV quasispecies characteristics between 8 pairs of mothers and infants at 7 months and 3 years old. Dynamic change of serum markers and HBV quasispecies characteristics at full-length HBV genome level (A). The quasispecies complexities (B), mutation frequencies (C) and genetic distances (D) of Core, PreS2, P, RT, S and NTCP-BD regions at nucleotide level at three time points. The quasispecies complexities (E), mutation frequencies (F) and genetic distances (G) of Core, PreS2, P and RT regions at amino acid level at three time points. The dynamic change of synonymous substitution rates (dS) (H) and the non-synonymous substitution rates (dN) (I) of Core, PreS2, P, RT, S and NTCP-BD regions at three time points. The phylogenetic trees of 8 mothers and their paired infants at 7 months and 3 years old. (J) Black point represents the clones from mothers, light blue point for 7-month-old infants and dark blue point for 3-year-old infants. *represents P < 0.05. “nt” represents nucleotide level. “aa” represents amino acid level. NTCP-BD: sodium taurocholate cotransporting polypetide binding domain; Quasispecies complexity was measured using normalized Shannon entropy (Sn). Genetic distance was calculated at the nucleotide level under Tamura 3-parameter method and at the amino acid level under Jones–Taylor–Thornton matrix-based method. dS and dN were calculated under modified Nei–Gojobori model with Jukes–Cantor correction

At nucleotide acid level, the quasispecies complexities, mutation frequencies and genetic distances of Core, PreS2 and other regions increased to near the maternal level at 3 years of age (Additional file 1: Table S4 and Fig. 2B–D). For amino acid level, the phenomenon was noticed in Core and PreS2 regions (Figs. 2E–G and Additional file 1: Table S5). The synonymous substitution rates of Core, P and NTCP-BD regions (Fig. 2H and Additional file 1: Table S6) and the non-synonymous substitution rates of Core and PreS2 regions also significantly increased to near the maternal level at 3 years old (Fig. 2I and Additional file 1: Table S6).

The phylogenetic tree of each pair revealed that most sequences of infants at 7 months and 3 years old partially mixed (Fig. 2J). The mutation rates of single nucleotide site were analyzed, and 19 mutants with mutation rates significantly changed have been found in 5 infants (Table 3). Among them, the mutation rates of 3 sites significantly decreased from 7 months to 3 years old, and caused amino acid substitutions; While in 15 nucleotide sites with the mutation rates significantly increased, 13 of them caused amino acid substitution, including 10 located at Core region. Core region is the most diverse region during HBV quasispecies evolution in the early stage of infection. It is worth noticing that these significant changing mutants in Core region were found in 4 of 5 male infants at 3 years old.

Table 3 The detailed nucleotide sites with mutation rate changed significantly from 7 months to 3 years of age

Dynamics of the potential NAs-resistant mutations in RT region and the mutation rates of single amino acid site in core and PreS2 regions

The potential NAs-resistant mutants defined by previous reports [24, 25] were searched from the clones of 8 mothers (118 clones) and their paired infants at 7 months (133 clones) and 3 years old (141 clones). The deletion rate in RT region was significantly higher in mothers (8.47%, 10/118) than that in 7-month-old infants (0, 0/133) (P < 0.001) (Fig. 3A). Totally, 21 potential NAs-resistant mutants were found and listed in Additional file 1: Table S7. The cumulative rate of NAs-resistant mutations (the number of clone with NAs-resistant mutants/total clone number, clones contain multiple NAs-resistant mutants were counted multiple times) was significantly lower in 7-month-old infants (26.32%, 35/133) than that in mothers (49.15%, 58/118) (P = 0.001) and 3-year-old infants (39.01%, 55/141) (P = 0.025), whereas there was no significant difference between mothers and 3-year-old infants (Fig. 3B). Similarly, the ratio of clones with NAs-resistant mutants (clones contain multiple NAs-resistant mutations were counted only once) was also significantly lower in 7-month-old infants (23.31%, 31/133) than that in mothers (43.22%, 51/118) (P = 0.001) and 3-year-old infants (35.46%, 50/141) (P = 0.028), whereas there was no significant difference between mothers and 3-year-old infants (Fig. 3C).

Fig. 3
figure 3

Comparative analysis of the potential NAs-resistant mutations in RT region and the mutation rates of single amino acid site in Core and PreS2 regions. Dynamic changes of the indel mutation rates in Core and RT regions of HBV genome from mothers and infants at 7 months and 3 years old. (A) The cumulative rate of NAs-resistant mutations in RT region of HBV genome from mothers and infants at 7 months and 3 years old. (B) The ratio of clones with NAs-resistant mutations in RT region of HBV genome from mothers and infants at 7 months and 3 years old. (C) The single amino acid site mutation rate in Core and PreS2 region of all clones from mothers and infants at 7 months and 3 years old. (D) Black line represents the data from mothers, light blue line for 7-month-old infants and dark blue line for 3-year-old infants. Sites with mutation rate > 1% in Core region and sites with mutation rate > 0.5% in PreS2 region were noted. All the mutations were defined based on a same consensus sequence synthesized by all clones from mothers

Due to the dramatic change of diversity in Core and PreS2 regions during MTCT and early infancy, the substitution rates of single amino acid site were calculated in these regions. More indels were found in mothers (3.40%, 4/118) than that in 7-month-old infants (0, 0/133) (P = 0.048), and there were also more indels in 3-year-old infants (5.67%, 8/141) than that in 7-month-old infants (P = 0.007), whereas there was no significant difference between mothers and 3-year-old infants (Fig. 3A). As shown in Fig. 3D, more substitutions in B cell and CD4+ T cell epitopes in Core and PreS2 regions were found in the clones from mothers and 3-year-old infants than that in 7-month-old infants. For 184 amino acid sites in Core region, the ratio of sites with substitution rate over 1% was lower in 7-month-old infants (2.72%, 5/184) than that in mothers (8.15%, 15/184) (P = 0.021) and 3-year-old infants (13.59%, 25/184) (P < 0.001) (Fig. 3D). For 55 amino acid sites in PreS2 region, the ratio of sites with substitution rate over 0.5% was lower in 7-month-old infants (18.18%, 10/55) than that in mothers (40%, 22/55) (P = 0.012) and 3-year-old infants (38.18%, 21/55) (P = 0.02) (Fig. 3D).

Discussion

In our previous study, the results of full-length HBV genome clone-based sequencing showed superiority in comprehensively depicting the feature of HBV quasispecies from whole genome landscape [19]. In this study, the infants and mothers both were antiviral treatment-naive, thus these subjects were suitable for exploring HBV quasispecies’ natural dynamic changes and its interaction with host.

MTCT is considered as a bottleneck event for HBV proliferation, therefore the viral population is expected to have competent replication capacity and low mutation rate after MTCT [19]. Therefore, HBV quasispecies characteristic values were significantly decreased after MTCT, especially for the regions related with immune escape (PreS2 and Core regions) and replication capacity (RT, ENI, SPI and XP regions). Further, the decline of many substitutions in B cell and CD4+ T cell epitopes of PreS2 and Core regions after MTCT, might due to the immature immunity of infants. Since the virus strains of immune escape commonly accompanied with replication capacity undermined, the virus with these mutants might be outnumbered by wild type which have stronger replication capacity in a new environment with feeble immunity. Meanwhile, the different human leukocyte antigen (HLA) types might also play an important role during this period as the individualized mutants found in 84.62% (11/13) mother-infant pairs (Additional file 1: Table S4). There were 18 mutations with an increased mutation rate after MTCT, and 12 of them lead to amino acid substitution, including 4 substitutions in S region (P62L, W74L, W165S and S193L), 3 substitutions in PreS2 region (A11T, V17E and A39V), 5 substitutions in PreS1 (V90A), P (Y95F), PreC (H5S), Core (P156T) and X region (V5M), respectively. However, these substitutions are individualized, and none of them was found in two or more infants. Since these substitutions are mainly distributed in PreS1, PreS2 and S region (66.67%, 8/12), it indicates the select pressure on HBV surface antigen. It’s worth noticing that no substitutions were found in “a” determinant region which was considered as the major region occurring vaccine escape mutants [26]. This result was consistent with our previous report that a more complex mutant spectrum in “a” determinant region might be more vulnerable to extinct through MTCT, and the vaccine-escape mutations was not a significant factor of immunoprophylaxis failure [19].

In this study, we found that the characteristics of HBV quasispecies in Core, PreS2, P, NTCP-BD and RT regions, which were strongly correlated with the host-immunity, virus infection and replication capacity, increased to near the maternal level at 3 years old, suggested that the complexity and diversity of HBV quasispecies increased along with age and reached at adult’s level at 3 years old. Several studies have reported a negative correlation between antiviral efficacy and quasispecies complexities in certain regions, such as RT region for HBV and HVR-1 region for hepatitis C virus (HCV) [15,16,17,18, 27, 28]. In this study, both the quasispecies complexity at nucleotide level and the ratios of potential NAs-resistant mutation in RT region were lower in the 7-month-old infants than that in mothers, and increased significantly to near the maternal levels at 3 years old. Therefore, during HBV evolution in infancy, the ability of HBV quasispecies against antiviral treatment increased along with age, especially for NAs treatment, which could explain the phenomenon that more efficacious outcomes were observed in infants younger than 1 year old and NAs treatment was more effective than IFN for infants [14]. Further, more mutations in B cell and CD4+ T cell epitopes in Core and PreS2 regions were found in older infants, suggesting that it might be easier for HBV quasispecies to achieve immune clearance in infants at 7 months old than that in infants at 3 years old in term of virological factors. Combined with the fast division of liver cells in infants, the virus might be diluted, which would ultimately accelerate HBV elimination [29]. Thus, these findings could explain the phenomenon that the younger infants are more responsive to antiviral treatment [13].

Conclusions

As presented in Fig. 4, our results demonstrate a relatively simple and pure viral population with low level of potential NAs-resistant and immune-escape mutants in HBV immunoprophylaxis failure infants at 7 months old, and the viral population will grow diverse to reach the maternal level at 3 years old. This study uncovered the evolution of HBV quasispecies in infancy after mother-to-child transmission, which may provide the virological evidence for explaning that younger children are more responsive to antiviral therapy.

Fig. 4
figure 4

Graphic abstract for the dynamics of hepatitis B virus quasispecies after MTCT and evolution in infancy

Availability of data and materials

Data are available upon reasonable request.

Abbreviations

HBV:

Hepatitis B virus

MTCT:

Mother-to-child transmission

HCC:

Hepatocellular carcinoma

CHB:

Chronic hepatitis B

HBIG:

Hepatitis B immunoglobin

NTCP-BD:

Sodium taurocholate cotransporting polypetide binding domain

MHR:

Major hydrophilic region

NAs:

Nucleotide analogues

HLA:

Human leukocyte antigen

References

  1. Dienstag JL. Hepatitis B virus infection. N Engl J Med. 2008;359(14):1486–500.

    Article  CAS  Google Scholar 

  2. Chen CH, et al. Hepatitis B virus transmission and hepatocarcinogenesis: a 9 year retrospective cohort of 13676 relatives with hepatocellular carcinoma. J Hepatol. 2004;40(4):653–9.

    Article  Google Scholar 

  3. Dienstag JL, et al. Lamivudine as initial treatment for chronic hepatitis B in the United States. N Engl J Med. 1999;341(17):1256–63.

    Article  CAS  Google Scholar 

  4. Jonas MM, et al. Clinical trial of lamivudine in children with chronic hepatitis B. N Engl J Med. 2002;346(22):1706–13.

    Article  CAS  Google Scholar 

  5. Kobak GE, et al. Interferon treatment for chronic hepatitis B: enhanced response in children 5 years old or younger. J Pediatr. 2004;145(3):340–5.

    Article  CAS  Google Scholar 

  6. Wirth S, et al. Efficacy and safety of peginterferon Alfa-2a (40KD) in children with chronic hepatitis B: the PEG-B-ACTIVE study. Hepatology. 2018;68(5):1681–94.

    Article  CAS  Google Scholar 

  7. Choe BH, et al. Long-term therapeutic efficacy of lamivudine compared with interferon-alpha in children with chronic hepatitis B: the younger the better. J Pediatr Gastroenterol Nutr. 2007;44(1):92–8.

    Article  CAS  Google Scholar 

  8. Fan H, et al. Interferon alpha treatment leads to a high rate of hepatitis B surface antigen seroconversion in Chinese children with chronic hepatitis B. J Viral Hepat. 2019;26(Suppl 1):77–84.

    Article  CAS  Google Scholar 

  9. Liu Y, et al. Long-term efficacy and safety of peginterferon in the treatment of children with HBeAg-positive chronic hepatitis B. J Viral Hepat. 2019;26(Suppl 1):69–76.

    Article  CAS  Google Scholar 

  10. D’Antiga L, et al. Combined lamivudine/interferon-alpha treatment in “immunotolerant” children perinatally infected with hepatitis B: a pilot study. J Pediatr. 2006;148(2):228–33.

    Article  CAS  Google Scholar 

  11. Poddar U, et al. Cure for immune-tolerant hepatitis B in children: is it an achievable target with sequential combo therapy with lamivudine and interferon? J Viral Hepat. 2013;20(5):311–6.

    Article  CAS  Google Scholar 

  12. Zhu S, et al. Antiviral therapy in hepatitis B virus-infected children with immune-tolerant characteristics: a pilot open-label randomized study. J Hepatol. 2018;68(6):1123–8.

    Article  CAS  Google Scholar 

  13. Zhu SS, et al. A randomized controlled study on factors influencing the curative effect of sequential combined interferon and lamivudine therapy in children with immune-tolerant phase chronic hepatitis B. Chin J Hepatol. 2019;27(8):604–9.

    CAS  Google Scholar 

  14. Zhu S, et al. Early initiation of antiviral therapy contributes to a rapid and significant loss of serum HBsAg in infantile-onset hepatitis B. J Hepatol. 2019. https://doi.org/10.1016/j.jhep.2019.07.022.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Chambers TJ, et al. Quasispecies heterogeneity within the E1/E2 region as a pretreatment variable during pegylated interferon therapy of chronic hepatitis C virus infection. J Virol. 2005;79(5):3071–83.

    Article  CAS  Google Scholar 

  16. Chen L, et al. Early changes of hepatitis B virus quasispecies during lamivudine treatment and the correlation with antiviral efficacy. J Hepatol. 2009;50(5):895–905.

    Article  CAS  Google Scholar 

  17. Liu F, et al. Evolutionary patterns of hepatitis B virus quasispecies under different selective pressures: correlation with antiviral efficacy. Gut. 2011;60(9):1269–77.

    Article  CAS  Google Scholar 

  18. Shire NJ, et al. HCV kinetics, quasispecies, and clearance in treated HCV-infected and HCV/HIV-1-coinfected patients with hemophilia. Hepatology. 2006;44(5):1146–57.

    Article  CAS  Google Scholar 

  19. Xiao Y, et al. Quasispecies characteristic in “a” determinant region is a potential predictor for the risk of immunoprophylaxis failure of mother-to-child-transmission of sub-genotype C2 hepatitis B virus: a prospective nested case-control study. Gut. 2020;69(5):933–41.

    Article  CAS  Google Scholar 

  20. Lu Y, et al. The maternal viral threshold for antiviral prophylaxis of perinatal hepatitis B virus transmission in settings with limited resources: a large prospective cohort study in China. Vaccine. 2017;35(48 Pt B):6627–33.

    Article  CAS  Google Scholar 

  21. Chen BF, et al. High prevalence and mapping of pre-S deletion in hepatitis B virus carriers with progressive liver diseases. Gastroenterology. 2006;130(4):1153–68.

    Article  CAS  Google Scholar 

  22. Li F, et al. Whole genome characterization of hepatitis B virus quasispecies with massively parallel pyrosequencing. Clin Microbiol Infect. 2015;21(3):280–7.

    Article  Google Scholar 

  23. Yang G, et al. Quasispecies characteristics in mother-to-child transmission of hepatitis B virus by next-generation sequencing. J Infect. 2017;75(1):48–58.

    Article  Google Scholar 

  24. Choi YM, Lee SY, Kim BJ. Naturally occurring hepatitis B virus reverse transcriptase mutations related to potential antiviral drug resistance and liver disease progression. World J Gastroenterol. 2018;24(16):1708–24.

    Article  CAS  Google Scholar 

  25. Zhang Q, et al. Incidence of natural resistance mutations in naive chronic hepatitis B patients: a systematic review and meta-analysis. J Gastroenterol Hepatol. 2015;30(2):252–61.

    Article  Google Scholar 

  26. Carman WF, et al. Fulminant reactivation of hepatitis B due to envelope protein mutant that escaped detection by monoclonal HBsAg ELISA. Lancet. 1995;345(8962):1406–7.

    Article  CAS  Google Scholar 

  27. Koizumi K, et al. Diversity of quasispecies in various disease stages of chronic hepatitis C virus infection and its significance in interferon treatment. Hepatology. 1995;22(1):30–5.

    CAS  PubMed  Google Scholar 

  28. Abbate I, et al. HVR-1 quasispecies modifications occur early and are correlated to initial but not sustained response in HCV-infected patients treated with pegylated- or standard-interferon and ribavirin. J Hepatol. 2004;40(5):831–6.

    Article  CAS  Google Scholar 

  29. Allweiss L, et al. Proliferation of primary human hepatocytes and prevention of hepatitis B virus reinfection efficiently deplete nuclear cccDNA in vivo. Gut. 2018;67(3):542–52.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank the study investigators, coordinators, nurses, patients and their families for their contributions.

Funding

This study was supported by grants from the National Major Scientific and Technological Special Project during the Thirteenth Five-year Plan Period (2017ZX10201201003), and the National Natural Science Foundation of China (81571949).

Author information

Authors and Affiliations

Authors

Contributions

JL obtained funding. JL and JW designed the study. YL, YX, LL, YS, LY, FD, JW and JL performed the study. XZ, JL, LZ, JJ, LL, CL and JL were in charge of patients enrolled in this study. YL, YX, LL, YS, JL, ZD, HZ, JW and JL analyzed data. YL drafted the manuscript. JL and JW contributed to the interpretation of the results and critical revision of the manuscript for important intellectual content and approved the final version of the manuscript. JL and JW are the study guarantors. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Jie Wang or Jie Li.

Ethics declarations

Ethics approval and consent to participate

The experimental protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki. This study was approved by the Ethics Committee of Peking University Health Science Center. Written informed consent was obtained from each subject.

Consent to 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

. Detailed materials and supplementary Tables and Figures.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Xiao, Y., Li, L. et al. The dynamic changes of HBV quasispecies diversity in infancy after immunoprophylaxis failure: a prospective cohort study. Virol J 18, 236 (2021). https://doi.org/10.1186/s12985-021-01707-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-021-01707-9

Keywords