Evolutionary patterns of hepatitis B virus quasispecies under different selective pressures: correlation with antiviral efficacy

Genetic variability of hepatitis B virus in acute and in different phases of chronic infection in Brazil

The selection pressure imposed by the host immune system impacts on hepatitis B virus (HBV) variability. This study evaluates HBV genetic diversity, nucleos(t)ide analogs resistance and HBsAg escape mutations in HBV patients under distinct selective pressures. One hundred and thirteen individuals in different phases of HBV infection were included: 13 HBeAg-positive chronic infection, 9 HBeAg-positive chronic hepatitis, 47 HBeAg-negative chronic infection (ENI), 29 HBeAg-negative chronic hepatitis (ENH) and 15 acute infected individuals. Samples were PCR amplified, sequenced and genetically analyzed for the overlapping POL/S genes. Most HBV carriers presented genotype A (84/113; 74.3%), subgenotype A1 (67/84; 79.7%), irrespective of group, followed by genotypes D (20/113; 17.7%), F (8/113; 7.1%) and E (1/113; 0.9%). Clinically relevant mutations in polymerase (tL180M/M204V) and in the Major Hydrophilic Region of HBsAg (sY100C, T118A/M, sM133T, sD144A and sG145R) were observed. Our findings, however, indicated that most polymorphic sites were located in the cytosolic loops (CYL1-2) and transmembrane domain 4 (TMD4) of HBsAg. Lower viral loads and higher HBV genetic diversity were observed in ENI and ENH groups (p < 0.001), suggesting that these groups are subjected to a higher selective pressure. Our results provide information on the molecular characteristics of HBV in a diverse clinical setting, and may guide future studies on the balance of HBV quasispecies at different stages of infection.

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.

Collateral sensitivity: An evolutionary trade-off between antibiotic resistance mechanisms, attractive for dealing with drug-resistance crisis

Background

The discovery and development of antimicrobial drugs were one of the most significant advances in medicine, but the evolution of microbial resistance limited the efficiency of these drugs.

Aim

This paper reviews the collateral sensitivity in bacteria and its potential and limitation as a new target for treating infections.

Results and Discussion

Knowledge mechanisms of resistance to antimicrobial agents are useful to trace a practical approach to treat and control of resistant pathogens. The effect of a resistance mechanism to certain antibiotics on the susceptibility or resistance to other drugs is a key point that may be helpful for applying a strategy to control resistance challenges. In an evolutionary trade-off known as collateral sensitivity, the resistance mechanism to a certain drug may be mediated by the hypersensitivity to other drugs. Collateral sensitivity has been described for different drugs in various bacteria, but the molecular mechanisms affecting susceptibility are not well demonstrated. Collateral sensitivity could be studied to detect its potential in the battle against resistance crisis as well as in the treatment of pathogens adapting to antibiotics. Collateral sensitivity-based antimicrobial therapy may have the potential to limit the emergence of antibiotic resistance.

Evidence of Residual Ongoing Viral Replication in Chronic Hepatitis B Patients Successfully Treated With Nucleos(t)ide Analogues

BACKGROUND

Chronic hepatitis B is usually treated with nucleos(t)ide analogues (NAs). However, a cure is rarely achieved, even with years of treatment. Here, we investigated whether viral replication is completely halted and how long covalently closed circular DNA (cccDNA) persists in patients successfully treated with NAs.

METHODS

A series of longitudinal serum samples and a collection of cross-sectional liver biopsies were obtained from patients successfully treated with NAs. Viral variants in serum HBV RNA were enumerated by deep sequencing. Viral replication intermediates in hepatocytes were directly visualized by in situ hybridization. The apparent half-life of each cccDNA was estimated.

RESULTS

Three of 6 successfully treated patients demonstrated clear evidence of a small proportion of virus evolution, although the overwhelming proportion of variants were identical or possessed a similar degree of divergence through time. The apparent half-life of variants was estimated to be from approximately 7.42 weeks to infinite. Hepatocytes remained positive for cytoplasmic nucleocapsids-associated relaxed circular DNA in 4 of 7 liver needle biopsies.

CONCLUSIONS

We conclude that even after prolonged treatment, a small proportion of the cccDNA reservoir is constantly replenished by continued low-level HBV replication, whereas a large proportion of the cccDNA reservoir persists over time.

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

Hepatitis B virus (HBV) remains a global public health concern, with over 250 million individuals living with chronic HBV infection (CHB) and no curative therapy currently available. Viral diversity is associated with CHB pathogenesis and immunological control of infection. Improved methods to characterize the viral genome at both the population and intra-host level could aid drug development efforts. Conventionally, HBV sequencing data are aligned to a linear reference genome and only sequences capable of aligning to the reference are captured for analysis. Reference selection has additional consequences, including sample-specific 'consensus' sequence construction. It remains unclear how to select a reference from available sequences and whether a single reference is sufficient for genetic analyses. Using simulated short-read sequencing data generated from full-length publicly available HBV genome sequences and HBV sequencing data from a longitudinally sampled individual with CHB, we investigate alternative graph-based alignment approaches. We demonstrate that using a phylogenetically representative 'genome graph' for alignment, rather than linear reference sequences, avoids issues of reference ambiguity, improves alignment, and facilitates the construction of sample-specific consensus sequences genetically similar to an individual's infection. Graph-based methods can therefore improve efforts to characterize the genetics of viral pathogens, including HBV, and may have broad implications in host pathogen research.

Compartmentalisation of Hepatitis B virus X gene evolution in hepatocellular carcinoma microenvironment and the genotype-phenotype correlation of tumorigenicity in HBV-related patients with hepatocellular carcinoma

Hepatitis B virus (HBV) exists as quasispecies (QS). However, the evolutionary characteristics of haplotypes of HBV X gene in the hepatocellular carcinoma (HCC) microenvironment remain unclear. Mutations across X gene are essential for the tumorigenicity of HBV X protein (HBx). However, the functional phenotypes of many mutant HBx remain unknown. This study aims to compare the characteristics of X gene evolution between tumour and non-tumour tissues in HCC patients and investigate the tumorigenic phenotype of HBx harbouring mutation T81P/S101P/L123S. This study included 24 HCC patients. Molecular cloning of X gene was performed to analyse characteristics of haplotypes in liver tissues. HCC cell lines stably expressing wild-type or mutant HBx and subcutaneous tumour xenograft mouse model were used to assess HBx-T81P/S101P/L123S tumorigenicity. The mean heterogeneity of HBV QS across X gene in tumour tissues was lower than that in non-tumour tissues. A location bias was observed in X gene clones with genotype C or D in tumour tissues compared to those with genotype B. Mutations in genotype-C or - D clones were mainly clustered in the dimerization region and aa110-aa140 within the transactivation region. A novel mutation combination at residues 81, 101 and 123 was identified in tumour tissues. Further, HBx-T81P/S101P/L123S promotes cell proliferation and increases genomic instability, which was mediated by MYC. This study elucidates the compartmentalized evolution patterns of HBV X gene between intra tumour and non-tumour tissues in HCC patients and provides a new mechanism underlying HBV-driven hepatocarcinogenesis, suggesting a potential viral marker for monitoring HCC.

Liver Histology of Treatment-Naïve Children with Chronic Hepatitis B Virus Infection in Shanghai China

BACKGROUND & AIMS

Chronic hepatitis B (CHB) is associated with high morbidity and mortality. We aimed to investigate associations between hepatic histology and clinical characteristics in treatment-naïve children with CHB in Shanghai, China.

METHODS

The liver biopsy specimens of 278 treatment-naïve children with CHB virus infection were scored for inflammation and fibrosis, and correlations with clinical and laboratory data were determined.

RESULTS

CHB clinical, virologic and pathologic features were studied in 278 treatment-naïve children (177 males (63.7%)) from Shanghai, China. Maternal sera were hepatitis B surface antigen (HBsAg) positive for 277 children. At biopsy, 87.4% of patients were hepatitis B e antigen positive. The median age at biopsy was 5.1 years (interquartile range 2.8-8.4 years). Hepatitis B virus (HBV) DNA levels were generally high (mean 7.4 log₁₀ IU/ml), as were levels of serum alanine aminotransferase (ALT, median 105 U/L). Using the Metavir histology activity index scoring system, no, mild, moderate and severe inflammation was seen in 2.9%, 22.3%, 73.4%, and 1.4% of patients, respectively. No fibrosis, mild fibrosis, moderate fibrosis, and cirrhosis was seen in 11.5%, 32.7%, 47.5%, and 8.3% of patients, respectively. When the serum ALT level was ≤80 (2× the upper limit of normal) and >80 U/L, the inflammation score (P<0.0001) was significantly different. And the fibrosis score was significantly different (P<0.0001), either. Inflammation and fibrosis were aggravated with increasing ALT levels. Fibrosis scores were significantly higher in children aged ≤3 than aged >3 years (P<0.0001). The rates of moderate fibrosis and cirrhosis were higher in children aged ≤3 years at biopsy. No correlations were found between histologic changes and sex, HBV genotype or HBV DNA level.

CONCLUSION

Substantial heterogeneity in inflammatory and fibrotic levels was observed in treatment-naïve children with CHB in Shanghai, China. Serum ALT levels >80 U/L may be a strong indicator of the degree of hepatic inflammation and fibrosis severity. Moderate fibrosis and cirrhosis can appear in children aged <3 years or younger.

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.

Early dynamic changes of quasispecies in the reverse transcriptase region of hepatitis B virus in telbivudine treatment

BACKGROUND

Telbivudine (LdT) - a synthetic thymidine β-L-nucleoside analogue (NA) - is an effective inhibitor for hepatitis B virus (HBV) replication. The quasispecies spectra in the reverse transcriptase (RT) region of the HBV genome and their dynamic changes associated with LdT treatment remains largely unknown.

METHODS

We prospectively recruited a total of 21 treatment-naive patients with chronic HBV infection and collected sequential serum samples at five time points (baseline, weeks 1, 3, 12, and 24 after LdT treatment). The HBV RT region was amplified and shotgun-sequenced by the Ion Torrent Personal Genome Machine (PGM)® system. We reconstructed full-length haplotypes of the RT region using an integrated bioinformatics framework, including de novo contig assembly and full-length haplotype reconstruction. In addition, we investigated the quasispecies' dynamic changes and evolution history and characterized potential NAs resistant mutations over the treatment course.

RESULTS

Viral quasispecies differed obviously between patients with complete (n = 8) and incomplete/no response (n = 13) at 12 weeks after LdT treatment. A reduced d_N/d_S ratio in quasispecies demonstrated a selective constraint resulting from antiviral therapy. The temporal clustering of sequential quasispecies showed different patterns along with a 24-week observation, although its statistic did not differ significantly. Several patients harboring pre-existing resistant mutations showed different clinical responses, while NAs resistant mutations were rare within a short-term treatment.

CONCLUSION

A complete profile of quasispecies reconstructed from in-depth shotgun sequencing may has important implications for enhancing clinical decision in adjusting antiviral therapy timely.

Compartmentalized evolution of hepatitis B virus contributes differently to the prognosis of hepatocellular carcinoma

Serum hepatitis B virus (HBV) mutations can predict hepatocellular carcinoma (HCC) occurrence. We aimed to clarify if HBV evolves synchronously in the sera, adjacent liver and tumors and predict HCC prognosis equally. A total of 203 HBV-positive HCC patients with radical hepatectomy in Shanghai, China, during 2011-15 were enrolled in this prospective study. Quasispecies complexity (QC) in HBV core promoter region was assessed using clone-based sequencing. We performed RNA sequencing on tumors and paired adjacent tissues of another 15 HCC patients and analyzed it with three public data sets containing 127 samples. HBV QC was positively correlated to APOBEC3s' expression level (r = 0.28, P < 0.001), higher in the adjacent tissues than in the tumors (P = 6.50e-3), and higher in early tumors than in advanced tumors (P = 0.039). The evolutionary distance between the sera-derived HBV strains and the tumor-derived ones was significantly longer than that between the sera-derived ones and the adjacent tissue-derived ones (P < 0.001). Multivariate Cox regression analyses indicated that high HBV QC in the sera predicted an unfavorable overall survival (P = 0.002) and recurrence-free survival (RFS; P = 0.004) in HCC, whereas, in the tumors, it predicted a favorable RFS (P < 0.001). APOBECs-related HBV mutations, including G1764A, were more frequent in the sera than in the adjacent tissues. High-frequent A1762T/G1764A in the sera predicted an unfavorable RFS (P < 0.001), whereas, in the tumors, it predicted a favorable RFS (P = 0.035). In conclusion, HBV evolves more advanced in the sera than in the tumors. HBV QC and A1762T/G1764A in the sera and tumors have contrary prognostic effects in HCC.

Mutational characterization of HBV reverse transcriptase gene and the genotype-phenotype correlation of antiviral resistance among Chinese chronic hepatitis B patients

Background and Aims: The drug resistance of hepatitis B virus (HBV) originates from mutations within HBV reverse transcriptase (RT) region during the prolonged antiviral therapy. So far, the characteristics of how these mutations distribute and evolve in the process of therapy have not been clarified yet. Thus we aimed to investigate these characteristics and discuss their contributing factors.

Methods: HBV RT region was direct-sequenced in 285 treatment-naive and 214 post-treatment patients. Mutational frequency and Shannon entropy were calculated to identify the specific mutations differing between genotypes or treatment status. A typical putative resistance mutation rtL229V was further studied using in-vitro susceptibility assays and molecular modeling.

Results: The classical resistance mutations were rarely detected among treatment-naive individuals, while the putative resistance mutations were observed at 8 AA sites. rtV191I and rtA181T/V were the only resistance mutations identified as genotype-specific mutation. Selective pressure of drug usage not only contributed to the classical resistance mutations, but also induced the changes at a putative resistance mutation site rt229. rtL229V was the major substitution at the site of rt229. It contributed to the most potent suppression of viral replication and reduced the in-vitro drug susceptibility to entecavir (ETV) when coexisting with rtM204V, consistent with the hypothesis based on the molecular modeling and clinical data analysis.

Conclusions: The analysis of mutations in RT region under the different circumstances of genotypes and therapy status might pave the way for a better understanding of resistance evolution, thus providing the basis for a rational administration of antiviral therapy.

Dynamic changes of HCV genomes under selective pressure from DAAs therapy in relapsed patients

Currently, direct-acting antiviral drugs (DAAs) are widely used as therapeutic methods for hepatitis C virus (HCV)-positive patients, however, patients may experience treatment failure, and the dynamic changes of HCV genomes in these patients are unknown. In this study, three real-world DAAs cohorts were enrolled to observe clinical efficacy. In addition, serum samples from treatment failure patients at baseline and relapse were used to analyze changes of the HCV genomes at near full-length genome level, including resistance-associated variants (RAVs), viral quasispecies diversity and selection analysis. Next-generation sequencing was used as the detection method. The overall sustained virological response at 12 w after the end of treatment was achieved in 91.9% (57/62) of HCV patients, and 3 paired samples obtained from relapsed patients. All the 3 patients harbored baseline NS5A RAVs, the frequency of NS5A RAVs increased in 2 patients and a new NS5A RAV emerged in 1 patient at relapse, and almost all the viral strains existed with NS5A RAVs at relapse. The results of the viral quasispecies diversity analysis revealed that viral quasispecies diversity decreased at relapse compared to baseline, and the results of selection analysis indicated that the virus population experienced a bottleneck phenomenon, recent selective sweep and population expansion or was under purification selection after DAAs treatment. This study indicated that the clinical efficacy was excellent in real-world DAAs cohorts, and the viral strains existed at relapse were selective by DAAs therapy.

Viral quasispecies quantitative analysis: a novel approach for appraising the immune tolerant phase of chronic hepatitis B virus infection

Few non-invasive models were established for precisely identifying the immune tolerant (IT) phase from chronic hepatitis B (CHB). This study aimed to develop a novel approach that combined next-generation sequencing (NGS) and machine learning algorithms using our recently published viral quasispecies (QS) analysis package. 290 HBeAg positive patients from whom liver biopsies were taken were enrolled and divided into a training group (n = 148) and a validation group (n = 142). HBV DNA was extracted and QS sequences were obtained by NGS. Hierarchical clustering analysis (HCA) and principal component analysis (PCA) based on viral operational taxonomic units (OTUs) were performed to explore the correlations among QS and clinical phenotypes. Three machine learning algorithms, including K-nearest neighbour, support vector machine, and random forest algorithm, were used to construct diagnostic models for IT phase classification. Based on histopathology, 90 IT patients and 200 CHB patients were diagnosed. HBsAg titres for IT patients were higher than those of CHB patients (p < 0.001). HCA and PCA analysis grouped IT and CHB patients into two distinct clusters. The relative abundance of viral OTUs differed mainly within the BCP/precore/core region and was significantly correlated with liver inflammation and fibrosis. For the IT phase classification, all machine-learning models showed higher AUC values compared to models based on HBsAg, APRI, and FIB-4. The relative abundance of viral OTUs reflects the severity of liver inflammation and fibrosis. The novel QS quantitative analysis approach could be used to diagnose IT patients more precisely and reduce the need for liver biopsy.

New therapeutic options for persistent low-level viremia in patients with chronic hepatitis B virus infection: Increase of entecavir dosage

Chronic hepatitis B virus (HBV) infection (CHB) is a public health concern worldwide. Current therapies utilizing nucleos(t)ide analogs (NA) have not resulted in a complete cure for CHB. Furthermore, patients on long-term NA treatment often develop low-level viremia (LLV). Persistent LLV, in addition to causing the progression of liver disease or hepatocellular carcinoma, may shed light on the current plight of NA therapy. Here, we review the literature on LLV, NA treatment, and various doses of entecavir to find a strategy for improving the efficacy of this antiviral agent. For LLV patients, three therapeutic options are available, switching to another antiviral monotherapy, interferon-α switching therapy, and continuing monotherapy. In real-world clinical practice, entecavir overdose has been used in antiviral therapy for CHB patients with NA refractory and persistent LLV, which encouraged us to conduct further in-depth literature survey on dosage and duration related entecavir studies. The studies of pharmacodynamics and pharmacokinetics show that entecavir has the maximal selected index for safety, and has great potential in inhibiting HBV replication, in all of the NAs. In the particular section of the drug approval package published by the United States Food and Drug Administration, entecavir doses 2.5-20 mg/d do not increase adverse events, and entecavir doses higher than 1.0 mg/d might improve the antiviral efficacy. The literature survey led us to two suggestions: (1) Increasing entecavir dose to 1.0 mg/d for the treatment of NA naïve patients with HBV DNA >2 × 10⁶ IU/mL is feasible and would provide better prognosis; and (2) Further research is needed to assess the long-term toxic effects of higher entecavir doses (2.5 and 5.0 mg/d), which may prove beneficial in treating patients with prior NA treatment, partial virological response, or LLV state.

The association of the heterogeneity of HBV reverse transcriptase quasispecies with antiviral efficacy after treatment with nucleos(t)ide analogues for 10 years

To assess the heterogeneity of HBV reverse transcriptase (RT) quasispecies during 10 years of antiviral therapy and their association with antiviral efficacy. Nineteen patients with chronic hepatitis B (CHB) infection receiving nucleos(t)ide analogues (NAs) were enrolled. Based on the antiviral efficacy after 1 year of treatment, 5 patients were grouped into an early virologic response (EVR) group, while 8 patients were grouped into a late virologic response (LVR) group. Furthermore, 6 CHB patients that had undergone antiviral treatment for 10 years were grouped into a virologic breakthrough (VBT) group. The HBV RT from each patient were amplified, cloned, and sequenced. The complexity of the RT gene in the EVR group was significantly higher than that in the LVR (P = 0.0393) and VBT groups (P = 0.0141). Phylogenetic tree analysis showed that the average branch length of the EVR and LVR groups were significantly greater than that of VBT group (P < 0.001). The complexity (at the nucleotide level) of the RT quasispecies was negatively correlated with the corresponding HBV DNA load (P = 0.0163) at one year post-antiviral treatment. Moreover, both the LVR and VBT groups accumulated more deleterious mutations than the EVR group. After 1 year of NAs treatment, the increased HBV quasispecies complexity and evolutionary topologies, coupled with less deleterious mutations, are likely associated with a favorable efficacy during long-term antiviral treatment.

Comparison and Correlation of Genetic Variability of the HBV Pre-S Region in HIV/HBV Co-Infected Patients: Quasispecies Perspective

Background

Human immunodeficiency virus (HIV)/hepatitis B virus (HBV) co-infection can accelerate HBV-induced liver disease. A previous study showed that variation in the HBV pre-S region and quasispecies heterogeneity (Sn, mean genetic distance, dS, dN, and dS/dN) are both related to HBV-induced terminal liver disease in HBV mono-infection. Currently, data are lacking on quasispecies variation of the HBV pre-S region in HIV/HBV co-infection. Investigating the quasispecies variation of the HBV pre-S region and its related factors in HIV/HBV co-infection will help to better explore the pathogenic mechanism of HIV/HBV co-infection.

Methods

According to the HIV antibody results obtained before treatment, chronic HBV-infected patients were divided into HIV/HBV co-infected and HBV mono-infected groups. The clinical characteristics of all patients were collected, and DNA was extracted from the serum. The HBV pre-S region was amplified by nested PCR and was further TA cloned. BioEdit software 7.0 was used for sequence alignment with reference to the standard sequence of the matched HBV genotype. We used 1:1 propensity score matching (PSM) to control for baseline confounding factors between the two groups.

Results

After 1:1 PSM, we identified 100 patients with similar propensities: 50 HIV/HBV co-infected patients and 50 HBV mono-infected patients. HBV quasispecies indices were lower in the HIV/HBV co-infected group than those in the HBV mono-infected group. A significant correlation was observed between all quasispecies indices and soluble cluster of differentiation 163 (sCD163) and interleukin-18 (IL-18) in the HIV/HBV co-infected group; however, this phenomenon was not found in the HBV mono-infected group.

Conclusion

Combined HIV infection reduces quasispecies heterogeneity in the HBV pre-S region, and the quasispecies heterogeneity is related to the sCD163 and IL-18 levels.

Comparison of Serum Hepatitis B Virus RNA Levels and Quasispecies Evolution Patterns between Entecavir and Pegylated-Interferon Mono-treatment in Chronic Hepatitis B Patients

Hepatitis B virus (HBV) RNA may independently predict virological and serological response. This study aimed to compare dynamic changes in serum HBV RNA levels and HBV quasispecies evolution patterns between entecavir and pegylated-interferon mono-treatment in chronic hepatitis B patients and to determine the clinical significance during treatment. TaqMan real-time PCR was used for quantitative analysis. HBV RNA levels were retrospectively determined in serial serum samples from 178 chronic hepatitis B patients who received either entecavir or pegylated-interferon treatment.

Hepatitis B virus (HBV) RNA may independently predict virological and serological response. This study aimed to compare dynamic changes in serum HBV RNA levels and HBV quasispecies evolution patterns between entecavir and pegylated-interferon mono-treatment in chronic hepatitis B patients and to determine the clinical significance during treatment. TaqMan real-time PCR was used for quantitative analysis. HBV RNA levels were retrospectively determined in serial serum samples from 178 chronic hepatitis B patients who received either entecavir or pegylated-interferon treatment. Both serum HBV DNA and RNA quasispecies were analyzed via next-generation sequencing. Receiver operating characteristics (ROC) analysis was performed to evaluate the prediction value of individual biomarkers for hepatitis B e antigen (HBeAg) seroconversion. Patients who received pegylated-interferon treatment showed stronger declines in HBV RNA levels than did those who received entecavir treatment. Serum HBV RNA levels were lower in patients with subsequent HBeAg seroconversion. At baseline, the level of HBV RNA was better than other indicators in predicting HBeAg seroconversion. Moreover, the predictive value of serum HBV RNA levels was better in the entecavir group. Baseline HBV RNA exhibited a significantly higher genetic diversity than HBV DNA and had a significant decline after 4 weeks of entecavir treatment. Higher baseline genetic diversity may result in a better outcome in pegylated-interferon-treated patients. Serum HBV RNA levels showed different decline kinetics, and HBV RNA quasispecies showed different evolution patterns in entecavir and pegylated-interferon mono-treatment. Taken together, serum HBV RNA may serve as a promising biomarker of HBeAg seroconversion in patients during antiviral treatment.

Pre-S Deletions are Predominant Quasispecies in HIV/HBV Infection: Quasispecies Perspective

Background

Combined HIV infection can accelerate HBV-induced liver disease. It is known that HBV Pre-S deletion is closely related to HBV-associated terminal liver disease in HBV mono-infection. Currently, data on HBV Pre-S quasispecies feature deletion in HIV/HBV co-infected patients are lacking.

Methods

The characteristics and blood samples of patients with chronic HBV infection were collected and classified into an HIV/HBV co-infection group and an HBV mono-infection group according to HIV antibody results before treatment. HBV DNA in serum was extracted. The HBV Pre-S region was amplified by nested-PCR and was further T-A cloned. Using the standard sequence of the matched genotype HBV as a reference, BioEdit 7.0 software was employed for sequence alignment.

Results

HBV Pre-S regions were successfully amplified from 147 patients, including 71 cases in the HIV/HBV co-infected group and 76 cases in the HBV mono-infected group. The proportion of the HIV/HBV co-infected group with Pre-S quasispecies deletion was lower than that of the HBV mono-infected group. By analyzing the frequency of Pre-S quasispecies in the two groups, the frequency of Pre-S quasispecies in HIV/HBV co-infected patients with Pre-S quasispecies was higher than HBV mono-infected patients. The frequency of Pre-S quasispecies deletion of the S protein promoter region in the HIV/HBV co-infected group was significantly higher than that in the HBV mono-infected group.

Conclusion

High-frequency Pre-S quasispecies deletions are predominant in HIV/HBV co-infected patients; however, low-frequency Pre-S deletions are predominant in HBV mono-infected patients, providing a reference for the pathogenesis of the accelerated progression of liver disease in HIV/HBV co-infection.

An integrated software for virus community sequencing data analysis

BACKGROUND

A virus community is the spectrum of viral strains populating an infected host, which plays a key role in pathogenesis and therapy response in viral infectious diseases. However automatic and dedicated pipeline for interpreting virus community sequencing data has not been developed yet.

RESULTS

We developed Quasispecies Analysis Package (QAP), an integrated software platform to address the problems associated with making biological interpretations from massive viral population sequencing data. QAP provides quantitative insight into virus ecology by first introducing the definition "virus OTU" and supports a wide range of viral community analyses and results visualizations. Various forms of QAP were developed in consideration of broader users, including a command line, a graphical user interface and a web server. Utilities of QAP were thoroughly evaluated with high-throughput sequencing data from hepatitis B virus, hepatitis C virus, influenza virus and human immunodeficiency virus, and the results showed highly accurate viral quasispecies characteristics related to biological phenotypes.

CONCLUSIONS

QAP provides a complete solution for virus community high throughput sequencing data analysis, and it would facilitate the easy analysis of virus quasispecies in clinical applications.

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

OBJECTIVE

This study was performed to explore the correlation between the characteristics of hepatitis B virus (HBV) quasispecies in HBV-infected pregnant women and the risk of immunoprophylaxis failure for their infants.

DESIGN

In this prospective nested case-control study, the characteristics of HBV quasispecies in mothers whose infants were immunoprophylaxis success (control group) and those whose infants were immunoprophylaxis failure (case group) were analysed by the clone-based sequencing of full-length HBV genome and next-generation sequencing (NGS) of "a" determinant region, and were compared between the two groups.

RESULTS

The quasispecies characteristics including mutant frequency, Shannon entropy and mean genetic distance at amino acid level of "a" determinant region were significantly lower in case group than that in control group, using the full-length HBV genome clone-based sequencing assay. These results were confirmed by NGS assay. Notably, we discovered that the differences were also significant at nucleotide level by NGS assay. Furthermore, the risk of immunoprophylaxis failure could be predicted by analysing the three HBV quasispecies characteristics either at nucleotide level or at amino acid level of "a" determinant region, and the corresponding predictive values were tentatively set up.

CONCLUSIONS

HBV quasispecies with a more complex mutant spectrum in "a" determinant region might be more vulnerable to extinct through mother-to-child-transmission (MTCT). More importantly, analysing HBV quasispecies characteristics in pregnant women with high HBV DNA load might be helpful to predict the high-risk population of immunoprophylaxis failure, and consequently provide accurate intervention against MTCT of HBV.

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.

Analysis of HBV X gene quasispecies characteristics by next-generation sequencing and cloning-based sequencing and its association with hepatocellular carcinoma progression

OBJECTIVES

This study aimed to describe the differences between next-generation sequencing (NGS) and cloning-based sequencing (CBS) in HBX quasispecies research and primitively investigate the relationship between the dominant HBX quasispecies and hepatocellular carcinoma (HCC).

METHODS

A total of 12 serum samples were collected. Serum hepatitis B virus (HBV) DNA was extracted, and the HBV X-region (HBX) was amplified by nested polymerase chain reaction (PCR). The PCR products were simultaneously tested with NGS and CBS to detect quasispecies of the HBX.

RESULTS

A total of 9348 eligible quasispecies sequences were obtained by NGS, which were much larger than the 98 of that by CBS. By the phylogenetic tree, the dominant quasispecies sequence of each sample could be found, although they had several nucleotides differences between the dominant quasispecies sequences found by CBS and NGS. By comparing the quasispecies heterogeneity, it was found that the quasispecies complexity value of HBV X-region obtained by NGS was higher than CBS (P < 0.05). The diversity values, including d, dS, dN, an d d N/ dS obtained by NGS were lower than by CBS (all of P < 0.01). The relativity of Spearman(rs) in d, dS, and dN were statistically significant (rs_ d = 0.865, P = 0.001; rs_ dS = 0.722, P = 0.014; and rs_ dN = 0.738, P = 0.011, respectively). There were 21 different bases between the HBX quasispecies of case A and control B.

CONCLUSION

The results of this can be used as guidance when researchers plan to choose a suitable method to study quasispecies, especially the HBV X gene quasispecies. Some high-risk mutations of HBX quasispecies were also found in this study and their relationship with HCC need deeper exploration.

Viral quasispecies of hepatitis B virus in patients with YMDD mutation and lamivudine resistance may not predict the efficacy of lamivudine/adefovir rescue therapy

The association between hepatitis B virus (HBV) quasispecies (QS) and the efficacy of nucleos(t)ide analog therapy is currently not well defined, particularly in the case of lamivudine (LAM)/adefovir (ADV) combination rescue therapy for patients with chronic HBV infection (CHB) presenting with LAM resistance. In the present study, 16 CHB patients with the rtM204I/V mutation in the tyrosine-methionine-aspartate-aspartate motif of the C domain of the polymerase gene who switched to LAM/ADV treatment due to LAM resistance were assessed. HBV DNA was isolated from these patients and the reverse transcriptase (RT) region was sequenced. The QS heterogeneity and distribution was analyzed, the mutation sites were recorded and the phylogenetic trees were constructed. The results indicated that QS heterogeneity and distribution in the RT and S regions were not significantly different between responders (RS) and non-RS (NRS) at baseline (P>0.05), except for the higher frequency of a dominant strain in the RT region at the nucleotide level in the RS group (P=0.039). In addition, in NRS, no significant difference in QS heterogeneity or distribution in these regions was identified at six months vs. the baseline. Furthermore, although in the non-responder group the frequency of the LAM resistance-associated mutations (rtM204V/I) decreased at 6 months compared with the baseline, it did not disappear in any of the patients after six months of treatment. Analysis of individual patients did not indicate any consistent selection of specific HBV mutants during LAM/ADV rescue therapy. In conclusion, the baseline HBV QS within the RT and S regions may not be a valid predictor of the response to LAM/ADV rescue treatment in CHB patients with LAM resistance.

Clinical Manifestations and Laboratory Tests of AECHB and Severe Hepatitis (Liver Failure)

This chapter describes the clinical symptoms and signs of AECHB and HBV ACLF, classification, grading of HBV ACLF and their features, diagnostic principles and standards in liver pathology, biochemistry, and virology of HBV ACLF.

Liver failure is defined as serious damage to the liver cause by a variety of etiologies, leading to liver function disorder or even decompensation, and clinical syndromes with coagulopathy, jaundice, hepatic encephalopathy, and ascites.

Severe hepatitis B can be indicated pathologically by apparent hepatocellular necrosis, including extensive multifocal, confluent, bridging, sub-massive or massive necrosis.

Laboratory tests during the course of severe exacerbation of chronic hepatitis B can reflect pathological changes and liver function in a timely manner, providing objective and informative reference data for evaluation of disease severity and treatment efficacy. Among the most important laboratory tests are those for prothrombin activity, international normalized ratio, and increases in total bilirubin concentration.

Severe hepatitis B is associated with interactions between the virus and host factors. Detection of HBV DNA, HBV genotype, quasispecies and HBV mutation can provide important theoretical bases for the prevention, control or mitigation of the progress of severe hepatitis B.

Noninvasive imaging modalities can be used to visualize the entire liver and parts of it. Measuring liver volume to evaluate liver size and liver reserve capacity is regarded as important in diagnosis, surgical approach and prognostic evaluation of patients with severe exacerbation of chronic hepatitis B and liver failure.

Model for End-Stage Liver Disease (MELD) is the first quantitative method developed to assess whether a patient with liver failure requires a liver transplant. The predictive value of the MELD model has been improved by the MELD-Na, iMELD, and MESO models. Several other valuable prognostic models have been developed. For example, for patients with HBV-ACLF, the established TPPM scoring system was found to be more predictive than MELD score.

How viral genetic variants and genotypes influence disease and treatment outcome of chronic hepatitis B. Time for an individualised approach?

Chronic hepatitis B virus (HBV) infection remains a global problem. Several HBV genotypes exist with different biology and geographical prevalence. Whilst the future aim of HBV treatment remains viral eradication, current treatment strategies aim to suppress the virus and prevent the progression of liver disease. Current strategies also involve identification of patients for treatment, namely those at risk of progressive liver disease. Identification of HBV genotype, HBV mutants and other predictive factors allow for tailoured treatments, and risk-surveillance pathways, such as hepatocellular cancer screening. In the future, these factors may enable stratification not only of treatment decisions, but also of patients at risk of higher relapse rates when current therapies are discontinued. Newer technologies, such as next-generation sequencing, to assess drug-resistant or immune escape variants and quasi-species heterogeneity in patients, may allow for more information-based treatment decisions between the clinician and the patient. This article serves to discuss how HBV genotypes and genetic variants impact not only upon the disease course and outcomes, but also current treatment strategies. Adopting a personalised genotypic approach may play a role in future strategies to combat the disease. Herein, we discuss new technologies that may allow more informed decision-making for response guided therapy in the battle against HBV.

A comparative study on the characterization of hepatitis B virus quasispecies by clone-based sequencing and third-generation sequencing

Hepatitis B virus (HBV) has a high mutation rate due to the extremely high replication rate and the proofreading deficiency during reverse transcription. The generated variants with genetic heterogeneity are described as viral quasispecies (QS). Clone-based sequencing (CBS) is thought to be the ‘gold standard’ for assessing QS complexity and diversity of HBV, but an important issue about CBS is cost-effectiveness and laborious. In this study, we investigated the utility of the third-generation sequencing (TGS) DNA sequencing to characterize genetic heterogeneity of HBV QS and assessed the possible contribution of TGS technology in HBV QS studies. Parallel experiments including 3 control samples, which consisted of HBV full gene genotype B and genotype C plasmids, and 10 patients samples were performed by using CBS and TGS to analyze HBV whole-genome QS. Characterization of QS heterogeneity was conducted by using comprehensive statistical analysis. The results showed that TGS had a high consistency with CBS when measuring the complexity and diversity of QS. In addition, to detect rare variants, there were strong advantages conferred by TGS. In summary, TGS was considered to be practicable in HBV QS studies and it might have a relevant role in the clinical management of HBV infection in the future.

Relationship between serum HBV-RNA levels and intrahepatic viral as well as histologic activity markers in entecavir-treated patients

BACKGROUND & AIMS

In diagnostics, serum hepatitis B virus (HBV)-RNA levels are valuable when the HBV-DNA load in circulation is effectively suppressed by nucleos(t)ide analogue (NUC) therapy. This study aimed to determine the intrahepatic viral replication activity reflected in serum HBV-RNA and whether HBV-RNA contributes to liver histological changes in patients treated with NUC.

METHODS

A cross-sectional set of serum and liver biopsy samples was obtained from patients treated with entecavir, who had undetectable levels of serum HBV-DNA. The correlations between serum HBV-RNA concentration and levels of peripheral and intrahepatic viral replicative forms, as well as histological scores, were analyzed. Quasispecies of serum HBV-RNA and intrahepatic viral replicative forms were examined by deep sequencing. HBV-RNA-positive hepatocytes were visualized by in situ hybridization.

RESULTS

Serum HBV-RNA was detected in 35 of 47 patients (74.47%, 2.33-4.80log10copies/ml). These levels correlated not only with the intrahepatic HBV-RNA level and the ratio of intrahepatic HBV-RNA to covalently closed circular DNA (cccDNA), but also with the histological scores for grading and staging. Regarding quasispecies, serum HBV-RNA was dynamic and more genetically homogenous to simultaneously sampled intrahepatic HBV-RNA than to the cccDNA pool. In situ histology revealed that HBV-RNA-positive hepatocytes were clustered in foci, sporadically distributed across the lobules, and co-localized with hepatitis B surface antigen.

CONCLUSION

Serum HBV-RNA levels reflect intrahepatic viral transcriptional activity and are associated with liver histopathology in patients receiving NUC therapy. Our study sheds light on the nature of HBV-RNA in the pathogenesis of chronic HBV infection and has implications for the management of chronic hepatitis B during NUC therapy.

LAY SUMMARY

Serum HBV-RNA levels are indicative of the intrahepatic transcriptional activity of covalently closed circular DNA and are associated with liver histological changes in patients with chronic B hepatitis who are receiving nucleos(t)ide analogue therapy.

Quantifying perinatal transmission of Hepatitis B viral quasispecies by tag linkage deep sequencing

Despite full immunoprophylaxis, mother-to-child transmission (MTCT) of Hepatitis B Virus still occurs in approximately 2–5% of HBsAg positive mothers. Little is known about the bottleneck of HBV transmission and the evolution of viral quasispecies in the context of MTCT. Here we adopted a newly developed tag linkage deep sequencing method and analyzed the quasispecies of four MTCT pairs that broke through immunoprophylaxis. By assigning unique tags to individual viral sequences, we accurately reconstructed HBV haplotypes in a region of 836 bp, which contains the major immune epitopes and drug resistance mutations. The detection limit of minor viral haplotypes reached 0.1% for individual patient sample. Dominance of “a determinant” polymorphisms were observed in two children, which pre-existed as minor quasispecies in maternal samples. In all four pairs of MTCT samples, we consistently observed a significant overlap of viral haplotypes shared between mother and child. We also demonstrate that the data can be potentially useful to estimate the bottleneck effect during HBV MTCT, which provides information to optimize treatment for reducing the frequency of MTCT.

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

OBJECTIVES

To identify within-host quasispecies characteristics of hepatitis B virus (HBV) in mothers and children infected via mother-to-child transmission (MTCT).

METHODS

Using next-generation sequencing (NGS), we analyzed sequences within the non-overlapping pre-core/core (pre-C/C) gene in 37 mother-child pairs.

RESULTS

Phylogenetic and Highlighter analyses suggested that both a single strain and multiple distinct strains may be transmitted in MTCT of HBV. However, analysis of reassembled viral sequences revealed a relatively narrow distribution of variants in children, which was confirmed by a lower viral diversity in children than that in mothers. New closely related variants with combinations of two to five high-frequency mutations were observed in seven children with elevated ALT levels; the new variants out-competed the transmitted maternal variants to become the dominant strains in five of them. Furthermore, 30 mutations with a frequency >1% of all viruses within-host were present in those children; the mutations caused 19 amino-acid substitutions. Interestingly, almost all were located within the well-known T-cell or B-cell epitopes.

CONCLUSIONS

There are restrictive changes that occur in the early stages of chronic HBV infection through MTCT with different clinical consequences. These data might have important implications for future investigations of interrelated immunopathogenesis and therapeutic strategies.

Clinical features and viral quasispecies characteristics associated with infection by the hepatitis B virus G145R immune escape mutant

Coexistence of the hepatitis B surface antigen (HBsAg) and hepatitis B surface antibody (anti-HBs) is an uncommon phenomenon, and the underlying mechanisms remain largely unknown. Amino-acid (aa) substitution from glycine to arginine at aa 145 (G145R), in the major hydrophilic region, has been reported in patients with HBsAg and anti-HBs coexistence. However, there is limited knowledge about the clinical features and viral quasispecies characteristics associated with G145R mutant hepatitis B virus (HBV) infection. We herein describe the dynamic changes in the serological and virological markers in a case of hepatitis B with coexisting HBsAg and anti-HBs, caused by a G145R immune escape mutant (genotype C). Entecavir was administered during the 4th week after admission. Alanine aminotransferase peaked in the 16th week, while both the HBsAg and HBeAg declined rapidly. HBsAg clearance and hepatitis B e antigen (HBeAg)/hepatitis B e antibody (anti-HBe) seroconversion were achieved in the 36th week, and then entecavir was withdrawn. A follow-up of 96 weeks showed that HBV DNA remained undetectable and that anti-HBs was maintained above 100 mIU/mL. The quasispecies characteristics of the G145R mutant HBV were investigated via ultra-deep sequencing. The complexity and genetic distance of the S and RT regions were much higher in the 8th week than at baseline or in the 4th week. Moreover, the frequencies of mutations (L173P, Q181R and A184V) in cytotoxic T lymphocyte epitopes increased before entecavir treatment. These findings extend understanding of the evolution of HBV under host immune pressure and of the clinical outcomes of affected patients.

HBV quasispecies composition in Lamivudine-failed chronic hepatitis B patients and its influence on virological response to Tenofovir-based rescue therapy

The present study sought to evaluate the structure of HBV quasispecies in Lamivudine (LMV)-failed chronic hepatitis B (CHB) patients and its impact in defining the subsequent virological responses to Tenofovir (TDF)-based rescue-therapy. By analyzing HBV clones encompassing reverse transcriptase (RT) and surface (S) region from LMV-failed and treatment-naïve CHB patients, we identified 5 classical and 12 novel substitutions in HBV/RT and 9 substitutions in immune-epitopes of HBV/S that were significantly associated with LMV failure. In silico analysis showed spatial proximity of some of the newly-identified, mutated RT residues to the RT catalytic centre while most S-substitutions caused alteration in epitope hydrophobicity. TDF administration resulted in virological response in 60% of LMV-failed patients at 24-week but non-response in 40% of patients even after 48-weeks. Significantly high frequencies of 6 S-substitutions and one novel RT-substitution, rtH124N with 6.5-fold-reduced susceptibility to TDF in vitro, were noted at baseline in TDF non-responders than responders. Follow-up studies depicted greater evolutionary drift of HBV quasispecies and significant decline in frequencies of 3 RT and 6 S-substitutions in responder-subgroup after 24-week TDF-therapy while most variants persisted in non-responders. Thus, we identified the HBV-RT/S variants that could potentially predict unfavorable response to LMV/TDF-therapy and impede immune-mediated viral clearance.

Characteristics of CpG Islands and their quasispecies of full-length hepatitis B virus genomes from patients at different phases of infection

BACKGROUND

CpG islands in hepatitis B virus (HBV) genome are potential targets for methylation mediated gene silencing, and may be involved in the pathogenesis of HBV infection. To date, their characteristics in HBV quasispecies (QS) remain largely unknown. The purpose of this study was to investigate the characteristics of CpG islands in HBV QS.

METHODS

Forty patients diagnosed as acute hepatitis B (AHB, n = 10), immune-tolerant HBV carriers (IT, n = 9), chronic hepatitis B (CHB, n = 11), or acute on chronic liver failure (ACLF, n = 10), were enrolled in this case-control study. A total of 599 clones were isolated, and full-length HBV genomes were sequenced.

RESULTS

CpG island II (CGII) in AHB group was shorter in length and its QS heterogeneity was lower than that in the chronic infection group. Among the chronic infection subgroups, CGII and CpG island III (CGIII) in IT group were longer and their heterogeneity was lower compared to CHB and ACLF groups. Length of CGII correlated with HBV DNA levels positively while the complexity and diversity of CGII correlated with HBV DNA levels negatively. Moreover, CGII and CGIII were shorter in genotype B than those in genotype C, while QS complexity and diversity of either CGII or CGIII had no significant difference between genotype B and C.

CONCLUSIONS

Overall, our results suggest that the distribution, length and QS heterogeneity of CpG islands in full-length HBV genome differ across clinical phases of infection, of which the mechanism warrants further study.

A strong conservative tendency in HBV transcriptase (RT): a majority of natural RT mutations derived from the S gene

BACKGROUND & AIMS

Little is known about natural mutations in the HBV reverse transcriptase (RT) region. Our study aimed to characterize the natural RT mutation along the natural course of chronic Hepatitis B (CHB).

METHODS

Sixty CHB patients (immune-tolerant phase, IT, n = 20; immune-active phase, IA, n = 20 and inactive carriers phase, IC, n = 20) were selected from the Focal study, including 25 subjects with median 18 months follow-up. Mutations were evaluated at both RT and main S protein encoding region by clone-based sequencing.

RESULTS

The HBV RT quasispecies had significant lower heterogeneity in IT than IA and IC phases (P < 0.05), but not between IA and IC phases (P > 0.05). Limited heterogeneity over time was further confirmed in a longitudinal study. Locations of RT mutations were primarily located in the interdomians and the lowest in functional domains in each phase. Mutations in human leukocyte antigen (HLA) I epitopes (IT, 0.95%; IA, 1.31%; IC, 1.28%, P < 0.05) and HLA II epitopes (IT, 0.70%; IA, 0.90%; IC, 1.45%, P < 0.01) varied significantly over time. More frequent mutations were detected in the ORF of S gene from the same clones (HBsAg vs. RT: IT, 75 vs. 45; IA, 83 vs. 64; IC, 80 vs. 65). The majority of RT mutations were shared with genetic changes in the main S gene.

CONCLUSIONS

Our findings suggested that HBV RT showed a strong conservative tendency and a majority of their natural mutations were derived from the same genetic changes in the S gene.

Efficacy of interferon for chronic hepatitis B in patients with nucleoside and nucleotide combination therapy failure

BACKGROUND AND AIM

In China, inappropriate therapies with nucleos(t)ide analogues (NA) have induced hepatitis B virus resistance, combination therapy with nucleoside and nucleotide (ComTNsNt) failure, or multi-drug resistant mutations. However, the efficacy of combination therapy with entecavir plus tenofovir for ComTNsNt failure is limited. In the current study, the regimens of interferon-α (IFN-α) therapy, switching from NAs to IFN-α, and subsequent re-treatment with IFN-α were applied to treat ComTNsNt failure. We further evaluated the efficacy of this therapy.

METHODS

Eleven patients with ComTNsNt failure were enrolled in this study. Nine subjects (9/11) received IFN-α switching therapy. Combination therapy with IFN-α and ComTNsNt was administered in the first 4 weeks. Then, ComTNsNt was discontinued at the end of Week 4, and IFN-α monotherapy was continued for 6 months. Two (2/11) patients discontinued ComTNsNt without receiving IFN-α treatment. All 11 patients received the first re-treatment of IFN-α when they experienced hepatitis relapses after the withdrawal of IFN-α or ComTNsNt. Six (6/11) patients received a second re-treatment of IFN-α. Follow up was conducted after IFN-α therapy in all 11 patients.

RESULTS

Two patients (2/9) receiving IFN-α switching therapy experienced alanine aminotransferase (ALT) flare. In contrast, the two patients without IFN-α switching therapy experienced ALT flare. Multiple re-treatments with IFN-α resulted in a sustained response.

CONCLUSIONS

Interferon-α switching therapy and IFN-α re-treatment might be applied for treatment of ComTNsNt failure. IFN-α switching therapy resulted in safe ComTNsNt cessation, and IFN-α re-treatment induced a sustained response of IFN-α in all patients. This IFN-α treatment is an optional treatment for ComTNsNt failure.

Dynamics of Genotypic Mutations of the Hepatitis B Virus Associated With Long-Term Entecavir Treatment Determined With Ultradeep Pyrosequencing: A Retrospective Observational Study

The aim of the study is to explore the evolution of genotypic mutations within the reverse transcriptase region in partial virological responders (PVRs) receiving long-term entecavir (ETV) treatment. A total of 32 patients were classified as completely virological responders (CVRs) (n = 12) or PVRs (n = 20). Five partial responders were hepatitis B virus (HBV)-DNA positive after long-term therapy, which lasted for >3 years. A total of 71 serum samples from these 32 patients were assayed by ultra-deep pyrosequencing (UDPS): 32 samples were from all patients at baseline, and 39 were from PVRs with sequential inter-treatment. Approximately 84,708 sequences were generated per sample. At baseline, the quasispecies heterogeneity did not significantly differ between the 2 groups. The frequencies of substitutions indicating pre-existence of nucleos(t)ide analog resistant (NAr) mutants ranged from 0.10% to 6.70%, which did not statistically differ between groups either. However, the substitutions associated with the NAr mutants were significantly different from those associated with the non-NAr mutants in 13 patients; 6 of these patients were PVRs and the others were CVRs. Five patients were HBV DNA positive after regular ETV monotherapy for >3 years, and 4 of these patients underwent mild NAr substitution fluctuations (<20%). One patient developed virological breakthrough while bearing single, double, and triple (rtL180 M, rtM204 V, rtS202G) substitutions. In addition to the common substitutions, unknown amino acid substitutions, such as rtL145 M/S, rtF151Y/L, rtR153Q, rtI224 V, rtN248H, rtS223A, rtS256C, need to be further verified. NAr substitutions are observed at frequencies of 0.10% to 6.7% before therapy. Long-term ETV therapy generally results in virological responses, as long as the proportion of resistance mutations remains at a relatively low level. Genotypic resistance to ETV is detected in all PVRs receiving long-term ETV therapy.

Evidence of an Exponential Decay Pattern of the Hepatitis Delta Virus Evolution Rate and Fluctuations in Quasispecies Complexity in Long-Term Studies of Chronic Delta Infection

Chronic HDV infection can cause a severe form of viral hepatitis for which there is no specific treatment. Characterization of the hepatitis B or C viral quasispecies has provided insight into treatment failure and disease recurrence following liver transplantation, has proven useful to understand hepatitis B e antigen seroconversion, and has helped to predict whether hepatitis C infection will resolve or become chronic. It is likely that characterization of the hepatitis delta virus (HDV) quasispecies will ultimately have similar value for the management of this infection. This study sought to determine the RNA evolution rates in serum of chronic hepatitis delta (CHD) treatment-naïve patients, using next-generation sequencing methods. The region selected for study encompassed nucleotide positions 910 to 1270 of the genome and included the amber/W codon. Amber/W is a substrate of the editing process by the ADAR1 host enzyme and is essential for encoding the 2 delta antigens (HDAg). The amber codon encodes the small (unedited) HDAg form and the W codon the large (edited) HDAg form. The evolution rate was analyzed taking into account the time elapsed between samples, the percentage of unedited and edited genomes, and the complexity of the viral population. The longitudinal studies included 29 sequential samples from CHD patients followed up for a mean of 11.5 years. In total, 121,116 sequences were analyzed. The HDV evolution rate ranged from 9.5x10-3 to 1.2x10-3 substitutions/site/year and showed a negative correlation with the time elapsed between samples (p<0.05). An accumulation of transition-type changes was found to be responsible for higher evolution rates. The percentages of unedited and edited genomes and the quasispecies complexity showed no relationships with the evolution rate, but the fluctuations in the percentages of genomes and in complexity suggest continuous adaptation of HDV to the host conditions.

Composition and Interactions of Hepatitis B Virus Quasispecies Defined the Virological Response During Telbivudine Therapy

Reverse transcriptase (RT) mutations contribute to hepatitis B virus resistance during antiviral therapy with nucleos(t)ide analogs. However, the composition of the RT quasispecies and their interactions during antiviral treatment have not yet been thoroughly defined. In this report, 10 patients from each of 3 different virological response groups, i.e., complete virological response, partial virological response and virological breakthrough, were selected from a multicenter trial of Telbivudine treatment. Variations in the drug resistance-related critical RT regions in 107 serial serum samples from the 30 patients were examined by ultra-deep sequencing. A total of 496,577 sequence reads were obtained, with an average sequencing coverage of 4,641X per sample. The phylogenies of the quasispecies revealed the independent origins of two critical quasispecies, i.e., the rtA181T and rtM204I mutants. Data analyses and theoretical modeling showed a cooperative-competitive interplay among the quasispecies. In particular, rtM204I mutants compete against other quasispecies, which eventually leads to virological breakthrough. However, in the absence of rtM204I mutants, synergistic growth of the drug-resistant rtA181T mutants with the wild-type quasispecies could drive the composition of the viral population into a state of partial virological response. Furthermore, we demonstrated that the frequency of drug-resistant mutations in the early phase of treatment is important for predicting the virological response to antiviral therapy.

Increased intrahepatic quasispecies heterogeneity correlates with off-treatment sustained response to nucleos(t)ide analogues in e antigen-positive chronic hepatitis B patients

Finite treatment with nucleos(t)ide analogues (NAs) remains a great challenge for chronic hepatitis B in the clinic. This study aimed to investigate the relationship between intrahepatic quasispecies heterogeneity and the NAs off-treatment outcomes in a prospective cohort. Eighteen HBeAg-positive patients with chronic hepatitis B who achieved the cessation criteria underwent liver biopsy, and stopped treatment thereafter. Patients were followed up prospectively for 1 year. The reverse transcriptase (RT) gene of intrahepatic hepatitis B virus (HBV) was cloned and sequenced. Intrahepatic quasispecies heterogeneity and specific gene mutations were analysed using bioinformatic methods. Ten patients achieved sustained response, and eight patients developed viral relapse. The intrahepatic quasispecies Shannon entropy and nucleotide diversity within either RT or the surface (S) region of patients with sustained response were significantly higher (p < 0.05) than those of patients who had a viral relapse. Intrahepatic quasispecies Shannon entropy at the nucleotide level predicted the sustained off-treatment response (area under receiver operating characteristics curve 0.925; 95% CI 0.807-1.000; p 0.003). More positive selection sites and N-glycosylation mutations within the S region were found in patients with sustained response than in the patients with viral relapse (p < 0.01). Most of the positive selection sites in patients with sustained response were located in reported HLA-I-restricted or HLA-II-restricted epitopes. Intrahepatic quasispecies heterogeneity at the end of treatment was correlated with off-treatment outcomes in HBeAg-positive patients with chronic hepatitis B. More immune escape mutations were found within the S region in patients with sustained response. The higher intrahepatic quasispecies heterogeneity indicated a more robust immune control over HBV, which in turn maintained a sustained response after withdrawal of NAs.

Prediction of virological response by pretreatment hepatitis B virus reverse transcriptase quasispecies heterogeneity: the advantage of using next-generation sequencing

Prediction of antiviral efficacy prior to treatment remains largely unavailable. We have previously demonstrated the clinical value of on-treatment hepatitis B virus (HBV) reverse transcriptase (RT) quasispecies (QS) evolution patterns. In this study, we aimed to elucidate the relevance for prediction of pretreatment HBV RT QS characteristics by comparing the performance of next-generation sequencing (NGS) and clone-based Sanger sequencing (CBS). Thirty-six lamivudine-treated patients were retrospectively studied, including 18 responders and 18 non-responders. CBS and NGS data of pretreatment serum HBV were used to generate RT QS genetic complexity and diversity scores, according to our previous studies. The ability of both methods to predict responsiveness was evaluated with receiver operating characteristic (ROC) curves. A cut-off value was generated on the basis of prediction ability. Responders had significantly higher pretreatment RT QS genetic complexity and diversity (in the first two parts, which overlapped with the S gene, at both the nucleotide and amino acid levels) than non-responders by NGS-based testing. NGS-based algorithms predicted response better than CBS in the ROC curve analysis. The mean distance of the second contig had the highest area under the curve (AUC) value. When the cut-off value was set to 0.007186, the difference between survival curves was significant (p 0.0090). Pretreatment HBV RT QS heterogeneity in the overlapping region of the RT and S genes could be a predictor of antiviral efficacy. NGS improves the predictions of virological outcomes relative to CBS algorithms. This may have important implications for the clinical management of subjects chronically infected with HBV.

Prolonged entecavir therapy is not effective for HBeAg seroconversion in treatment-naive chronic hepatitis B patients with a partial virological response

The aims of this study were to investigate the efficacy of prolonged entecavir (ETV) therapy in treatment-naive chronic hepatitis B (CHB) patients and to determine whether continuous ETV therapy is feasible to achieve HBeAg seroconversion, particularly in patients with partial virological response (PVR). A total of 142 treatment-naive patients with CHB were enrolled. The mean duration of treatment was 65 (range, 26 to 90) months, and 86 patients (60.6%) were HBeAg positive. PVR was defined as detectable hepatitis B virus (HBV) DNA (>116 copies/ml) at year 1. The cumulative incidence of virological response (VR) increased from 54.9% at year 1 to 98.2% at year 7. HBeAg positivity (odds ratio [OR], 4.146; P = 0.001) and initial alanine aminotransferase (ALT) (OR, 0.997; P = 0.004) were independent risk factors for PVR. Among the 64 patients with PVR, 47 patients (73.4%) achieved VR within 4 years after prolonged ETV therapy without treatment adaptation. Three patients (2.1%) experienced virological breakthrough and HBV variants with genotypic resistance. The cumulative rate of HBeAg seroconversion was significantly higher in the patients with VR than in the patients with PVR (P = 0.018). None of the PVR patients with HBV DNA at ≥5,000 copies/ml at year 1 ever experienced HBeAg seroconversion. Multivariate analysis identified VR at year 1 as the only determinant of HBeAg seroconversion (hazard ratio [HR], 3.009; P = 0.010). In conclusion, although there were patients with PVR, prolonged ETV therapy showed excellent VR, with only 2.1% emergence of viral resistance during a 7-year follow-up. However, to achieve HBeAg seroconversion, drug modification is needed for HBeAg-positive patients with PVR (especially those with HBV DNA at ≥5,000 copies/ml at year 1).

Characterization of Full-Length Genomes of Hepatitis B Virus Quasispecies in Sera of Patients at Different Phases of Infection

Hepatitis B virus (HBV) infection results in different clinical presentation due to different levels of immune response. Our study aimed to characterize HBV full-length genome quasispecies (QS) in patients with different phases of infection to better understand its pathogenesis. Forty treatment-naive HBV-infected patients were enrolled, including 10 cases of acute hepatitis B (AHB), 9 cases of immunotolerant (IT) HBV carriers, 11 cases of chronic hepatitis B (CHB), and 10 cases of acute-on-chronic liver failure (ACLF). The present study was conducted by clone-based sequencing. QS heterogeneity within each open reading frame was calculated. The mutation frequency index (MFI) and amino acid variations within the large HBsAg, HBcAg, and HBxAg regions were analyzed based on the different infection phases. In total, 606 HBV full-length sequences were obtained. HBV QS had higher heterogeneity in ACLF and CHB than that in IT among chronically infected individuals. AHB patients had the lower QS heterogeneity at onset than those with chronic infection. ACLF patients had the highest frequency of mutations in the core promoter and precore region. A triple mutation (A1762T/G1764A/G1896A) was observed more frequently in genotype C than in genotype B. The MFI indicated that specific peptides of the studied regions had more frequent mutations in ACLF. Furthermore, several amino acid variations, known as T- and B-cell epitopes, were potentially associated with the immunoactive phase of infection. More HBV genome mutations and deletions were observed in patients with more severe diseases, particularly in specific regions of the core and preS regions, the clinical significance and mechanism of which need to be further investigated.

Hepatitis B virus (HBV) variants fluctuate in paired plasma and peripheral blood mononuclear cells among patient cohorts during different chronic hepatitis B (CHB) disease phases

Hepatitis B virus is classically considered a hepatotropic virus but also infects peripheral blood mononuclear cells. Chronic hepatitis B has different disease phases modulated by host immunity. We compared HBV variability, drug resistance and immune escape mutations in the overlapping HBV polymerase/surface gene in plasma and peripheral blood mononuclear cells in different disease phases. Plasma and peripheral blood mononuclear cells were isolated from 22 treatment naïve patient cohorts (five inactive, six immune-active, nine HBeAg negative and two immune-tolerant). HBV was genotyped via line probe assay, hepatitis B surface antigen titres were determined by an in-house immunoassay, and HBV DNA was quantified by kinetic PCR. The HBV polymerase/surface region, including full genome in some, was PCR-amplified and cloned, and ~20 clones/sample were sequenced. The sequences were subjected to various mutational and phylogenetic analyses. Clonal sequencing showed that only three of 22 patients had identical HBV genotype profiles in both sites. In immune-active chronic hepatitis B, viral diversity in plasma was higher compared with peripheral blood mononuclear cells. Mutations at residues, in a minority of clones, associated with drug resistance, and/or immune escape were found in both compartments but were more common in plasma. Immune escape mutations were more often observed in the peripheral blood mononuclear cells of immune-active CHB carriers, compared with other disease phases. During all CHB disease phases, differences exist between HBV variants found in peripheral blood mononuclear cells and plasma. Moreover, these data indicate that HBV evolution occurs in a compartment and disease phase-specific fashion.

Application of a hybrid method combining grey model and back propagation artificial neural networks to forecast hepatitis B in china

Accurate incidence forecasting of infectious disease provides potentially valuable insights in its own right. It is critical for early prevention and may contribute to health services management and syndrome surveillance. This study investigates the use of a hybrid algorithm combining grey model (GM) and back propagation artificial neural networks (BP-ANN) to forecast hepatitis B in China based on the yearly numbers of hepatitis B and to evaluate the method's feasibility. The results showed that the proposal method has advantages over GM (1, 1) and GM (2, 1) in all the evaluation indexes.

Effect of hepatitis B virus reverse transcriptase variations on entecavir treatment response

BACKGROUND

Entecavir therapy often reduces hepatitis B virus (HBV) DNA to an undetectable level, but HBV DNA remain detectable in some patients. We investigated whether baseline HBV reverse transcriptase (rt) polymorphism and quasispecies complexity and diversity were associated with treatment response.

METHODS

Pretreatment HBV DNA levels, HBV rt sequence, serology, and quasispecies complexity and diversity from 305 entecavir-treated patients were determined. These data were tested for their association with year 1 virological outcome, defined by optimal response (undetectable HBV DNA; lower limit of detection, ≤12 IU/mL) or partial response (detectable HBV DNA).

RESULTS

Four rt variants were more frequently detected in the 64 partial responders than in the 241 optimal responders (all P < .05). Multivariate analysis revealed that high baseline HBV DNA level (P < .0001; odds ratio [OR], 2.32), HBV e antigen (HBeAg) positivity (P < .001; OR, 3.70), and rt124N (P = .002; OR, 3.06) were associated with a partial entecavir response. Compared with the optimal responders, the partial responders had a lower quasispecies complexity and diversity.

CONCLUSIONS

Apart from the known factors (high baseline HBV DNA level and HBeAg positivity), a novel single nucleotide polymorphism (rt124N) and lower quasispecies complexity and diversity were associated with partial entecavir response at year 1.

Efficacy and safety of tenofovir in a kidney transplant patient with chronic hepatitis B and nucleos(t)ide multidrug resistance: a case report

INTRODUCTION

Five nucleos(t)ide analogs are used to treat chronic hepatitis B. Ideal nucleos(t)ide analog therapy in chronic hepatitis B patients with kidney transplantation must ensure virological suppression and minimize renal injury. However, resistance to nucleos(t)ide analogs frequently results in virological breakthrough, hepatitis flare, and complicated deterioration of the transplanted kidney. Inappropriate rescue therapy for drug resistance may subsequently cause hepatitis B virus multidrug resistance. Currently, tenofovir is used to treat chronic hepatitis B patients with kidney transplantation. In the field, we first reported combination therapy with tenofovir plus entecavir in a kidney transplant chronic hepatitis B patient with nucleos(t)ide analog multidrug resistance.

CASE PRESENTATION

A 50-year-old Chinese man with chronic hepatitis B and kidney transplantation received nucleos(t)ide analog therapy with sequential monotherapy and combination therapy. Virological parameters, hepatic enzymology and renal function were monitored. Drug-resistance mutations were detected by sequence analysis. Our patient received sequential nucleos(t)ide analog monotherapy and inappropriate combination therapy during 132 months, which caused multidrug resistance and renal functional injury. Entecavir plus adefovir was administered in month 106, resulting in decreased hepatitis B virus load, normal hepatic function, and stabilized creatinine clearance. As a result of rebounded viral load and significantly declining creatinine clearance, tenofovir plus entecavir was administered in month 133. After eight weeks, undetectable hepatitis B virus DNA, normal hepatic function and improved creatinine clearance were present. Compared with combination therapy with adefovir plus entecavir, tenofovir plus entecavir showed a potent antiviral effect for multidrug resistance and minimized renal injury.

CONCLUSIONS

In chronic hepatitis B patients with kidney transplantation, sequential monotherapy with antiviral agents with low barriers to resistance should be avoided, and initial therapy with entecavir is a better option. Combination therapy with tenofovir plus entecavir in this setting with multidrug resistance is safe and effective.

An improved method for simple and efficient hepatitis B virus genome cloning

The genetic variation of hepatitis B virus (HBV) is associated with a natural history of infection and with a response to antiviral therapy. Current attempts to investigate the genetic variation within patients infected with HBV have characterized only the subgenomic region or a limited number of full-length genomes. In this study, a simple and efficient full-length HBV cloning method, which is independent of restriction digestion and ligation, from a clinical sample is presented. The full-length HBV amplified from patients' sera serves as two megaprimers and contains two ends that are homologous to the insertion site of a newly constructed plasmid. This method could enable the extension of the plasmid backbone by DNA polymerase. To improve the cloning efficiency, a LacZ gene fragment was incorporated into the middle of the insertion site in the plasmid and was used for blue-white screening. Through optimization and evaluation, a high success rate of positive clones (90%) from an individual patient was obtained. Thus, this study provides a simple, efficient, and convenient method for full-length HBV cloning from many clinical samples, and this method could greatly facilitate studies concerning the genetic variation of HBV.

Inference with viral quasispecies diversity indices: clonal and NGS approaches

Given the inherent dynamics of a viral quasispecies, we are often interested in the comparison of diversity indices of sequential samples of a patient, or in the comparison of diversity indices of virus in groups of patients in a treated versus control design. It is then important to make sure that the diversity measures from each sample may be compared with no bias and within a consistent statistical framework. In the present report, we review some indices often used as measures for viral quasispecies complexity and provide means for statistical inference, applying procedures taken from the ecology field. In particular, we examine the Shannon entropy and the mutation frequency, and we discuss the appropriateness of different normalization methods of the Shannon entropy found in the literature. By taking amplicons ultra-deep pyrosequencing (UDPS) raw data as a surrogate of a real hepatitis C virus viral population, we study through in-silico sampling the statistical properties of these indices under two methods of viral quasispecies sampling, classical cloning followed by Sanger sequencing (CCSS) and next-generation sequencing (NGS) such as UDPS. We propose solutions specific to each of the two sampling methods-CCSS and NGS-to guarantee statistically conforming conclusions as free of bias as possible.

CONTACT

josep.gregori@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.