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

Viral Haplotypes in COVID-19 Patients Associated With Prolonged Viral Shedding

Background

Recently, more patients who recovered from the novel coronavirus disease 2019 (COVID-19) may later test positive for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) again using reverse transcription-polymerase chain reaction (RT-PCR) testing. Even though it is still controversial about the possible explanation for clinical cases of long-term viral shedding, it remains unclear whether the persistent viral shedding means re-infection or recurrence.

Methods

Specimens were collected from three COVID-19-confirmed patients, and whole-genome sequencing was performed on these clinical specimens during their first hospital admission with a high viral load of SARS-CoV-2. Laboratory tests were examined and analyzed throughout the whole course of the disease. Phylogenetic analysis was carried out for SARS-CoV-2 haplotypes.

Results

We found haplotypes of SARS-CoV-2 co-infection in two COVID-19 patients (YW01 and YW03) with a long period of hospitalization. However, only one haplotype was observed in the other patient with chronic lymphocytic leukemia (YW02), which was verified as one kind of viral haplotype. Patients YW01 and YW02 were admitted to the hospital after being infected with COVID-19 as members of a family cluster, but they had different haplotype characteristics in the early stage of infection; YW01 and YW03 were from different infection sources; however, similar haplotypes were found together.

Conclusion

These findings show that haplotype diversity of SARS-CoV-2 may result in viral adaptation for persistent shedding in multiple recurrences of COVID-19 patients, who met the discharge requirement. However, the correlation between haplotype diversity of SARS-CoV-2 virus and immune status is not absolute. It showed important implications for the clinical management strategies for COVID-19 patients with long-term hospitalization or cases of recurrence.

Next-generation sequencing for the diagnosis of hepatitis B: current status and future prospects

INTRODUCTION

Hepatitis B virus (HBV) causes a complex and persistent infection with a major impact on patients health. Viral-genome sequencing can provide valuable information for characterizing virus genotype, infection dynamics and drug and vaccine resistance.

AREAS COVERED

This article reviews the current literature to describe the next-generation sequencing progress that facilitated a more comprehensive study of HBV quasispecies in diagnosis and clinical monitoring.

EXPERT OPINION

HBV variability plays a key role in liver disease progression and treatment efficacy. Second-generation sequencing improved the sensitivity for detecting and quantifying mutations, mixed genotypes and viral recombination. Third-generation sequencing enables the analysis of the entire HBV genome, although the high error rate limits its use in clinical practice.

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 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.

Next generation sequencing identifies baseline viral mutants associated with treatment response to pegylated interferon in HBeAg-positive chronic hepatitis B

Current data of hepatitis B virus (HBV) variants associated with treatment outcome identified by next generation sequencing (NGS) are limited. This study was aimed at determining the role of baseline sequence variations in the enhancer II (EnhII), basal core promotor (BCP) and pre-core (PC) regions of HBV genotype C in patients treated with pegylated interferon (PEG-IFN). Patients with HBeAg-positive chronic hepatitis B (CHB) treated with 48-week PEG-IFN were enrolled. Combined response (CR) at week 96 was defined by HBeAg seroconversion plus HBV DNA < 2000 IU/mL and HBsAg < 1000 IU/mL. Pre-treatment viral mutations were characterized by Sanger sequencing and NGS (Miseq Illumina platform). Among 47 patients (32 male, mean age 32.4 years), CR was achieved in 12 (25.5%) individuals. Overall, NGS was superior to Sanger sequencing in detecting mutations (61.7% vs. 38.3%, P < 0.001). Based on NGS, the prevalence of T1753V (T1753C/A/G) and A1762T/G1764A variants were significantly lower in responders compared to non-responders (8.3% vs. 51.4%, P = 0.009 and 33.3% vs. 68.6%, P = 0.032, respectively). No significant difference between groups was found regarding C1653T and G1896A mutants. The absence of T1753V and A1762T/G1764A mutations were factors associated with CR (OR 11.65, 95%CI 1.36-100.16, P = 0.025, and OR 4.36, 95%CI 1.08-17.63, P = 0.039, respectively). The existence of pre-treatment T1753V, A1762T/G1764A mutations and their combination yielded negative predictive values of 94.7%, 85.7% and 93.8%, respectively. The presence of HBV mutants in the BCP region determined by NGS at baseline was associated with poor treatment outcome in patients with HBeAg-positive CHB receiving PEG-IFN.

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.

Insights From Deep Sequencing of the HBV Genome-Unique, Tiny, and Misunderstood

Hepatitis B virus (HBV) is a unique, tiny, partially double-stranded, reverse-transcribing DNA virus with proteins encoded by multiple overlapping reading frames. The substitution rate is surprisingly high for a DNA virus, but lower than that of other reverse transcribing organisms. More than 260 million people worldwide have chronic HBV infection, which causes 0.8 million deaths a year. Because of the high burden of disease, international health agencies have set the goal of eliminating HBV infection by 2030. Nonetheless, the intriguing HBV genome has not been well characterized. We summarize data on the HBV genome structure and replication cycle, explain and quantify diversity within and among infected individuals, and discuss advances that can be offered by application of next-generation sequencing technology. In-depth HBV genome analyses could increase our understanding of disease pathogenesis and allow us to better predict patient outcomes, optimize treatment, and develop new therapeutics.

Deep sequencing of HBV pre-S region reveals high heterogeneity of HBV genotypes and associations of word pattern frequencies with HCC

Hepatitis B virus (HBV) infection is a common problem in the world, especially in China. More than 60–80% of hepatocellular carcinoma (HCC) cases can be attributed to HBV infection in high HBV prevalent regions. Although traditional Sanger sequencing has been extensively used to investigate HBV sequences, NGS is becoming more commonly used. Further, it is unknown whether word pattern frequencies of HBV reads by Next Generation Sequencing (NGS) can be used to investigate HBV genotypes and predict HCC status. In this study, we used NGS to sequence the pre-S region of the HBV sequence of 94 HCC patients and 45 chronic HBV (CHB) infected individuals. Word pattern frequencies among the sequence data of all individuals were calculated and compared using the Manhattan distance. The individuals were grouped using principal coordinate analysis (PCoA) and hierarchical clustering. Word pattern frequencies were also used to build prediction models for HCC status using both K-nearest neighbors (KNN) and support vector machine (SVM). We showed the extremely high power of analyzing HBV sequences using word patterns. Our key findings include that the first principal coordinate of the PCoA analysis was highly associated with the fraction of genotype B (or C) sequences and the second principal coordinate was significantly associated with the probability of having HCC. Hierarchical clustering first groups the individuals according to their major genotypes followed by their HCC status. Using cross-validation, high area under the receiver operational characteristic curve (AUC) of around 0.88 for KNN and 0.92 for SVM were obtained. In the independent data set of 46 HCC patients and 31 CHB individuals, a good AUC score of 0.77 was obtained using SVM. It was further shown that 3000 reads for each individual can yield stable prediction results for SVM. Thus, another key finding is that word patterns can be used to predict HCC status with high accuracy. Therefore, our study shows clearly that word pattern frequencies of HBV sequences contain much information about the composition of different HBV genotypes and the HCC status of an individual.

HBV infection can lead to many liver complications including hepatocellular carcinoma (HCC), one of the most common liver cancers in China. High-throughput sequencing technologies have recently been used to study the genotype sequence compositions of HBV infected individuals and to distinguish chronic HBV (CHB) infection from HCC. We used NGS to sequence the pre-S region of a large number of CHB and HCC individuals and designed novel word pattern based approaches to analyze the data. We have several surprising key findings. First, most HBV infected individuals contained mixtures of genotypes B and C sequences. Second, multi-dimensional scaling (MDS) analysis of the data showed that the first principal coordinate was closely associated with the fraction of genotype B (or C) sequences and the second principal coordinate was highly associated with the probability of HCC. Third, we also designed K-nearest neighbors (KNN) and support vector machine (SVM) based classifiers for CHB and HCC with high prediction accuracy. The results were validated in an independent data set.

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.

Hepatitis B virus core antigen mutations predict post-operative prognosis of patients with primary hepatocellular carcinoma

The aim of this study was to explore the relationship between hepatitis B virus (HBV) core antigen (HBc) mutations and the post-operative prognosis of HBV-related hepatocellular carcinoma (HCC). In total, 98 patients suffering from HBV-related HCC and treated with surgery were enrolled, with a 48 month follow-up. The preCore/Core region of the HBV genome from tumour tissue (TT) and paired adjacent non-tumour tissue (ANTT) of these patients was sequenced, and a phylogenetic tree was reconstructed. The correlations between the viral features and evolutionary divergence of preCore/Core amino acid sequences from 67 paired TTs and ANTTs were analysed. Cox proportional hazard model analysis was applied for post-operative hazard risk evaluation. Phylogenetic analysis revealed that all of the sequences were ascribed to genotype C. The evolutionary divergence of amino acid sequences from matched TTs and ANTTs was significantly negatively correlated with serum and intrahepatic HBV DNA levels. Multivariate analysis showed that the HBc E77 mutation was associated with shorter overall survival, and HBc S87 and P156 mutations were independent risk factors for relapse. Furthermore, in contrast to with patients without the S87 mutation, no correlation was observed between serum HBV DNA and intrahepatic HBV DNA in HCC patients with the S87 mutation. Analysis of the intrahepatic sequence may advance our understanding of viral status; thus, it is useful for prognosis prediction for HBV-related HCC.

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.