Increased viral quasispecies evolution in HBeAg seroconverter patients treated with oral nucleoside therapy

Nonproductive Hepatitis B Virus Covalently Closed Circular DNA Generates HBx-Related Transcripts from the HBx/Enhancer I Region and Acquires Reactivation by Superinfection in Single Cells

Hepatitis B virus (HBV) infection remains a public health problem worldwide. Persistent HBV infection relies on active transcription of the covalently closed circular DNA (cccDNA) in hepatocytes, which is less understood at the single-cell level. In this study, we isolated primary human hepatocytes from liver-humanized FRG mice infected with HBV and examined cccDNA transcripts in single cells based on 5' end sequencing. Our 5' transcriptome sequencing (RNA-seq) analysis unambiguously assigns different viral transcripts with overlapping 3' sequences and quantitatively measures viral transcripts for structural genes (3.5 kb, 2.4 kb, and 2.1 kb) and the nonstructural X gene (0.7 kb and related) in single cells. We found that an infected cell either can generate all viral transcripts, signifying active transcription, or presents only transcripts from the X gene and its associated enhancer I domain and no structural gene transcripts. Results from cell infection assays with recombinant HBV show that nonproductive transcription of cccDNA can be activated by incoming virus through superinfection. Moreover, upon HBV infection, cccDNA apparently can be transcribed in the absence of HBx and produces HBx, needed for productive transcription of other viral genes. These results shed new light on cccDNA transcription at the single-cell level and provide insights useful for improving the treatment strategy against chronic HBV infection. IMPORTANCE Hepatitis B virus (HBV) infection can be effectively suppressed but rarely cured by available drugs. Chronic HBV infection is based on persistence of covalently closed circular DNA (cccDNA) and continuous infection and reinfection with HBV in the liver. Understanding transcriptional regulation of cccDNA will help to achieve permanent transcriptional silencing, i.e., functional cure of HBV. In our study, we found that an infected cell either can generate all viral transcripts, signifying active transcription, or presents only transcripts from the X gene and its associated enhancer I domain and no structural gene transcripts. The nonproductive transcription of cccDNA can be activated by incoming virus through superinfection. Upon an infection, cccDNA apparently can be transcribed in the absence of HBx to produce HBx, necessary for subsequent transcription of other HBV genes. Our studies shed new light on the mechanism of HBV infection and may have implications for a functional cure regimen for HBV.

Combined viral quasispecies diversity and hepatitis B core-related antigen predict off-nucleos(t)ide analog durability in HBeAg-negative patients

BACKGROUND

Viral quasispecies dynamics between pre- and post-nucleos(t)ide analog (NA) therapy remains unclear.

AIM

This study aimed to investigate the HBV quasispecies evolution and its relationship with durability of off-therapy responses in HBeAg-negative chronic hepatitis B (CHB) patients who stopped NA therapy.

METHODS

Fifty-four HBeAg-negative CHB patients who stopped NAs, including 19 virological controllers (VC) who maintained serum HBV DNA < 2000 IU/mL beyond 1-year off-therapy, and 35 virological relapsers (VR) experiencing virological relapse within 1-year off-therapy were recruited. Viral quasispecies was analyzed by deep sequencing. Hepatitis B core-related antigen (HBcrAg) and HBsAg were also measured.

RESULTS

VC had significantly higher baseline viral quasispecies diversity of the precore/core gene, measured by nucleotide diversity, than VR. Low baseline viral nucleotide diversity (< 0.01) and high HBcrAg (≧ 2.0 KU/mL), but not HBsAg, at end of treatment (EOT) were significantly associated with higher risk of 1-year virological relapse (hazard ratio [HR] 6.09 and 3.31, respectively). Combination of low baseline viral nucleotide diversity and high HBcrAg at EOT could identify patients at high risk (HR 15.82). Further analysis of the evolution of HBV whole genome showed that HBV nucleotide diversity negatively correlated with serum HBV DNA levels. Notably, the viral quasispecies diversity between pre- and post-NA treatment remained relatively unchanged.

CONCLUSION

Higher baseline HBV quasispecies diversity associates with more durable off-therapy viral suppression in HBeAg-negative CHB patients. Combination of baseline viral nucleotide diversity and HBcrAg at EOT can identify patients at high risk for virological relapse after stopping NAs.

Impact of HBeAg status on immune recovery in HIV/HBV co-infected patients after highly active antiretroviral therapy

A retrospective analysis was conducted on HIV-infected patients whose continuously HAART strategy was lamivudine +tenofovir+ efavirenz. Propensity matching for 35 HBeAg-positive/HIV co-infected patients, 35 HBeAg-negative/HIV co-infected patients, and 70 HIV mono-infected patients. Immune recovery (including CD4 cells count, CD4/CD8 ratio, CD4 count multiples and CD4/CD8 multiples) of HBeAg-negative/HIV co-infected group are continuously lower than HBeAg-positive/HIV co-infected group and HIV mono-infected group. The result indicated that the mechanisms associated with HBeAg-negative may be involved in the regulation of immune recovery after HAART.

Case Report: Application of hepatitis B virus (HBV) deep sequencing to distinguish between acute and chronic infection

Deep sequencing of the full-length hepatitis B virus (HBV) genome provides the opportunity to determine the extent to which viral diversity, genotype, polymorphisms, insertions and deletions may influence presentation and outcomes of disease. Increasing experience with analysis of HBV genomic data opens up the potential for using these data to inform insights into pathophysiology of infection and to underpin decision making in clinical practice. We here set out to undertake whole genome HBV sequencing from an adult who presented acutely unwell with a new diagnosis of HBV infection, and tested positive for both HBV anti-core IgM and IgG, possibly representing either acute hepatitis B infection (AHB) or chronic hepatitis B with an acute reactivation (CHB-AR). The distinction between these two scenarios may be important in predicting prognosis and underpinning treatment decisions, but can be challenging based on routine laboratory tests. Through application of deep whole-genome sequencing we typed the isolate as genotype-D1, and identified several minority variants including G1764A and G1986A substitutions in the pre-core promoter and pre-core regions, which support CHB-AR rather than AHB. In the longer term, enhanced deep sequencing data for HBV may provide improved evidence to distinguish between acute and chronic infection, to predict outcomes and to stratify treatment.

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.

Case Report: Application of hepatitis B virus (HBV) deep sequencing to distinguish between acute and chronic infection

Deep sequencing of the full-length hepatitis B virus (HBV) genome provides the opportunity to determine the extent to which viral diversity, genotype, polymorphisms, insertions and deletions may influence presentation and outcomes of disease. Increasing experience with analysis of HBV genomic data opens up the potential for using these data to inform insights into pathophysiology of infection and to underpin decision making in clinical practice. We here set out to undertake whole genome HBV sequencing from an adult who presented acutely unwell with a new diagnosis of HBV infection, and tested positive for both HBV anti-core IgM and IgG, possibly representing either acute hepatitis B infection (AHB) or chronic hepatitis B with an acute reactivation (CHB-AR). The distinction between these two scenarios may be important in predicting prognosis and underpinning treatment decisions, but can be challenging based on routine laboratory tests. Through application of deep whole-genome sequencing we typed the isolate as genotype-D1, and identified several minority variants including G1764A and G1986A substitutions in the pre-core promoter and pre-core regions, which support CHB-AR rather than AHB. In the longer term, enhanced deep sequencing data for HBV may provide improved evidence to distinguish between acute and chronic infection, to predict outcomes and to stratify treatment.

The evolution and clinical impact of hepatitis B virus genome diversity

The global burden of hepatitis B virus (HBV) is enormous, with 257 million persons chronically infected, resulting in more than 880,000 deaths per year worldwide. HBV exists as nine different genotypes, which differ in disease progression, natural history and response to therapy. HBV is an ancient virus, with the latest reports greatly expanding the host range of the Hepadnaviridae (to include fish and reptiles) and casting new light on the origins and evolution of this viral family. Although there is an effective preventive vaccine, there is no cure for chronic hepatitis B, largely owing to the persistence of a viral minichromosome that is not targeted by current therapies. HBV persistence is also facilitated through aberrant host immune responses, possibly due to the diverse intra-host viral populations that can respond to host-mounted and therapeutic selection pressures. This Review summarizes current knowledge on the influence of HBV diversity on disease progression and treatment response and the potential effect on new HBV therapies in the pipeline. The mechanisms by which HBV diversity can occur both within the individual host and at a population level are also discussed.

Illumina and Nanopore methods for whole genome sequencing of hepatitis B virus (HBV)

Advancing interventions to tackle the huge global burden of hepatitis B virus (HBV) infection depends on improved insights into virus epidemiology, transmission, within-host diversity, drug resistance and pathogenesis, all of which can be advanced through the large-scale generation of full-length virus genome data. Here we describe advances to a protocol that exploits the circular HBV genome structure, using isothermal rolling-circle amplification to enrich HBV DNA, generating concatemeric amplicons containing multiple successive copies of the same genome. We show that this product is suitable for Nanopore sequencing as single reads, as well as for generating short-read Illumina sequences. Nanopore reads can be used to implement a straightforward method for error correction that reduces the per-read error rate, by comparing multiple genome copies combined into a single concatemer and by analysing reads generated from plus and minus strands. With this approach, we can achieve an improved consensus sequencing accuracy of 99.7% and resolve intra-sample sequence variants to form whole-genome haplotypes. Thus while Illumina sequencing may still be the most accurate way to capture within-sample diversity, Nanopore data can contribute to an understanding of linkage between polymorphisms within individual virions. The combination of isothermal amplification and Nanopore sequencing also offers appealing potential to develop point-of-care tests for HBV, and for other viruses.

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.

Long-Term Follow-Up of Acute Hepatitis B: New Insights in Its Natural History and Implications for Antiviral Treatment

Acute hepatitis B infection (AHB) is still a common viral acute hepatitis worldwide. As vaccination, antiviral treatment, and immigration are bound to affect the epidemiological landscape of HBV infections, and some of its aspects need to be investigated: (1) the circulation of vaccine escape mutants and of primary drug resistant strains; (2) the change in HBV genotype prevalence; and (3) the clinical implications of AHB and the probability of chronification. The serological, virological, and clinical parameters of 75 patients, acutely infected by HBV, were gathered for a retrospective study. Long-term follow up, either to complete seroconversion or for up to five years, was possible for 44 patients. Sequence analysis of the reverse transcriptase/HBsAg and precore regions was performed to investigate the molecular epidemiology and pathogenesis of recent infections by HBV. Genotype distribution in AHB in Italian patients was radically different from that of chronic infections, with a dramatic increase of extra-European genotypes (A1, F), suggesting that a proportion of AHBs are currently related to imported strains. None of the documented infections occurred in vaccinated individuals, while HBsAg variants (potentially vaccine escape variants) were rare and less prevalent than in chronic infections. No drug resistant strains were observed. Spontaneous viral clearance occurred in all but three cases. Time to viral clearance was inversely proportional to liver damage, but HBsAg titer on day 28 and, better still, HBsAg decay from day 0 to day 28 after admission, were the best predictors of chronification. They are, thus, potentially useful to guide antiviral treatment to prevent chronic evolution.

Single-virion sequencing of lamivudine-treated HBV populations reveal population evolution dynamics and demographic history

BACKGROUND

Viral populations are complex, dynamic, and fast evolving. The evolution of groups of closely related viruses in a competitive environment is termed quasispecies. To fully understand the role that quasispecies play in viral evolution, characterizing the trajectories of viral genotypes in an evolving population is the key. In particular, long-range haplotype information for thousands of individual viruses is critical; yet generating this information is non-trivial. Popular deep sequencing methods generate relatively short reads that do not preserve linkage information, while third generation sequencing methods have higher error rates that make detection of low frequency mutations a bioinformatics challenge. Here we applied BAsE-Seq, an Illumina-based single-virion sequencing technology, to eight samples from four chronic hepatitis B (CHB) patients - once before antiviral treatment and once after viral rebound due to resistance.

RESULTS

With single-virion sequencing, we obtained 248-8796 single-virion sequences per sample, which allowed us to find evidence for both hard and soft selective sweeps. We were able to reconstruct population demographic history that was independently verified by clinically collected data. We further verified four of the samples independently through PacBio SMRT and Illumina Pooled deep sequencing.

CONCLUSIONS

Overall, we showed that single-virion sequencing yields insight into viral evolution and population dynamics in an efficient and high throughput manner. We believe that single-virion sequencing is widely applicable to the study of viral evolution in the context of drug resistance and host adaptation, allows differentiation between soft or hard selective sweeps, and may be useful in the reconstruction of intra-host viral population demographic history.

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.

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.

A model with combined viral and metabolic factors effectively predicts HBeAg status under long term entecavir therapy: a prospective cohort study

Background & Aim

The aim was to extract factors from virologic and biochemical profiles at baseline and 24 weeks of treatment to predict HBeAg seroconversion in patients treated with ETV.

Methods

HBeAg positive chronic hepatitis B patients receiving ETV naïve-treatment were enrolled. HBV DNA, ALT, and serological markers were prospectively monitored every 6 months for 240 weeks. The cumulative rates of virologic response (VR), biochemical response (BR), and HBeAg seroconversion were determined, and potential predictors for HBeAg seroconversion were identified through uni/multivariate analysis.

Result

Two hundred twenty nine patients were eligible for this study. The cumulative rates of VR, BR, and HBeAg seroconversion at 240 weeks were 88.4 %, 100 %, and 36.7 %, respectively. Multivariate analysis showed that HBV DNA (OR, 2.8, p = 0.003), ALT (OR, 2.6, p = 0.005) at baseline, undetectable HBV DNA within 24 weeks (OR = 3.2, p < 0.001), and body mass index (BMI) ≥24kg/m² (OR = 0.038, p = 0.013) were associated with HBeAg seroconversion. A prediction model for probability of HBeAg seroconversion was constructed. Patients can be classified into high (>40 %), intermediate (20–40 %), or low (≤20 %) groups based on the calculated probability of HBeAg seroconversion. The cumulative rates of HBeAg seroconversion were different among the three groups (p < 0.001). About 58 % patients in the high probability group achieved HBeAg seroconversion while almost 90 % patients within the low group remained HBeAg positive.

Conclusion

A combination of HBV DNA, ALT and BMI values at baseline, and undetectable HBV DNA level within 24 weeks can predict HBeAg seroconversion. Both viral and metabolic factors likely determine HBeAg status with ETV treatment.

Trial registration

CTR20132358

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.

Baseline interpatient hepatitis B viral diversity differentiates HBsAg outcomes in patients treated with tenofovir disoproxil fumarate

BACKGROUND AND AIMS

HBsAg loss is a desired, but rare, treatment-induced clinical endpoint in chronic hepatitis B (CHB). Few studies have evaluated viral factors contributing to HBsAg loss.

METHODS

This study evaluated baseline interpatient sequence diversity across the HBV genome in tenofovir disoproxil fumarate-treated patients who lost HBsAg and compared it to that of control patients with high HBsAg levels throughout therapy. Twenty-one HBeAg+ patients (14 genotype (GT) A and 7 GT D) who achieved HBsAg loss and 27 controls (17 GT A and 10 GT D), were analyzed. Population sequencing was performed on baseline samples and pairwise genetic distances were calculated for 17 overlapping regions across the HBV genome as a measure of interpatient viral diversity.

RESULTS

Overall, viral diversity was up to 10-fold higher across GT D patients compared to GT A patients throughout the HBV genome. Within the pol/RT and HBs genes, interpatient viral diversity was significantly lower among HBsAg loss patients for both GT A and D, with the difference driven largely by a reduction in diversity in the small S gene. Conversely, interpatient viral diversity was generally higher in HBsAg loss patients across the HBx gene regulatory elements and precore region.

CONCLUSION

In HBsAg loss patients, less interpatient viral diversity was observed within structural-coding regions while specific regions across the HBx and precore genes encoding nonstructural regulatory elements generally displayed higher interpatient viral diversity. These distinct patterns may reflect different responses to adaptive pressure for HBV genomic structural and nonstructural elements.

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.

Clinical application of estimating hepatitis B virus quasispecies complexity by massive sequencing: correlation between natural evolution and on-treatment evolution

AIM

To evaluate HBV quasispecies (QA) complexity in the preCore/Core regions in relation to HBeAg status, and explore QA changes under natural evolution and nucleoside analogue (NUC) treatment.

METHODS

Ultra-deep pyrosequencing of HBV preCore/Core regions in 30 sequential samples (baseline [diagnosis], treatment-free, and treatment-nonresponse) from 10 retrospectively selected patients grouped according to HBeAg status over time: HBeAg+ (N = 4), HBeAg- (N = 2), and fluctuating HBeAg (transient seroreversion/seroconversion pattern) (N = 4). QA complexity was defined by Shannon entropy, mutation frequency, nucleotide diversity, and mutation frequency of amino acids (MfAA) in preCore and Core.

RESULTS

The QA was less complex in HBeAg+ than in HBeAg- or fluctuating HBeAg. High complexity in preCore was associated with decreased viral replication (preCore MfAA negatively correlated with HBV-DNA, p = 0.005). QA complexity in the treatment-free period negatively correlated with values seen during treatment. Specific variants were mainly selected in the Core region in HBeAg- and fluctuating HBeAg patients, suggesting higher immune pressure than in HBeAg+.

CONCLUSIONS

The negative correlation between QA natural evolution and on-treatment evolution indicates the importance of pre-treatment QA study to predict QA changes in NUC nonresponders. Study of QA complexity could be useful for managing HBV infection.

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.

Clinical implications of evolutionary patterns of homologous, full-length hepatitis B virus quasispecies in different hosts after perinatal infection

Chronic hepatitis B virus (HBV) infection via perinatal transmission is common in the Asia-Pacific region, but related quasispecies (QS) characteristics are not yet defined. To investigate the homologous, full-length HBV QS after perinatal infection and their clinical implications, five pairs of mother-daughter patients with chronic HBV infection (one patient with liver cirrhosis, one with immune tolerance, and eight with chronic hepatitis) were included. Full-length HBV were amplified by PCR from serum samples before antiviral treatment and cloned; an average of 17 clones per sample were sequenced, and the QS complexities, diversities, and evolution patterns were analyzed. Full-length HBV sequence similarities within mother-daughter pairs were 91.3 to 98.3%. The QS complexities of full-length HBV were similar between mothers and daughters (median of 0.9734 compared to 0.9688, P>0.05), as were the diversities (median of 3.396×10(-3) compared to 4.617×10(-3) substitutions/site, P>0.05). However, the distribution patterns of QS complexities and diversities within specific genes were different, and QS genetic distances of the mothers were higher than those of daughters, both more evident in pairs with different antiviral responses and different immune phases or stages. The nucleotide substitution rate of full-length HBV was 14.388×10(-5) substitutions/site/year, whereas the preC/C gene rate was the highest. Mutations and indels were more common in clones from the mothers, which decreased the affinity of epitopes by 6- to 89-fold. The various genes from homologous HBV genomes evolved in different patterns despite numerically comparable full-length QS complexities and diversities. The different patterns may correlate with the immune stages of chronic HBV infection, which warrants further study.

Quasispecies structure, cornerstone of hepatitis B virus infection: mass sequencing approach

Hepatitis B virus (HBV) is a DNA virus with complex replication, and high replication and mutation rates, leading to a heterogeneous viral population. The population is comprised of genomes that are closely related, but not identical; hence, HBV is considered a viral quasispecies. Quasispecies variability may be somewhat limited by the high degree of overlapping between the HBV coding regions, which is especially important in the P and S gene overlapping regions, but is less significant in the X and preCore/Core genes. Despite this restriction, several clinically and pathologically relevant variants have been characterized along the viral genome. Next-generation sequencing (NGS) approaches enable high-throughput analysis of thousands of clonally amplified regions and are powerful tools for characterizing genetic diversity in viral strains. In the present review, we update the information regarding HBV variability and present a summary of the various NGS approaches available for research in this virus. In addition, we provide an analysis of the clinical implications of HBV variants and their study by NGS.