Quasispecies characteristics in mother-to-child transmission of hepatitis B virus by next-generation sequencing

A hepatitis B virus (HBV) sequence variation graph improves alignment and sample-specific consensus sequence construction

Nearly 300 million individuals live with chronic hepatitis B virus (HBV) infection (CHB), for which no curative therapy is available. As viral diversity is associated with pathogenesis and immunological control of infection, improved methods to characterize this diversity could aid drug development efforts. Conventionally, viral sequencing data are mapped/aligned to a reference genome, and only the aligned sequences are retained for analysis. Thus, reference selection is critical, yet selecting the most representative reference a priori remains difficult. We investigate an alternative pangenome approach which can combine multiple reference sequences into a graph which can be used during alignment. Using simulated short-read sequencing data generated from publicly available HBV genomes and real sequencing data from an individual living with CHB, we demonstrate alignment to a phylogenetically representative 'genome graph' can improve alignment, avoid issues of reference ambiguity, and facilitate the construction of sample-specific consensus sequences more genetically similar to the individual's infection. Graph-based methods can, therefore, improve efforts to characterize the genetics of viral pathogens, including HBV, and have broader implications in host-pathogen research.

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

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.

Intra-host diversity of hepatitis B virus during mother-to-child transmission: the X gene may play a key role in virus survival in children after transmission

Mother-to-child transmission of hepatitis B virus (HBV) is the main route of transmission in Asia, and characterization of HBV quasispecies is needed to further understand virus evolution and adaptation. To understand changes in HBV during mother-to-child transmission, we enrolled nine pairs of mothers and children in the study, including a set of twins. Three groups were infected with HBV genotype C, and six groups were infected with HBV genotype B. The full-length HBV genome was amplified by PCR from serum samples before antiviral treatment, the whole viral genomes from each pair were sequenced, and the complexity and diversity of the quasispecies were analyzed. The entropy of transmitted HBV in children was found to be lower than their mothers, suggesting that there was a bottleneck effect during HBV transmission from the mother to the child. Selective evolution was shown by calculating πN and πS in the whole genomes, and the highest values were obtained for the X gene, which plays a role in viral replication and immune escape. All genotype C patients and only one genotype B pair had a πN/πS greater than 1 ratio, indicating that positive selection had occurred. In addition, quasispecies were found to be different between the twin children despite having the same mother, indicating that virus evolution is host-specific.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

Establishment of Coamplification at Lower Denaturation Temperature PCR/Fluorescence Melting Curve Analysis for Quantitative Detection of Hepatitis B Virus DNA, Genotype, and Reverse Transcriptase Mutation and Its Application in Diagnosis of Chronic Hepatitis B

Dynamic and real-time hepatitis B virus (HBV) DNA, genotype, and reverse transcriptase mutation analysis plays an important role in diagnosing and monitoring chronic hepatitis B (CHB) and in assessing the therapeutic response. We established a highly sensitive coamplification at lower denaturation temperature PCR (COLD-PCR) coupled with probe-based fluorescence melting curve analysis (FMCA) for precision diagnosis of CHB patients. The imprecision with %CV and detection limit of HBV DNA detected by COLD-PCR/FMCA were 2.58% to 4.42% and 500 IU/mL, respectively. For mutation, the imprecision and detection limit were 3.35% to 6.49% and 1%, respectively. Compared with Sanger sequencing, the coincidence rates of genotype and mutation were 96.0% and 82.5%, respectively, whereas the inconsistent data resulted from a low proportion (<20%) of mixed genotypes or mixed mutations. The mutation ratio in HBV infection patients was as follows: hepatitis B e antigen (HBeAg)-positive infection (0/0.0%) < HBeAg-negative infection (16/4.5%) < HBeAg-positive hepatitis (30/5.5%) < HBeAg-negative hepatitis (36/6.5%). In patients with entecavir therapy, the proportion of mutation at baseline or week 4 in virologic response (VR) group was <4%, whereas in the partial VR group, it was mostly ≥4%. COLD-PCR/FMCA provides a novel tool with high sensitivity, convenience, and practicability for the simultaneous quantification of HBV DNA, genotype, and mutation. It might be used for distinguishing the different phases of HBV infection and predicting VR of CHB patients.

Characterization and Clinical Significance of Natural Variability in Hepatitis B Virus Reverse Transcriptase in Treatment-Naive Chinese Patients by Sanger Sequencing and Next-Generation Sequencing

Mutations in hepatitis B virus (HBV) reverse transcriptase (RT) are associated with nucleos(t)ide analogue (NA) resistance during long-term antiviral treatment. However, the characterization of mutations in HBV RT in untreated patients has not yet been well illustrated. The objective of this study was to investigate the characterization and clinical significance of natural variability in HBV RT in treatment-naive patients. HBV RT sequences were analyzed in 427 patients by Sanger sequencing and in 66 patients by next-generation sequencing. Primary or secondary NA resistance (NA^r) mutations were not found, except A181T in RT (rtA181T) by Sanger sequencing, but they were detected by next-generation sequencing. Mutations were found in 56 RT amino acid (aa) sites by Sanger sequencing, 36 of which had mutations that could lead to changes in B or T cell epitopes in the RT or S protein. The distribution of mutations was diverse in different sections within the RT region. Multiple mutations showed significant association with HBV DNA, HBsAg, HBeAg, age, and severity of liver fibrosis. Mutations at rt251, rt266, rt274, rt280, rt283, rt284, and rt286 were found most in the advanced liver disease (ALD) group by next-generation sequencing. The present study demonstrates that next-generation sequencing (NGS) was more suitable than Sanger sequencing to monitor NA^r mutations at a low rate in the treatment-naive patients, and that mutations in the RT region might be involved in the progression to ALD.

Horizontal Transmission of Hepatitis B Virus From Mother to Child Due to Immune Escape Despite Immunoprophylaxis

Hepatitis B virus (HBV) vaccination starting at birth is approximately 95% effective in preventing mother-to-child transmission to infants born to HBV-infected mothers. A higher risk of transmission is associated with birth to a highly viremic mother, often due to transplacental exposure, while later horizontal transmission is much less common, particularly following complete vaccination. This study reports a case of infection in an older child despite appropriate immunoprophylaxis starting at birth and an apparent protective immune response post-vaccination. Two immune escape mutations within the antigenic determinant of the surface antigen-coding region were observed in the child's dominant HBV sequence, whereas the maternal HBV variant lacked mutations at both sites. Ultra-deep sequencing confirmed the presence of 1 mutation at low levels within the maternal HBV quasispecies population, suggesting early exposure to the child followed by viral evolution resulting in immunoprophylaxis escape and chronic infection.

Evaluation of performance characteristics of hepatitis B e antigen serologic assays

BACKGROUND

Hepatitis B e antigen (HBeAg) is considered an indicator of high hepatitis B virus (HBV) replication. Performance characteristics of commercially available HBeAg assays have not been determined, thus it is unknown whether lack of HBeAg detection is because of test sensitivity or HBV basal core promoter and precore mutations.

OBJECTIVES

We studied the correlation between HBeAg reactivity with HBV DNA levels in three commercially available HBeAg assays using 335 HBsAg and HBV DNA positive serum/plasma samples.

STUDY DESIGN

Diagnostic sensitivity was determined by serial dilutions of a WHO HBeAg standard. The limit of HBeAg detection estimated through regression was 1 IU/mL (Centaur), 97 IU/mL (DiaSorin) and 129 IU/mL (Vitros). Of these 335 samples, enough sample volume remained in 253 samples for head-to-head comparison of the assays.

RESULTS

81 (32%), 41 (16%) and 36 (14%) of the samples were HBeAg positive by the Centaur, DiaSorin and Vitros assays, respectively. Compared to the FDA-approved Centaur assay the specificity of the other two assays was 98%, while sensitivity was 47% for the DiaSorin assay and 41% for the Vitros assay. Significant association was found between HBeAg positive samples and HBV DNA levels >20,000 IU/mL; 31% of HBeAg negative samples (Centaur) had HBV DNA levels >20,000 IU/mL, 26% of HBeAg positive samples had HBV DNA levels <20,000 IU/mL and 5 HBeAg positive samples had HBV DNA levels <2000 IU/mL.

CONCLUSION

Discordance was seen between these HBeAg assays, indicating reliance on HBeAg alone as a marker of high HBV replication can be misleading. Detection and quantification of HBV DNA remains the accurate and reliable marker of HBV replication.

Association of characteristics of HBV quasispecies with hepatitis B surface antigen seroconversion after pegylated interferon-α-2a treatment in child patients

BACKGROUND

The correlation between hepatitis B surface antigen (HBsAg) seroconversion and the characteristics of HBV quasispecies (QS) before and during pegylated interferon-α-2a (PEG-IFN-α-2a) treatment in hepatitis B e antigen (HBeAg)-positive chronic hepatitis B (CHB) children has not yet been reported.

METHODS

35 patients, including 18 HBsAg seroconverters (SS) and 17 non-seroconverters (SN), were enrolled. Serum samples were collected before treatment and at weeks 12 and 24 of treatment. Sequences within the basal core promoter/pre-core (BCP/PC) and S/reverse transcriptase (S/RT) region were analysed by next-generation sequencing.

RESULTS

There was no significant difference in the baseline diversity of HBV QS (Shannon entropy [Sn]; Hamming distance [HD]) in either region between the two groups. The baseline mutations A1762T/G1764A, C1913A, and T2003A/G or C2004T were correlated with non-response to therapy (P=0.025, P=0.036, P=0.032, respectively). After 24 weeks of therapy, HBV diversity within the BCP/PC region in the SS group notably declined (Sn: P=0.002; HD: P=0.011), while that of the SN group was nearly unchanged. As for the S/RT region, 24 weeks of treatment made no significant difference on QS diversity in either group.

CONCLUSIONS

Our data demonstrated that the baseline viral mutations and dynamic changes in HBV QS diversity within the BCP/PC region were closely related to HBsAg seroconversion in HBeAg-positive CHB children treated with PEG-IFN-α-2a.

Clinical implication and viral mutation in basal core promoter/pre-core of hepatitis B virus C/D recombinant

BACKGROUND AND AIM

Hepatitis B virus (HBV) C/D recombinant is predominant in Tibet in Western China. Although the geographical and ethnic distributions of the C/D recombinant have been described, the clinical implication and the characteristics of viral mutation in the basal core promoter (BCP)/pre-core (PC) region remain unclear.

METHODS

A total of 174 chronic HBV carriers, including 115 with chronic hepatitis B, 45 with liver cirrhosis, and 14 with hepatocellular carcinoma, were enrolled. Using next-generation sequencing, the S and BCP/PC genes were determined and analyzed.

RESULTS

Genotypes B, C2, D, and C/D recombinant were detected in 1.1% (2/174), 19.5% (34/174), 0.6% (1/174) and 78.7% (137/174) of the patients, respectively. The clinical parameters and viral mutation frequency in the BCP/PC region were compared between C2- and C/D recombinant-infected patients. The distribution of C2 and C/D did not differ by disease status or liver function. Significantly higher levels of HBV DNA (6.7 ± 1.6 vs. 5.9 ± 1.5, p = 0.014), HBeAg (263.5 vs. 20.0, p = 0.013) and A1762T/G1764A double-mutations (81.0 vs. 61.8%, p = 0.018), but a lower frequency of G1896A stop mutation (33.6 vs. 76.5%, p < 0.001) was observed in patients with the C/D recombinant than in patients with genotype C2. The clonal frequencies of A1762T, G1764A, G1896A and A1846T were lower in patients with C/D than C2.

CONCLUSION

The C/D recombinant has different mutation pattern in the BCP/PC region compared with genotype C2. The lower clonal frequencies of BCP/PC mutations may explain the higher levels of HBV DNA and HBeAg in C/D-infected patients.

Hepatitis B virus serosurvey and awareness of mother-to-child transmission among pregnant women in Shenyang, China: An observational study

Preventing hepatitis B virus (HBV) mother-to-child transmission (MTCT) is the key to controlling the prevalence of chronic HBV infection. Adequate awareness of hepatitis B in hepatitis B s antigen (HBsAg) positive pregnant women may be helpful to reduce HBV MTCT.The aim of this study was to explore HBV seroprevalence among pregnant women and investigate the level of hepatitis B awareness among HBsAg positive pregnant women.HBV serum biomarkers were tested among pregnant women visiting Shengjing Hospital of China Medical University. HBsAg-positive pregnant women received a HBV DNA test and completed a questionnaire. The different HBV DNA loads were interpreted as follows: 20 to  < 2 × 10 IU/mL was low viral load, 2 × 10 to  < 2 × 10 IU/mL was intermediate viral load and ≥2 × 10 IU/mL was high viral load. The pregnant women with high viral load were treated with telbivudine (LdT). HBV DNA at different times was tested. The rate of HBV MTCT was confirmed at 28 weeks postpartum.HBsAg prevalence among pregnant women was 3.1% (441/14314). There was significant difference in comparing HBsAg prevalence in different age groups (χ = 13.86, P < .01). Among 441 HBsAg-positive pregnant women, 151 (34.2%) were hepatitis B e antigen (HBeAg) positive and 112 (25.4%) had high viral load. After 4 weeks of treatment, the average HBV DNA load of 66 cases with high viral load was (5.0 ± 0.8) log10 IU/mL. The average HBV DNA load at 4 weeks postpartum rebounded to (7.9 ± 1.0) log10 IU/mL, which was not significantly different from that at baseline (t = 1.23, P = .22). At 28 weeks postpartum, the rate of HBV MTCT in the treatment group was significantly lower than that in the observation group (0% vs 12.2%; P = .02). Only 23.4% of pregnant women knew their HBV status before gestation and 17.7% of pregnant women knew the HBV status before delivery. However, only 21.3% of pregnant women realized to need antiviral treatment to prevent MTCT.The pregnant women in Shenyang had a low HBsAg prevalence. Antiviral treatment for pregnant women with high viral load can effectively reduce the rate of HBV MTCT. HBV screening and education among HBsAg-positive pregnant women should be strengthened.

Accurate quantification of within- and between-host HBV evolutionary rates requires explicit transmission chain modelling

Analyses of virus evolution in known transmission chains have the potential to elucidate the impact of transmission dynamics on the viral evolutionary rate and its difference within and between hosts. Lin et al. (2015, Journal of Virology, 89/7: 3512–22) recently investigated the evolutionary history of hepatitis B virus in a transmission chain and postulated that the ‘colonization–adaptation–transmission’ model can explain the differential impact of transmission on synonymous and non-synonymous substitution rates. Here, we revisit this dataset using a full probabilistic Bayesian phylogenetic framework that adequately accounts for the non-independence of sequence data when estimating evolutionary parameters. Examination of the transmission chain data under a flexible coalescent prior reveals a general inconsistency between the estimated timings and clustering patterns and the known transmission history, highlighting the need to incorporate host transmission information in the analysis. Using an explicit genealogical transmission chain model, we find strong support for a transmission-associated decrease of the overall evolutionary rate. However, in contrast to the initially reported larger transmission effect on non-synonymous substitution rate, we find a similar decrease in both non-synonymous and synonymous substitution rates that cannot be adequately explained by the colonization-adaptation-transmission model. An alternative explanation may involve a transmission/establishment advantage of hepatitis B virus variants that have accumulated fewer within-host substitutions, perhaps by spending more time in the covalently closed circular DNA state between each round of viral replication. More generally, this study illustrates that ignoring phylogenetic relationships can lead to misleading evolutionary estimates.